Hybrid Density Functional Theory Research Articles

Strain and Electric Field Engineering of G-ZnO/SnXY (X, Y = S, Se) S-Scheme Heterostructures for Photocatalyst and Electronic Device Applications: A Hybrid DFT Calculation.

Using hybrid density functional theory calculations, we systematically study the biaxial strain and electric field modulated electronic properties of g-ZnO/SnS2, g-ZnO/SnSe2, and g-ZnO/SnSSe S-scheme van der Waals heterostructures (vdWHs). g-ZnO/SnS2 and g-ZnO/SnSSe are found to be promising photocatalysts for water splitting with high solar-to-hydrogen efficiencies, even under acidic, alkaline, and high-stress conditions. The strain effect on the bandgaps of g-ZnO/SnXY is explained in detail according to the correlation between geometry structure and orbital hybridization of SnXY, which could help understand the strain-induced band structure evolutions in other SnXY (X, Y = S, Se)-based vdWHs. It is surprising that under an external electric field, g-ZnO/SnS2, g-ZnO/SnSe2, and g-ZnO/SnSSe can offer the occupied nearly free-electron (NFE) states. In many materials, NFE states are usually unoccupied and is not conducive to the charge transport. The NFE state in g-ZnO/SnSe2 is the most sensitive to the electric field and might be promising electron transport channel in nanoelectronic devices. g-ZnO/SnSe2 might also have application potential in gas sensors and high-temperature superconductors.

Benchmark Investigation of SCC-DFTB against Standard and Hybrid DFT to Model Electronic Properties in Two-Dimensional MOFs for Thermoelectric Applications.

Recent studies have shown that metal-organic frameworks (MOFs) have potential as thermoelectric materials, and the topic has received increasing attention. The main motivation for this project is to further our knowledge of thermoelectric properties in MOFs and find which available self-consistent-charge density functional tight binding (SCC-DFTB) method can best predict (at least trends in) the electronic properties of MOFs at a lower computational cost than standard density functional theory (DFT). In this work, the electronic properties of monolayer, serrated, AA-stacked, and/or AB-stacked Zn3C6O6, Cd3C6O6, Zn-NH-MOF─for which no previous calculations of thermoelectric performance exist─and Ni3(HITP)2 MOFs are modeled with DFT-PBE, DFT-HSE06, GFN1-xTB, GFN2-xTB, and DFTB-3ob/mio. The band structures, density of states, and their relative orbital contributions, as well as the electrical conductivity, Seebeck coefficient, and power factor, are compared across methods and geometries. Our results suggest that GFN-xTB is adequate to predict the MOFs' band structure shape and density of states but not band gap. Our calculations further indicate that Zn3C6O6, Cd3C6O6, and Zn-NH-MOF have higher power factor values than Ni3(HITP)2, one of the highest performing synthesized MOFs, and are therefore promising for thermoelectric applications.

Topological phase transition and tunable surface states in YBi.

A unique co-existence of extremely large magnetoresistance (XMR) and topological characteristics in non-magnetic rare-earth monopnictides has stimulated intensive research on these materials. Yttrium monobismuthide (YBi) has been reported to exhibit XMR up to 105% but its Topological properties still need clarification. Here we use the hybrid density functional theory to probe the structural, electronic, and topological properties of YBi in detail. We observe that YBi is topologically trivial semimetal at ambient pressure which is in accordance with reported experimental results. The topological phase transitions i.e., trivial to non-trivial are obtained with volumetric pressure of 6.5 GPa and 3% of epitaxial strain. This topological phase transitions are well within the structural phase transition of YBi (24.5 GPa). The topological non-trivial state is characterized by band inversions among Y-d band and Bi-p band at Γ- and X-point which is further verified with the help of surface band structure along (001) plane. The Z2 topological invariants are calculated with the help of product of parities and evolution of Wannier charge centers. The occurrence of non-trivial phase in YBi with a relatively small epitaxial strain, which a thin film geometry can naturally has, make it an ideal candidate to probe inter-relationship between XMR and non-trivial topology.

Stable, while Still Active? A DFT Study of Cu, Ag, and Au Single Atoms at the C3N4/TiO2 Interface.

