Theoretical Investigation of the Lattice Thermal Conductivities of II-IV-V2 Pnictide Semiconductors.

Ternary pnictide semiconductors with II-IV-V2 stoichiometry hold potential as cost-effective thermoelectric materials with suitable electronic transport properties, but their lattice thermal conductivities (κ) are typically too high. Insights into their vibrational properties are therefore crucial to finding strategies to reduce κ and achieve improved thermoelectric performance. We present a theoretical exploration of the lattice thermal conductivities for a set of pnictide semiconductors with ABX2 composition (A = Zn, Cd; B = Si, Ge, Sn; and X = P, As) using machine-learning-based regression algorithms to extract force constants from a reduced number of density functional theory simulations and then solving the Boltzmann transport equation for phonons. Our results align well with available experimental data, decreasing the mean absolute error by ∼3 W m-1 K-1 with respect to the best previous set of theoretical predictions. Zn-based ternary pnictides have, on average, more than double the thermal conductivity of the Cd-based compounds. Anisotropic behavior increases with the mass difference between A and B cations, but while the nature of the anion does not affect the structural anisotropy, the thermal conductivity anisotropy is typically higher for arsenides than for phosphides. We identify compounds such as CdGeAs2, for which nanostructuring to an affordable range of particle sizes could lead to κ values low enough for thermoelectric applications.

Design Principles Guided by DFT Calculations and High-Throughput Frameworks for the Discovery of New Diamond-like Chalcogenide Thermoelectric Materials.

Rational design principles are one pathway to discovering new materials. However, technological breakthroughs rarely occur in this way because these design principles are usually based on incremental advances that seldom lead to disruptive applications. The emergence of machine-learning (ML) and high-throughput (HT) techniques has changed the paradigm, opening up new possibilities for efficiently screening large chemical spaces and creating on-the-fly design principles for the discovery of novel materials with desired properties. In this work, the approach is used to discover novel thermoelectric (TE) materials based on quaternary diamond-like chalcogenides. A HT framework that integrates density functional theory calculations, ML, and the solution of the Boltzmann transport equation is used to efficiently rationalize the transport properties of these compounds and identify those with potential as TE materials, achieving ZT values above 2.

High-Throughput Prediction of the Thermal and Electronic Transport Properties of Large Physical and Chemical Spaces Accelerated by Machine Learning: Charting the ZT of Binary Skutterudites.

Thermal and electronic transport properties are the keys to many technological applications of materials. Thermoelectric, TE, materials can be considered a singular case in which not only one but three different transport properties are combined to describe their performance through their TE figure of merit, ZT. Despite the availability of high-throughput experimental techniques, synthesizing, characterizing, and measuring the properties of samples with numerous variables affecting ZT are not a cost- or time-efficient approach to lead this strategy. The significance of computational materials science in discovering new TE materials has been running in parallel to the development of new frameworks and methodologies to compute the electron and thermal transport properties linked to ZT. Nevertheless, the trade-off between computational cost and accuracy has hindered the reliable prediction of TE performance for large chemical spaces. In this work, we present for the first time the combination of new ab initio methodologies to predict transport properties with machine learning and a high-throughput framework to establish a solid foundation for the accurate prediction of thermal and electron transport properties. This strategy is applied to a whole family of materials, binary skutterudites, which are well-known as good TE candidates. Following this methodology, it is possible not only to connect ZT with the experimental synthetic (carrier concentration and grain size) and operando (temperature) variables but also to understand the physical and chemical phenomena that act as driving forces in the maximization of ZT for p-type and n-type binary skutterudites.

High-Throughput Screening of the Thermoelastic Properties of Ultrahigh-Temperature Ceramics.

