Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy density functional”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 559 records · Page 31

Datasets for Custom-trained Machine-learning Interatomic Potentials: Nitric Acid Aqueous Solution

This dataset was generated using an iterative active learning strategy with the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials (MLIPs) for aqueous nitric acid. Each active-learning cycle consisted of three stages: (1) training, (2) exploration, and (3) labeling. The initial training set comprised approximately 800 randomly selected configurations from a previous study by Lewis et al. (https://doi.org/10.1021/jp205510q), which investigated nitric acid solutions at 2, 3, 4, and 5 mol/L. For all configurations, single-point calculations of atomic forces and total energies were performed at the quantum density functional theory BLYP-D2 and PBE-D3 levels of theory using the CP2K Quickstep module. Valence electrons were treated explicitly, while core electrons on all atoms were represented by norm-conserving Goedecker–Teter–Hutter (GTH) pseudopotentials. Long-range dispersion interactions were accounted for using Grimme dispersion corrections. Wave functions were expanded in a mixed Gaussian-and-plane-wave scheme using TZV2P-MOLOPT basis sets for all elements and an 800 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent field convergence was accelerated using orbital transformation and Direct Inversion in the Iterative Subspace, with a convergence threshold of 10^{-6}. All single-point calculations were carried out in periodic orthorhombic cells whose dimensions match those of the molecular configurations sampled from earlier trajectories. The CELL_REF keyword in CP2K was used to define a fixed reference cell, ensuring consistency in the reference data used for MLIP training, particularly when cell fluctuations are present in NpT simulations. The resulting high-fidelity energies and forces constitute the ground-truth labels used to train the MLIPs contained in this dataset.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗

Ni coarsening under humid atmosphere in the electrode of solid oxide cells: A combined study of density-functional theory and phase-field modeling

Ni coarsening in solid oxide cell (SOC) electrodes is known to be significantly faster under a humid atmosphere. The underlying mechanisms, though, have not been fully understood. In this work, we examine the surface diffusion of Ni(OH) x by a combination of density-functional theory and phase-field modeling. Here, density-functional theory is used to evaluate the adsorption and diffusion of the Ni(OH) x species on Ni (111) surface, and the results are used to obtain the effective surface diffusivity of Ni versus steam and hydrogen gas partial pressures. This diffusivity is used as an input for a phase-field model to investigate the Ni coarsening in SOC electrodes. It is found that Ni(OH) x formed through chemical reaction cannot accelerate coarsening unless an extremely high steam to hydrogen ratio is reached. However, if Ni(OH) x is formed through electrochemical reactions near triple phase boundaries (TPBs), surface diffusion of Ni(OH) x may cause faster Ni coarsening in fuel cell mode under a large overpotential. Specifically, surface diffusion of Ni(OH) may be comparable to or faster than that of Ni under an overpotential that is large but still possible under fuel cell operating conditions.

25 ENERGY STORAGE↗

Cobalt–Nickel Exchange in Exfoliated Battery Layered Cathode Sheets with Application to Recycling

With the emerging dominance of electric vehicles (EV) in the transportation sector, recycling or upcycling spent battery materials will be required to reduce EV costs, lessen waste, and ease critical material supply chain issues for EV batteries. Here, motivated by work in the literature describing the exfoliation of layered oxides, first‐principles calculations are performed to show that Li x CoO 2 , if exfoliated into nanosheets, can readily undergo transition metal cation exchange in aqueous media. The substitution of Co 3+ or Co 4+ cations inside the sheet by Ni 2+ is associated with modest reaction barriers (Δ G* ≈ 0.3–0.7 eV) and is at most mildly endothermic (Δ G ≈ 0.2–0.3 eV). In contrast, previous battery degradation studies have shown that Co 3+ diffusion is strongly inhibited inside bulk layered oxides. This suggests that processing spent layered oxides as nanosheets can provide a potentially low‐energy‐cost pathway to altering the transition metal and/or dopant stoichiometry, which can be used toward developing new room‐temperature upcycling routes for cathodes from end‐of‐life batteries.

battery recycling↗

Pressure-Dependent Magnetic Properties of Quasi-2D Cr 2 Si 2 Te 6 and Mn 3 Si 2 Te 6

Recently, pressure has been used to induce structural and magnetic phase transitions in many layered quantum materials whose layers are linked by van der Waals forces. Materials with such weakly held layers allow for relatively easy manipulation of the superexchange mechanism, which can give rise to novel magnetic behavior. Here, using hydrostatic pressure as a disorderless means to manipulate the interlayer coupling, we applied pressure on two quasi-2D sister compounds, namely, Cr 2 Si 2 Te 6 (CST) and Mn 3 Si 2 Te 6 (MST), up to ~1 GPa. Magnetic property measurements with the application of pressure revealed that the ferromagnetic transition temperature decreases in CST, while the opposite occurs for the ferrimagnetic MST. In MST, magnetization decreases with the increase in pressure, and such a trend is not clearly observed within the pressure range studied for CST. The overall pressure effect on magnetic characteristics such as exchange couplings and magnetic anisotropy energies is also examined theoretically using density functional theory. Exchange coupling in MST is strongly frustrated, and the first nearest neighbor interaction is the most dominant of the components with the strongest pressure dependence. In CST, the exchange coupling parameters exhibit very little dependence on pressure. This combined experimental and theoretical work has the potential to expand to other relevant quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Machine learning to accelerate screening for Marcus reorganization energies

Understanding and predicting the charge transport properties of π-conjugated materials is an important challenge for designing new organic electronic devices, such as solar cells, plastic transistors, light-emitting devices, and chemical sensors. A key component of the hopping mechanism of charge transfer in these materials is the Marcus reorganization energy which serves as an activation barrier to hole or electron transfer. While modern density functional methods have proven to accurately predict trends in intramolecular reorganization energy, such calculations are computationally expensive. In this work, we outline active machine learning methods to predict computed intramolecular reorganization energies of a wide range of polythiophenes and their use toward screening new compounds with low internal reorganization energies. Our models have an overall root mean square error (RMSE) of ±0.113 eV, but a much smaller RMSE of only ±0.036 eV on the new screening set. Since the larger error derives from high-reorganization energy compounds, the new method is highly effective to screen for compounds with potentially efficient charge transport parameters.

14 SOLAR ENERGY↗

Electrochemical CO 2 Reduction Reaction over Cu Nanoparticles with Tunable Activity and Selectivity Mediated by Functional Groups in Polymeric Binder

Electrochemical carbon dioxide reduction reaction (CO 2 RR) using copper (Cu)-based catalysts has received significant attention mainly because Cu is an element capable of producing hydrocarbons and oxygenates. One possible way to control the CO 2 RR performance at the electrode interface is by modifying catalysts with specific functional groups of different polymeric binders, which are necessary components in the process of electrode fabrication. However, the modification effect of the key functional groups on the CO 2 RR activity and selectivity is poorly understood over Cu-based catalysts. In this work, the role of functional groups (e.g., -COOH and -CF 2 groups) in hydrophilic and hydrophobic polymeric binders on the CO 2 RR of Cu-based catalysts is investigated using a combination of electrochemical measurements, in-situ characterization and density functional theory (DFT) calculations. DFT results reveal that functional groups influence the binding energies of key intermediates involved in both CO 2 RR and the competing hydrogen evolution reaction, consistent with experimental observation of binder-dependent product distributions among formic acid, CO, CH 4 , and H 2 . This study provides a fundamental understanding that the selection of desired polymeric binders is a useful strategy for tuning the CO 2 RR activity and selectivity.

25 ENERGY STORAGE↗

A first-principles machine-learning force field for heterogeneous ice nucleation on microcline feldspar

The formation of ice in the atmosphere affects precipitation and cloud properties, and plays a key role in the climate of our planet. Although ice can form directly from liquid water at deeply supercooled conditions, the presence of foreign particles can aid ice formation at much warmer temperatures. Over the past decade, experiments have highlighted the remarkable efficiency of feldspar minerals as ice nuclei compared to other particles present in the atmosphere. However, the exact mechanism of ice formation on feldspar surfaces has yet to be fully understood. Here, we develop a first-principles machine-learning model for the potential energy surface aimed at studying ice nucleation at microcline feldspar surfaces. The model is able to reproduce with high-fidelity the energies and forces derived from density-functional theory (DFT) based on the SCAN exchange and correlation functional. Our training set includes configurations of bulk supercooled water, hexagonal and cubic ice, microcline, and fully-hydroxylated feldspar surfaces exposed to vacuum, liquid water, and ice. In this work, we apply the machine-learning force field to study different fully-hydroxylated terminations of the (100), (010), and (001) surfaces of microcline exposed to vacuum. Our calculations suggest that terminations that do not minimize the number of broken bonds are preferred in vacuum. We also study the structure of supercooled liquid water in contact with microcline surfaces, and find that water density correlations extend up to around 10 Å from the surfaces. Finally, we show that the force field maintains a high accuracy during the simulation of ice formation at microcline surfaces, even for large systems of around 30,000 atoms. Future work will be directed towards the calculation of nucleation free energy barriers and rates using the force field developed herein, and understanding the role of different microcline surfaces on ice nucleation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The role of H–H interactions and impurities on the structure and energetics of H/Pd(111)

Understanding hydrogen incorporation into palladium requires detailed knowledge of surface and subsurface structure and atomic interactions as surface hydrogen is being embedded. Using density functional theory (DFT), we examine the energies of hydrogen layers of varying coverage adsorbed on Pd(111). Here we find that H–H and H–Pd interactions promote the formation of the well-known ($\sqrt{3}$ x $\sqrt{3}$) phases but also favor an unreported (3 × 3) phase at high H coverages for which we present experimental evidence. We relate the stability of isolated H vacancies of the (3 × 3) phase to the need of H 2 molecules to access bare Pd before they can dissociate. Following higher hydrogen dosage, we observe initial steps of hydride formation, starting with small clusters of subsurface hydrogen. The interaction between H and Pd is complicated by the persistent presence of carbon at the surface. X-ray photoelectron spectroscopy experiments show that trace amounts of carbon, emerging from the Pd bulk despite many surface cleaning cycles, become mobile enough to repopulate the C-depleted surface at temperatures above 200 K. When exposed to hydrogen, these surface carbon atoms react to form benzene, as evidenced by scanning tunneling microscopy observations interpreted with DFT.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spectroscopic Probe Molecule Selection Using Quantum Theory, First-Principles Calculations, and Machine Learning

Probe molecule vibrational spectra have a long history of being used to characterize materials including metals, oxides, metal-organic frameworks, and even human proteins. Furthermore, recent advances in machine learning have enabled computationally generated spectra to aid in detailed characterization of complex surfaces with probe molecules. Despite widespread use of probe molecules, the science of probe molecule selection is underdeveloped. Here, we develop physical concepts, including orbital interaction energy and the energy overlap integral, to explain and predict the ability of probe molecules to discriminate structural descriptors. We resolve the crystal orbital overlap population (COOP) to specific molecular orbitals and quantify their bonding character, which directly influences vibrational frequencies. Using only a single adsorbate calculation from density function theory (DFT), we compute the interaction energy of individual adsorbate molecular orbitals with adsorption site atomic orbitals across many different sites. Combining the molecular orbital resolved COOP and changes in orbital interaction energy enables probe molecule selection for improved discrimination of various sites. We demonstrate these concepts by comparing the predicted effectiveness of carbon monoxide (CO), nitric oxide (NO), and ethylene (C 2 H 4 ) to probe Pt adsorption sites. Finally, using a previously developed machine learning framework, we show that models trained on hundreds of thousand C 2 H 4 spectra, computed from DFT, which regress surface binding-type and generalized coordination number (GCN), outperform those trained using CO and NO spectra. Lastly, a python package, pDOS_overlap, for implementing the electron density based analysis on any combination of adsorbates and materials, is also made available.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Geodyn Material Library: Pseudocap models for dry porous tocks

This report describes the second edition of the Pseudocap Strength models for porous rocks implemented in GEODYN material library. The first model was developed in 2007 and calibrated for concrete. Then, the model parameters were calibrated based on triaxial tests reported for limestones and sandstones of various porosities. In these models some key parameters were chosen as functions of the reference porosity,Φ. Since then multiple modifications were implemented in the model, therefore, it has been recalibrated for some common porous materials such as limestones, sandstones, alluvium, tuffs and granite. Two types of models are described in this repot. The first type (called Pseudocap Model or PM) is for rocks from a specific location. Parameters were calibrated for several specific geologic materials. The second type (called Generic Pseudocap Model or GPM) is useful for the sites where only basic information (rock type, porosity) is available. Generic models include built-in correlations between porosities and other mechanical properties observed for certain rock types. The models of both types were validated by comparing not only to quasi-static triaxial tests for these materials but also to shock Hugoniot data and spherical explosion data for some materials. All models were derived in the frame of isotropic plasticity. They are designed to be used in explicit finite element/difference codes. Tangent stiffness tensor is not provided but can be calculated numerically for the model to be used in implicit finite element codes. For an isotropic material the stress can be decomposed into volumetric and deviatoric parts. The volumetric part is modeled using an Equation of state (EOS) which calculates the pressure and the bulk sound speed as functions of the internal specific energy and density. Here a simple, Mie-Gruneisen EOS is presented, but tabulated EOS (LEOS) provided by the library can be used as well. The stress is limited by the yield surface which depends on three invariants of the stress tensor and specific internal energy. In addition, to capture the strain-rate dependence a simple multiplier is used for the yield surface which depends on the equivalent plastic strain rate. The failure surface (the ultimate yield, Y f , defined later) is chosen in the Hoek-Brown form, commonly used in rock mechanics. It includes measurable parameters such as Unconfined Compressive Strength (UCS) as well as scale parameters characterizing the quality of the rock such as GSI (Geologic Strength Index). Thus, even though the model is calibrated for small samples it offers a way to extrapolate the strength to the field scale using geological characterization of the rock mass. The model captures effects of brittle-ductile transition in rocks by introducing a cap multiplier to the yield function. The rate of dilatancy (bulking) is proportional to the slope of the yield surface affected by the cap. Therefore, it takes place only at low confinements when the pressure is less than the brittle-ductile transition pressure, P BD . On the contrary, the porous compaction takes place at pressures higher than P BD . The cap moves as the porosity is compacted or new porosity is generated due to dilatancy. The porous compaction is modeled using an evolution equation which includes deviatoric stress so that the onset of compaction corresponds to the cap surface. The model captures effects of shear-enhanced compaction which is an important for porous rocks. Section 2 describes the modeling framework and Section 3 presents the model calibration procedure. Section 4 compares experimental data for various rocks versus model predictions. The model parameters used for this comparison are given in Appendix. The files with material constants are available with the latest GEODYN material library distribution.

58 GEOSCIENCES↗

Machine Learning-Assisted Recovery of Delicate Kinetic Information from Transient Reactor Experiments

Identifying active sites and their roles in chemical reaction steps remains a vital challenge in heterogeneous catalysis. Transient experiments offer a unique way to probe active sites and distinguish subtle kinetic features. Although physics-based analysis methods may be well-developed, they can be highly susceptible to experimental noise, and smoothing methods may erase or even distort important features; a smooth curve is not always the best curve. We demonstrate a new workflow for the direct interpretation of intrinsic kinetic information from exit flux curves measured in transient reactor experiments. This workflow contains three artificial neural networks (ANNs), including a noise reducer, a concentration predictor, and a rate predictor to analyze experimental data, followed by the virtual TAP (VTAP) physics-based reactor model and density functional theory (DFT) calculations of adsorption energies on specific sites. We use this workflow to analyze the data from experiments titrating Pt/Al 2 O 3 and Pt/SiO 2 catalysts with carbon monoxide (CO) in the temporal analysis of products (TAP) reactor. Our workflow separates the time-evolving chemical reaction and mass transfer information contained in the TAP pulse response. The existence of strong- and weak-binding sites on the Pt/Al 2 O 3 catalyst is observed in the catalyst titration experiment in the transient reactor. The structures of the strong- and weak-binding sites are then identified by using DFT calculations. We find that the Pt/SiO 2 catalyst has only strong-binding sites, which aligns with the inactive support effect of SiO 2 . We demonstrate how machine learning methods provide unique insights with high-resolution data analysis that cannot be achieved by using state-of-the-art physics-based methods.

Adsorption↗

Role of Pairwise Reactions on the Synthesis of Li 0.3 La 0.57 TiO 3 and the Resulting Structure–Property Correlations

The performance of single-ion conductors is highly sensitive to the material’s defect chemistry. Tuning these defects is limited for solid-state reactions as they occur at particle–particle interfaces, which provide a complex evolving energy landscape for atomic rearrangement and product formation. In this report, we investigate the (1) order of addition and (2) lithium precursor decomposition temperature and their effect on the synthesis and grain boundary conductivity of the perovskite lithium lanthanum titanium oxide (LLTO). We use an intimately mixed sol–gel, a solid-state reaction of Li precursor + La 2 O 3 + TiO 2 , and Li precursor + amorphous La 0.57 TiO x as different chemical routes to change the way in which the elements are brought together. The results show that the perovskite can accommodate a wide range of Li deficiencies (upward of 50%) while maintaining the tetragonal LLTO structure, indicating that X-ray diffraction (XRD) is insufficient to fully characterize the chemical nature of the product (i.e., Li-deficient LLTO may behave differently than stoichiometric LLTO). Variations in the relative intensities of different reflections in XRD suggest variations in the La ordering within the crystal structure between synthesis methods. Furthermore, the choice of the precursor and the order of addition of the reactants lower the time required to form a pure phase. Density functional theory calculations of the formation energy of possible reaction intermediates support the hypothesis that a greater thermodynamic driving force to form LLTO leads to a greater LLTO yield. Here, the retention of lithium is correlated with the thermal decomposition temperature of the Li precursor and the starting material mixing strategy. Taking the results together suggests that cations that share a site with Li should be mixed early to avoid ordering. Such cation ordering inhibits Li motion, leading to higher Li ion resistance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A numerical Poisson solver with improved radial solutions for a self-consistent locally scaled self-interaction correction method

Abstract The universal applicability of density functional approximations is limited by self-interaction error made by these functionals. Recently, a novel one-electron self-interaction-correction (SIC) method that uses an iso-orbital indicator to apply the SIC at each point in space by scaling the exchange-correlation and Coulomb energy densities was proposed. The locally scaled SIC (LSIC) method is exact for the one-electron densities, and unlike the well-known Perdew–Zunger SIC (PZSIC) method recovers the uniform electron gas limit of the uncorrected density functional approximation, and reduces to PZSIC method as a special case when isoorbital indicator is set to the unity. Here, we present a numerical scheme that we have adopted to evaluate the Coulomb potential of the electron density scaled by the iso-orbital indicator required for the self-consistent LSIC calculations. After analyzing the behavior of the finite difference method (FDM) and the green function solution to the radial part of the Poisson equation, we adopt a hybrid approach that uses the FDM for the Coulomb potential due to the monopole and the GF for all higher-order terms. The performance of the resultant hybrid method is assessed using a variety of systems. The results show improved accuracy than earlier numerical schemes. We also find that, even with a generic set of radial grid parameters, accurate energy differences can be obtained using a numerical Coulomb solver in standard density functional studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Hot and dense matter equation of state probability distributions for astrophysical simulations

We add an ensemble of nuclei to the equation of state for homogeneous nucleonic matter to generate a new set of models suitable for astrophysical simulations of core-collapse supernovae and neutron star mergers. We implement empirical constraints from (i) nuclear mass measurements, (ii) proton-proton scattering phase shifts, and (iii) neutron star observations. Our model is also guided by microscopic many-body theory calculations based on realistic nuclear forces, including the zero-temperature neutron matter equation of state from quantum Monte Carlo simulations and thermal contributions to the free energy from finite-temperature many-body perturbation theory. We ensure that the parameters of our model can be varied while preserving thermodynamic consistency and the connection to experimental or observational data, thus providing a probability distribution of the astrophysical hot and dense matter equation of state. Furthermore, we compare our results with those obtained from other available equations of state. While our probability distributions indeed represent a large number of possible equations of state, we cannot yet claim to have fully explored all of the uncertainties, especially with regard to the structure of nuclei in the hot and dense medium.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Filament Formation in TaO x Thin Films for Memristor Device Application: Modeling Electron Energy Loss Spectra and Electron Transport

Abstract Although understanding filament formation in oxide‐based memristive devices by theory has emerged, there are still fundamental unanswered questions. Importantly, for practical application of thin films the material in its amorphous state is to be considered, but mostly lacking so far, and details on sub‐stoichiometry are also scarce. To gain insight into the optical and electronic properties of sub‐stoichiometric amorphous tantalum oxide (TaO x ), the electron energy loss spectrum (EELS) of model systems is characterized theoretically and electron transport characteristics are analyzed in detail. Calculated blue‐shifts by increasing sub‐stoichiometry explained the measurements, potentially suggesting estimation of oxygen vacancy concentrations through EEL spectra. Electron transport results based on TaO x material models validated by EELS measurements show that oxygen vacancy filamentary paths are initiated at low bias upon increasing sub‐stoichiometry yet noting an interplay with the local amorphous structure. Contact resistances at interfaces of the TaO x switching layer and a tantalum scavenging layer or titanium nitride electrode are quantified, indicating the possibility for either oxygen vacancy‐ or metal cluster‐based conduction mechanisms at the interface. The computational work, combined with experimental characterization for validation, provides a basis for investigating effects of sub‐stoichiometry on filament formation in TaO x thin film memristive devices.

36 MATERIALS SCIENCE↗

Real-Space Local Dynamics in 1,2,3-Triazole Using Inelastic Neutron Scattering

Enhancing proton transport in polymer electrolytes is crucial for advancing next-generation solid-state batteries, yet our understanding of proton conductivity in nonaqueous environments remains limited due to a lack of atomic-scale insights. Here, in this study, we investigated the atomic-scale dynamics of 1,2,3-triazole, a small molecule capable of dynamic hydrogen bonding, as a model system for proton hopping in nonaqueous environments. Using the real-space correlation function determined by the double Fourier transformation of inelastic neutron scattering spectra, we identified that the self-motion of protons and intermolecular dynamics occur on comparable time scales. Furthermore, we observed that the activation energy associated with the intermolecular dynamics matches the energy barrier for molecular rotations determined through Density Functional Theory calculations. These findings underscore the importance of controlling molecular dynamics at the atomic scale to control proton transport. Additionally, we demonstrated that intermolecular dynamics in systems involving protons can be studied using inelastic neutron scattering even without deuteration, thereby providing a broader avenue for studying atomic-scale dynamics in soft matter systems.

Shinohara, Yuya [Oak Ridge National Laboratory (OR↗

Understanding Cation Selectivity in Carbon Nanopores with Hybrid First-Principles/Continuum Simulations: Implications for Water Desalination and Separation Technologies

Understanding ion adsorption in porous carbons is crucial for a range of technologies, including water desalination and energy storage. In this work, we combined density functional theory with a continuum solvation model to investigate thermodynamics and kinetics of the adsorption process of alkali metal ions from aqueous solutions into carbon nanopores with different sizes and geometries. We found that cations with a larger ionic radius are more favorable to enter the nanopores because of a lower energy penalty of dehydration. In addition, the pore size and geometry were found to have a significant impact on the ion–pore interaction under confinement and cation selectivity. Our study highlights a complex interplay among nanopore geometry, ion size, and hydration on the cation adsorption selectivity, suggesting that tuning the porosity could represent a general strategy for improving ion separations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine learning defect properties in Cd-based chalcogenides

Impurity energy levels in the band gap can have serious consequences for a semiconductor's performance as a photovoltaic absorber. Data-driven approaches can help accelerate the prediction of point defect properties in common semiconductors, and thus lead to the identification of potential deep lying impurity states. In this work, we use density functional theory (DFT) to compute defect formation energies and charge transition levels of hundreds of impurities in CdX chalcogenide compounds, where X = Te, Se or S. We apply machine learning techniques on the DFT data and develop on-demand predictive models for the formation energy and relevant transition levels of any impurity atom in any site. The trained ML models are general and accurate enough to predict the properties of any possible point defects in any Cd-based chalcogenide, as we prove by testing on a few selected defects in mixed chalcogen compounds CdTe 0.5 Se 0.5 and CdSe 0.5 S 0.5 . The ML framework used in this work can be extended to any class of semiconductors.

density functional theory↗