Hybrid DFT calculations are employed to compare the adsorption and stabilization of Cu, Ag, and Au atoms on graphitic C3N4 and on the heterojunction formed by g- C3N4 and TiO2. While Cu and Ag can be strongly chemisorbed in form of cations on g- C3N4, Au is only weakly physisorbed. On g- C3N4/TiO2, all coinage metal adatoms can be strongly chemisorbed, but, while Cu and Ag forms cations, Au form an Au- species. Ab Initio Molecular Dynamics simulations confirm that the metal adatoms on g-C3N4 are highly mobile at room temperature, while they remain confined in the interfacial spacing between C3N4 and TiO2 on the heterojunction, being both stably bound and reachable for the reactants in a catalytic cycle. Doping g-C3N4/TiO2 with metal single atoms permits thus to generate catalytic systems with tunable charge and chemical properties and improved stability with respect to bare C3N4. Moreover, the changes in the electronic structure of g-C3N4/TiO2 induced by the presence of the metal single atoms are beneficial also for photocatalytic applications.

Electronic Structure of Mg-, Si-, and Zn-Doped SnO2 Nanowires: Predictions from First Principles

We investigated the electronic structure of Mg-, Si-, and Zn-doped four-faceted [001]- and [110]-oriented SnO2 nanowires using first-principles calculations based on the linear combination of atomic orbitals (LCAO) method. This approach, employing atomic-centered Gaussian-type functions as a basis set, was combined with hybrid density functional theory (DFT). Our results show qualitative agreement in predicting the formation of stable point defects due to atom substitutions on the surface of the SnO2 nanowire. Doping induces substantial atomic relaxation in the nanowires, changes in the covalency of the dopant–oxygen bond, and additional charge redistribution between the dopant and nanowire. Furthermore, our calculations reveal a narrowing of the band gap resulting from the emergence of midgap states induced by the incorporated defects. This study provides insights into the altered electronic properties caused by Mg, Si, and Zn doping, contributing to the further design of SnO2 nanowires for advanced electronic, optoelectronic, photovoltaic, and photocatalytic applications.

Sparse-Stochastic Fragmented Exchange for Large-Scale Hybrid Time-Dependent Density Functional Theory Calculations.

We extend our recently developed sparse-stochastic fragmented exchange formalism for ground-state near-gap hybrid DFT to calculate absorption spectra within linear-response time-dependent generalized Kohn-Sham DFT (LR-GKS-TDDFT) for systems consisting of thousands of valence electrons within a grid-based/plane-wave representation. A mixed deterministic/fragmented-stochastic compression of the exchange kernel, here using long-range explicit exchange functionals, provides an efficient method for accurate optical spectra. Both real-time propagation as well as frequency-resolved Casida-equation-type approaches for spectra are presented, and the method is applied to large molecular dyes.

Hybrid Density Functional Theory Calculations for the Crystal Structure and Electronic Properties of Al3+ Doped KDP Crystals

Intentionally adding select ions such as Al3+ could be helpful in controlling the crystal habit of KDP crystal for high yield of optics. The study of how Al3+ ions affect crystal quality can provide a basis for selecting an appropriate doping level without negatively affecting the optical properties of crystals. Here, the influence of Al3+ ions on the crystal structure and properties of KDP crystals have been investigated by using first-principles calculations. Theoretical calculations show that Al3+ ions mainly replace K sites in KDP crystals and could complex with intrinsic VH− point defects to form AlK2+ + 2VH− cluster defects. The linear absorption spectra indicate that the presence of Al3+ ions has minimal impact on the linear absorption of KDP crystals, aligning well with the experimental findings. And Al3+ ions could cause a slight shortening of the band gap of KDP crystals. However, these ions could bring significant deformations of O-H bonds. As the concentration of Al3+ ions increase, more O-H bonds linking to PO4 groups are distorted in KDP crystals. As a result, the structural instability could be fast enhanced with increasing the defect concentration. Therefore, high concentrations of Al3+ ions could cause the instability of the crystal structure, which finally affects the laser-induced damage resistance of the KDP crystals. This manuscript contributes to a more comprehensive understanding of the physical mechanisms by which different impurity ions affect the optical properties of KDP crystals.

Deep polaronic acceptors in LiGa5O8