Ultrahigh-temperature ceramics (UHTCs) are a group of materials with high technological interest because of their applications in extreme environments. However, their characterization at high temperatures represents the main obstacle for their fast development. Obstacles are found from an experimental point of view, where only few laboratories around the world have the resources to test these materials under extreme conditions, and also from a theoretical point of view, where actual methods are expensive and difficult to apply to large sets of materials. Here, a new theoretical high-throughput framework for the prediction of the thermoelastic properties of materials is introduced. This approach can be systematically applied to any kind of crystalline material, drastically reducing the computational cost of previous methodologies up to 80% approximately. This new approach combines Taylor expansion and density functional theory calculations to predict the vibrational free energy of any arbitrary strained configuration, which represents the bottleneck in other methods. Using this framework, elastic constants for UHTCs have been calculated in a wide range of temperatures with excellent agreement with experimental values, when available. Using the elastic constants as the starting point, other mechanical properties such a bulk modulus, shear modulus, or Poisson ratio have been also explored, including upper and lower limits for polycrystalline materials. Finally, this work goes beyond the isotropic mechanical properties and represents one of the most comprehensive and exhaustive studies of some of the most important UHTCs, charting their anisotropy and thermal and thermodynamical properties.

Comprehensive Experimental and Theoretical Study of the CO + NO Reaction Catalyzed by Au/Ni Nanoparticles.

The catalytic and structural properties of five different nanoparticle catalysts with varying Au/Ni composition were studied by six different methods, including in situ X-ray absorption spectroscopy and density functional theory (DFT) calculations. The as-prepared materials contained substantial amounts of residual capping agent arising from the commonly used synthetic procedure. Thorough removal of this material by oxidation was essential for the acquisition of valid catalytic data. All catalysts were highly selective toward N2 formation, with 50-50 Au:Ni material being best of all. In situ X-ray absorption near edge structure spectroscopy showed that although Au acted to moderate the oxidation state of Ni, there was no clear correlation between catalytic activity and nickel oxidation state. However, in situ extended X-ray absorption fine structure spectroscopy showed a good correlation between Au-Ni coordination number (highest for Ni50Au50) and catalytic activity. Importantly, these measurements also demonstrated substantial and reversible Au/Ni intermixing as a function of temperature between 550 °C (reaction temperature) and 150 °C, underlining the importance of in situ methods to the correct interpretation of reaction data. DFT calculations on smooth, stepped, monometallic and bimetallic surfaces showed that N + N recombination rather than NO dissociation was always rate-determining and that the activation barrier to recombination reaction decreased with increased Au content, thus accounting for the experimental observations. Across the entire composition range, the oxidation state of Ni did not correlate with activity, in disagreement with earlier work, and theory showed that NiO itself should be catalytically inert. Au-Ni interactions were of paramount importance in promoting N + N recombination, the rate-limiting step.

Improving the activity of gold nanoparticles for the water-gas shift reaction using TiO2-Y2O3: an example of catalyst design.

In the last ten years, there has been an acceleration in the pace at which new catalysts for the water-gas shift reaction are designed and synthesized. Pt-based catalysts remain the best solution when only activity is considered. However, cost, operation temperature, and deactivation phenomena are important variables when these catalysts are scaled in industry. Here, a new catalyst, Au/TiO2-Y2O3, is presented as an alternative to the less selective Pt/oxide systems. Experimental and theoretical techniques are combined to design, synthesize, characterize and analyze the performance of this system. The mixed oxide demonstrates a synergistic effect, improving the activity of the catalyst not only at large-to-medium temperatures but also at low temperatures. This effect is related to the homogeneous dispersion of the vacancies that act both as nucleation centers for smaller and more active gold nanoparticles and as dissociation sites for water molecules. The calculated reaction path points to carboxyl formation as the rate-limiting step with an activation energy of 6.9 kcal mol-1, which is in quantitative agreement with experimental measurements and, to the best of our knowledge, it is the lowest activation energy reported for the water-gas shift reaction. This discovery demonstrates the importance of combining experimental and theoretical techniques to model and understand catalytic processes and opens the door to new improvements to reduce the operating temperature and the deactivation of the catalyst.

Adsorption of Prototypical Asphaltenes on Silica: First-Principles DFT Simulations Including Dispersion Corrections.

