Source author record

Ricardo Grau-Crespo

Ricardo Grau-Crespo appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2022arXiv

Co-substituted BiFeO3: electronic, ferroelectric, and thermodynamic properties from first principles

Bismuth ferrite, BiFeO3, is a multiferroic solid that is attracting increasing attention as a potential photocatalytic material, because the ferroelectric polarisation enhances the separation of photogenerated carriers. With the motivation of finding routes to engineer the band gap and the band alignment, while conserving or enhancing the ferroelectric properties, we have investigated the thermodynamic, electronic and ferroelectric properties of BiCoxFe1 xO3 solid solutions, with 0 < x < 0.13, using density functional theory. We show that the band gap can be reduced from 2.9 eV to 2.1 eV by cobalt substitution, while simultaneously increasing the spontaneous polarisation, which is associated with a notably larger Born effective charge of Co compared to Fe cations. We discuss the interaction between Co impurities, which is strongly attractive and would drive the aggregation of Co, as evidenced by Monte Carlo simulations. Phase separation into a Co-rich phase is therefore predicted to be thermodynamically preferred, and the homogeneous solid solution can only exist in metastable form, protected by slow cation diffusion kinetics. Finally, we discuss the band alignment of pure and Co-substituted BiFeO3 with relevant redox potentials, in the context of its applicability in photocatalysis.

preprint2022arXiv

Spinel nitride solid solutions: charting properties in the configurational space with explainable machine learning

Ab initio prediction of the variation of properties in the configurational space of solid solutions is computationally very demanding. We present an approach to accelerate these predictions via a combination of density functional theory and machine learning, using the cubic spinel nitride GeSn$_2$N$_4$ as a case study, exploring how formation energy and electronic bandgap are affected by configurational variations. Furthermore, we demonstrate the utility of applying explainable machine learning to understand the crystal chemistry origins of the trends that we observe. Different configuration descriptors (Coulomb matrix eigenspectrum, many-body tensor representation, and cluster correlation function vectors) are combined with different models (linear regression, gradient-boosted decision tree, and multi-layer perceptron) to extrapolate the calculation of ab initio properties from a small set of configurations to the full space with thousands of configurations. We discuss the performance of different descriptors and models. SHAP (SHapley Additive exPlanations) analysis of the machine learning models highlights how values of formation energy are dominated by variations in local crystal structure (single polyhedral environments), while values of electronic bandgap are dominated by variations in more extended structural motifs. Finally, we demonstrate the usefulness of this approach by constructing structure-property maps, identifying important configurations of GeSn$_2$N$_4$ with extremal properties, as well as by calculating accurate equilibrium properties using configurational averaging.

preprint2021arXiv

Mixed-anion mixed-cation perovskite (FAPbI$_3$)$_{0.875}$(MAPbBr$_3$)$_{0.125}$: an ab-initio molecular dynamics study

Mixed-anion mixed-cation perovskites with (FAPbI$_3$)$_{1-x}$(MAPbBr$_3$)$_x$ composition have allowed record efficiencies in photovoltaic solar cells, but their atomic-scale behaviour is not well understood yet, in part because their theoretical modelling requires consideration of complex and interrelated dynamic and disordering effects. We present here an ab initio molecular dynamics investigation of the structural, thermodynamic, and electronic properties of the (FAPbI$_3$)$_{0.875}$(MAPbBr$_3$)$_{0.125}$ perovskite. A special quasi-random structure is proposed to mimic the disorder of both the molecular cations and the halide anions, in a stoichiometry that is close to that of one of today's most efficient perovskite solar cells. We show that the rotation of the organic cations is more strongly hindered in the mixed structure in comparison with the pure compounds. Our analysis suggests that this mixed perovskite is thermodynamically stable against phase separation despite the endothermic mixing enthalpy, due to the large configurational entropy. The electronic properties are investigated by hybrid density functional calculations including spin-orbit coupling in carefully selected representative configurations extracted from the molecular dynamics. Our model, that is validated here against experimental information, provides a more sophisticated understanding of the interplay between dynamic and disordering effects in this important family of photovoltaic materials.

preprint2016arXiv

Electron and phonon transport in shandite-structured Ni$_3$Sn$_2$S$_2$