Recently, LiGa5O8 was claimed to be a p-type dopable ultrawide-bandgap oxide, based on measurements of undoped material. Here, the electronic properties of potential acceptor dopant impurities in LiGa5O8 are calculated using hybrid density functional theory to evaluate their potential for causing p-type conductivity. As with the related compound LiGaO2, the heavy oxygen-derived valence bands lead to stable self-trapped holes in LiGa5O8. Acceptor defects and dopants also bind trapped holes (or small polarons), which lead to large acceptor ionization energies. The calculations here indicate that neither native acceptor defects (such as cation vacancies or antisites) nor impurity dopants can give rise to p-type conductivity in LiGa5O8. Optical transitions associated with these defects are also calculated, in order to allow for possible experimental verification of their behavior.

High-throughput screening of dual atom catalysts for oxygen reduction and evolution reactions and rechargeable zinc-air battery

We provide the rational design of dual atom catalysts (DACs) supported on nitrogen-doped graphene for the oxygen reduction reaction (ORR) and oxygen evolution reaction (OER) through high-throughput computational screening of M1M2N6-DAC systems, where M1 and M2 represent Fe, Co, Ni, Ru, Rh, Pd, Os, Ir, or Pt metals. We predict that FeRuN6-DAC at the summit of the volcano plot exhibits a low theoretical ORR overpotential (ηORR) of 0.24 V and a low theoretical OER overpotential (ηOER) of 0.19 V. The low ηORR and ηOER result from the catalytic performance of the Fe site being tuned to electronic properties that facilitate adsorption and desorption of the OH* intermediate. Inspired by these hybrid density functional theory (DFT) computational and machine learning (ML) results, we synthesized FeRuN6-DAC, FeN4-SAC, and RuN4-SAC and characterized them using X-ray photoelectron spectroscopy (XPS), X-ray absorption spectroscopy (XAS), scanning transmission electron microscopy (STEM), and in-situ electron spin resonance (ESR). Our in-situ ESR spectroscopy signifies that the spin of the Fe active site increases with increasing applied potential due to the increase in the concentration of OH* intermediate on Fe. We verified experimentally the predicted catalytic performances, finding that FeRuN6-DAC leads to an experimental ORR overpotential of 0.29 V with a Tafel slope of 104 mVdec−1 and an OER overpotential of 0.27 V with a Tafel slope of 124 mVdec−1. The rechargeable Zinc-air battery setup was fabricated with FeRuN6-DAC in place of the cathode, showing a maximum power density of 0.45 W/cm2 at the current density of 0.44 A/cm2 and good stability after 120 cycles. According to our findings, we demonstrate that DFT-guided strategies are useful for designing advanced DACs applicable to ORR, OER, and Zinc-air battery applications.

Prediction of Feasibility of Polaronic OER on (110) Surface of Rutile TiO2.

The polaronic effects at the atomic level hold paramount significance for advancing the efficacy of transition metal oxides in applications pertinent to renewable energy. The lattice-distortion mediated localization of photoexcited carriers in the form of polarons plays a pivotal role in the photocatalysis. This investigation focuses on rutile TiO2, an important material extensively explored for solar energy conversion in artificial photosynthesis, specifically targeting the generation of green H2 through photoelectrochemical (PEC) H2O splitting. By employing Hubbard-U corrected and hybrid density functional theory (DFT) methods, we systematically probe the polaronic effects in the catalysis of oxygen evolution reaction (OER) on the (110) surface of rutile TiO2. Theoretical understanding of polarons within the surface, coupled with simulations of OER at distinct titanium (Ti) and oxygen (O) active sites, reveals diverse polaron formation energies within the lattice sites with strong preference for bulk and surface bridge (Ob) oxygen sites. Moreover, we provide the evidence for the facilitative role of polarons in OER. We find that hole polarons situated at the equatorial oxygen sites near the Ti-active site, along with bridge site hole polarons distal from the Ob active site yield a small reduction in OER overpotential by ~0.06 eV and ~0.12 eV, respectively. However, subsurface, equatorial, and bridge site hole polarons significantly reduce the Ti-active site OER overpotential by ~0.4 eV through the peroxo-type oxygen pathway. We also observe that the presence of hole polarons stabilizes the *OH, *O, and *OOH intermediate species compared to the scenario without hole polarons. Overall, this study provides a detailed mechanistic insight into polaron-mediated OER, offering a promising avenue for improving the catalytic activity of transition metal oxide-based photocatalysts catering to renewable energy requisites.