In this work, we explore the interaction between some prototypical asphaltene and porphyrin molecules with a fully hydroxylated (0001) surface of α-quartz by means of theoretical calculations based on the density functional theory (DFT) under periodic boundary conditions. The influence of dispersion forces, adsorption geometries, and size of the side chain is analyzed. The inclusion of London dispersion forces is overriding as they increase the interaction by about 1 order of magnitude. All of the considered molecules strongly interact with the hydroxylated surface and prefer to adsorb in a parallel position instead of vertically. It is also found that adsorption energy always increases with larger side chains because dispersion interactions also augment. Interestingly, in the case of porphyrin, the less stable isomer in the gas phase is the preferred one after adsorption, which is substantiated by a differential stabilization induced by the surface. Finally, we present a comparative study of the adsorption of these model molecules in terms of energy per area unit and energy per interacting π electron.

Ag2S Quantum Dot-Sensitized Solar Cells by First Principles: The Effect of Capping Ligands and Linkers.

Quantum dots solar cells, QDSCs, are one of the candidates for being a reliable alternative to fossil fuels. However, the well-studied CdSe and CdTe-based QDSCs present a variety of issues for their use in consumer-goods applications. Silver sulfide, Ag2S, is a promising material, but poor efficiency has been reported for QDSCs based on this compound. The potential influence of each component of QDSCs is critical and key for the development of more efficient devices based on Ag2S. In this work, density functional theory calculations were performed to study the nature of the optoelectronic properties for an anatase-TiO2(101) surface sensitized with different silver sulfide nanoclusters. We demonstrated how it is possible to deeply tune of its electronic properties by modifying the capping ligands and linkers to the surface. Finally, an analysis of the electron injection mechanism for this system is presented.

Effects of the capping ligands, linkers and oxide surface on the electron injection mechanism of copper sulfide quantum dot-sensitized solar cells.

Quantum dot-sensitized solar cells, QDSCs, are a clean and effective alternative to fossil fuels to reduce CO2 emissions. However, the different components that constitute the QDSCs and the difficulty of isolating experimentally their effects on the performance of the whole system slow down the development of more efficient devices. In this work, DFT calculations are combined with a bottom-up approach to differentiate the effect of each component on the electronic structure and absorption spectra. First, Cu2S QDs were built including a U parameter to effectively describe the localization of electrons. The effect of capping agents is addressed using ligands with different electron-donating/withdrawing groups. The role of linkers and their adsorption on the oxide surface are also examined. Finally, we propose a main indirect electron injection mechanism based on the position of the peaks of the spectra.

Host-guest interactions between cyclodextrins and surfactants with functional groups at the end of the hydrophobic tail.

The aim of this work was to investigate the influence of the incorporation of substituents at the end of the hydrophobic tail on the binding of cationic surfactants to α-, ß-, and γ-cyclodextrins. The equilibrium binding constants of the 1:1 inclusion complexes formed follow the trend K1(α-CD)>K1(ß-CD)≫K1(γ-CD), which can be explained by considering the influence of the CD cavity volume on the host-guest interactions. From the comparison of the K1 values obtained for dodecyltriethylammonium bromide, DTEAB, to those estimated for the surfactants with the substituents, it was found that the incorporation of a phenoxy group at the end of the hydrocarbon tail does not affect K1, and the inclusion of a naphthoxy group has some influence on the association process, slightly diminishing K1. This makes evident the importance of the contribution of hydrophobic interactions to the binding, the length of the hydrophobic chain being the key factor determining K1. However, the presence of the aromatic rings does influence the location of the host and the guest in the inclusion complexes. The observed NOE interactions between the aromatic protons and the CD protons indicate that the aromatic rings are partially inserted within the host cavity, with the cyclodextrin remaining close to the aromatic rings, which could be partially intercalated in the host cavity. To the authors' knowledge this is the first study on the association of cyclodextrins with monomeric surfactants incorporating substituents at the end of the hydrophobic tail.

Communication: improving the density functional theory+U description of CeO2 by including the contribution of the O 2p electrons.

