Abdesslem
Jedidi
*a,
Shatha M.
Alamri
a,
Norah O.
Alotaibi
a,
Souraya
Goumri-Said
*b and
Mohammed Benali
Kanoun
*c
aChemistry Department, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia. E-mail: ajedidi@kau.edu.sa
bPhysics Department, College of Science and General Studies, Alfaisal University, P.O. Box 50927, Riyadh 11533, Saudi Arabia. E-mail: sosaid@alfaisal.edu
cDepartment of Mathematics and Sciences, College of Humanities and Sciences, Prince Sultan University, P.O. Box 66833, Riyadh 11586, Saudi Arabia. E-mail: mkanoun@psu.edu.sa
First published on 1st August 2024
The solar cell and light-emitting device research community is currently focusing on investigating two-dimensional (2D) hybrid perovskite materials owing to their remarkable stability and intriguing optoelectronic characteristics, which hold significant promise for various applications. In general, the introduction of chirality in hybrid perovskites arises from symmetry breaking within their inorganic frameworks. Nevertheless, despite this understanding, the specific factors driving the observed increase in splitting remain obscure due to a lack of comprehensive investigations. Our research delves into the electronic properties of 2D layered hybrid perovskites, considering their behavior with and without spin–orbit coupling. We specifically focus on effect of Rashba splitting and the impact of electronic structure variation across a range of chiral perovskites by introducing chiral organic cations with differing degrees of π-conjugation, resulting in significant changes in spin-splitting magnitude. Systematic first principles investigations confirm that the distortion of the cage and d-spacing of chiral perovskites are crucial design parameters for achieving strong spin-splitting in 2D layered perovskites. Furthermore, our investigation reveals that these systems exhibit remarkable absorption capabilities in the visible light spectrum, as demonstrated by their computed optoelectronic characteristics. The chiral perovskites described in this study, which exhibit substantial spin-splitting, present a distinctive prototype with potential implications for spintronics and photovoltaics.
Numerous studies14–19 have enriched the field of chiral metal halides, with the majority centered on Pb systems and employing specific chiral organic linkers.16–21 The necessity of eliminating toxic lead has driven the shift towards utilizing alternative metals such as Cu, Bi, Sb, and Sn in the synthesis of lead-free chiral metal halides,13,15,19–25 unveiling a variety of intriguing properties. Substituting tin for lead presents a common strategy aimed at creating less toxic materials, simultaneously narrowing the band gap towards the optimum25 in accordance with the Shockley–Queisser limit.26 Chiral tin halide perovskites, notably (R/S-MBA)2SnI4, have garnered significant interest, initially investigated by Lu et al.27 Their study introduced layered tin iodine perovskites characterized by highly distorted octahedra, influenced by specific hydrogen bonding arrangements. Furthermore, they demonstrated that the presence of chiral MBA cations induces circularly polarized absorption from the inorganic Sn–I sublattice, showcasing chiroptical activity.
Recent studies have employed chiral molecules in Sn-based 2D halide perovskites to fabricate field-effect transistor28 and solar cell devices,29 showcasing enhanced charge injection and greater stability of the inorganic perovskite layers. Very recently, F. De Angelis and colleagues utilized the chiral Cl-MBA to design chiral pairs, Cl-MBA2SnI4.30 Their research unveiled a novel low-bandgap Sn-based 2D chiral perovskite capable of tuning chiroptical properties through Cl-MBA ligand inclusion. This innovation induced Rashba band splitting, resulting in significant spin splitting values observed for the R- and S-enantiomers. These studies have convincingly illustrated the promising capabilities of 2D tin halide perovskite materials.
The advancement of novel lead-free chiral metal halide systems serves as a platform for elucidating the intricate structure–property relationships, particularly regarding Rashba spin splitting, aimed at facilitating efficient spin manipulation. It is noteworthy that the 2D tin halide perovskites have received comparatively less attention than their lead-based counterparts. It is established that organic molecules, serving to facilitate charge transport, can be incorporated into the semiconducting inorganic framework.31 Moreover, chiral organic cations serve dual roles as encapsulating layers, offering potential to enhance the stability of perovskite materials.31 Additionally, by varying the organic cation sites, it is possible to manipulate the characteristics of the Sn iodide framework, thereby allowing for the modulation of optical, spin, and electronic properties. Recent experimental observations suggest that chiral perovskites containing a large π-conjugated organic spacer exhibit the most substantial circular dichroism signals ever reported in chiral perovskites.32,33 However, the influence of chiral organic size and the presence of π bonds remains unexplored in relation to the spintronic properties of chiral perovskites.
To address the concerns raised earlier, we present an in-depth theoretical analysis delving into the stability, crystal structures, and electronic properties of two distinct 2D tin halide perovskite systems: S-cyclohexylethylamine (CYHEA) and S-methylbenzylamine (MBA), as depicted in Scheme 1. Our investigation aims to shed light on their versatile tunable optoelectronic and spintronic attributes. We have harnessed density functional theory calculations to offer a deeper understanding of impact of organic spacers on the intricate interplay between structure and property within 2D hybrid tin halide perovskites. These insights hold the potential to govern and finely tune the optoelectronic properties of these materials. Furthermore, we explore the effect of spacer organic cation on the electronic structures and optical properties of the materials. Through first-principles computations incorporating spin–orbit coupling and considering broken inversion symmetry, we elucidate the potential role of Rashba band splitting in 2D hybrid perovskites. This phenomenon has implications for tailoring their electronic and optical absorption properties, making them suitable for applications in optoelectronics, spintronics. Our findings highlight the significant role of structural distortions in driving the observed Rashba effect.
To further quantify the distortion of inorganic [SnI6]2− octahedra, we calculated the bond angle variances and bond length distortion
![]() | (1) |
![]() | (2) |
To examine the energetic stability of (S-MBA)2SnI4 and (S-CYHEA)2SnI4, the formation energy, Ef, is calculated, defined as follows:
| ET((MBA/CYHEA)2SnI4) = ET((MBA/CYHEA)2SnI4) − 2ET(MBA/CYHEAI) − ET(SnI2) | (3) |
To investigate the optical properties of the 2D tin halide perovskite systems, the optical absorption coefficient, α(ω), are computed using the following equation:
![]() | (4) |
Each CYHEA spacer cation is hydrogen-bonded to the distorted SnI6 octahedrons through three hydrogen bonds, adopting an axial configuration. The hydrogen bond distances to the bridging halides measure longer at 2.85 Å and 2.69 Å compared to the axial halides at 2.45 Å and 2.60 Å, and 2.57 Å and 2.41 Å, respectively, for two adjacent cations. Similarly, in CYHEA-based 2D perovskite, each ordered MBA spacer cation is linked to the [SnI6]2− octahedron by three hydrogen bonds with an axial configuration, as illustrated in Fig. 2 and summarized in Table S3 (ESI†). Furthermore, chiral organic cations in both 2D chiral halide perovskites engage in van der Waals interactions with adjacent cations. These non-covalent interactions contribute to the stabilization of the crystal packing in the 2D systems. The energetically stable configuration was investigated by computing the formation energy, Ef, of the two structures, using eqn (3). The computed values for the formation energies of (S-MBA)2SnI4 and (S-CYHEA)2SnI4 are −9.55 eV and −14.92 eV, respectively. Additionally, the greater negative formation energy observed for (S-CYHEA)2SnI4 suggests enhanced stability.
![]() | ||
| Fig. 2 Hydrogen bonding interaction between iodide atoms and the hydrogen on the organic cation. (a) (S-MBA)2SnI4, and (b) (S-CYHEA)2SnI4. | ||
Next, the electronic properties of the (S-MBA)2SnI4 and (S-CYHEA)2SnI4 are explored by calculating the electronic band structure at GGA-PBE levels of approximation in Fig. 3. Both materials manifest a direct band gap characteristic, observed precisely at the high symmetry Γ-point within the Brillouin zone (BZ). The calculated band gap values for this feature are 1.63 eV for (S-MBA)2SnI4 and 1.56 eV for (S-CYHEA)2SnI4.
![]() | ||
| Fig. 3 Calculated band structure of (a) (S-MBA)2SnI4 and (b) (S-CYHEA)2SnI4 by GGA GGA-PBE without SOC. | ||
To enhance the accuracy of band gap determination while managing computational resources efficiently, we focused our HSE06 functional calculations on the total and partial densities of states (DOS) for both studied materials. The resulting DOS features are presented in Fig. 4a and b, offering a detailed view of the electronic state distributions and a more reliable estimation of the band gap.
The computed band gap values have increased to 2.1 eV and 2.0 eV for (S-MBA)2SnI4 and (S-CYHEA)2SnI4, respectively. Notably, there were no significant alterations observed in the contributions of atomic orbitals at the top of the valence band (VB) and the bottom of the conduction band (CB). However, the calculated band gap value of (S-MBA)2SnI4 aligns well with the values reported in the literature of 2.02 eV.27 Analysis of the DOS (refer to Fig. 4a and b) for (S-MBA)2SnI4 and (S-CYHEA)2SnI4, reveals that in both systems, the occupied band of edge of VB predominantly arises from contributions of I-p states, with a supplementary minor contribution from C p and Sn p states whereas the first unoccupied band is mainly influenced by Sn-p and I p states. In Fig. 4c and d, we present the isosurface of the electron charge density at the boundary between the valence and conduction bands of both 2D layered halide perovskites. Corresponding to the DOS characteristics, the lower of CB of the 2D chiral perovskite mainly consists of Sn p states with a small contribution of I p states, while the higher of VB is mostly composed of I-p, demonstrating excellent concordance with theoretical studies.22,24,27 These findings suggest that both materials show significant promise for electron transport, particularly since the lower of CB involves the p-states of the Sn-site atoms.
Furthermore, our investigation delved into the effective masses of charge carriers in 2D layered halide perovskite systems, shedding light on how the choice of organic spacer cations influences carrier mobilities.41 Specifically, we conducted calculations to determine the electron (me) and hole (mh) effective masses along the Γ direction of the Brillouin zone. For (S-MBA)2SnI4 and (S-CYHEA)2SnI4, we found the electron effective masses to be me = 0.241 m0 and me = 0.195 m0, respectively, while the hole effective masses are mh = 0.251 m0 and mh = 0.206 m0, respectively (m0 is a free electron mass). Moreover, we calculate the corresponding reduced mass of charges carriers using the following equation;
. The estimated reduced masses are found to be 0.123 and 0.100 for (S-MBA)2SnI4 and (S-CYHEA)2SnI4, respectively. The exciton stands as one of the most critical physical entities influencing the performance of optoelectronic and photovoltaic devices. Utilizing the effective-mass approximation and the hydrogen-like model, the exciton binding energy can be computed using the expression41,42;
, where the symbol RH (equal to 13.6 eV) represents the Rydberg constant, which signifies the binding energy of an electron in the ground-state 1s orbital of the hydrogen (H) atom and εr denotes the relative dielectric function. The εr represents the value of the real part of the dielectric function when it reaches its maximum, which is 2.1 for S-(MBA)2SnI4 and 2.07 for (S-CYHEA)2SnI4 (see Fig. S2, ESI†). The estimated exciton binding energy values for the (S-MBA)2SnI4 and (S-CYHEA)2SnI4 are found to be 379 and 317 meV, respectively. It is found that the Eb of (S-MBA)2SnI4 is greater than that of the 2D perovskite (S-CYHEA)2SnI4. Based on the findings, 2D layered perovskites exhibit high exciton binding energies, making them promising candidates for optoelectronic applications.
To investigate the optical absorption characteristics, we computed absorption spectra for (S-MBA)2SnI4 and (S-CYHEA)2SnI4 using the HSE06 functional across the photon energy that varies from 0 to 5 eV, as depicted in the Fig. 5. Interestingly, both 2D layered perovskites exhibit excellent visible light absorption capabilities, with a prominent absorption peak exceeding 105 cm−1 located around 2.9 eV. Notably, the absorption coefficient of (S-CYHEA)2SnI4 in the visible region and ultraviolet surpasses that of (S-MBA)2SnI4 perovskite. These findings underscore the potential application of the studied systems in photovoltaic cells.
![]() | ||
| Fig. 5 Calculated absorption coefficient spectra of (S-MBA)2SnI4 and (b) (S-CYHEA)2SnI4 systems using the HSE06 functional without SOC. | ||
The inclusion of heavy elements such as tin within the crystal structure triggers a significant SOC effect, potentially resulting in the Rashba effect in noncentrosymmetric systems. This effect could render the fundamental bandgap indirect, impacting the behavior of charge carriers. Consequently, it may lead to reduced recombination rates and an extended lifetime for charge carriers.31 To investigate the impact of chirality-induced inversion asymmetry, π-conjugated and the size of organic cations on the electronic band structures of chiral (S-MBA)2SnI4 and (S-CYHEA)2SnI4, we conducted electronic band structure calculations using DFT within GGA-PBE with SOC contribution, as plotted in Fig. 6. The inclusion of SOC significantly underestimates the band gap, yielding values of 1.48 eV and 1.39 eV for (S-MBA)2SnI4 and (S-CYHEA)2SnI4, respectively. This reduction can be attributed to the influence of heavy elements, present in the systems. In the band structure analysis of (S-MBA)2SnI4, the Rashba effect manifests by breaking the spin degeneracy present in both the valence and conduction bands. This leads to the formation of unique parabolic band minima, distinctly observed in the plane the Γ–Z path high-symmetry points of the Brillouin zone.
![]() | ||
| Fig. 6 Band structure of chiral (a) (S-MBA)2SnI4 (b) (S-CYHEA)2SnI4 along the high symmetry directions of Brillouin zone. Light purple and black curves for CB and VB band, respectively. | ||
The expression can be identified in the Rashba Hamiltonian43,44 that describes the linear-in-k space with effective SOC interaction in quasi-two-dimensional (2D) systems:
![]() | (5) |
, where, ER represents the energy splitting, and kR denotes the distance in k-space between the band extremum and the band crossing point. Based on our calculations, for the conduction band, we observe a splitting with an energy ER = 2.6 meV and a corresponding distance in k-space of ER = 0.0265 Å−1, resulting in a Rashba coupling parameter magnitude of 0.196 eV Å. For the valence band, we find an energy splitting of ER = 3 meV and a distance of kR = 0.029 Å−1, leading to a Rashba coupling parameter magnitude of 0.207 eV Å. Notably, (S-MBA)2SnI4 demonstrates a small value of α, and consequently, a minor band splitting.39 In the spacer cation-based system (S-CYHEA)2SnI4, characterized by the absence of π bonds within its molecular framework, an intriguing discovery emerges: the spin splitting values at both the valence and conduction band edges along the Γ–A points exhibit significant enhancements. Remarkably, they are measured at 0.388 eV Å for the VB and 0.339 eV Å for the CB, respectively, showcasing a notable increase when compared to those observed in MBA.
Our findings reveal a fascinating trend: the dramatic increase in spin-splitting observed from (S-MBA)2SnI4 to (S-CYHEA)2SnI4 cannot solely be ascribed to alterations in π-conjugation. Hence, it suggests that besides π-conjugation, another factor may play a significant role in influencing the splitting of chiral hybrid halide perovskites. As it is known, a large helical distortion of the octahedral cages generates a significant Rashba splitting. A comprehensive understanding of helical distortion is gained by comparing the cage distortions of (S-MBA)2SnI4 and S-(CYHEA)2SnI4 with that of (S-NEA = naphthylethylamine)2SnI4.12 As shown in Fig. S3 (ESI†), the cage of (S-NEA)2SnI4 significantly deviates from the ideal 180°, characteristic of an undistorted perovskite layer. This deviation may be attributed to the substantial π-conjugation in NEA cations, leading to strong interactions with neighboring organic cations. Additionally, the cage displays two distinct Sn–I–Sn angles: 149.88° on one side of the inorganic layer and 138.88° on the opposite side, resulting in distortions that propagate helically parallel to the y-axis. In contrast, the perovskites (S-MBA)2SnI4 and (S-CYHEA)2SnI4 demonstrate two very similar Sn–I–Sn angles, yielding nearly flat perovskite layers, as depicted in Fig. S3 (ESI†). Given that a larger helical distortion correlates with increased spin-splitting, these findings confirm that (S-NEA)2SnI4 cation exhibits (αR ≈ 0.707 eV Å) the largest Rashba splitting compared to perovskites based on CYHEA and MBA. However, despite the similar distortion of the octahedral cage between (S-MBA)2SnI4 and (S-CYHEA)2SnI4, the difference in splitting between these two 2D layered perovskites cannot be solely explained by cage distortion. Additionally, we explored another factor that could influence spin splitting in 2D layered perovskites: the d-spacing between inorganic layers. Surprisingly, we found that the estimated d-spacing of the chiral hybrid halide perovskites follows the order of 19.2 Å (S-NEA)2SnI4 > 16.7 Å (S-CYHEA)2SnI4 > 14.4 Å (S-MBA)2SnI4, which aligns with the magnitude of splitting, as depicted in Fig. S4 (ESI†). Remarkably, the observed spin-splitting is also linked to the d-spacing of chiral perovskites, where increasing distortion and d-spacing lead to a larger splitting size.
Footnote |
| † Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4cp02426e |
| This journal is © the Owner Societies 2024 |