Phase Separation of Br-Doped CsPbI3: A Combined Cluster Expansion, Monte Carlo, and DFT Study

Cluster expansion, which is a method that describes the concentration-dependent thermodynamic properties of materials while maintaining density functional theory accuracy, was used to predict new (CsPbIxBr1−x) structures. The cluster-expansion method generated 42 new stable (CsPb)xIyBrz (where x = 1 to 3 and y and z = 1 to 8) structures and these were ranked as meta-stable structures based on their enthalpies of formation. Monte Carlo calculations showed that CsPbI0.5Br0.5 composition separates into different phases at 300 K, but changes to a homogeneous phase at 700 K, suggesting that a different phase of CsPbI3 may exist at higher temperatures. Among the 42 predicted structures, randomly selected structures around iodide-rich, 50:50, and bromine-rich sites were studied further by determining their electronic, optical, mechanical, and thermodynamic properties using first-principle density functional theory. The materials possess similar properties as cubic Br-doped CsPbI3 perovskites. The mechanical properties of these compounds revealed that they are ductile in nature and mechanically stable. This work suggests that the introduction of impurities into CsPbI3 perovskite materials, as well as compositional engineering, can alter the electronic and optical properties, making them potential candidates for solar cell applications.


Introduction
The depletion of fossil fuels and the rapid growth in world population are the main drivers of research interest to find alternative renewable energy sources.Perovskite solar cells (PSCs) have been largely explored as a prospective source of clean and renewable energy.They have shown remarkable progress, with rapid increases in power conversion efficiency from early reports of approximately 3% in 2009 to over 25% today [1].Despite their excellent optoelectronic properties, including tunable bandgap, high absorption coefficients, high carrier mobility, long diffusion lengths for electrons and holes, small effective masses, and a simple reproducible process, PSC devices must retain high efficiencies, while exhibiting decent stability and acceptable degradation for practical applications.Research on material-constituent mixing often concentrates on examining the stability of various phases, such as ordered intermediate ground states and end members, as well as investigating the miscibility gaps among them [2][3][4][5].Certain intermediate phases can remain stable even at high temperatures and may prove advantageous for specific applications.For some purposes, disordered solid solutions, particularly mixed halides achieved through element mixing, are preferred, as they offer properties that can be adjusted easily.
The perovskite structure tends to show enhanced stability in relation to phase transitions.The impact of halides on halide perovskites and their optical band gap is widely recognized.These halides provide an effective means of adjusting the band gap of these perovskites.For instance, the energy band gap can be changed from 1.5-2.0eV by varying the types of halide ions present in the material [6][7][8].By controlling the energy band gap, the absorption edge of the perovskite may be adjusted, which allows optimization of the efficiency and performance of the device.The X-site anions in perovskite materials directly affect the band edge states and determine the opto-electronic properties of the materials [9].The most common X anions used in perovskites are halide anions F, Cl, Br, and I.These halides have eight valence electrons, and their respective ionic radii are 2.20, 1.96, 1.81, and 1.33 A for I, Br, Cl, and F, respectively.The decrease in ionic radius of the halide anion leads to a reduction in the bond length between the Pb atom and the X anion in a perovskite material.Consequently, the coupling of atomic orbitals between the Pb and the anion becomes stronger, resulting in increased overlap of the wavefunctions.This enhanced overlap has the effect of enlarging the energy band gap.In other words, the energy band gap of perovskite compounds becomes larger as the halide anion's ionic radius decreases.For example, CsPbI 3 has a reported band gap of approximately 1.73 eV [10], and that of CsPbBr 3 is approximately 2.30 eV [11].CsPbBr 3 is regarded as less suitable for solar cell applications because of its large band gap.However, the thermal and phase stability of CsPbBr 3 is better than that of CsPbI 3, which generated significant research interest.CsPbBr 3 has been utilized as an active absorber in various perovskites solar cell (PSC) configurations and has power-conversion efficiency (PCE) up to 9% [9].
Most recently, significant attention has been focused on CsPbX 3 , where X represents I and Br, and they are often combined in various proportions to form CsPbI 3−x Br x (0 < x < 3).By adjusting the halogen ratio, compounds such as CsPbIBr 2 , CsPbI 2 Br, or non-integer ratios are formed.Various mixtures of Cs-based perovskites containing bromide and iodide, with improved thermal stability and reduced band gap, have been tested in perovskite solar cells [12,13].This resulted in a notable increase in PCE, reaching approximately 16%.Beal et al. conducted the first photoluminescence (PL) experimental study on the light stability of mixed halide perovskite CsPbI 3−x Br x with a range of halide compositions [14].The PL spectra of the films were discovered to be stable for 0 ≤ x ≤ 1 and demonstrated spectral redshift for 1 ≤ x ≤ 3 during the illumination procedure (intensity of 100 mW cm 2 ).Chen et al. used fluorescence lifetime imaging to evaluate and obtain a better understanding of ion transport in CsPbI 2 Br films.The study was able to spot Br vacancies moving into dark regions, which led to the existence of I-phase photoluminescence elsewhere [15].Phase segregation has the potential, under certain conditions, to have a significant impact on the performance of mixed halide perovskite photovoltaic devices.A thermodynamic model developed by Brivio et al. showed that MAPb(I 1−x Br x ) 3 ranging from x = 0.3 and 0.6 can function as an extra energy source with light illumination at 300 K to help get beyond the kinetic energy barrier of halide segregation [16].Lin et al. analyzed the electronic properties of CsPbI 3−x Br x to better understand the band gap structure of the mixed halogen perovskite crystal [17].The calculated band gap of CsPbI 3−x Br x at atmospheric pressure showed that the crystal valence band maximum (VBM) and the conduction band minimum (CBM) are located at the Z-point, with high symmetry.This results in direct band-gap characteristics for various halogen ratios.The electron band gap demonstrated an increase and a blue-shift absorption edge when Br was incorporated, resulting in the varied electronegativities of the elements found in the CsPbI 3−x Br x perovskite materials.The presence of I-I and Pb-I bonding results in a Coulombic attraction.When the I -is substituted with the Br − , there is transfer of electrons occurring between the two atoms because of their different electronegativities, causing a change in the total Coulomb energy and the resulting electrontransport characteristics.
Saliba et al. reported that replacing methylammonium with Cs + and halides mixing can improve the electronic coupling across quantum dots, which, in turn, increases carrier mobility [18].The stability of perovskites has been shown to improve with chemical intervention at grain boundaries; this may cause the ingress of water into the perovskite layer to be more difficult.The mixed-halide CsPbI 2 Br showed an enhanced Goldschmidt tolerance factor (0.84), as well as a better stability of the cubic phase, which paves the way for an incredibly stable all-inorganic PSC.The substitution of bromide in CsPbI 3 reduces the temperature required for the formation of the CsPbI 2 Br black phase from 623 K to 523 K, and the CsPbI 2 Br exhibits improved stability compared to that of CsPbI 3 .By introducing Br ions into CsPbI 3 perovskite using mixed halide composition via a two-step approach, Lee et al. reported experimental results with greatly improved ambient thermodynamic stability, exceeding that of hybrid perovskite [19].
The perovskite configuration offers outstanding tunability of the crystal, the electronic structure, the dimensionality, and the chemical content.Previous studies reported unconventional A x B y X 3z perovskites with impressive optoelectronic properties and stability in ambient conditions.Studies suggested that the perovskite formula can be altered into A 2 B 2 X 6 , A 3 B 3 X 9 , and A 3 B 2 X 9 [20][21][22].Saparov et al. characterized a unique lead-free Cs 3 Sb 2 I 9 perovskite derivative as a potential candidate for high-band gap photovoltaics with an optical band gap of 2.05 eV [22].Furthermore, attempts have been made to develop unique lead-free perovskites such as Cs 3 Bi 2 I 3 [23], Cs 3 Cu 2 I 5 [24], and Cs 4 SnBr 6 [25].Dong et al. reported on the geometrical, electronic, and optical properties of doubleperovskite material Cs 2 NaMX 6 (M = In, Tl, Sb, Bi; X = Cl, Br, I) using the first-principles calculations to investigate the halogen anions and metal anions affecting the properties [26].The structural, opto-electronic, and mechanical properties of all-inorganic vacancy-ordered double-perovskite A 2 Sn 1−x Ti x Y 6 (A = K, Rb, Cs; Y = Cl, Br, I) were investigated using DFT, as reported by Chen et al. [27].The band structure and optical absorption results demonstrated that all A 2 Sn 1−x Ti x Br 6 (A = K, Rb, Cs) materials are semiconductors with band gaps between 0.9-1.6 eV, and A 2 TiBr 6 (A = K, Rb, Cs) compounds have the highest optical absorption both in the visible and ultraviolet regions, making them excellent materials for photovoltaic applications.Urmi et al. computed the opto-electronic properties of Cs 2 TiI 6 (Cs 2 TiI x Br 6-x ) halide perovskite and the findings showed that this material has a considerable potential as a low-band gap alternative material for photovoltaic applications [28].In this work, the cluster expansion (CE) method based on a genetic algorithm was used to predict new structures of Br-doped CsPbI 3 and to further predict their structural, electronic, optical, mechanical, and thermodynamic properties using the first-principle calculations.As far as we are aware, theoretical investigation of the three structures with chemical compositions Cs 2 Pb 2 I 3 Br 3 , Cs 2 Pb 2 I 5 Br, and Cs 3 Pb 3 IBr 8 have not been reported before.

