cometadata/datacite-titles-descriptions-related-identifiers huggingface dataset:
DOI: 10.24435/materialscloud:rm-cx
title: Temperature dependence of the electronic structure of semiconductors and insulators

description: The renormalization of electronic eigenenergies due to electron-phonon coupling (temperature dependence and zero-point motion effect) is sizable in many materials with light atoms. This effect, often neglected in ab initio calculations, can be computed using the perturbation-based Allen-Heine-Cardona theory in the adiabatic or non-adiabatic harmonic approximation. After a short description of the recent progresses in this field and a brief overview of the theory, we focus on the issue of phonon wavevector sampling convergence, until now poorly understood. Indeed, the renormalization is obtained numerically through a slowly converging q-point integration. For non-zero Born effective charges, we show that a divergence appears in the electron-phonon matrix elements at q → Γ, leading to a divergence of the adiabatic renormalization at band extrema. This problem is exacerbated by the slow convergence of Born effective charges with electronic wavevector sampling, which leaves residual Born effective charges in ab initio calculations on materials that are physically devoid of such charges. Here, we propose a solution that improves this convergence. However, for materials where Born effective charges are physically non-zero, the divergence of the renormalization indicates a breakdown of the adiabatic harmonic approximation, which we assess here by switching to the non-adiabatic harmonic approximation. Also, we study the convergence behavior of the renormalization and develop reliable extrapolation schemes to obtain the converged results. Finally, the adiabatic and non-adiabatic theories, with corrections for the slow Born effective charge convergence problem (and the associated divergence) are applied to the study of five semiconductors and insulators: α-AlN, β-AlN, BN, diamond, and silicon. For these five materials, we present the zero-point  renormalization, temperature dependence, phonon-induced lifetime broadening, and the renormalized electronic band structure.
DOI: 10.24435/materialscloud:9v-pt
title: Hidden spontaneous polarisation in the chalcohalide photovoltaic Sn2SbS2I3

description: Enormous research efforts are currently devoted to the discovery of 'perovskite-inspired materials', aiming to replicate the astonishing optoelectronic performance of lead-halide perovskites (LHPs). Recently, chalco halides of group IV/V elements have attracted attention due to the stability provided by stronger metal-chalcogen bonds, alongside compositional flexibility and ns2 cations — a performance-defining feature of LHPs. Following the experimental report of stable, solution-grown tin-antimony sulfoiodide (Sn2SbS2I3) solar cells, with power conversion efficiencies above 4%, we comprehensively characterise the structural and electronic properties of this emerging material. We find that the experimentally-reported centrosymmetric Cmcm crystal structure represents an average over multiple polar Cmc2_1 configurations. This dynamic crystal structure and ferroelectric behaviour could benefit photovoltaic performance. Using state-of-the-art ab initio methods, we assess the efficiency limits of this material, finding maximal solar-conversion efficiencies η_max > 30 % with film thicknesses t > 0.5μm, at the radiative limit.

Open-access Materials Horizons paper: https://doi.org/10.1039/D1MH00764E
Talks on this and other works at: https://www.youtube.com/channel/UCoVGnBeZeWmKzv8_-PzCKCw
DOI: 10.24435/materialscloud:az-s2
title: Yttrium tantalum oxynitride multiphases as photoanodes for water oxidation

description: Perovskite yttrium tantalum oxynitride is theoretically proposed as a promising semiconductor for solar water splitting because of the predicted bandgap and energy positions of band edges. In experiment, however, we show here that depending on processing parameters, yttrium tantalum oxynitrides exist in multiphases, including the desired perovskite YTaON2, defect fluorite YTa(O,N,o)4, and N-doped YTaO4. These multiphases have bandgaps ranging between 2.13 and 2.31 eV, all responsive to visible light. The N-doped YTaO4, perovskite main phase, and fluorite main phase derived from crystalline fergusonite oxide precursors exhibit interesting photoelectrochemical performances for water oxidation, while the defect fluorite derived from low crystallized scheelite-type oxide precursors show negligible activity. Preliminarily measurements show that loading IrOx cocatalyst on N-doped YTaO4 significantly improves its photoelectrochemical performance encouraging further studies to optimize this new material for solar fuel production.
DOI: 10.24435/materialscloud:zm-09
title: Yttrium tantalum oxynitride multiphases as photoanodes for water oxidation

