Modeling of the Lattice Dynamics in Strontium Titanate Films of Various Thicknesses: Raman Scattering Studies.

While the bulk strontium titanate (STO) crystal characteristics are relatively well known, ultrathin perovskites' nanostructure, chemical composition, and crystallinity are quite complex and challenging to understand in detail. In our study, the DFT methods were used for modelling the Raman spectra of the STO bulk (space group I4/mcm) and 5-21-layer thin films (layer group p4/mbm) in tetragonal phase with different thicknesses ranging from ~0.8 to 3.9 nm. Our calculations revealed features in the Raman spectra of the films that were absent in the bulk spectra. Out of the seven Raman-active modes associated with bulk STO, the frequencies of five modes (2Eg, A1g, B2g, and B1g) decreased as the film thickness increased, while the low-frequency B2g and higher-frequency Eg modes frequencies increased. The modes in the films exhibited vibrations with different amplitudes in the central or surface parts of the films compared to the bulk, resulting in frequency shifts. Some peaks related to bulk vibrations were too weak (compared to the new modes related to films) to distinguish in the Raman spectra. However, as the film thickness increased, the Raman modes approached the frequencies of the bulk, and their intensities became higher, making them more noticeable in the Raman spectrum. Our results could help to explain inconsistencies in the experimental data for thin STO films, providing insights into the behavior of Raman modes and their relationship with film thickness.

Water Splitting on Multifaceted SrTiO3 Nanocrystals: Calculations of Raman Vibrational Spectrum.

Various photocatalysts are being currently studied with the aim of increasing the photocatalytic efficiency of water splitting for production of hydrogen as a fuel and oxygen as a medical gas. A noticeable increase of hydrogen production was found recently experimentally on the anisotropic faces (facets) of strontium titanate (SrTiO3, STO) nanoparticles. In order to identify optimal sites for water splitting, the first principles calculations of the Raman vibrational spectrum of the bulk and stepped (facet) surface of a thin STO film with adsorbed water derivatives were performed. According to our calculations, the Raman spectrum of a stepped STO surface differs from the bulk spectrum, which agrees with the experimental data. The characteristic vibrational frequencies for the chemisorption of water derivatives on the surface were identified. Moreover, it is also possible to distinguish between differently adsorbed hydrogen atoms of a split water molecule. Our approach helps to select the most efficient (size and shape) perovskite nanoparticles for efficient hydrogen/oxygen photocatalytic production.

Oxygen Vacancy Formation and Migration within the Antiphase Boundaries in Lanthanum Scandate-Based Oxides: Computational Study.

The atomic structure of antiphase boundaries in Sr-doped lanthanum scandate (La1-xSrxScO3-δ) perovskite, promising as the proton conductor, was modelled by means of DFT method. Two structural types of interfaces formed by structural octahedral coupling were constructed: edge- and face-shared. The energetic stability of these two interfaces was investigated. The mechanisms of oxygen vacancy formation and migration in both types of interfaces were modelled. It was shown that both interfaces are structurally stable and facilitate oxygen ionic migration. Oxygen vacancy formation energy in interfaces is lower than that in the regular structure, which favours the oxygen vacancy segregation within such interfaces. The calculated energy profile suggests that both types of interfaces are advantageous for oxygen ion migration in the material.

BaCoO3 monoclinic structure and chemical bonding analysis: hybrid DFT calculations.

Cobalt based perovskites have great potential for numerous applications. Contrary to a generally assumed hexagonal space group (SG P63/mmc) model as the ground state of BaCoO3 (BCO), our hybrid DFT calculations with B1WC density functional and the symmetry group-subgroup derived crystal structure model support the ground state of BCO to be indeed monoclinic, in agreement with recent experimental predictions [Chin et al., Phys. Rev. B, 2019, 100, 205139]. We found for the monoclinic BCO that the C-type anti-ferromagnetic low-spin (AFM LS) state (SG P2/c) is energetically only slightly more preferential at 0 K than the ferromagnetic (FM) LS state (SG C2/c). In turn, these monoclinic structures are energetically more favourable than the hexagonal ones, due to slight z-axis tilting. The analysis of density of states (DOS) and crystal orbital overlap population (COOP) shows a significant (almost 2 eV) separation between occupied and empty t2g states (in the spin-down channel and corresponding anti-bonding states) induced by the z-axis tilting.