Ground State Search
To begin the search for the ground states of the CsPb(I/Br) 3 system using CE, the DFT enthalpies of formation of CsPb(I/Br) 3 were initially calculated.A total of 42 new structures suggested by the genetic algorithm (GA) were predicted.Using these structures, the GA was started and a 2 meV cross-validation score per atom was obtained.In the first run of CE, a new input structure was selected by the CE, which predicted it as a ground state.The total energy of the new structure was then calculated using VASP and added to the input set.Cluster-expansion predictions were modified and fitted in the procedure.The system was considered to have converged when no new ground-state structures were discovered.Figure 1 depicts the final ground-state line related to the FCC ground states of atoms and the predicted enthalpies of formation.

Miscibility Gap
From Table 1, it is apparent that there would be variations within the calculations because of the stochastic method used in the cluster-expansion method.However, based on the final line of the table (the 4th iteration), it can be inferred that the final CE includes 42 structures from the training set, with a CSV score of 0.652 meV/pos for the system.The next column remained unpopulated until the 2nd iteration, implying that the optimized scheme transitioned from the miscible constituents to the miscibility gap mode.There-fore, the system with a miscibility gap has no stable structures with respect to existing pure phases.

Miscibility Gap
From Table 1, it is apparent that there would be variations within the calculations because of the stochastic method used in the cluster-expansion method.However, based on the final line of the table (the 4th iteration), it can be inferred that the final CE includes 42 structures from the training set, with a CSV score of 0.652 meV/pos for the system.The next column remained unpopulated until the 2nd iteration, implying that the optimized scheme transitioned from the miscible constituents to the miscibility gap mode.Therefore, the system with a miscibility gap has no stable structures with respect to existing pure phases.