description: Perovskite yttrium tantalum oxynitride is theoretically proposed as a promising semiconductor for solar water splitting because of the predicted bandgap and energy positions of band edges. In experiment, however, we show here that depending on processing parameters, yttrium tantalum oxynitrides exist in multiphases, including the desired perovskite YTaON2, defect fluorite YTa(O,N,o)4, and N-doped YTaO4. These multiphases have bandgaps ranging between 2.13 and 2.31 eV, all responsive to visible light. The N-doped YTaO4, perovskite main phase, and fluorite main phase derived from crystalline fergusonite oxide precursors exhibit interesting photoelectrochemical performances for water oxidation, while the defect fluorite derived from low crystallized scheelite-type oxide precursors show negligible activity. Preliminarily measurements show that loading IrOx cocatalyst on N-doped YTaO4 significantly improves its photoelectrochemical performance encouraging further studies to optimize this new material for solar fuel production.
DOI: 10.24435/materialscloud:zw-be
title: Predicting polymerization reactions via transfer learning using chemical language models

description: Polymers are candidate materials for a wide range of sustainability applications such as carbon capture and energy storage. However, computational polymer discovery lacks automated analysis of reaction pathways and stability assessment through retro-synthesis. Here, we report the first extension of transformer-based language models to polymerization reactions for both forward and retrosynthesis tasks. We curated a polymerization dataset for vinyl polymers covering reactions and retrosynthesis for representative homo-polymers and co-polymers. Overall, we report a forward model accuracy of 80% and a backward model accuracy of 60%. We further analyse the model performance on a set of case studies by providing polymerization and retro-synthesis examples and evaluating the model's predictions quality from a materials science perspective.
DOI: 10.24435/materialscloud:2019.0028/v1
title: Intrinsic defects in amorphous TiO2

description: This entry provides the atomic structures of three bulk amorphous TiO2 models generated through the melt-and-quench method with different cooling rates and of ten O-O peroxy linkages obtained by adding two holes to the bulk model constructed with the lowest cooling rate.
DOI: 10.24435/materialscloud:fn-r7
title: Predicting polymerization reactions via transfer learning using chemical language models

description: Polymers are candidate materials for a wide range of sustainability applications such as carbon capture and energy storage. However, computational polymer discovery lacks automated analysis of reaction pathways and stability assessment through retro-synthesis. Here, we report the first extension of transformer-based language models to polymerization reactions for both forward and retrosynthesis tasks. We curated a polymerization dataset for vinyl polymers covering reactions and retrosynthesis for representative homo-polymers and co-polymers. Overall, we report a forward model accuracy of 80% and a backward model accuracy of 60%. We further analyse the model performance on a set of case studies by providing polymerization and retro-synthesis examples and evaluating the model's predictions quality from a materials science perspective.
DOI: 10.24435/materialscloud:6d-a0
title: Intrinsic defects in amorphous TiO2

description: This entry provides the atomic structures of three bulk amorphous TiO2 models generated through the melt-and-quench method with different cooling rates and of ten O-O peroxy linkages obtained by adding two holes to the bulk model constructed with the lowest cooling rate.
DOI: 10.24435/materialscloud:2019.0053/v1
title: On-surface synthesis of antiaromatic and open-shell indeno[2,1-b]fluorene polymers and their lateral fusion into porous ribbons

description: This record contains the experimental and computational data needed to support the work done on on-surface synthesis of conjugated polymers consisting of indeno[2,1-b]fluorene units, which are antiaromatic and open-shell biradicaloids. The observed reaction products have been characterized via low-temperature scanning tunneling microscopy/spectroscopy and noncontact atomic force microscopy, complemented by density-functional theory calculations. These polymers present a low band gap when adsorbed on Au(111). Moreover, their pronounced antiaromaticity and radical character, elucidated by ab initio calculations, make them promising candidates for applications in electronics and spintronics. Further, they provide a rich playground to explore magnetism in low-dimensional organic nanomaterials.
DOI: 10.24435/materialscloud:2019.0032/v1
title: Free-Energy Surface Prediction by Flying Gaussian Method

