Effect of Exchange-Correlation Functional on the Structural, Mechanical, and Optoelectronic Properties of Orthorhombic RbSrBr3 Perovskite ()
1. Introduction
Over the past decade, metal-halide perovskites (MHP) have generated interest for their exceptional efficiency as both a photovoltaic material and a light-emitting device (LED) [1] [2] [3]. The perovskite ideal structure follows a specific pattern, with certain elements occupying specific positions. In this structure, a monovalent cation takes up the A position, a divalent cation (such as Sr2+ or Pb2+) occupies the B position, and a monovalent halogen anion (such as Cl−, Br−, or I−) fills the X position [4]. In 1978, hybrid organic-inorganic perovskites (HOIP) were reported [5] [6]. Hybrids use organic cations like methylammonium (CH3NH3+), formamidinium (HC(NH2)2+), or ethylammonium ((C2H5)NH3)+ in the A-site [7].
Perovskite possesses an extraordinary variety of physical and chemical properties, making it highly suitable for a diverse range of advanced applications. These include ferroelectricity, piezoelectricity, extremely high temperatures superconductivity, enormous magnetoresistance, charge ordering, thermoelectricity, as well as magneto-transport features [8] [9] [10] [11] [12]. Using the DFT FP-LAPW approach, Rai et al. utilized the GGA and mBJ-based functionals to calculate the electronic properties of the RbMF3 (M = Be, Mg, Ca, Sr, Ba) compound. They stated that the mBJ approach was more accurate when compared to the GGA and LDA approaches [13]. In 2018, a study was conducted by AH Larbi et al. on the electronic and optical properties of orthorhombic RbSrCl3 materials. They utilized the TB-mBJ functional, which yielded superior outcomes compared to the GGA-PBE approximation. The study found that the materials had a bandgap of 7.132 eV and 5.452 eV, respectively [14]. In 2020, a study by HM Ghaithan et al. examined the structural, electronic, and optical properties of CsPbBr3 perovskite in cubic, tetragonal, and orthorhombic structures. The researchers used various exchange-correlation functionals (GGA-PBE, GGA-PBEsol, EV-GGA, mBJ-GGA, nmBJ-GGA, umBJ-GGA) to calculate these properties. They acquired the bands from mBJ-GGA and umBJ-GGA approaches, which demonstrated excellent agreement with experimental data [15]. Multiple researchers, such as K.E. Babu et al. [16] and Zhen-Li et al. [17], have explored the properties of cubic perovskites and their optical characteristics. In these investigations, the FP-LAPW method was employed. More specifically, the generalized gradient approximation (GGA) and LDA approaches were utilized.
A work that was very comparable to this one on orthorhombic RbSrBr3 was not reported previously, as far as we are aware. We developed and explored orthorhombic RbSrBr3 perovskite’s geometrical, electrical, optical, and elastic properties using DFT simulations. For DFT calculations, we utilized four exchange-correlation functionals to analyze the different properties of RbSrBr3.
2. Computational Details
The structural, mechanical, electronic, and optical properties of orthorhombic RbSrBr3 perovskite have been studied employing DFT based on the Cambridge Serial Total Energy Package (CASTEP) code [18] [19] [20] [21] [22]. The Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm is applied to geometrical optimization [23]. To get all the properties at the lowest ground state energy, the structure has been optimized through generalized gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) [24] [25], revised PBE functional for solids (PBEsol) [26], Perdew-Wang (PW91) and local density approximation (LDA) functional Ceperley-Alder-Perdew-Zunger (CAPZ) using ultrasoft pseudopotential [27] with a plane-wave basis set with cut-off energy of 400 eV to get a comprehensive solution along with medium self-consistent field (SCF) tolerance in the program [28]. These functionals were chosen to study electronics and optical properties due to their broad theoretical foundations and good performance in similar experiments. Benchmarking standard GGA-PBE is accurate and efficient. The structure is improved by solid-state material optimization. Even if PBE replaced GGA-PW91, consistency is needed. LDA-CAPZ, a simplified local density approximation, shows advanced GGAs’ effects. We receive stable, non-functional effects [29]. The structure has been optimized using 4 × 1 × 1 Monkhorst-Pack [30] grid in the Brillouin zone of the unit cell. The lattice parameters and the atomic positions within the cell were optimized until the forces and total energy converged within 0.05 eV/Å, 2 × 10−5 eV/atom, and atomic displacement during geometry optimization was less than 0.002Å, respectively.
3. Result and Discussion
3.1. Structural Properties
In the present work, the orthorhombic structure of the RbSrBr3 perovskite that is crystallized with cm cm space group (No. 63) has been optimized using GGA-PBE, GGA-PBEsol, GGA-PW91, and LDA-CAPZ functionals to find out the details of the structure and to judge the functional variation on the properties of the studied compound. The optimized geometrical structure of orthorhombic RbSrBr3 perovskite obtained using various functionals has been depicted in Figure 1. The lattice constants (Å), cell volume (Å3), ground state energy (eV), and bulk modulus have been presented in Table 1. The minimum values of lattice constants, volume, and ground state energy have been observed by utilizing the LDA-CAPZ functional. The GGA-PBE, GGA-PBEsol, GGA-PW91, and LDA-CAPZ functionals have been revealed in the orthorhombic phase of RbSrBr3 perovskite. Compared to other methods, PBEsol-GGA lattice constants were comparable to experimental values.
3.2. Mechanical Properties
The elastic constants Cij describes the crystal’s mechanical and dynamical nature and shows how pressure deforms matter and subsequently returns it to its original shape [31]. Elastic constants impact qualities such as anisotropy, ductility, and hardness and are governed by stress-strain relationships. In an orthorhombic structure, the nine independent elastic constants Cij (C11, C12, C13, C22, C23, C33, C44, C55, C66) can satisfy the Born stability conditions, which are provided in Table 2 and elastic stability criteria are C11 > 0, C44 > 0, C55 > 0, C66 > 0, C11C22 > C212, C11C22C33 + 2C12C13C23 - C11C223 - C22C2 13 - C33C212 > 0 [32].
Figure 1. (a)-(d) The optimized geometrical structure of RbSrBr3 perovskite using the GGA-PBE, GGA-PBEsol, GGA-PW91, and LDA-CAPZ approximation.
Table 1. Summarized lattice parameters, angle, volume, and ground state energy of RbSrBr3 perovskite for different functional.
Functional |
Lattice Parameters (Å) |
Angle (Degree) |
Volume (Å3) |
Ground state energy (eV) |
a |
b |
c |
α |
β |
γ |
GGA-PBE |
4.53 |
15.41 |
11.22 |
90 |
90 |
90 |
782.38 |
−10408.96 |
GGA-PBEsol |
4.43 |
14.85 |
11.07 |
90 |
90 |
90 |
727.95 |
−10408.83 |
GGA-PW91 |
4.50 |
15.27 |
11.28 |
90 |
90 |
90 |
774.93 |
−10417.19 |
LDA-CAPZ |
4.35 |
13.89 |
11.03 |
90 |
90 |
90 |
665.95 |
−10439.52 |
Table 2. Calculated elastic constants (GPa) of orthorhombic RbSrBr3 structure for different functionals.
Functional |
C11 |
C12 |
C13 |
C22 |
C23 |
C33 |
C44 |
C55 |
C66 |
GGA-PBE |
28.06 |
6.85 |
6.38 |
14.57 |
12.23 |
30.18 |
5.27 |
5.19 |
10.35 |
GGA-PBEsol |
33.94 |
10.92 |
8.68 |
17.12 |
16.85 |
45.37 |
6.67 |
5.49 |
14.48 |
GGA-PW91 |
29.57 |
10.91 |
9.13 |
23.05 |
13.84 |
25.27 |
5.01 |
4.74 |
4.82 |
LDA-CAPZ |
40.05 |
16.50 |
13.00 |
21.78 |
21.45 |
44.30 |
7.81 |
7.86 |
17.05 |
The elastic constants presented in Table 2 indicate that orthorhombic RbSrBr3 has met all stability criteria, indicating mechanical stability for all the functionals. The elastic constants found for RbSrBr3 are consistent with the theoretical results given in the literature [33]. The stiffness matrix [34] is positive definite using various functionals, as shown in the eigenvalues listed in Table 3. Furthermore, the parameters, including Bulk modulus (B), shear modulus (G), Young’s modulus (Y), Pugh’s ratio (B/G), Poisson ratio (υ), and anisotropy (A) were measured. The nature of the bonding can be described by Cauchy pressure (Cp), according to Pettifor [35]. If CP is positive, it means that the structure is flexible, and if it is negative, the combination will be brittle [36]. According to the condition above, the structure RbSrBr3 is flexible using various functionals.
Table 3. Calculated eigenvalues of the stiffness matrix in GPa of RbSrBr3 structure using various functionals
Functional |
λ1 |
λ2 |
λ3 |
λ4 |
λ5 |
λ6 |
GGA-PBE |
5.19 |
5.27 |
7.39 |
10.35 |
22.96 |
42.46 |
GGA-PBEsol |
5.49 |
6.67 |
7.57 |
14.48 |
29.35 |
59.52 |
GGA-PW91 |
4.74 |
4.82 |
5.01 |
10.00 |
19.32 |
48.57 |
LDA-CAPZ |
6.69 |
7.81 |
7.86 |
17.05 |
29.06 |
70.39 |
The shear modulus (G) shows how resistant the material is to changes in shape (bond angle) when pressure is applied from the outside. The bulk modulus (B) indicates the material’s resistance to volumetric variations (bond length) [37]. ‘E’ represents Young’s modulus, a method for determining the strength of a material. The value equals the stress-to-strain ratio [38]. Several methods, such as the Voigt, Reuss, and Hill techniques, can be used to derive estimates. Voigt suggested that the strain distribution in composite materials is homogenous [39]. The Reuss approximation is a commonly employed micromechanical model for studying composite materials. This model assumes that the composite phases experience equal stress, specifically normal stress [40]. The Hill convention is a numerical method that averages Reuss and Voigt values [41]. The parameters that determine the strength of the structure have been explained in Table 4.
(1)
(2)
These symbols BH, BV, BR, GH, GV, and GR correspond to the Hill bulk modulus, Voigt bulk modulus, Reuss bulk modulus, Hill shear modulus, Voigt shear modulus, and Reuss shear modulus.
Again, the following equation gives Young’s modulus (E) and Poisson’s ratio (υ) [42]
(3)
(4)
The Poisson ratio is utilized in finite element simulations to design axially stressed components [43]. With the Hill G and B values, we can determine Young’s modulus and Poisson’s ratio. For all functionals, the RbSrBr3 perovskite has been predicted to be ductile, where υ ≥ 1/3 is comparable to the reference [44]. Pugh’s ratio, Poisson’s ratio, and Cauchy’s pressure correlate with the material’s stability and fragility. Pugh’s ratio [45] suggested that the B/G ratio can indicate brittle or ductile material properties. The material is considered flexible when the ratio B to G is greater than 1.75; otherwise, it is classified as brittle. The orthorhombic of RbSrBr3 structure has been showing flexible behavior for all the functionals.
Table 4. Calculated Bulk modulus B (GPa), Shear modulus G (GPa), Young modulus E (GPa), Pugh’s ratio, Poisson’s ratio υ, Compressibility 1/β (GPa−1), Anisotropic factor A, Cauchy’s pressure Cp = C23 - C44 (GPa) of RbSrBr3 using various functionals.
Functional |
B (GPa) |
G (GPa) |
B/G |
E (GPa) |
υ |
1/β (GPa−1) |
A |
Cp |
GGA-PBE |
13.05 |
6.67 |
1.96 |
17.10 |
0.28 |
0.08 |
1.19 |
6.96 |
GGA-PBEsol |
17.30 |
8.11 |
2.13 |
21.05 |
0.30 |
0.06 |
1.95 |
10.18 |
GGA-PW91 |
16.18 |
5.64 |
2.87 |
15.15 |
0.34 |
0.06 |
0.39 |
8.83 |
LDA-CAPZ |
21.94 |
8.80 |
2.49 |
23.28 |
0.32 |
0.05 |
2.05 |
13.64 |
Ranganathan and Ostoja-Starzewski established the universal anisotropic factor [46] to measure crystal anisotropy. The anisotropic factor is essential for understanding defect dynamics, structure stability, and compound elastic anisotropy [47].
(5)
The equation shows that the universal anisotropic A is greater than 1, which means that the crystalline structure of RbSrBr3 has elastic anisotropy [48].
Young’s modulus, shear modulus, compressibility (the opposite of bulk modulus), and Poisson’s ratio should be spherical for isotropic solids. If they don’t stay in the shape of a sphere, it means that anisotropy is present. Figure 2 revealed the directional dependency of Young’s modulus, shear modulus, compressibility, and Poisson’s ratio for the ELATE-generated RbSrBr3 system [49].
The Debye temperature (θD), which correlates to the crystal’s most extensive normal mode of vibration, is one of the most fundamental properties of solids. This parameter establishes a relationship between a solid’s elastic and thermodynamic properties, containing phonons, thermal expansion, conductivity, specific heat, and lattice enthalpy. The average sound velocity, Vm and ρ, the mass density of the solid has to be used to calculate the Debye temperature (θD) The Debye temperature (θD) is given by [50].
(6)
Figure 2. 3D directional dependences of RbSrBr3’s (a) Young modulus, (b) compressibility, (c) shear modulus, and (d) Poisson’s ratio using various functionals.
In Equation (6), h represents Planck’s constant, kB represents Boltzmann’s constant, NA represents Avogadro’s number, ρ represents mass density, M represents molecular weight, and n represents the number of atoms in the cell. The LDA-CAPZ functional has been calculated to be the highest, and the GGA-PW91 functional has shown the lowest Debye Temperature, (θD) for RbSrBr3, which is listed in Table 5. The cut-off debye temperature determines the lattice stability [51]. At temperatures below the Debye temperature, the atomic motion within the crystal lattice is negligible, enabling the free movement of electrons through the lattice planes with minimal dispersion due to the low electron-phonon coupling [52]. This solid becomes a weak heat conductor when the scattering level increases above the Debye temperature.
Table 5. Calculated values of transverse sound wave velocity (VT) and longitudinal sound wave velocity (VL), average sound velocity (Vm), melting temperature (Tm), and Debye temperature (θD) of RbSrBr3 using various functionals.
Functional |
VT (m/s) |
VL (m/s) |
Average sound velocity, Vm (m/s) |
Melting Temperature, Tm (K) |
Debye Temperature, θD (K) |
GGA-PBE |
2751.33 |
4991.11 |
3066.16 |
483.45 |
131.90 |
GGA-PBEsol |
2925.28 |
5446.45 |
3359.95 |
523.88 |
142.37 |
GGA-PW91 |
2518.62 |
5166.56 |
2830.79 |
480.62 |
123.67 |
LDA-CAPZ |
2916.43 |
5704.42 |
3267.57 |
540.60 |
147.72 |
From the obtained bulk modulus (B) and shear modulus (G), the calculations of the average sound velocity (Vm), transverse sound velocity (VT), and longitudinal sound wave velocity (VL) are derived using Navier’s equation [53] which is given below:
(7)
(8)
(9)
Table 5 shows that the speed of the longitudinal wave is twice that of the transverse wave, which is comparable to the literature [36].
Fine et al. put together an empirical formula to obtain the melting temperature Tm of this structure via elastic constants [54]
(10)
The melting temperature of crystalline materials indicates cohesive energy, atomic bond strength, and lattice anharmonicity. This material’s melting temperature was found to be low, which means it is soft. The observed Tm values for all the functionals are lower than those of PbMO3 (M = Sb, Bi) [55].
3.3. Electronic Structure
3.3.1. Band Structure
The electronic band structure of RbSrBr3 perovskite has been calculated using various functionals (GGA-PBE, GGA-PBEsol, GGA-PW91, LDA-CAPZ) through the symmetry points of the Brillouin zone (G→Z→T→Y→S→X→U→R) where the valence band maximum (VBM) was fixed at 0 eV. The band structure of RbSrBr3 perovskite for all the functionals has been presented in Figure 3. The values of bandgap for four different functions lie between 4.293 eV to 4.494 eV. The lowest band gap has been found using GGA-PBE functional, whereas the highest band gap has been obtained using GGA-PBEsol functional. The obtained values of band gap of the present studied perovskite signify that they are insulating-like. The size, form, symmetry, and geometry of the Brillouin zone have been modified by the band structure [56]. Direct band gap nature has been found in RbSrBr3 structure for all functions because the valence band maximum (VBM) and conduction band minimum (CBM) have been found at the “G” of the K-point.
Figure 3. Calculated electronic band structure of RbSrBr3 found using (a) GGA-PBE, (b) GGA-PBEsol, (c) GGA-PW91, (d) LDA-CAPZ functionals.
3.3.2. Density of State
The total density of state (TDOS) and partial density of state (PDOS) of orthorhombic RbSrBr3 structure calculated using various functionals have been presented in Figure 4(a)-(d). Figure 4(e) represents the compared total density of states (TDOS) of RbSrBr3 structure among four different functions. The atomic bonding in the RbSrBr3 compound has been investigated through PDOS. The valence band maximum for all functionals has been formed from the p orbital of Br, with a small contribution from the d orbital of Sr. On the other hand, the
(a) (b)
(c) (d)
(e)
Figure 4. Calculated Partial Density of State (PDOS) of RbSrBr3 utilizing (a) GGA-PBE, (b) GGA-PBEsol, (c) GGA-PW91, (d) LDA-CAPZ functionals, and (e) Compared Total Density of State (TDOS) of RbSrBr3 structure for four different functions.
conduction band minimum has been derived from the d orbital of Sr with a minor contribution from the s orbital of Sr using various functionals [15]. The total density of states (TDOS) plotted in Figure 4 revealed that the halogen atom (Br) at the Fermi level edges contributes significantly to DOS in VB, resulting in flat bands at the band edges [57].
3.4. Optical Properties
The optical characteristics of a material have an inherent connection to its ground-state electrical configuration. Different optical properties of the orthorhombic RbSrBr3 structure with the variation of energy (eV) for different functions have been presented in Figure 5 and Figure 6.
Figure 5. Variations of absorption coefficient, optical conductivity, and reflectivity for change in energy (eV) for RbSrBr3 structure using different functionals.
The absorption coefficient is a measure of the photon absorption ability of materials [58]. The absorption coefficient of orthorhombic RbSrBr3 perovskite has been presented in Figure 5(a). The major absorption peaks have been obtained at 9.21 eV, 9.21 eV, 9.30 eV, and 8.88 eV for GGA-PBE, GGA-PBEsol, GGA-PW91, and LDA-CAPZ, respectively. The higher value of absorption peak (>104 cm−1) has been found almost at the same energy for all functions and all major peaks lie in the UV region.
The optical conductivity spectrum represents the quantity of free charge carriers generated through bond breakage during electron-photon interaction [59]. The optical conductivity is shown in Figure 5(b). The highest peaks have been obtained at 7.52 eV, 7.73 eV, 7.44 eV, and 7.69 eV for GGA-PBE, GGA-PBEsol, GGA-PW91, and LDA-CAPZ, respectively. In this structure, zero optical conductivity (OC) and higher optical conductivity (OC) have been found in the UV, and Visible regions for all the functions.
Figure 6. The variation of (a) dielectric function, (b) refractive index, and (c) loss function concerning change in energy (eV) for RbSrBr3 for different functionals.
Reflectivity measures the fraction of energy that is reflected from the crystal structure [60]. It has been varied from zero to unit. The absorption of light and reflectivity are closely related to each other. The variation of reflectivity with the change of energy in the RbSrBr3 structure for all functions is presented in Figure 5(c). The lowest value of reflectivity has been found in the visible and IR regions for all the functions.
The dielectric function can be expressed as [61]
(11)
Here, ε1(ω) represents a real part of the dielectric function. The higher value of the real part of the dielectric function indicates a stronger ability of the material to polarize in response to the electric field. The imaginary part ε2(ω) signifies the molecular polarization loss due to fluctuations in the external electric field [62]. The dielectric function of orthorhombic RbSrBr3 perovskite has been illustrated in Figure 6(a). The imaginary part and the real part are connected through the Kramers-Kronig relationships [63]. The absorption coefficient has been found in the imaginary part of the dielectric function. The threshold energy for the imaginary dielectric function has been found at 3.90 eV for all the functionals that were correlated with the energy band gap and worked for the UV-region. The principal peak has been located at 6.87 eV, 7.52 eV, 6.87 eV and 6.24 eV for GGA-PBE, PBEsol, PW91 and LDA-CAPZ functionals, respectively. The static dielectric function has been obtained at 5.67 eV for all functionals. The Penn relation (ε1(0) ≈ 1 + (ħω/Eg)2) represents that the static dielectric function is inversely proportional to the energy band gap. The dielectric function becomes negative with the increment of energy, indicating that the medium reflects electromagnetic waves, exhibiting its metallic properties.
The complex refractive function has been calculated by the relation:
(12)
where n(ω) and k(ω) are the refractive index and the extinction index, respectively, which could be found from a real and imaginary portion of the dielectric function employing the previously mentioned relation [64]. The refractive function versus energy curve has been found in Figure 6(b). The value of the refractive index of RbSrBr3 structure was 1.78, 1.79, 1.70, and 1.87 for GGA-PBE, PBEsol, PW91, and LDA-CAPZ functions, respectively. The refractive index peak is found in the UV region for all the functions, and it decreases with the increase of energy. Extinction coefficient k(ω) signifies the light absorption ability of the compound at a certain range. The profiles of k(ω) and ε2(ω) were closely connected in their respective functionals. The loss function L(ω) has been displaced in Figure 6(c). The function L(ω) is essential for representing the energy that fast electrons lose as they move through a solid [65]. The peaks observed in the L(ω) spectra have been suggestive of the properties related to the plasma resonance. The resonant energy loss was obtained at 12.04 eV.
4. Conclusion
In summary, the structural, mechanical, and optoelectronic properties of inorganic orthorhombic RbSrBr3 perovskite have been studied employing GGA-PBE, GGA-PBEsol, GGA-PW91, and LDA-CAPZ functions using density functional theory (DFT). The optimized lattice parameters of the RbSrBr3 perovskite fulfilled the criteria of orthorhombic crystal structure using all the functions. The obtained electronic band gaps of orthorhombic RbSrBr3 perovskite are 4.296 eV, 4.494 eV, 4.339 eV, and 4.309 eV using GGA-PBE, GGA-PBEsol, GGA-PW91, and LDA-CAPZ functionals, respectively. Moreover, every band gap has been found direct band gap in nature. The occupied valence band (VB) and unoccupied conduction band (CB) have occurred at the ‘G’ point for every function. The GGA-PBEsol functional has the highest band gap, and the GGA-PBE functional has the lowest band gap. The total and partial density of states were discussed in detail, which is comparable to the reference. The orthorhombic RbSrBr3 perovskite has been found mechanically stable at all functionals. According to Pugh’s ratio and Poisson’s ratio, the structure has shown ductile behavior. The value of the anisotropic factor was greater than 1, and the 3D structure of Young’s modulus, Shear modulus, Compressibility, and Poisson ratio are an anisotropic nature. The higher absorption > 104 cm−1) makes the studied structure as potential candidate for the absorber layer. Refractive index and reflectivity measurement indicate significant photon energy loss in the RbSrBr3 structure. The present study reflects that the RbSrBr3 structure is a promising candidate for scintillator applications and an anti-reflection coating material because of having a wide band gap, optical isotropy, and structural anisotropy nature in the structure.
Data Availability Statement
The data that support the findings of this study are available from the corresponding author upon reasonable request.