Monte Carlo
To visually examine the behavior of phase segregation in CsPbI 1−x Br x compounds, simulation cells with a fixed halides content (50.0%concentration) underwent annealing at temperatures that varied from 0 K to 3200 K.This annealing process was carried out using the canonical ensemble simulation technique called Monte Carlo.Observations were made of the spatial distributions of halides (I-Br) atoms within annealing simulation units.Figure 2 presents the spatial distribution patterns of I-Br atoms at four representative temperatures (300 K, 500 K, 600 K, and 700 K) for CsPbI 1−x Br x compounds, as these temperatures exhibited similar spatial morphologies.To enhance the visibility of the spatial distribution morphology, the figures include the display of Cs and Pb atoms in CsPbI 1−x Br x compounds.From the observations depicted in Figure 3, it is evident that there is a segregation of I and Br halide atoms in the considered CsPbI 1−x Br x compounds between 300 K and 500 K.This indicates that phase segregation can occur in CsPbI 1−x Br x compounds, leading to the separation of CsPbI 3 (or I-rich CsPbI 1−x Br x ) and CsPbBr 3 (or Br-rich CsPbI 1−x Br x ) at a low temperature.As the temperature increases, the agglomerated clusters gradually disperse at 600 K and the spatial distribution of the halides (I and Br) atoms becomes homogenous at 700 K.The homogeneous spatial distribution at elevated temperatures could be attributed to the increased thermal energy allowing the constituent atoms and ions (cesium, lead, and halides) to move more freely, filling the available space more uniformly.tial distribution morphology, the figures include the display of Cs and Pb atoms in CsPbI1−xBrx compounds.From the observations depicted in Figure 3, it is evident that there is a segregation of I and Br halide atoms in the considered CsPbI1−xBrx compounds between 300 K and 500 K.This indicates that phase segregation can occur in CsPbI1−xBrx compounds, leading to the separation of CsPbI3 (or I-rich CsPbI1−xBrx) and CsPbBr3 (or Br-rich CsPbI1−xBrx) at a low temperature.As the temperature increases, the agglomerated clusters gradually disperse at 600 K and the spatial distribution of the halides (I and Br) atoms becomes homogenous at 700 K.The homogeneous spatial distribution at elevated temperatures could be a ributed to the increased thermal energy allowing the constituent atoms and ions (cesium, lead, and halides) to move more freely, filling the available space more uniformly.

Structural Properties
From the 42 structures predicted by cluster expansion, we randomly selected three structures to investigate their structural, electronic, optical, elastic, and thermodynamic properties to predict their opto-electronic behavior.The random selection was mainly