The shandite family of solids, with hexagonal structure and composition A3M2X2 (A = Ni,Co,Rh,Pd; M = Pb,In,Sn,Tl; X = S,Se), has attracted recent research attention due to promising applications as thermoelectric materials. Herein we discuss the electron and phonon transport properties of shandite-structured Ni3Sn2S2, based on a combination of density functional theory (DFT), Boltzmann transport theory, and experimental measurements. Ni3Sn2S2 exhibits a metallic and non-magnetic groundstate with Ni$^0$ oxidation state and very low charge on Sn and S atoms. Seebeck coefficients obtained from theoretical calculations are in excellent agreement with those measured experimentally between 100 and 600 K. From the calculation of the ratio $σ$/$τ$ between the electronic conductivity and relaxation time, and the experimental determination of electron conductivity, we extract the variation of the scattering rate (1/$τ$) with temperature between 300 and 600 K, which turns out to be almost linear, thus implying that the dominant electron scattering mechanism in this temperature range is via phonons. The electronic thermal conductivity, which deviates only slightly from the Wiedemann-Franz law, provides the main contribution to thermal transport. The small lattice contribution to the thermal conductivity is calculated from the phonon structure and third-order force constants, and is only ~2 Wm$^{-1}$K$^{-1}$ at 300 K (less than 10% of the total thermal conductivity), which is confirmed by experimental measurements. Overall, Ni$_3$Sn$_2$S$_2$ is a poor thermoelectric material (ZT~0.01 at 300 K), principally due to the low absolute value of the Seebeck coefficient. However, the understanding of its transport properties will be useful for the rationalization of the thermoelectric behavior of other, more promising members of the shandite family.

preprint2016arXiv

The symmetry-adapted configurational ensemble approach to the computer simulation of site-disordered solids

Site-occupancy disorder, defined as the non-periodic occupation of lattice sites in a crystal structure, is a ubiquitous phenomenon in solid-state physics and chemistry. Examples are mineral solid solutions, synthetic non-stoichiometric compounds and metal alloys. The experimental investigation of these materials using diffraction techniques only provides averaged information of their structure. However, many properties of interest in these solids are determined by the local geometry and degree of disorder, which escape an "average crystal" description, either from experiments or from theory. In this paper, we describe a methodology for the computer simulation of site-disordered solids, based on the consideration of configurational ensembles and statistical mechanics, where the number of occupancy configurations is reduced by taking advantage of the crystal symmetry of the lattice. Thermodynamics and non-thermodynamic properties are then defined from the statistics in the symmetry-adapted configurational ensemble. We will briefly summarize and discuss some recent applications of this methodology to problems in materials science.

preprint2015arXiv

Electronic structure of porphyrin-based metal-organic frameworks and their suitability for solar fuel production photocatalysis

Metal-organic frameworks (MOFs) can be exceptionally good catalytic materials thanks to the presence of active metal centres and a porous structure that is advantageous for molecular adsorption and confinement. We present here a first-principles investigation of the electronic structure of a family of MOFs based on porphyrins connected through phenyl-carboxyl ligands and AlOH species, in order to assess their suitability for the photocatalysis of fuel production reactions using sunlight.

preprint2015arXiv

Engineering the electronic bandgaps and band edge positions in carbon-substituted 2D boron nitride: a first-principles investigation

Modification of graphene to open a robust gap in its electronic spectrum is essential for its use in field effect transistors and photochemistry applications. Inspired by recent experimental success in the preparation of homogeneous alloys of graphene and boron nitride (BN), we consider here engineering the electronic structure and bandgap of C$_{2x}$B$_{1-x}$N$_{1-x}$ alloys via both compositional and configurational modification. We start from the BN end-member, which already has a large bandgap, and then show that (a) the bandgap can in principle be reduced to about 2 eV with moderate substitution of C $(x<0.25)$; and (b) the electronic structure of C$_{2x}$B$_{1-x}$N$_{1-x}$ can be further tuned not only with composition $x$, but also with the configuration adopted by C substituents in the BN matrix. Our analysis, based on accurate screened hybrid functional calculations, provides a clear understanding of the correlation found between the bandgap and the level of aggregation of C atoms: the bandgap decreases most when the C atoms are maximally isolated, and increases with aggregation of C atoms due to the formation of bonding and anti-bonding bands associated with hybridization of occupied and empty defect states. We determine the location of valence and conduction band edges relative to vacuum and discuss the implications on the potential use of 2D C$_{2x}$B$_{1-x}$N$_{1-x}$ alloys in photocatalytic applications. Finally, we assess the thermodynamic limitations on the formation of these alloys using a cluster expansion model derived from first-principles.

preprint2015arXiv

The importance of anisotropic Coulomb interaction in LaMnO$_{3}$