New Charge-Transfer Complexes of Organochalcogenide Compound Based on Aryl Acetamide Group with Quinones: Synthesis, Characterization, Antioxidant, and Computational Study

This study aims to prepare charge transfer complexes derived from organochalcogenide of arylamide derivatives with different quinones. A new charge-transfer complexes have been developed through a direct reaction between (PhNHCOCH2)2Se, (o-CH3PhNHCOCH2)2Se, and (PhCH2NHCOCH2)2E, where E = S, Se, and Te are electron donors and different quinones are electron acceptors. The quinones used in the reaction were 2,3-dichloro-5,6-dicyanobenzoquinones (DDQ), 7,7’,8,8’-tetracyanoquinodimethane, and tetracyanoethane. The electron donors and electron acceptor mol were 1:1, and the reaction was conducted in acetonitrile. Infrared, 1H and 13C-NMR spectroscopic data characterized all complexes. The complexes’ antioxidant activity was evaluated through α,α-diphenyl-β-picrylhydrazyl at 10–0.312 mg/mL. The results showed that all complexes exhibited promising antioxidant activities. Among them, (PhCH2NHCOCH2)2S·DDQ compound had the least IC50 value of 6.725 mg/mL, indicating a potent scavenging property compared to other compounds. The molecular structures of charge-transfer complexes were investigated using hybrid density functional theory (B3LYP) and basis set 3-21G. We obtained geometrical structures' highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) surfaces and energy gaps through geometric optimization. We also investigated the molecular shapes and contours of the prepared compounds through geometrical optimization and compared the HOMO energy of the CT compounds to investigate donor and acceptor properties.

Temperature-transferable tight-binding model using a hybrid-orbital basis.

Finite-temperature calculations are relevant for rationalizing material properties, yet they are computationally expensive because large system sizes or long simulation times are typically required. Circumventing the need for performing many explicit first-principles calculations, tight-binding and machine-learning models for the electronic structure emerged as promising alternatives, but transferability of such methods to elevated temperatures in a data-efficient way remains a great challenge. In this work, we suggest a tight-binding model for efficient and accurate calculations of temperature-dependent properties of semiconductors. Our approach utilizes physics-informed modeling of the electronic structure in the form of hybrid-orbital basis functions and numerically integrating atomic orbitals for the distance dependence of matrix elements. We show that these design choices lead to a tight-binding model with a minimal amount of parameters that are straightforwardly optimized using density functional theory or alternative electronic-structure methods. The temperature transferability of our model is tested by applying it to existing molecular-dynamics trajectories without explicitly fitting temperature-dependent data and comparison with density functional theory. We utilize it together with machine-learning molecular dynamics and hybrid density functional theory for the prototypical semiconductor gallium arsenide. We find that including the effects of thermal expansion on the onsite terms of the tight-binding model is important in order to accurately describe electronic properties at elevated temperatures in comparison with experiment.

Hybrid density functional theory simulation of sodium impurity and impurity–vacancy defect complexes in germanium: perspectives of defect engineering for activation of shallow donors

At present, germanium (Ge) as a high-mobility material is considered as a possible replacement for Si in microelectronics. The main obstacle in the method is the large vacancy concentration obtained after the implantation of shallow donors. This prevents the formation of n+ regions and has a negative impact on device performance. One way to eliminate the detrimental effects of such defects is through the codoping approach. In this study, the behavior of Na impurities in germanium together with the formation of defect complexes with vacancies, divacancies and shallow donors (P and As) were investigated using hybrid density functional theory calculations. The site preference, diffusion barriers, charge states and binding energies of various complexes, as well as the activation of donor–vacancy complexes with both F and Na codoping were investigated. For this purpose, two extreme cases were considered. Firstly, we calculate the concentrations of substitutional and interstitial Na together with monovacancies in the case of quasi-equilibrium conditions realized for doping from melt. Another extreme case is Na codoping during ion implantation. Our calculations show that Na passivates vacancies, divacancies and donor–vacancy complexes (E-centers) with a large energy gain, which prevents rapid donor diffusion similar to the case of doping with fluorine. We have shown that the passivation of E-centers with both F and Na impurities leads to the back transformation of such complexes into shallow donors. This means that codoping with Na may neutralize the harmful effects of vacancies and is a perspective for defect engineering.