Structural Properties
From the 42 structures predicted by cluster expansion, we randomly selected three structures to investigate their structural, electronic, optical, elastic, and thermodynamic properties to predict their opto-electronic behavior.The random selection was mainly guided by previous studies on unconventional perovskites, indicating that the perovskite formula can be altered into A 2 B 2 X 6 , A 3 B 3 X 9 , and A 3 B 2 X 9 phases, which we reported to be stable.The calculated structural parameters of the selected structures are shown in Figure 4 and Table 2  The calculated electronic band structures for Cs2Pb2I3Br3, Cs2Pb2I5Br, and Cs3Pb3IBr8 using the GGA-PBE functional are shown in Figure 5.The Fermi level represented with a horizontal red line was set at 0 eV to coincide with the VBM and CBM.For Cs2Pb2I3Br3 and Cs2Pb2I5Br, the VBM and CBM were observed at Z-point and at R point for Cs3Pb3IBr8, indicating direct band gap materials, as illustrated in Figure 5.The calculated band gaps of Cs2Pb2I3Br3, Cs2Pb2I5Br, and Cs3Pb3IBr8 were 1.692, 1.520, and 1.901 eV, as seen in Table 3. Again, there are no reported experimental or theoretical values to compare with.However, comparing these structures with the experimental band gap of CsPbI3 (1.73 eV) and CsPbBr3 (2.30 eV) demonstrated that mixing halogens can reduce the band gap of perovskites materials, making them suitable for solar cell applications [4,10].3. Again, there are no reported experimental or theoretical values to compare with.However, comparing these structures with the experimental band gap of CsPbI 3 (1.73 eV) and CsPbBr 3 (2.30eV) demonstrated that mixing halogens can reduce the band gap of perovskites materials, making them suitable for solar cell applications [4,10].

Density of States
The total and partial densities of states of the predicted structures Cs2Pb2I3Br3, Cs2Pb2I5Br, and Cs3Pb3IBr8, which determine the orbital contribution of atoms to the electronic states in the valence and conduction bands, are shown in Figure 6.The Fermi level represented with a vertical red line was set at 0 eV to coincide with the VB and CB.The valence band between −15.0 eV and −5 eV is made up of halides (I-5s and Br-4s orbitals) and a small contribution from Cs-5p and Pb-(4s and 4p) orbitals, as depicted in Figure 6.The valence upper band close to the Fermi level located from −5.0 eV to 0 eV shows contributions of Pb-(4s and 4p) and halides (I-5p and Br-4p) states.The conduction band is mainly comprised of Pb-4p orbitals near the band edge with a small contribution from the halide atoms (I-5p and Br-4p).These halogen atoms contribute to both valence and conduction bands, and this positively contributes to the electronic properties of the perovskite materials.These predicted structures have densities of states that are similar to those of CsPbI3 and CsPbBr3.

Density of States
The total and partial densities of states of the predicted structures Cs 2 Pb 2 I 3 Br 3 , Cs 2 Pb 2 I 5 Br, and Cs 3 Pb 3 IBr 8 , which determine the orbital contribution of atoms to the electronic states in the valence and conduction bands, are shown in Figure 6.The Fermi level represented with a vertical red line was set at 0 eV to coincide with the VB and CB.The valence band between −15.0 eV and −5 eV is made up of halides (I-5s and Br-4s orbitals) and a small contribution from Cs-5p and Pb-(4s and 4p) orbitals, as depicted in Figure 6.The valence upper band close to the Fermi level located from −5.0 eV to 0 eV shows contributions of Pb-(4s and 4p) and halides (I-5p and Br-4p) states.The conduction band is mainly comprised of Pb-4p orbitals near the band edge with a small contribution from the halide atoms (I-5p and Br-4p).These halogen atoms contribute to both valence and conduction bands, and this positively contributes to the electronic properties of the perovskite materials.These predicted structures have densities of states that are similar to those of CsPbI 3 and CsPbBr 3.  2.6.Optical Properties 2.6.1.Dielectric Function Optical characteristics are necessary to determine the effectiveness of the photovoltaic properties of solar cell devices.These computed optical parameters are very important for perovskite materials because they determine the type of response these materials will show when the photons or light are incident upon them.The optical properties of the cubic perovskites CsPbI 3 , CsPbI 2 Br, CsPbBr 2 I, and CsPbBr 3 were calculated using the frequency-dependent dielectric function, as follows: where ε 1 (ω) and ε 2 (ω) are real and imaginary parts, respectively, of the dielectric function, ε(ω).Using the dielectric function (ε 1 (ω) and ε 2 (ω)), other optical parameters, such as absorption coefficient α(ω), reflectivity R(ω), conductivity σ(ω), extinction coefficient K(ω), and refractive index n(ω) [29], can be determined using the following expressions: (2) (3) Figure 7 displays the real and imaginary parts of the dielectric function plotted against the frequency of the predicted structures Cs 2 Pb 2 I 3 Br 3 , Cs 2 Pb 2 I 5 Br, and Cs 3 Pb 3 IBr 8 .The static value of the ε 1 (0) defines the amount of carrier recombination and performance of the opto-electronic device [20].The static values of the dielectric constant ε 1 (0) are 4.82, 5.88, and 4.80 for Cs 2 Pb 2 I 3 Br 3 , Cs 2 Pb 2 I 5 Br, and Cs 3 Pb 3 IBr 8 , respectively.These compounds display the most intense peaks at low energy (visible) regions.The calculated real part maximum values were 6.95 at 2.41 eV, 8.88 at 2.26 eV, and 4.41 at 2.96 eV for Cs 2 Pb 2 I 3 Br 3 , Cs 2 Pb 2 I 5 Br, and Cs 3 Pb 3 IBr 8 respectively.The real-part dielectric function of these compounds has a higher peak in the visible region (low energy) but decreases in the ultra-violet region, which demonstrates that they are suitable for solar cell applications.It was also observed that these predicted structures are red shifted, compared to CsPbI 3 and CsPbBr 3 .
The imaginary dielectric function in Figure 7b, ε 2 , begins at approximately 1.63 eV, 1.47 eV, and 1.82 eV for Cs 2 Pb 2 I 3 Br 3 , Cs 2 Pb 2 I 5 Br, and Cs 3 Pb 3 IBr 8 , respectively, which is related to the band gap of these compounds.As photon energy increases, the calculated imaginary part also increases, and reaches its highest peaks of 5.87 at 3.70 eV, 8.29 at 3.03 eV, and 5.53 at 3.86 eV for Cs 2 Pb 2 I 3 Br 3 , Cs 2 Pb 2 I 5 Br, and Cs 3 Pb 3 IBr 8 , respectively.Once these maximum values are reached, the ε 2 of dielectric functions begins to drop, due to an increase in the energy of the photon.Cs 2 Pb 2 I 5 Br is the second most-intense peak, compared to others, indicating that mixed halides can also improve the intensity of the imaginary dielectric function.