First principles calculations of the vibrational properties of single and dimer F-type centers in corundum crystals.

The present paper investigates the F-type centers in α-Al2O3 through their electronic and vibrational properties from first principle calculations using a periodic supercell approach, a hybrid functional, and all-electron Gaussian basis sets as implemented in the CRYSTAL17 code. Single F-type and dimer F2-type centers related to oxygen vacancies in various charge states were considered. The defect-induced vibrational modes were identified and found to appear mainly in the low (up to 300 cm-1) and high (above 700 cm-1) frequency regions, depending on the defect charge. The perturbation introduced by the defects to the thermal nuclear motion in the crystal lattice is discussed in terms of atomic anisotropic displacement parameters. The calculated Raman spectra are discussed for the first time for such defects in α-Al2O3, suggesting important information for future experimental and theoretical studies and revealing deeper insight into their behavior.

Hybrid density functional theoretical study of NASICON-type NaxTi2(PO4)3 (x = 1-4).

Sodium Super Ionic Conductor (NASICON) structured phosphate framework compounds represent a very attractive class of materials for their use as Na-ion battery electrodes. A series of NASICON-structured NaxTi2(PO4)3 compounds corresponding to varying degrees of sodiation (x = 1-4) have been investigated using high-level hybrid density functional theory calculations using the Linear Combination of Atomic Orbitals and Gaussian-type basis set formalism together with hybrid B1WC and HSE06 exchange-correlation functionals. Using primitive cells of NaxTi2(PO4)3 compounds with different stoichiometry, sodium sublattice structure and titanium oxidation states are constructed and analyzed using group theoretical symmetry considerations. The existence of mixed titanium oxidation states for x = 4 (Ti2+/Ti3+) and x = 2 (Ti3+/Ti4+) and a single oxidation state for x = 1 (Ti4+) and x = 3 (Ti3+) has been demonstrated. The results show a necessary set of symmetry reductions taking place due to the highest possible sodium/vacancy and titanium charge ordering with changing x. For each composition, an electroneutrality condition for the oxidation states of all atoms was applied which led to the discovery of several energy minima corresponding to different electronic configurations as identified by different Ti magnetic moments. An interesting relation between the bulk electronic properties of NaxTi2(PO4)3 compounds and the variation of sodium content was also found. In addition to sodium and titanium oxidation state charge ordering, the existence of large differences between the origin and the size of the band gap is shown. The band gap changes from the 4.05 eV 2p-3d gap in Na1Ti2(PO4)3 to the 0.59 eV 3d-3d gap in Na4Ti(PO4)3 with extra states due to mixed titanium valence. These results serve as an important electronic structure benchmark for further studies of such polyanion materials and help to explain some important properties of these systems relevant to battery applications.

Interface-induced enhancement of piezoelectricity in the (SrTiO3)m/(BaTiO3)M-m superlattice for energy harvesting applications.