Development of analytic gradients for the Huzinaga quantum embedding method and its applications to large-scale hybrid and double hybrid DFT forces.

The theory of analytic gradients is presented for the projector-based density functional theory (DFT) embedding approach utilizing the Huzinaga-equation. The advantages of the Huzinaga-equation-based formulation are demonstrated. In particular, it is shown that the projector employed does not appear in the Lagrangian, and the potential risk of numerical problems is avoided at the evaluation of the gradients. The efficient implementation of the analytic gradient theory is presented for approaches where hybrid DFT, second-order Møller-Plesset perturbation theory, or double hybrid DFT are embedded in lower-level DFT environments. To demonstrate the applicability of the method and to gain insight into its accuracy, it is applied to equilibrium geometry optimizations, transition state searches, and potential energy surface scans. Our results show that bond lengths and angles converge rapidly with the size of the embedded system. While providing structural parameters close to high-level quality for the embedded atoms, the embedding approach has the potential to relax the coordinates of the environment as well. Our demonstrations on a 171-atom zeolite and a 570-atom protein system show that the Huzinaga-equation-based embedding can accelerate (double) hybrid gradient computations by an order of magnitude with sufficient active regions and enables affordable force evaluations or geometry optimizations for molecules of hundreds of atoms.

Reduction of molecular oxygen in flavodiiron proteins - Catalytic mechanism and comparison to heme-copper oxidases

The family of flavodiiron proteins (FDPs) plays an important role in the scavenging and detoxification of both molecular oxygen and nitric oxide. Using electrons from a flavin mononucleotide cofactor molecular oxygen is reduced to water and nitric oxide is reduced to nitrous oxide and water. While the mechanism for NO reduction in FDPs has been studied extensively, there is very little information available about O2 reduction. Here we use hybrid density functional theory (DFT) to study the mechanism for O2 reduction in FDPs. An important finding is that a proton coupled reduction is needed after the O2 molecule has bound to the diferrous diiron active site and before the OO bond can be cleaved. This is in contrast to the mechanism for NO reduction, where both NN bond formation and NO bond cleavage occurs from the same starting structure without any further reduction, according to both experimental and computational results. This computational result for the O2 reduction mechanism should be possible to evaluate experimentally. Another difference between the two substrates is that the actual OO bond cleavage barrier is low, and not involved in rate-limiting the reduction process, while the barrier connected with bond cleavage/formation in the NO reduction process is of similar height as the rate-limiting steps. We suggest that these results may be part of the explanation for the generally higher activity for O2 reduction as compared to NO reduction in most FDPs. Comparisons are also made to the O2 reduction reaction in the family of heme‑copper oxidases.

Comment on “High‐Efficient Blue Emission and Bandgap Engineering from Jahn‒Teller Distorted Halide Double Perovskites”

AbstractIn a recent publication, Liu et al. investigated the doping strategy of pentavalent antimony (Sb5+) in the benchmark material Cs2NaInCl6, achieving a bright blue emission with a high photoluminescence quantum yield of ≈90% at low doping concentrations and bandgaps as low as 1.24 eV at high doping concentrations (Adv. Opt. Mater. 2024, 12, 2301576). In this comment, by reanalyzing Liu et al.' data and conducting hybrid density functional theory calculations, the study provides evidence supporting the incorporation of Sb3+ ions rather than “Sb5+ ions”, as the cause of the efficient blue emission observed in the lightly‐doped Cs2NaInCl6 samples. Furthermore, it is proposed that the narrowing of the bandgap in the “heavily‐doped Cs2NaInCl6” samples is a result of the formation of a new phase related to a mixed‐valence salt, Cs4SbIIISbVCl12.

Computational Prediction of an Antimony-Based n-Type Transparent Conducting Oxide: F-Doped Sb2O5.

Transparent conducting oxides (TCOs) possess a unique combination of optical transparency and electrical conductivity, making them indispensable in optoelectronic applications. However, their heavy dependence on a small number of established materials limits the range of devices that they can support. The discovery and development of additional wide bandgap oxides that can be doped to exhibit metallic-like conductivity are therefore necessary. In this work, we use hybrid density functional theory to identify a binary Sb(V) system, Sb2O5, as a promising TCO with high conductivity and transparency when doped with fluorine. We conducted a full point defect analysis, finding F-doped Sb2O5 to exhibit degenerate n-type transparent conducting behavior. The inherently large electron affinity found in antimony oxides also widens their application in organic solar cells. Following our previous work on zinc antimonate, this work provides additional support for designing Sb(V)-based oxides as cost-effective TCOs for a broader range of applications.