Refractive Index and Extinction Coefficient
The refractive index is defined by the amount of light bent or refracted when entering a material.Figure 8a shows the calculated refractive index as a function of photon energy for the predicted cluster-expansion structures Cs2Pb2I3Br3, Cs2Pb2I5Br, and Cs3Pb3IBr8.The static refractive indices n(0) are 2.20 for Cs2Pb2I3Br3, 2.43 for Cs2Pb2I5Br, and 2.19 for Cs3Pb3IBr8.The refractive indices increase and reach the maximum values of 2.67 at 2.51 eV for Cs2Pb2I3Br3, 3.06 at 2.36 eV for Cs2Pb2I5Br, and 2.69 at 2.72 eV for Cs3Pb3IBr8.The degree of light refraction is determined by the refractive index, which is an essential parameter that is particularly useful in photoelectric applications.The refractive index n(ω) is greater than 1 for all the compounds because photons interact with electrons, slowing down when entering a substance.The higher the refractive index, the greater the amount of slowing down that occurs when photons travel through the material.Typically, any mechanism that increases the electron density in a substance will also increase its refractive index.Cs2Pb2I5Br has the most intense peak, of around 2 eV, when compared to those CsPbI3 and CsPbBr3, indicating that mixed halides can also improve the intensity of the refractive index.
The ability of a material to reduce the amount of light that passes through it is referred to as the extinction coefficient, K(ω). Figure 8b demonstrates the calculated extinction coefficient for the predicted structures.The imaginary part of the refractive index, (n(ω)), is the extinction coefficient (K(ω)) of the materials.It is indicative of the radiation absorbed and has similar behavior an imaginary part, (ε2(ω)), of the dielectric function.The extinction coefficients, K(max), of the highest peaks were 1.37 at 6.71 eV for Cs2Pb2I3Br3, 1.73 at 3.24 eV for Cs2Pb2I5Br, and 1.34 at 4.12 eV for Cs3Pb3IBr3.Among the predicted structures, Cs2Pb2I5Br had the most intense peak in the visible region, compared to those of CsPbI3 and CsPbBr3.

Refractive Index and Extinction Coefficient
The refractive index is defined by the amount of light bent or refracted when entering a material.Figure 8a  The degree of light refraction is determined by the refractive index, which is an essential parameter that is particularly useful in photoelectric applications.The refractive index n(ω) is greater than 1 for all the compounds because photons interact with electrons, slowing down when entering a substance.The higher the refractive index, the greater the amount of slowing down that occurs when photons travel through the material.Typically, any mechanism that increases the electron density in a substance will also increase its refractive index.Cs 2 Pb 2 I 5 Br has the most intense peak, of around 2 eV, when compared to those CsPbI 3 and CsPbBr 3 , indicating that mixed halides can also improve the intensity of the refractive index.

Reflectivity and Conductivity
For a photon of light incident on a material, we can determine how much light is reflected using reflectivity spectra, as shown in Figure 9a, which displays the plots of R(ω) against frequencies from 0 eV to 25 eV.The values of static reflectivity R(0) were observed to be 0.141 for Cs2Pb2I3Br3, 0.174 for Cs2Pb2I5Br, and 0.137 for Cs2Pb2I5Br.As the photon energy increases, the reflectivity also increases and reaches maximum, R(max), values of 0.308 at 14.699 eV for (Cs2Pb2I3Br3)3, 0.339 at 3.191 eV for Cs2Pb2I5Br, and 0.353 at 14.80 eV The ability of a material to reduce the amount of light that passes through it is referred to as the extinction coefficient, K(ω). Figure 8b demonstrates the calculated extinction coefficient for the predicted structures.The imaginary part of the refractive index, (n(ω)), is the extinction coefficient (K(ω)) of the materials.It is indicative of the radiation absorbed and has similar behavior an imaginary part, (ε 2 (ω)), of the dielectric function.The extinction coefficients, K(max), of the highest peaks were 1.37 at 6.71 eV for Cs 2 Pb 2 I 3 Br 3 , 1.73 at 3.24 eV for Cs 2 Pb 2 I 5 Br, and 1.34 at 4.12 eV for Cs 3 Pb 3 IBr 3 .Among the predicted structures, Cs 2 Pb 2 I 5 Br had the most intense peak in the visible region, compared to those of CsPbI 3 and CsPbBr 3 .

