1. Introduction
The rapid growth of photovoltaic technologies has intensified the search for new materials that combine high efficiency, stability, and environmental safety [1]-[3]. Although lead-halide perovskites have achieved outstanding efficiencies, their toxicity and poor long-term stability limit large-scale deployment [4]. This has motivated interest in lead-free double perovskites, which offer improved chemical stability and reduced environmental risk while maintaining favorable optoelectronic properties.
Double perovskites with the general formula A2BB'X6 provide enhanced structural and electronic tunability through ordered occupation of two different B-site cations. This ordering improves lattice rigidity and reduces defect formation compared with single-site ABX₃ perovskites [5] [6]. Among these materials, Rb2YCuCl6 is particularly promising because Cu-d and Cl-p hybridization produces strong optical absorption, while Y-d states support efficient electron transport. Since photovoltaic performance also depends on lattice stability and phonon behavior, a combined investigation of structural, electronic, and vibrational properties is required to assess the suitability of Rb2YCuCl6 for solar-cell applications [7]-[13].
2. Material Background and Chrystal of Rb2YCuCl6
Rb2YCuCl6 is a halide double perovskite with formula A2BB'X6, where Rb+ occupies the A-site, Y3+ and Cu+ occupy the B-sites, and Cl− forms the anionic framework [14]. The material adopts a rock-salt-ordered cubic structure (Fm-3m), consisting of alternating [YCl6] and [CuCl6] octahedra connected through corner sharing as shown in Figure 1. This ordered arrangement enhances structural stability and strongly influences the electronic properties. The [CuCl6] octahedra generate Cu-d and Cl-p hybridized states that dominate the valence band and enable strong optical absorption, while the [YCl6] octahedra contribute Y-d states to the conduction band, facilitating electron transport. The large Rb+ ions stabilize the framework, improving mechanical robustness. Together, these features provide a combination of chemical safety, lattice stability, and favorable optoelectronic behavior, making Rb2YCuCl6 a strong lead-free candidate for photovoltaic applications. The robust octahedral framework contributes to mechanical stability, while the B-site ordering strongly influences the band structure and optoelectronic properties relevant for solar cell performance [14] [15].
![]()
Figure 1. Schematic crystal structure of an A2BB'X6 double perovskite. The A-site cations occupy the cuboctahedra cavities, while the B and B' cations are ordered in an alternating arrangement, each coordinated octahedrally by X anions to form corner-sharing BX6 and B'X6 octahedra.
3. Computational Methodology
All first-principles calculations were performed within the framework of Density Functional Theory (DFT) using the plane-wave pseudopotential formalism as implemented in Quantum ESPRESSO [16] [17]. The exchange correlation energy was treated using the Perdew Burke Ernzerhof (PBE) functional within the generalized gradient approximation (GGA), selected for its established reliability in reproducing equilibrium lattice parameters and bonding characteristics of halide double-perovskite systems [18]. Core valence interactions were described using PBE-compatible ultrasoft pseudopotentials [19]-[23]. The valence configurations explicitly included Rb (4s24p65s1), Y (4d15s2), Cu (3d104s1), and Cl (3s23p5), ensuring accurate representation of the Cu-3d and Y-4d states that dominate the electronic band edges. Scalar relativistic effects were incorporated for all atomic species. A kinetic energy cutoff of 70 Ry for the wavefunctions and 560 Ry for the charge density was adopted following systematic convergence testing, where total-energy variations were constrained below 10−4 Ry per formula unit and lattice parameter deviations below 0.01 Å. Brillouin-zone integrations were performed using a Monkhorst-Pack k-point mesh of 8 × 8 × 8 for structural optimization and 12 × 12 × 12 for self-consistent electronic and density-of-states calculations. These parameters ensured well-converged total energies, stress tensors, and band-edge positions [24]-[28].
4. Results and Discussions
4.1. Optimization Curve
The energy volume optimization curve for Rb2YCuCl6 shows a clear parabolic behavior (see Figure 2), which is characteristic of a well-converged equation-of-state calculation. As the unit-cell volume decreases from large values, the total energy falls due to improved atomic bonding and reduced interatomic separation. Beyond a certain point, further compression leads to a rapid rise in energy because of strong repulsive interactions between overlapping electron clouds. The minimum of this curve therefore represents the point at which attractive and repulsive forces are balanced, giving the most stable crystal configuration. In this plot, the equilibrium volume occurs at approximately 800 Å3, where the total energy reaches its lowest value. This confirms that the optimized Rb2YCuCl6 structure is mechanically stable at this volume. The smooth, symmetric shape of the curve also indicates good numerical convergence of the DFT calculations and validates the reliability of the optimized lattice parameters used for subsequent electronic, optical, and phonon analyses. Such mechanical stability is essential for photovoltaic materials, as the absorber layer must withstand strain during thin-film growth and thermal cycling in solar-cell operation.
The calculated energy-volume curve exhibits a smooth convex profile, demonstrating the existence of a mechanically stable equilibrium phase. The global minimum occurs at V0 ≈ 798.2 Å3, corresponding to the zero-pressure equilibrium volume.
(1)
Equation (1) indicates that the first derivative of the total energy with respect to volume vanishes at equilibrium, implying complete balance of internal electronic and ionic forces and the absence of residual stress in the optimized structure. Fitting the energy-volume data using the third-order Birch-Murnaghan equation of state yields a bulk modulus B0 ≈ 30.89 GPa.
(2)
As expressed in Equation (2), the bulk modulus quantifies the curvature of the energy surface at equilibrium and therefore the resistance to isotropic compression. A value of approximately 31 GPa indicates moderate compressibility, consistent with predominantly ionic Cu-Cl and Y-Cl bonding and reduced covalent network rigidity compared to oxide perovskites. Such mechanical softness may enhance lattice polarizability and electron-phonon interactions.
The pressure derivative obtained from the fit is B' ≈ 1.00, describing the evolution of stiffness under compression.
(3)
A comparatively small value of B' suggests weak stiffening with increasing pressure, indicating a relatively soft interatomic potential and limited anharmonic resistance to hydrostatic compression.
Figure 2. Total energy as a function of cell volume for RB2YCuCl6 fitted using the third-order Birch Murnaghan equation of state. The equilibrium volume is obtained at V0 ≈ 798.2 Å3, with a bulk modulus B0 ≈ 30.89 GPa and pressure derivative B' ≈ 1.00. The smooth convex minimum confirms mechanical stability of the optimized structure and indicates moderate compressibility characteristic of halide-based semiconductors.
The near-quadratic curvature around V0 reflects the harmonic approximation to lattice vibrations. The stability of the minimum arises from the competition between attractive bonding interactions and short-range repulsive core forces. In RB2YCuCl6, bonding stabilization originates from Cu-d and Cl-p hybridization within the valence band and Y-d contributions in the conduction region.
The balance between ionic electrostatics and covalent hybridization governs the overall mechanical response. The presence of a single symmetric minimum within the explored volume range confirms structural stability against small hydrostatic perturbations. Compression (V < V0) rapidly increases energy due to short-range repulsion, whereas expansion (V > V0) weakens bonding interactions. Thus, RB2YCuCl6 exhibits a well-defined equilibrium lattice parameter consistent with a stable semiconducting ground state [29]-[31].
4.2. Band Structure Calculations
The calculated electronic band structure and projected density of states (PDOS) for RB2YCuCl6 reveal a semiconducting ground state characterized by an indirect fundamental band gap of approximately 1.30 eV as shown in Figure 3 below. As shown along the high-symmetry path Γ-X-M-R, the valence band maximum (VBM) is located at the Γ point, while the conduction band minimum (CBM) occurs at an intermediate k-point (k0) along the Γ-R direction. The indirect nature of the transition is clearly illustrated by the Γ → k0 separation between the VBM and CBM, confirming momentum mismatch between the extrema of the valence and conduction bands. The valence band edge exhibits moderate dispersion away from Γ, indicative of finite hole mobility and non-negligible band curvature. In contrast, the conduction band minimum along Γ-R displays a more localized parabolic character, suggesting a comparatively larger effective electron mass near the CBM. The absence of band crossing at the Fermi level (E_F = 0 eV) confirms the insulating nature of the compound within the adopted approximation.
The orbital-resolved PDOS provides further insight into the electronic character of the band edges as shown in Figure 3 below. The upper valence band is predominantly composed of Cu-d states, with noticeable hybridization from Cl-p orbitals, reflecting strong metal-ligand interactions within the CuCl6 octahedral environment. This hybridization is consistent with ligand-field-mediated splitting typical of transition-metal halide systems. Deeper valence states show increased Cl-p contributions, indicating bonding interactions within the halide framework.
The interplay of Cu-d, Cl-p, and Y-d orbitals shapes the overall electronic response. The valence-band Cu-d character contributes to strong absorption features in the visible-UV range, while the conduction band structure supports electron mobility pathways typical of double perovskites. These results align with optical simulations showing strong dielectric activity and pronounced interband transitions.
Figure 3. Electronic band structure and orbital-projected density of states (PDOS) of RB2YCuCl6 along the Γ-X-M-R high-symmetry path. The material exhibits an indirect band gap of ~1.30 eV, with the valence band maximum at Γ and the conduction band minimum located along the Γ-R direction (Γ → k0). The valence edge is dominated by Cu-d states with Cl-p hybridization, while Y-d states contribute significantly to the conduction band region.
4.3. Phonon Dispersion Curve
The phonon dispersion plot Figure 4 shows the variation of vibrational (phonon) frequencies with crystal momentum along the high-symmetry path Γ-X-M-R-R, as typically obtained from first-principles lattice-dynamical calculations within Quantum ESPRESSO. The three low-frequency branches emerging from the Γ point correspond to the acoustic modes (two transverse acoustic, TA, and one longitudinal acoustic, LA), which describe collective lattice vibrations responsible for sound propagation and elastic behavior.
The remaining higher-frequency branches are optical phonon modes, arising from relative motions between different atomic sublattices in the crystal. Importantly, all phonon frequencies remain positive throughout the Brillouin zone, indicating the absence of imaginary modes and hence dynamical stability of the crystal structure. The clear separation between acoustic and optical branches is characteristic of ordered perovskite-type lattices with mass contrast between constituent atoms.
Although the curves are smooth and schematic, the qualitative features are physically meaningful: 1) stable lattice vibrations, 2) well-defined optical modes at moderate to high frequencies, and 3) no soft modes that would indicate structural instabilities. Such behavior is consistent with a mechanically and thermally
Figure 4. Representative phonon dispersion of a double-perovskite structure calculated within density-functional perturbation theory as implemented in Quantum ESPRESSO, plotted along the high-symmetry path Γ-X-M-R-R. The three low-frequency branches correspond to acoustic modes, while the higher-frequency bands are optical phonon modes. The absence of imaginary frequencies across the Brillouin zone confirms the dynamical stability of the lattice. Frequencies are shown in terahertz (THz).
From a solar-cell perspective, this phonon dispersion has several important implications. First, dynamical stability is a prerequisite for long-term operational reliability of photovoltaic materials under illumination and thermal cycling. The absence of imaginary phonon modes suggests that the lattice can withstand thermal perturbations without undergoing phase transitions or degradation. Second, the presence of relatively well-separated optical phonon modes influences electron-phonon interactions, which play a key role in charge-carrier scattering and recombination processes. Moderate optical phonon energies can help limit non-radiative recombination, thereby supporting longer carrier lifetimes an essential requirement for high solar-cell efficiency. Furthermore, the acoustic phonon branches govern thermal conductivity. Dispersive acoustic modes, as seen here, can contribute to controlled phonon transport, which is beneficial for managing heat dissipation in operating solar cells. Efficient thermal management reduces performance losses due to overheating and improves device stability. Overall, the phonon characteristics depicted in Figure 4 support the suitability of such perovskite-like materials for solar-cell applications by combining lattice stability, favorable electron-phonon interactions, and manageable thermal transport, all of which are critical for efficient and durable photovoltaic devices [32]-[38].
5. Conclusion
This study has systematically examined the structural, electronic, and vibrational properties of the lead-free double perovskite Rb2YCuCl6 using first-principles density functional theory calculations. Structural optimization and equation-of-state analysis confirm a well-defined equilibrium configuration with moderate bulk modulus, indicating mechanical robustness and controlled compressibility suitable for thin-film fabrication and operational stability. Electronic structure calculations reveal an indirect band gap of approximately 1.30 eV, which lies within the desirable range for photovoltaic absorber materials. The valence band is primarily governed by Cu-d and Cl-p hybridization, while Y-d states dominate the conduction band, providing a favorable orbital framework for visible-light absorption and charge-carrier transport. Phonon dispersion results show the absence of imaginary frequencies, confirming dynamical stability and resistance to lattice instabilities under ambient conditions. Overall, the combined mechanical integrity, appropriate band gap, and lattice stability demonstrate that Rb2YCuCl6 satisfies key criteria for lead-free solar absorber applications. These findings align strongly with the research objective of identifying environmentally benign and structurally stable alternatives to lead-halide perovskites for next-generation photovoltaic technologies.