In low-temperature anti-ferromagnetic LaMnO$_{3}$, strong and localized electronic interactions among Mn 3d electrons prevent a satisfactory description from standard local density and generalized gradient approximations in density functional theory calculations. Here we show that the strong on-site electronic interactions are described well only by using direct and exchange corrections to the intra-orbital Coulomb potential. Only DFT+U calculations with explicit exchange corrections produce a balanced picture of electronic, magnetic and structural observables in agreement with experiment. To understand the reason, a rewriting of the functional form of the +U corrections is presented that leads to a more physical and transparent understanding of the effect of these correction terms. The approach highlights the importance of Hund's coupling (intra-orbital exchange) in providing anisotropy across the occupation and energy eigenvalues of the Mn d states. This intra-orbital exchange is the key to fully activating the Jahn-Teller distortion, reproducing the experimental band gap and stabilizing the correct magnetic ground state in LaMnO$_{3}$. The best parameter values for LaMnO$_{3}$ within the DFT (PBEsol) +U framework are determined to be U = 8 eV and J = 1.9 eV.

preprint2014arXiv

Adsorption of organic molecules at the TiO2(110) surface: the effect of van der Waals interactions

Understanding the interaction of organic molecules with TiO2 surfaces is important for a wide range of technological applications. While density functional theory (DFT) calculations can provide valuable insight about these interactions, traditional DFT approaches with local exchange-correlation functionals suffer from a poor description of non-bonding van der Waals (vdW) interactions. We examine here the contribution of vdW forces to the interaction of small organic molecules (methane, methanol, formic acid and glycine) with the TiO2 (110) surface, based on DFT calculations with the optB88-vdW functional. The adsorption geometries and energies at different configurations were also obtained in the standard generalized gradient approximation (GGA-PBE) for comparison. We find that the optB88-vdW consistently gives shorter surface adsorbate-to-surface distances and slightly stronger interactions than PBE for the weak (physisorbed) modes of adsorption. In the case of strongly adsorbed (chemisorbed) molecules both functionals give similar results for the adsorption geometries, and also similar values of the relative energies between different chemisorption modes for each molecule. In particular both functionals predict that dissociative adsorption is more favourable than molecular adsorption for methanol, formic acid and glycine, in general agreement with experiment. The dissociation energies obtained from both functionals are also very similar, indicating that vdW interactions do not affect the thermodynamics of surface deprotonation. However, the optB88-vdW always predicts stronger adsorption than PBE. The comparison of the methanol adsorption energies with values obtained from a Redhead analysis of temperature programmed desorption data suggests that optB88-vdW significantly overestimates the adsorption strength, although we warn about the uncertainties involved in such comparisons.

preprint2012arXiv

Band gap control via tuning of inversion degree in CdIn$_2$S$_4$ spinel

Based on theoretical arguments we propose a possible route for controlling the band-gap in the promising photovoltaic material CdIn$_2$S$_4$. Our \textit{ab initio} calculations show that the experimental degree of inversion in this spinel (fraction of tetrahedral sites occupied by In) corresponds approximately to the equilibrium value given by the minimum of the theoretical inversion free energy at a typical synthesis temperature. Modification of this temperature, or of the cooling rate after synthesis, is then expected to change the inversion degree, which in turn sensitively tunes the electronic band-gap of the solid, as shown here by accurate screened hybrid functional calculations.

preprint2012arXiv

Density functional theory study of rutile VO2 surfaces

We present the results of a density functional theory (DFT) investigation of the surfaces of rutile-like vanadium dioxide, VO2(R). We calculate the surface energies of low Miller index planes, and find that the most stable surface orientation is the (110). The equilibrium morphology of a VO2(R) particle has an acicular shape, laterally confined by (110) planes and topped by (011) planes. The redox properties of the (110) surface are investigated by calculating the relative surface free energies of the non-stoichiometric compositions as a function of oxygen chemical potential. It is found that the VO2(110) surface is oxidized with respect to the stoichiometric composition, not only at ambient conditions but also at the more reducing conditions under which bulk VO2 is stable in comparison with bulk V2O5. The adsorbed oxygen forms surface vanadyl species much more favorably than surface peroxo species.

preprint2012arXiv

Mind the gap but also the spin: why the Heyd-Scuseria-Ernzerhof hybrid functional description of VO2 phases is not correct