Reflectivity and Conductivity
For a photon of light incident on a material, we can determine how much light is reflected using reflectivity spectra, as shown in Figure 9a

Reflectivity and Conductivity
For a photon of light incident on a material, we can determine how much light is reflected using reflectivity spectra, as shown in Figure 9a, which displays the plots of R(ω) against frequencies from 0 eV to 25 eV.The values of static reflectivity R(0) were observed to be 0.141 for Cs2Pb2I3Br3, 0.174 for Cs2Pb2I5Br, and 0.137 for Cs2Pb2I5Br.As the photon energy increases, the reflectivity also increases and reaches maximum, R(max), values of 0.308 at 14.699 eV for (Cs2Pb2I3Br3)3, 0.339 at 3.191 eV for Cs2Pb2I5Br, and 0.353 at 14.80 eV for Cs3Pb3IBr8.The optical conductivities, showing optically activated free charge carriers of the predicted structures as a function of incident photon energy, are depicted in Figure 9b.The maximum σ(ω) was observed at 6.55 eV for Cs2Pb2I3Br3, 13.51 eV for Cs2Pb2I5Br, and 14.184 eV for Cs3Pb3IBr8, and then it started to decrease with increasing photon energy.CsPbI3 had the most intense peak.Cs2Pb2I5Br had the second-most intense peak.).However, the predicted structures red-shifted optical absorption in the visible region suggested that these materials can be applied to solar cells.
Cs2Pb2I3Br3, Cs2PbI5Br, and Cs3Pb3IBr8, respectively.The primary absorption peak was attributed to charge-transfer excitation from the VB to the CB.Peak transitions in these compounds started from electrons of the empty orbital valence band, Cs 6 s and Pb 6p, to I 5p and Br 4p that dominated the deep conduction band.Mixing halogens (I and Br) depicted peaks that were similar to those of the two parent compounds (CsPbI3 and CsPbBr3).However, the predicted structures red-shifted optical absorption in the visible region suggested that these materials can be applied to solar cells.

Elastic Properties
The Cs2Pb2I5Br material belongs to the tetragonal system and has six independent elastic constants, C11, C12, C13, C33, C44, and C66, while Cs2Pb2I3Br3 and Cs3Pb3IBr8 belong to orthorhombic system with nine elastic constants, C11, C12, C13, C22, C23, C33, C44, C55, and C66, as shown in Table 3.These elastic constants determined the mechanical stability of compounds.Using Equations ( 7) and ( 8), the elastic constants of these structures were determined, and they suggested that these predicted structures were mechanically stable.The corresponding mechanical stability criterion for tetragonal crystal structure is defined as: The corresponding mechanical stability criterion for orthorhombic crystal structure is defined as:   3.These elastic constants determined the mechanical stability of compounds.Using Equations ( 7) and ( 8), the elastic constants of these structures were determined, and they suggested that these predicted structures were mechanically stable.The corresponding mechanical stability criterion for tetragonal crystal structure is defined as: The corresponding mechanical stability criterion for orthorhombic crystal structure is defined as: The bulk modulus and shear moduli are crucial characteristics of crystals that demonstrate their capability to resist compression and shear deformations.Table 4 presents the elastic constants, bulk modulus, shear modulus, Young's modulus, Pugh's ratio, and Poisson's ratio of the predicted structures Cs 2 Pb 2 I 3 Br 3 , Cs 2 Pb 2 I 5 Br, and Cs 3 Pb 3 IBr 8 .It was evident that Cs 2 Pb 2 I 3 Br 3 has a high bulk, shear, and Young's moduli.This indicated that the Cs 2 Pb 2 I 3 Br 3 compound had relatively good resistance to deformation, a good capability to resist shear strain, and good rigidity.All the B/G and Poisson ratios for these materials were higher than 1.75 and 0.26, indicating that the materials are ductile [19].