description: Molecular simulations are computationally expensive, especially in systems with multiple free energy minima. To address this problem many enhanced sampling techniques have been developed. Metadynamics uses a bias potential defined as a sum of Gaussian hills in space of few (one or two) collective variables. This bias potential disfavors states that have been visited since the beginning of the simulation. Multiple walker metadynamics simulates the system in multiple parallel replicas. The bias potential disfavors states that have been visited since the beginning of the simulation in any replica. Flying Gaussian method presented here also simulates the system in multiple parallel replicas. It disfavors states that are, at certain moment, similar in two or more replicas. It was demonstrated on Alanine Dipeptide in vacuum and water, cis/trans-isomerisation of Proline-containing peptides and Met-enkephalin.
DOI: 10.24435/materialscloud:fb-3g
title: On-surface synthesis of antiaromatic and open-shell indeno[2,1-b]fluorene polymers and their lateral fusion into porous ribbons

description: This record contains the experimental and computational data needed to support the work done on on-surface synthesis of conjugated polymers consisting of indeno[2,1-b]fluorene units, which are antiaromatic and open-shell biradicaloids. The observed reaction products have been characterized via low-temperature scanning tunneling microscopy/spectroscopy and noncontact atomic force microscopy, complemented by density-functional theory calculations. These polymers present a low band gap when adsorbed on Au(111). Moreover, their pronounced antiaromaticity and radical character, elucidated by ab initio calculations, make them promising candidates for applications in electronics and spintronics. Further, they provide a rich playground to explore magnetism in low-dimensional organic nanomaterials.
DOI: 10.24435/materialscloud:2019.0035/v1
title: Electron and Hole Polarons at the BiVO4–Water Interface

description: We determine the transition levels of electron and hole polarons at the BiVO4–water interface through thermodynamic integration within a hybrid functional scheme, thereby accounting for the liquid nature of the water component. The electron polaron is found to be less stable at the interface than in the bulk by 0.18 eV, while for the hole polaron the binding energy increases by 0.20 eV when the charge localizes in the surface layer of BiVO4. These results indicate that interfacial effects on the polaron binding energy and charge distribution are sizeable and cannot trivially be inferred from bulk calculations.
DOI: 10.24435/materialscloud:p7-js
title: Free-Energy Surface Prediction by Flying Gaussian Method

description: Molecular simulations are computationally expensive, especially in systems with multiple free energy minima. To address this problem many enhanced sampling techniques have been developed. Metadynamics uses a bias potential defined as a sum of Gaussian hills in space of few (one or two) collective variables. This bias potential disfavors states that have been visited since the beginning of the simulation. Multiple walker metadynamics simulates the system in multiple parallel replicas. The bias potential disfavors states that have been visited since the beginning of the simulation in any replica. Flying Gaussian method presented here also simulates the system in multiple parallel replicas. It disfavors states that are, at certain moment, similar in two or more replicas. It was demonstrated on Alanine Dipeptide in vacuum and water, cis/trans-isomerisation of Proline-containing peptides and Met-enkephalin.
DOI: 10.24435/materialscloud:rr-9b
title: Electron and Hole Polarons at the BiVO4–Water Interface

description: We determine the transition levels of electron and hole polarons at the BiVO4–water interface through thermodynamic integration within a hybrid functional scheme, thereby accounting for the liquid nature of the water component. The electron polaron is found to be less stable at the interface than in the bulk by 0.18 eV, while for the hole polaron the binding energy increases by 0.20 eV when the charge localizes in the surface layer of BiVO4. These results indicate that interfacial effects on the polaron binding energy and charge distribution are sizeable and cannot trivially be inferred from bulk calculations.
DOI: 10.24435/materialscloud:2020.0032/v1
title: Ni Nanoparticles on CeO2(111): Energetics, Electron Transfer and Structure by Ni Adsorption Calorimetry, Spectroscopies and DFT