In contrast with recent claims that the Heyd-Scuseria-Ernzerhof (HSE) screened hybrid functional can provide a good description of the electronic and magnetic structure of VO2 phases [V. Eyert, Phys. Rev. Lett. 107, 016401 (2011)], we show here that the HSE lowest-energy solutions for both the low-temperature monoclinic (M1) phase and the high-temperature rutile (R) phase, which are obtained upon inclusion of spin polarization, are at odds with experimental observations. For the M1 phase the groundstate is (but should not be) magnetic, while the groundstate of the R phase, which is also spin-polarized, is not (but should be) metallic. The energy difference between the low-temperature and high-temperature phases is also in strong discrepancy with the experimental latent heat.

preprint2011arXiv

The interplay between dopants and oxygen vacancies in the magnetism of V-doped TiO2

Density functional theory calculations indicate that the incorporation of V into Ti lattice positions of rutile TiO2 leads to magnetic V4+ species, but the extension and sign of the coupling between dopant moments confirm that ferromagnetic order cannot be reached via low-concentration doping in the non-defective oxide. Oxygen vacancies can introduce additional magnetic centres, and we show here that one of the effects of vanadium doping is to reduce the formation energies of these defects. In the presence of both V dopants and O vacancies all the spins tend to align with the same orientation. We conclude that V doping favours the ferromagnetic behaviour of TiO2 not only by introducing spins associated with the dopant centres but also by increasing the concentration of oxygen vacancies with respect to the pure oxide.

preprint2010arXiv

Dopant-vacancy binding effects in Li-doped magnesium hydride

We use a combination of ab initio calculations and statistical mechanics to investigate the substitution of Li+ for Mg2+ in magnesium hydride (MgH2) accompanied by the formation of hydrogen vacancies with positive charge (with respect to the original ion at the site). We show here that the binding energy between dopants and vacancy defects leads to a significant fraction of trapped vacancies and therefore a dramatic reduction of the number of free vacancies available for diffusion. The concentration of free vacancies initially increases with dopant concentration, but reaches a maximum at around 1 mol% Li doping and slowly decreases with further doping. At the optimal level of doping, the corresponding concentration of free vacancies is much higher than the equilibrium concentrations of charged and neutral vacancies in pure MgH2 at typical hydrogen storage conditions. We also show that Li-doped MgH2 is thermodynamically metastable with respect to phase separation into pure magnesium and lithium hydrides at any significant Li concentration, even after considering the stabilization provided by dopant-vacancy interactions and configurational entropic effects. Our results suggest that lithium doping may enhance hydrogen diffusion hydride, but only to a limited extent determined by an optimal dopant concentration and conditioned to the stability of the doped phase.

preprint2010arXiv

Phase separation and surface segregation in ceria-zirconia solid solutions

Using a combination of density functional theory calculations and statistical mechanics, we show that a wide range of intermediate compositions of ceria - zirconia solid solutions are thermodynamically metastable with respect to phase separation into Ce-rich and Zr-rich oxides. We estimate that the maximum equilibrium concentration of Zr in CeO2 at 1373 K is ~2%, and therefore equilibrated samples with higher Zr content are expected to exhibit heterogeneity at the atomic scale. We also demonstrate that in the vicinity of the (111) surface, cation redistribution at high temperatures will occur with significant Ce enrichment of the surface, which we attribute to the more covalent character of Zr-O bonds compared to Ce-O bonds. Although the kinetic barriers for cation diffusion normally prevent the decomposition/segregation of ceria-zirconia solid solutions in typical catalytic applications, the separation behaviour described here can be expected to occur in modern three-way catalytic converters, where very high temperatures are reached.

preprint2010arXiv

Vacancy ordering and electronic structure of gamma-Fe2O3 (maghemite): a theoretical investigation

The crystal structure of the iron oxide gamma-Fe2O3 is usually reported in either the cubic system (space group P4332) with partial Fe vacancy disorder or in the tetragonal system (space group P41212) with full site ordering and c/a\approx 3. Using a supercell of the cubic structure, we obtain the spectrum of energies of all the ordered configurations which contribute to the partially disordered P4332 cubic structure. Our results show that the configuration with space group P41212 is indeed much more stable than the others, and that this stability arises from a favourable electrostatic contribution, as this configuration exhibits the maximum possible homogeneity in the distribution of iron cations and vacancies. Maghemite is therefore expected to be fully ordered in equilibrium, and deviations from this behaviour should be associated with metastable growth, extended anti-site defects and surface effects in the case of small nanoparticles. The confirmation of the ordered tetragonal structure allows us to investigate the electronic structure of the material using density functional theory (DFT) calculations. The inclusion of a Hubbard (DFT+U) correction allows the calculation of a band gap in good agreement with experiment. The value of the gap is dependent on the electron spin, which is the basis for the spin-filtering properties of maghemite.