Thermodynamic Properties
The calculated entropy (S), free energy (F), enthalpy (H), and heat capacity against temperature are presented in Figure 11.Entropy refers to the level of disorder of a material, and the middle lines in Figure 11d-f demonstrate the product of entropy and temperature for these materials.The thermal motion of molecules increases as the temperature rises, causing a high degree of disorder and entropy.The top solid lines in Figure 11d-f demonstrate materials' enthalpy, which increases with increasing temperature as all three quantities (internal energy (U), volume (V), and pressure (P)) increase.The bottom three curves in Figure 11 show the free energy of the compound, expressed as F = U − TS.Both U and TS increase with increasing temperature; the increase in TS is greater, resulting in an increase in the value of TS − H.However, the value of H-TS decreases as the temperature increases and, hence, F = U − TS decreases with increasing temperature.The heat capacity of a material is composed of contributions from lattice vibration and electron motion, and Figure 11 depicts the variation in heat capacity with temperature.The melting point of these perovskite compounds increases rapidly below 200 K.At higher temperatures, the melting point increases gradually until it reaches the Dulong-Petit limit.

Cluster Expansion
To determine the structural properties and stability of CsPbI3−xBrx systems, the total energy for all the potential configurations must be calculated.Cluster expansion accounts for the energies of various atomic configurations by generating an Ising-like Hamiltonian

Cluster Expansion
To determine the structural properties and stability of CsPbI 3−x Br x systems, the total energy for all the potential configurations must be calculated.Cluster expansion accounts for the energies of various atomic configurations by generating an Ising-like Hamiltonian and was used to calculate the configurational total energy.The cluster-expansion method was implemented within the UNCLE code [30].The UNCLE code progressively increases the number of clusters considered in the cluster expansion until an ideal degree of precision is attained.This method was employed to determine the ground-state structure and the phase diagram of CsPbI 3−x Br x systems.The cross-validation score is provided as an indicator to assess the accuracy of the cluster-expansion fit.

Monte Carlo
The UNCLE code [30] was employed to perform the Monte Carlo simulations for CsPb(I/Br) 3 systems.To achieve the desired precision in the average concentration of the CsPb(I/Br) 3 compounds, the averaging times were set to 0.1%.The Monte Carlo moves involved exchanging the positions of I and Br atoms.The phase diagrams were calculated using the canonical ensemble (NVT) with specific parameters such as system volume (V), number of particles (N), and temperature (T), as implemented in the UNCLE code.The simulations used a supercell of 25 × 25 × 25 atoms.The temperature range explored in the simulations was set at 100 K to 3200 K.

Geometry Optimization
The first principle DFT was employed to optimize and calculate the properties of predicted structures implemented within VASP [31][32][33][34].GGA was used as an exchangecorrelation functionals [35].The kinetic energy cut-offs were 520, 550, and 500 eV with the k-points of 4 × 4 × 5, 5 × 4 × 5 and 5 × 5 × 6 for the cluster-expansion predicted structures of Cs 2 Pb 2 I 3 Br 3 , Cs 2 Pb 2 I 5 Br, and Cs 3 Pb 3 IBr 8 , respectively.The convergence threshold for the total energy was set at 10 −6 eV/atom.The compounds were considered fully optimized when forces were smaller than 10 −4 eV/A.

Conclusions
The combined methods, including cluster expansion, Monte Carlo, and density functional theory, were employed to study the ground state search and phase stability of (CsPb) x I y Br z (where x = 2 or 3 and y and z = 1 to 8) materials.To conduct Monte Carlo calculations, a set of effective cluster interactions, capable of providing an accurate formation energy for ordered CsPbI 1−x Br x configurations, was initially established within the cluster-expansion framework.There were 42 new ordered phases that were generated from a binary phase diagram.These newly generated phases demonstrated a miscibility gap because of thermodynamically unstable and phase separation.To further study the thermodynamic properties and phase separation at different temperatures for CsPbI 1−x Br x at 50% concentrations, the Monte Carlo method was used.The electronic properties showed that Cs 2 Pb 2 I 3 Br 3 , Cs 2 Pb 2 I 5 Br, and Cs 3 Pb 3 IBr 8 are narrow band gap semiconductors with energy band gap values between 1.60, 1.59, and 1.90 eV, respectively.The optical properties demonstrated a strong ability of photon absorption because of their narrow band gaps.In all these predicted perovskite materials, specific heat capacity at low temperature obeys Debye's law and at high temperature approaches the Dulong-Petit limit.The comparative analysis of predicted structures with CsPbI 3 and CsPbBr 3 showed that the opto-electronic properties of the perovskite compounds are enhanced with the mixed halides and can be used in solar cell applications.

Figure 2 .
Figure 2. Temperature profile cross section through 25 × 25 × 25 Monte Carlo simulation for CsPb(IBr)3.Blue and orange balls are I and Br atoms, respectively.

Figure 3 .
Figure 3. Monte Carlo simulation cells of CsPbI 0.5 Br 0.5 at different temperatures.Blue and orange balls are I and Br atoms, respectively.
. The lattice constants for Cs 2 Pb 2 I 3 Br 3 and Cs 3 Pb 3 IBr 8 are of an orthorhombic system, where a = b = c, whereas Cs 2 Pb 2 I 5 Br has tetragonal system lattice constants with a = b = c, as shown in Table 2.The Pb-I bond length is longer than the Pb-Br bond length for Cs 2 Pb 2 I 3 Br 3 , Cs 2 Pb 2 I 5 Br, and Cs 3 Pb 3 IBr 8 structures because of the difference in the ionic radii of Br and I.We could not find published experimental or theoretical data to compare our calculated structural properties for these predicted structures.Molecules 2024, 29, 92 7 of 18