description: The morphology, interfacial bonding energetics and charge transfer of Ni clusters and nanoparticles on slightly-reduced CeO2-x(111) surfaces at 100 to 300 K have been studied using single crystal adsorption calorimetry (SCAC), low-energy ion scattering spectroscopy (LEIS), X-ray photoelectron spectroscopy (XPS), low energy electron diffraction (LEED) and density functional theory (DFT).  The initial heat of adsorption of Ni vapor decreased with the extent of pre-reduction (x) of the CeO2-x(111), showing that stoichiometric ceria adsorbs Ni more strongly than oxygen vacancies.  On CeO1.95(111) at 300 K, the heat dropped quickly with coverage in the first 0.1 ML, attributed to nucleation of Ni clusters on stoichiometric steps, followed by the Ni particles spreading onto less favorable terrace sites. At 100 K, the clusters nucleate on terraces due to slower diffusion. Adsorbed Ni monomers are in the +2 oxidation state, and they bind by ~45 kJ/mol more strongly to step sites than terraces. The measured heat of adsorption versus average particle size on terraces is favorably compared to DFT calculations.  The Ce 3d XPS lineshape showed an increase in Ce3+/Ce4+ ratio with Ni coverage, providing the number of electrons donated to the ceria per Ni atom. The charge transferred per Ni is initially large but strongly decreases with increasing cluster size for both experiments and DFT, and shows large differences between clusters at steps versus terraces. This charge is localized on the interfacial Ni and Ce atoms in their atomic layers closest to the interface. This knowledge is crucial to understanding the nature of the active sites on the surface of Ni-CeO2 catalysts for which metal-oxide interactions play a very important role in the activation of O-H and C-H bonds.  The changes in these interactions with Ni particle size (metal loading) and the extent of reduction of the ceria help to explain how previously reported catalytic activity and selectivity change with these same structural details.
DOI: 10.24435/materialscloud:2020.0038/v1
title: Oxygen evolution reaction: Bifunctional mechanism breaking the linear scaling relationship

description: The bifunctional mechanism for the oxygen evolution reaction (OER) involving two distinct reaction sites is studied through the computational hydrogen electrode method for a set of catalyst materials including rutile TiO2(110), anatase TiO2(101), SnO2(110), RuO2(110), IrO2(110), Ni2P(0001), and BiVO4(001). The calculations are performed both at the semilocal level and at the hybrid functional level. Moreover, anodic conditions are modeled and their effect on the OER free energy steps is evaluated. The free energies of the reaction steps indicate that for specific combinations of catalysts, the limitations due to the linear scaling relationship can be overcome, leading to smaller overpotentials for the overall OER. At the same time, a detailed analysis of the results reveals a strong dependence on the adopted functional. For both functionals, it is shown that the energy level of the highest occupied electronic state can serve as a descriptor to guide the search for the optimal catalyst acting as a hydrogen acceptor. These results support the bifunctional mechanism as a means to break the linear scaling relationship and to further reduce the overpotential of the OER.
DOI: 10.24435/materialscloud:b7-cn
title: Ni Nanoparticles on CeO2(111): Energetics, Electron Transfer and Structure by Ni Adsorption Calorimetry, Spectroscopies and DFT

description: The morphology, interfacial bonding energetics and charge transfer of Ni clusters and nanoparticles on slightly-reduced CeO2-x(111) surfaces at 100 to 300 K have been studied using single crystal adsorption calorimetry (SCAC), low-energy ion scattering spectroscopy (LEIS), X-ray photoelectron spectroscopy (XPS), low energy electron diffraction (LEED) and density functional theory (DFT).  The initial heat of adsorption of Ni vapor decreased with the extent of pre-reduction (x) of the CeO2-x(111), showing that stoichiometric ceria adsorbs Ni more strongly than oxygen vacancies.  On CeO1.95(111) at 300 K, the heat dropped quickly with coverage in the first 0.1 ML, attributed to nucleation of Ni clusters on stoichiometric steps, followed by the Ni particles spreading onto less favorable terrace sites. At 100 K, the clusters nucleate on terraces due to slower diffusion. Adsorbed Ni monomers are in the +2 oxidation state, and they bind by ~45 kJ/mol more strongly to step sites than terraces. The measured heat of adsorption versus average particle size on terraces is favorably compared to DFT calculations.  The Ce 3d XPS lineshape showed an increase in Ce3+/Ce4+ ratio with Ni coverage, providing the number of electrons donated to the ceria per Ni atom. The charge transferred per Ni is initially large but strongly decreases with increasing cluster size for both experiments and DFT, and shows large differences between clusters at steps versus terraces. This charge is localized on the interfacial Ni and Ce atoms in their atomic layers closest to the interface. This knowledge is crucial to understanding the nature of the active sites on the surface of Ni-CeO2 catalysts for which metal-oxide interactions play a very important role in the activation of O-H and C-H bonds.  The changes in these interactions with Ni particle size (metal loading) and the extent of reduction of the ceria help to explain how previously reported catalytic activity and selectivity change with these same structural details.
DOI: 10.24435/materialscloud:n2-b8
title: Oxygen evolution reaction: Bifunctional mechanism breaking the linear scaling relationship