Density functional theory (DFT) based approaches within the local-density approximation or generalized gradient approximation frameworks fail to predict the correct electron localization in strongly correlated systems due to the lack of cancellation of the Coulomb self-interaction. This problem might be circumvented either by using hybrid functionals or by introducing a Hubbard-like term to account for the on site interactions. This latter DFT+U approach is less expensive and therefore more practical for extensive calculations in solid-state computational simulations. By and large, the U term only affects the metal electrons, in our case the Ce 4f ones. In the present work, we report a systematic analysis of the effect of adding such a U term also to the oxygen 2p electrons. We find that using a set of U(f) = 5 eV and U(p) = 5eV effective terms leads to improved description of the lattice parameters, band gaps, and formation and reduction energies of CeO(2).

Comparative Study on the Performance of Hybrid DFT Functionals in Highly Correlated Oxides: The Case of CeO2 and Ce2O3.

The outstanding catalytic properties of cerium oxides rely on the easy Ce(3+) ↔ Ce(4+) redox conversion, which however constitutes a challenge in density functional based theoretical chemistry due to the strongly correlated nature of the 4f electrons present in the reduced materials. In this work, we report an analysis of the performance of five exchange-correlation functionals (HH, HHLYP, PBE0, B3LYP, and B1-WC) implemented in the CRYSTAL06 code to describe three properties of ceria: crystal structure, band gaps, and reaction energies of the CeO2 → Ce2O3 process. All five functionals give values for cell parameters that are in fairly good agreement with experiment, although the PBE0 hybrid functional is found to be the most accurate. Band gaps, 2p-4f-5d in the case of CeO2 and 4f-5d in the case of Ce2O3, are found to be, in general, overestimated and drop off when the amount of Hartree-Fock exchange in the exchange-correlation functional decreases. In contrast, the reaction energies are found to be underestimated, and increase when the amount of HF exchange lowers. Overall, at its standard formulation, the B1-WC functional seems to be the best choice as it provides good band gaps and reaction energies, and very reasonable crystal parameters.

Synthesis, solid-state structure, and bonding analysis of the beryllocenes [Be(C5Me4H)2], [Be(C5Me5)2], and [Be(C5Me5)(C5Me4H)].

The beryllocenes [Be(C(5)Me(4)H)(2)] (1), [Be(C(5)Me(5))(2)] (2), and [Be(C(5)Me(5))(C(5)Me(4)H)] (3) have been prepared from BeCl(2) and the appropriate KCp' reagent in toluene/diethyl ether solvent mixtures. The synthesis of 1 is facile (20 degrees C, overnight), but generation of decamethylberyllocene 2 demands high temperatures (ca. 115 degrees C) and extended reaction times (3-4 days). The mixed-ring beryllocene 3 is obtained when the known [(eta(5)-C(5)Me(5))BeCl] is allowed to react with K[C(5)Me(4)H], once more under somewhat forcing conditions (115 degrees C, 36 h). The structures of the three metallocenes have been determined by low-temperature X-ray studies. Both 1 and 3 present eta5/eta1 geometries of the slip-sandwich type, whereas 2 exhibits an almost regular, ferrocene-like, sandwich structure. In the mixed-ring compound 3, C(5)Me(5) is centrally bound to beryllium and the eta(1)-C(5)Me(4)H ring bonds to the metal through the unique CH carbon atom. This is also the binding mode of the eta(1)-ring of 1. To analyze the nature of the bonding in these molecules, theoretical calculations at different levels of theory have been performed on compounds 2 and 3, and a comparison with the bonding in [Be(C(5)H(5))(2)] has been made. As for the latter molecule, energy differences between the eta5/eta5 and the eta5/eta1 structures of 2 are very small, being of the order of a few kcal mol(-1). Constrained space orbital variations (CSOV) calculations show that the covalent character in the bonding is larger for [Be(C(5)Me(5))(2)] than for [Be(C(5)H(5))(2)] due to larger charge delocalization and to increased polarizability of the C(5)Me(5) fragment.