We present the results of a detailed first principles study of the piezoelectric properties of the (SrTiO3)m/(BaTiO3)M-m heterostructure using the 3D STOm/BTOM-m superlattice model. The atomic basis set, hybrid functionals and slabs with different numbers of STO and BTO layers were used. The interplay between ferroelectric (FEz) and antiferrodistortive (AFDz) displacements is carefully analyzed. Based on the experimental data and group theoretical analysis, we deduce two possible space groups of tetragonal symmetry which allow us to reproduce the experimentally known pure STO and BTO bulk phases in the limiting cases, and to model the corresponding intermediate superlattices. The characteristic feature of the space group P4mm (#99) model is atomic displacements in the [001] direction, which allows us to simulate the FEz displacements, whereas the P4 (#75) model besides FEz displacements permits oxygen octahedra antiphase rotations around the [001] direction and thus AFDz displacements. Our calculations demonstrate that for m/M≤ 0.75 layer ratios both models show similar geometries and piezoelectric constants. Moreover, both models predict an approximately 6-fold increase of the piezoelectric constant e33 compared to the BaTiO3 bulk value, albeit at slightly different layer ratios. The obtained results clearly demonstrate that piezoelectricity arises due to the coordinated collective FEz displacements of atoms in both STO and BTO slabs and interfaces and reaches its maximum when the superlattice approaches the point where the tetragonal phase becomes unstable and transforms to a pseudo-cubic phase. We demonstrate that even a single or double layer of BTO is sufficient to trigger FEz displacements in the STO slab, in P4mm and P4 models, respectively.

First-principles calculations of iodine-related point defects in CsPbI3.

We present here first principles hybrid functional calculations of the atomic and electronic structure of several iodine-related point defects in CsPbI3, a material relevant for photovoltaic applications. We show that the presence of neutral interstitial I atoms or electron holes leads to the formation of di-halide dumbbells of I2- (analogous to the well-known situation in alkali halides). Their formation and one-electron energies in the band gap are determined. The formation energy of the Frenkel defect pair (I vacancies and neutral interstitial I atoms) is found to be ∼1 eV, and as such is smaller than the band gap. We conclude that both iodine dumbbells and iodine vacancies could be, in principle, easily produced by interband optical excitation.

Thermodynamic stability of non-stoichiometric SrFeO3-δ: a hybrid DFT study.

SrFeO3-δ is a mixed ionic-electronic conductor with a complex magnetic structure that reveals a colossal magnetoresistance effect. This material and its solid solutions are attractive for various spintronic, catalytic and electrochemical applications, including cathodes for solid oxide fuel cells and permeation membranes. Its properties strongly depend on oxygen non-stoichiometry. Ab initio hybrid functional approach was applied herein to study the thermodynamic stability of a series of SrFeO3-δ compositions with several non-stoichiometries δ, ranging from 0 to 0.5 (SrFeO3-SrFeO2.875-SrFeO2.75-SrFeO2.5) as a function of temperature and oxygen pressure. The results obtained by two approaches, in which either (i) all electrons at Fe atoms explicitly described or (ii) inner core electrons at Fe atoms are replaced by effective core potential, are compared. Based on our calculations, phase diagrams were constructed, allowing the determination of environmental conditions for the existence of stable phases. It is shown that (within an employed model) only the SrFeO2.5 phase appears to be stable. The stability region for this phase was re-drawn on the contour map of oxygen chemical potential, presented as a function of temperature and oxygen partial pressure. A similar analysis was also performed using experimental Gibbs energies of perovskite formation from the elements. The present modelling strongly suggest a significant attraction between neutral oxygen vacancies. These vacancies are created during a series of the abovementioned SrFeO3-δ mutual transformations accompanied by oxygen release.

Electromechanical Properties of Ba(1-x)SrxTiO3 Perovskite Solid Solutions from First-Principles Calculations.

An enhancement of the piezoelectric properties of lead-free materials, which allow conversion of mechanical energy into electricity, is a task of great importance and interest. Results of first-principles calculations of piezoelectric/electromechanical properties of the Ba(1-x)SrxTiO3 (BSTO) ferroelectric solid solution with a perovskite structure are presented and discussed. Calculations are performed within the linear combination of atomic orbitals (LCAO) approximation and periodic-boundary conditions, using the advanced hybrid functionals of density functional theory (DFT). A supercell model allows the investigation of multiple chemical compositions x. In particular, three BSTO solid solutions with x = 0, 0.125, 0.25 are considered within the experimental stability domain of the ferroelectric tetragonal phase of the solid solution (x < 0.3). The configurational disorder with x = 0.25 composition is also investigated explicitly considering the seven possible atomic configurations corresponding to this composition. It is predicted that Sr-doping of BaTiO3 makes it mechanically harder and enhances its electromechanical/piezoelectric properties, which are important for relevant applications.

First-principles calculations of oxygen interstitials in corundum: a site symmetry approach.

Using site symmetry analysis, four possible positions of interstitial oxygen atoms in the α-Al2O3 hexagonal structure have been identified and studied. First principles hybrid functional calculations of the relevant atomic and electronic structures for interstitial Oi atom insertion in these positions reveal differences in energies of ∼1.5 eV. This approach allows us to get the lowest energy configuration, avoiding time-consuming calculations. It is shown that the triplet oxygen atom is barrierless displaced towards the nearest regular oxygen ion, forming a singlet dumbbell (split interstitial) configuration with an energy gain of ∼2.5 eV. The charge and spatial structure of the dumbbell is discussed. Our results are important, in particular, for understanding the radiation properties and stability of α-Al2O3 and other oxide crystals.

Thermodynamic stability of stoichiometric LaFeO3 and BiFeO3: a hybrid DFT study.

BiFeO3 perovskite attracts great attention due to its multiferroic properties and potential use as a parent material for Bi1-xSrxFeO3-δ and Bi1-xSrxFe1-yCoyO3-δ solid solutions in intermediate temperature cathodes of oxide fuel cells. Another iron-based LaFeO3 perovskite is the end member for well-known solid solutions (La1-xSrxFe1-yCoyO3-δ) used for oxide fuel cells and other electrochemical devices. In this study an ab initio hybrid functional approach was used for the study of the thermodynamic stability of both LaFeO3 and BiFeO3 with respect to decompositions to binary oxides and to elements, as a function of temperature and oxygen pressure. The localized (LCAO) basis sets describing the crystalline electron wave functions were carefully re-optimized within the CRYSTAL09 computer code. The results obtained by considering Fe as an all-electron atom and within the effective core potential technique are compared in detail. Based on our calculations, the phase diagrams were constructed allowing us to predict the stability region of stoichiometric materials in terms of atomic chemical potentials. This permits determining the environmental conditions for the existence of stable BiFeO3 and LaFeO3. These conditions were presented as contour maps of oxygen atoms' chemical potential as a function of temperature and partial pressure of oxygen gas. A similar analysis was also performed using the experimental Gibbs energies of formation. The obtained phase diagrams and contour maps are compared with the calculated ones.

Formation and migration of oxygen vacancies in La(1-x)Sr(x)Co(1-y)Fe(y)O(3-δ) perovskites: insight from ab initio calculations and comparison with Ba(1-x)Sr(x)Co(1-y)Fe(y)O(3-δ).

The formation and migration of oxygen vacancies in the series of (La,Sr)(Co,Fe)O(3-δ) perovskites, which can be used as mixed conducting SOFC cathode materials and oxygen permeation membranes, are explored in detail by means of first principles density functional calculations. Structure distortions, charge redistributions and transition state energies during the oxygen ion migration are obtained and analyzed. Both the overall chemical composition and vacancy formation energy are found to have only a small impact on the migration barrier; it is rather the local cation configuration which affects the barrier. The electron charge transfer from the migrating O ion towards the transition metal ion in the transition state is much smaller in (La,Sr)(Co,Fe)O(3-δ) compared to (Ba,Sr)(Co,Fe)O(3-δ) perovskites where such a charge transfer makes a significant contribution to the low migration barriers observed (in particular for high Ba and Co content).

Adsorption of atomic and molecular oxygen on the LaMnO3(001) surface: ab initio supercell calculations and thermodynamics.

We present and discuss the results of ab initio DFT plane-wave supercell calculations of the atomic and molecular oxygen adsorption and diffusion on the LaMnO(3) (001) surface which serves as a model material for a cathode of solid oxide fuel cells. The dissociative adsorption of O(2) molecules from the gas phase is energetically favorable on surface Mn ions even on a defect-free surface. The surface migration energy for adsorbed O ions is found to be quite high, 2.0 eV. We predict that the adsorbed O atoms could penetrate the electrode first plane when much more mobile surface oxygen vacancies (migration energy of 0.69 eV) approach the O ions strongly bound to the surface Mn ions. The formation of the O vacancy near the O atom adsorbed atop surface Mn ion leads to an increase of the O-Mn binding energy by 0.74 eV whereas the drop of this adsorbed O atom into a vacancy possesses no energy barrier. Ab initio thermodynamics predicts that at typical SOFC operation temperatures (approximately 1200 K) the MnO(2) (001) surface with adsorbed O atoms is the most stable in a very wide range of oxygen gas pressures (above 10(-2) atm).

Evidence for interfacial-storage anomaly in nanocomposites for lithium batteries from first-principles simulations.

We present theoretical support for a mass storage anomaly proposed for nanocomposites in the context of lithium batteries which forms the transition between an electrostatic capacitive mechanism and an electrode mechanism. Ab initio atomic and electronic structure calculations, performed on the Ti(0001)/Li2O(111) model interface, indicate the validity of the phenomenological model of interfacial Li storage and provide a deeper insight into the local situation. Beyond the specific applicability to storage devices, the possibility of a two-phase effect on mass storage generally highlights the availability of novel degrees of freedom in materials research when dealing with nanocomposites.