Basis Set Requirements of σ-Functionals for Gaussian- and Slater-Type Basis Functions and Comparison with Range-Separated Hybrid and Double Hybrid Functionals.

σ-Functionals belong to the class of Kohn-Sham (KS) correlation functionals based on the adiabatic-connection fluctuation-dissipation theorem and are technically closely related to the random phase approximation (RPA). They have the same computational demand as the latter, with the computational effort of an energy evaluation for both methods being lower than that of a preceding hybrid DFT calculation for typical systems but yield much higher accuracy, reaching chemical accuracy of 1 kcal/mol for quantities such as reactions and transition energies in main group chemistry. In previous work on σ-functionals, rather large Gaussian basis sets have been used. Here, we investigate the actual basis set requirements of σ-functionals and present three setups that employ smaller Gaussian basis sets ranging from quadruple-ζ (QZ) to triple-ζ (TZ) quality and represent a good compromise between accuracy and computational efficiency. Furthermore, we introduce an implementation of σ-functionals based on Slater-type basis sets and present two setups of QZ and TZ quality for this implementation. We test the accuracy of these setups on a large database of various physical properties and types of reactions, as well as equilibrium geometries and vibrational frequencies. As expected, the accuracy of σ-functional calculations becomes somewhat lower with a decreasing basis set size. However, for all setups considered here, calculations with σ-functionals are clearly more accurate than those within the RPA and even more so than those of the conventional KS methods. For the smallest setup using Gaussian-type basis functions and Slater-type basis functions, we introduce a reparametrization that reduces the loss in accuracy due to the basis set error to some extent. A comparison with the range-separated hybrid ωB97X-V and the double hybrid DSD-BLYP-D3 shows that σ functionals outperform in accuracy both of these accurate and, for their class, representative functionals.

The performance and efficiency of twelve range-separated hybrid DFT functionals for calculation of the magnetic exchange coupling constants of di-nuclear first row transition metal complexes

The performance and efficiency of twelve range-separated hybrid functionals are compared with each other in calculating magnetic exchange coupling constant (J) data on eleven di-nuclear first-row transition metal complexes having di-copper (Cu) and di-vanadium (V) as their central atoms. The J values are then compared with the experimental findings. In addition to that we also computed the B3LYP results of these complexes. We compared them with the results obtained from range-separated functionals as well as with the experimental results to understand their performances in predicting the magnetic exchange coupling constant. When computing, we reoptimized all the structures of the complexes taken from the crystal data, and then the values of Js were determined. The performances are predicted in terms of four statistical error matrices, namely, the mean absolute error (MAE), mean fractional error (MFE), mean signed error (MSE), root mean square error (RMSE) in cm−1.

Thermoelectric properties of the bismuth oxychalcogenides Bi2SO2, Bi2SeO2 and Bi2TeO2

We present a detailed theoretical study of the thermoelectric properties of the bismuth oxychalcogenides Bi2ChO2 (Ch = S, Se, Te). The electrical transport is modelled using semi-classical Boltzmann transport theory with electronic structures from hybrid density-functional theory, including an approximate model for the electron lifetimes. The lattice thermal conductivity is calculated using first-principles phonon calculations with an explicit treatment of anharmonicity, yielding microscopic insight into how partial replacement of the chalcogen in the bismuth chalcogenides impacts the phonon transport. We find very good agreement between the predicted transport properties and a favourable cancellation of errors that allows for near-quantitative predictions of the thermoelectric figure of merit ZT. Our calculations suggest recent experiments on n-doped Bi2SeO2 have achieved close to the largest ZT possible in bulk materials, whereas the largest reported ZT for Bi2TeO2 could be improved sixfold by optimising the carrier concentration. We also predict that much larger ZT > 2.5, competitive with the benchmark thermoelectric SnSe, could be obtained for Bi2SO2 and Bi2SeO2 with heavy p-type doping. This study demonstrates the predictive power of this modelling approach for studying thermoelectrics and highlights several avenues for improving the performance of the Bi2ChO2.