Figure 4 .
Figure 4.The cluster-expansion predicted structure of Cs2Pb2I3Br3, Cs2Pb2I5Br, and Cs3Pb3IBr8.Blue, yellow and orange balls are Cs, Pb and halides (I and Br) atoms, respectively.

Figure 6 .
Figure 6.The calculated total and partial densities of states of the predicted (a) Cs 2 Pb 2 I 3 Br 3 , (b) Cs 2 Pb 2 I 5 Br, and (c) Cs 3 Pb 3 IBr 8 structures.
shows the calculated refractive index as a function of photon energy for the predicted cluster-expansion structures Cs 2 Pb 2 I 3 Br 3 , Cs 2 Pb 2 I 5 Br, and Cs 3 Pb 3 IBr 8 .The static refractive indices n(0) are 2.20 for Cs 2 Pb 2 I 3 Br 3 , 2.43 for Cs 2 Pb 2 I 5 Br, and 2.19 for Cs 3 Pb 3 IBr 8 .The refractive indices increase and reach the maximum values of 2.67 at 2.51 eV for Cs 2 Pb 2 I 3 Br 3 , 3.06 at 2.36 eV for Cs 2 Pb 2 I 5 Br, and 2.69 at 2.72 eV for Cs 3 Pb 3 IBr 8 .
, which displays the plots of R(ω) against frequencies from 0 eV to 25 eV.The values of static reflectivity R(0) were observed to be 0.141 for Cs 2 Pb 2 I 3 Br 3 , 0.174 for Cs 2 Pb 2 I 5 Br, and 0.137 for Cs 2 Pb 2 I 5 Br.As the photon energy increases, the reflectivity also increases and reaches maximum, R(max), values of 0.308 at 14.699 eV for (Cs 2 Pb 2 I 3 Br 3 ) 3 , 0.339 at 3.191 eV for Cs 2 Pb 2 I 5 Br, and 0.353 at 14.80 eV for Cs 3 Pb 3 IBr 8 .optical conductivities, showing optically activated free charge carriers of the predicted structures as a function of incident photon energy, are depicted in Figure 9b.The maximum σ(ω) was observed at 6.55 eV for Cs 2 Pb 2 I 3 Br 3 , 13.51 eV for Cs 2 Pb 2 I 5 Br, and 14.184 eV for Cs 3 Pb 3 IBr 8 , and then it started to decrease with increasing photon energy.CsPbI 3 had the most intense peak.Cs 2 Pb 2 I 5 Br had the second-most intense peak.

Figure 8 .
Figure 8.The calculated (a) refractive index and (b) extinction coefficient of the pure phase structures (CsPbI3 and CsPbBr3) and the predicted structures (Cs2Pb2I3Br3, Cs2Pb2I5Br, and Cs3Pb3IBr8).

2. 6
Figure 10 depicts the absorption coefficient α(ω) spectra of Cs 2 Pb 2 I 3 Br 3 , Cs 2 Pb 2 I 5 Br, and Cs 3 Pb 3 IBr 8 as a function of photon energy.The absorption spectra demonstrate the penetration capability of photons passing through a material before absorption.The absorption spectra start below 2 eV due to the direct energy band gap of the materials.The strongest peaks were observed in the ultraviolet region at 18.33, 18.28, and 12.16 eV for Cs 2 Pb 2 I 3 Br 3 , Cs 2 PbI 5 Br, and Cs 3 Pb 3 IBr 8, respectively.The primary absorption peak was attributed to charge-transfer excitation from the VB to the CB.Peak transitions in these compounds started from electrons of the empty orbital valence band, Cs 6 s and Pb 6p, to I 5p and Br 4p that dominated the deep conduction band.Mixing halogens (I and Br) depicted peaks that were similar to those of the two parent compounds (CsPbI 3 and CsPbBr 3 ).However, the predicted structures red-shifted optical absorption in the visible region suggested that these materials can be applied to solar cells.

2. 7 .
Elastic Properties The Cs 2 Pb 2 I 5 Br material belongs to the tetragonal system and has six independent elastic constants, C 11 , C 12 , C 13 , C 33 , C 44 , and C 66 , while Cs 2 Pb 2 I 3 Br 3 and Cs 3 Pb 3 IBr 8 belong to orthorhombic system with nine elastic constants, C 11 , C 12 , C 13 , C 22 , C 23 , C 33 , C 44 , C 55 , and C 66 , as shown in Table

Table 1 .
The development of an iterative procedure that lists stable structures constructed of numerous pure phases.

Table 1 .
The development of an iterative procedure that lists stable structures constructed of numerous pure phases.No.