description: The bifunctional mechanism for the oxygen evolution reaction (OER) involving two distinct reaction sites is studied through the computational hydrogen electrode method for a set of catalyst materials including rutile TiO2(110), anatase TiO2(101), SnO2(110), RuO2(110), IrO2(110), Ni2P(0001), and BiVO4(001). The calculations are performed both at the semilocal level and at the hybrid functional level. Moreover, anodic conditions are modeled and their effect on the OER free energy steps is evaluated. The free energies of the reaction steps indicate that for specific combinations of catalysts, the limitations due to the linear scaling relationship can be overcome, leading to smaller overpotentials for the overall OER. At the same time, a detailed analysis of the results reveals a strong dependence on the adopted functional. For both functionals, it is shown that the energy level of the highest occupied electronic state can serve as a descriptor to guide the search for the optimal catalyst acting as a hydrogen acceptor. These results support the bifunctional mechanism as a means to break the linear scaling relationship and to further reduce the overpotential of the OER.
DOI: 10.24435/materialscloud:2020.0013/v1
title: Molecular Dynamics Simulations of Crystal Nucleation from Solution at Constant Chemical Potential

description: A widespread method of crystal preparation is to precipitate it from a supersaturated solution. In such a process, control of solution concentration is of paramount importance. The nucleation process, polymorph selection, and crystal habits depend crucially on this thermodynamic parameter. When performing molecular dynamics simulations with a fixed number of molecules in the canonical ensemble, crystal growth is accompanied by a decrease in the solution concentration. This modification of the thermodynamic condition leads to significant artifacts. Inspired by the recent development of the constant chemical potential molecular dynamics simulation method by Perego et al. [J. Chem. Phys. 2015, 142, 144113], we develop a spherical variant of it to study nucleation from solution. Our method allows determining the crystal nucleus size and nucleation rates at constant supersaturation. As an example, we study the homogeneous nucleation of sodium chloride from its supersaturated aqueous solution.
DOI: 10.24435/materialscloud:2019.0031/v1
title: Reaction pathway of oxygen evolution on Pt(1 1 1) revealed through constant Fermi level molecular dynamics

description: The pathway of the oxygen evolution reaction at the Pt(1 1 1)/water interface is disclosed through constant Fermi level molecular dynamics. Upon the application of a positive bias potential H2O_ads and O_ads adsorbates are found to arrange in a hexagonal lattice with an irregular alternation. Increasing further the electrode potential then induces the oxygen evolution reaction, which is found to proceed through a hydrogen peroxide intermediate. Calculation of the associated overpotential shows a reduction of 0.2 eV compared to the associative mechanism. This result highlights the forcefullness of the applied scheme in exploring catalytic reactions in an unbiased way.
DOI: 10.24435/materialscloud:yp-kr
title: Molecular Dynamics Simulations of Crystal Nucleation from Solution at Constant Chemical Potential

description: A widespread method of crystal preparation is to precipitate it from a supersaturated solution. In such a process, control of solution concentration is of paramount importance. The nucleation process, polymorph selection, and crystal habits depend crucially on this thermodynamic parameter. When performing molecular dynamics simulations with a fixed number of molecules in the canonical ensemble, crystal growth is accompanied by a decrease in the solution concentration. This modification of the thermodynamic condition leads to significant artifacts. Inspired by the recent development of the constant chemical potential molecular dynamics simulation method by Perego et al. [J. Chem. Phys. 2015, 142, 144113], we develop a spherical variant of it to study nucleation from solution. Our method allows determining the crystal nucleus size and nucleation rates at constant supersaturation. As an example, we study the homogeneous nucleation of sodium chloride from its supersaturated aqueous solution.