Ethylammonium as an alternative cation for efficient perovskite solar cells from first-principles calculations.

Mixed-cation lead halide perovskites have emerged as a new class of promising photovoltaic materials for perovskite solar cells. Formamidinium (FA), methylammonium (MA), and Cs cations are widely studied in the field of mixed-cation hybrid halide perovskites. In this work, we have investigated ethylammonium (CH3CH2NH3, EA) as an alternative cation to explore the stabilities and electronic properties of mixed MA1-x EA x PbI3 perovskites. The results indicate that replacing MA with EA is a more effective way to improve the stabilities of the mixed MA1-x EA x PbI3 perovskites except for MA0.75EA0.25PbI3. The band gap of MA1-x EA x PbI3 slightly increases with x from 0.25 to 1.00, which is quite different from the MA-FA mixed-cation perovskites. The results indicate that the c axis distortion of the Pb-I-Pb bond angles can play a greater role in tuning the band gap. Moreover, the mixed MA1-x EA x PbI3 perovskites show comparable absorption abilities in the visible light region to the pure MAPbI3 structure. We hope that our study will be greatly helpful for further experiments to find more efficient perovskite materials in the future. This journal is © The Royal Society of Chemistry.

Introduction

Hybrid organic–inorganic lead halide perovskites with outstanding photovoltaic performance have been identified as desirable light absorption materials for perovskite solar cells.[1-3] The power conversion efficiency of perovskite solar cells has been improved from an intial 3.8%[4] to 23.7%.[5] In a typical ABX3 perovskite structure, A is a monovalent cation such as methylammonium (CH3NH3+, MA), formamidinium (CH(NH2)2+, FA), and Cs+. B is a metal ion (usually Pb2+, Sn2+, or Ge2+) and X is a halide anion (usually I−, Br−, or Cl−). Recently, composition engineering on ABX3 perovskites has been carried out on the cation and halide. This has proved to be a powerful approach to tune the stability and band gap.[6-11] The most popular cations for composition engineering are Cs, MA, and FA.[12-16] These works have demonstrated that mixed cations have advantages in improving the stabilities and optoelectronic properties of hybrid perovskites. Therefore, there is a strong desire to explore new stable perovskites with suitable photovoltaic performances.

In general, the mixed-cation composition engineering can be evaluated by the Goldschmidt tolerance factor (t), which is defined as follows:[17]where RA and RB are the ionic radii of A and B cations respectively, and RX is the ionic radius of the anion. The results suggest that the value of the tolerance factor t could be between 0.81 and 1.11.[18] To maintain a cubic structure, the ideal t should be close to 1. When t > 1.11, the A-site cation is too large and usually precludes the formation of a perovskite. When t < 0.81, the A-site cation is too small, and this often leads to other structures.

The radius of the ethylammonium cation (CH3CH2NH3, EA) is about 2.3 Å,[19] which is greater than that of MA (1.8 Å) but smaller than that of FA (2.6 Å).[20] Based on the eqn (1), the tolerance factor of MA1−EAPbI3 ranges from 0.83 to 0.94, as shown in Fig. 1. The tolerance factors of the MA1−EAPbI3 perovskites are consistent with the relevant conclusions.[18] This indicates that incorporating EA into MAPbI3 is a effective way to form stable mixed-cation MA1−EAPbI3 perovskites. However, the EAPbI3 structure was synthesized in 2012.[21] Its structure was described as a 2H perovskite type with a larger band gap of 2.2 eV.[21,22] The incorporation of EA cations with larger effective radius to improve material stability has been investigated experimentally.[19,23,24] A small amount of EA cations (x ≤ 0.3) have been successfully incorporated into MAPbI3 to form the stable mixed MA1−EAPbI3 perovskites.[25] Recently, a new 3D perovskite with a chemical structure of (1,3-Pr(NH3)2)0.5PbI3 was successfully synthesized based on a larger organic cation of propane-1,3-diammonium (1,3-Pr(NH3)22+).[26] The results showed that this new perovskite exhibited a tetragonal structure and enhanced stability.[26] These results indicate that EA as a alternative cation can form a stable 3D perovskite structure under proper experimental conditions. In this work, we have investigated the structural stabilities and electronic properties of MA–EA mixed-cation MA1−EAPbI3 perovskites, whose electronic properties are substantially different from the MA–FA mixed-cation perovskites.

Fig. 1

Correlations between the tolerance factor and the structures of the MA1−EAPbI3 perovskites.

Computational details

All of the first-principles calculations were performed with the Vienna Ab initio Simulation Package (VASP) based on DFT.[27] The projector-augmented wave method was used for the electron-ion interactions and the Perdew–Burke–Ernzerhof (PBE)[28] generalized gradient approximation (GGA)[29] was used to describe the exchange–correlation effect of the electrons. The cut-off energy for the plane wavefunctions was set to 550 eV. A Monkhorst–Pack 6 × 6 × 6 k-point mesh over the first Brillouin zone was used for structural optimization and a denser 8 × 8 × 8 k-point mesh was used for the electronic and optical properties. The total energy was converged to 10−5 eV. The atom coordinates were fully optimized until the residual forces were smaller than 0.01 eV Å−1. It is well known that van der Waals interactions play an important role in hybrid organic–inorganic halide perovskite materials.[30,31] Thus, the DFT-D3 functional, which includes dispersion interactions, was used for the structural optimization.[32]

It is well known that the PBE functional always underestimated the band gaps of semiconductors. In the previous DFT calculation, our calculated result shows that the band gap of MAPbI3 is 1.50 eV with PBE functional,[33] which is quite close experimental value of 1.55 eV.[34] Thus, the PBE functional could give the reasonable band gaps for hybrid lead iodine perovskites. Besides, the calculated optical properties of cation-doped MAPbI3 with PBE functional also can show the right results. Note that DFT-D3 method was adopted to account for van der Waals interactions in all calculations.

To evaluate the stabilities of the pure and mixed perovskites, we calculated the formation energies of the pure and EA-doped MAPbI3 perovskites. The crystal structures of MAI was taken from the ref. 35 The crystal structure of EAI was obtained from monoclinic (C2/m) CH3CH2NH2I[36] with full relaxation of the structural parameters. The stable structure of hexagonal PbI2 was used. The energies of PbI2, MAI, and EAI were calculated using the same computational parameters.

Results and discussion

Geometric structures

The simulated mixed MA1−EAPbI3 perovskites are performed by gradually reducing the percentage of MA from the starting I4cm tetragonal phase. The simulated unit cell contains four MAPbI3 units, so we can simulate three values of the MA–EA ratio, namely 3 : 1 (x = 0.25), 2 : 2 (x = 0.50), 1 : 3 (x = 0.75), and the pure phase of x = 0 and 1. For x = 0.25, 0.50, and 0.75, all of the possible substitution sites are considered. The optimized structures of the mixed MA1−EAPbI3 perovskites are shown in Fig. 2. For the structures with x = 0.50 and x = 0.75, all the possible substitution sites are approximately equivalent with the energy differences within 0.02 eV. However, the results show that different structures with x = 0.25 have different stabilities, with calculated maximum energy differences within 0.08 eV. The most stable structures of MA1−EAPbI3 with x = 0.25, 0.50, and 0.75 are chosen to investigate their structural stabilities, electronic and optical properties.

Fig. 2

Optimized stable geometries of the mixed MA1−EAPbI3 perovskites.

The optimized lattice constants of the mixed MA1−EAPbI3 perovskites are listed in Table 1. In the tetragonal MAPbI3 structure, the relaxed lattice constants a is 8.76 Å and c is 12.95 Å, which is good agreement with the experimental results.[37] The volume slightly increases with increasing EA percentage because of large-size cation doping.

Calculated lattice parameters for the MA1−EAPbI3 perovskites

x	a/Å	b/Å	c/Å	V/Å3
0.00	8.76	8.76	12.95	993.74
0.25	8.92	8.70	12.97	1005.76
0.50	9.05	8.67	12.98	1018.52
0.75	9.00	8.81	13.02	1031.69
1.00	8.97	8.94	13.04	1045.63

To further evaluate of the stabilities of the pure and mixed perovskites, we calculated the formation energies of the pure and EA-doped perovskites using the following equation:where E(MA1−EAPbI3), E(MAI), E(EAI), and E(PbI2) refer to the formation energies of MA1−EAPbI3, MAI, EAI, and PbI2 phases, respectively. The calculated results for the mixed MA1−EAPbI3 perovskites are shown in Fig. 3. For MAPbI3, our calculated formation energy is well consistent with the previous calculated values of −0.01 eV per f.u.[38] and −0.02 eV per f.u.[39] All of the ΔH values are negative for the MA1−EAPbI3 perovskites. This indicates that substituting MA with EA is energetically favourable. Interestingly, the stability of MA0.75EA0.25PbI3 is much lower than that of MAPbI3. The similar results were also found for the mixed MA1−AzPbI3 perovskites.[33] The stabilities of the mixed MA1−EAPbI3 perovskites are apparently improved when the EA percentage increases from 0.50 to 1.00.

Fig. 3

Calculated formation energies of the MA1−EAPbI3 perovskites.

Electronic properties

To obtain insight into the effect of EA-doped cation on the electronic properties, we performed band structure calculations for all of the structures. Fig. 4 shows the MA1−EAPbI3 band structures with x = 0.00, 0.25, 0.50, 0.75, 1.00. Notably, regardless of the EA percentage, the band structures of perovskites possess direct bandgaps. The band gap of MAPbI3 at the PBE level is 1.50 eV, which agrees well with a previous theoretical result of 1.53 eV (ref. 40) and the experimental value of 1.55 eV.[34] The band gaps are 1.57, 1.61, 1.64, 1.65 eV for MA1−EAPbI3 with x = 0.25, 0.50, 0.75, and 1.00, respectively. However, the band gap of EAPbI3 is much lower than its experimental result (2.2 eV).[21] In fact, an orthorhombic structure in the space group of Pmmn was found for EAPbI3.[21] EAPbI3 can form a stable 3D structure according to its radius. It indicates that the radius of EA may be measured. The band gaps of the mixed MA1−EAPbI3 perovskites exhibit a slightly increase with EA concentration, which is in good agreement with the experimental results.[25] However, the results show that the EA incorporation is significantly different from the FA incorporation although the two cations have similar radii. The (1,3-Pr(NH3)2)0.5PbI3 perovskite also showed a larger band gap (1.62 eV) than that of the MAPbI3 structure.[26] These suitable band gaps of the mixed MA1−EAPbI3 perovskites, from 1.50 to 1.65 eV, makes them highly promising for optoelectronic applications.

Fig. 4

Calculated band structures of the MA1−EAPbI3 perovskites.

To understand the band structure modification, we calculated the projected density of states (PDOS) of the mixed MA1−EAPbI3 perovskites. Fig. 5 shows the dominant PDOS near the band edgeds of MA1−EAPbI3 with x = 0.00, 0.25, 0.50, 0.75, and 1.00. The results indicate that the valence band maximum (VBM) of MA1−EAPbI3 originates mainly from the 5p orbital of I. The contribution of the conduction band minimum (CBM) is mainly composed of the 6p orbital of Pb. The organic cation does not contribute to the band edge states because the states of organic cation are far from the Fermi level.

Fig. 5

DOS structures of (a) MAPbI3, (b) MA0.75EA0.25PbI3, (c) MA0.50EA0.50PbI3, (d) MA0.25EA0.75PbI3, and (e) EAPbI3.

Previous report has indicated that the electronic structures of hybrid perovskites are affected by structural distortion, and the main factor is the Pb–I–Pb bond angle.[41] The Pb–I–Pb bond angles of the MA1−EAPbI3 perovskites are shown in Table 2. The larger distortion of the Pb–I–Pb bond angles is from ab plane. The Pb–I–Pb bond angle from ab plane has little change except for EAPbI3. MA0.75EA0.25PbI3 gives rise to an increased band gap of 1.57 eV compared with the pure MAPbI3 structure. This results indicate that the c axis distortion can play a greater role in tuning the band gap. It can be seen that the change of the Pb–I–Pb bond angles is responsible for the difference in band gaps of the MA1−EAPbI3 perovskites.

Pb–I–Pb bond angles and band gaps of the MA1−EAPbI3 perovskites

x	Pb–I–Pb bond angles (°)		E g (eV)
	In ab plane	Along c axis
0.00	150.0	175.6	1.50
0.25	149.9	174.2	1.57
0.50	149.6	171.2	1.61
0.75	151.2	171.0	1.64
1.00	155.8	171.5	1.65

To get more insight into the electronic properties, we next examine the Bader charge analysis for the pure and mixed perovskites. The results of Bader charge analysis are shown in Table 3. Charge transfer of organic cation and I ion gradually decrease with increasing x. However, charge transfer of Pb cation gradually increases with increasing x.

Bader charge analysis of the MA1−EAPbI3 perovskites (A represents the organic cation)

x	A	Pb	I
0.00	+0.727	+0.941	−0.556
0.25	+0.719	+0.944	−0.554
0.50	+0.714	+0.946	−0.553
0.75	+0.710	+0.947	−0.552
1.00	+0.707	+0.949	−0.552

Optical properties

To evaluate the impact of the EA incorporation on the optical performance of MA–EA mixed-cation perovskites, we have calculated the optical absorption coefficient of each structure. The absorption coefficient I(ω) was given as below:[42]where ε1(ω) and ε2(ω) represent the real and imaginary parts of the dielectric function depending on the light frequency ω. The absorption spectra of the perovskite compounds are shown in Fig. 6. The optical absorption of MAPbI3 agrees well with the previous result.[40] It can be seen that the MA1−EAPbI3 perovskites have slightly weaker absorption than that of MAPbI3 in the range of 300–400 nm, which are attributed to their band gaps. The MA1−EAPbI3 perovskites show the comparable absorption abilities in the visible light region compared with the pure MAPbI3 perovskite. Based on the stability of material, the MA1−EAPbI3 (0.50 ≤ x ≤ 1.00) perovskites can become the potential candidates for efficient perovskite solar cells.

Fig. 6

Calculated optical absorption spectra of the MA1−EAPbI3 perovskites.

Conclusions

In summary, first-principles calculations based on DFT were carried out on structural stabilities, electronic and optical properties of the mixed MA1−EAPbI3 perovskites. The results indicate that incorporating EA into MAPbI3 to form mixed-cation MA1−EAPbI3 perovskite is a potentially effective way to improve its stability. The band gap of MA1−EAPbI3 slightly increases with x from 0.25 to 1.00, which is quite different from the MA–FA mixed-cation perovskites. The results indicate that the c axis distortion of the Pb–I–Pb bond angles can play a greater role in tuning the band gap. The mixed MA1−EAPbI3 perovskites show the comparable absorption abilities in the visible light region compared with the MAPbI3 structure. Based on the above results, we hope that our study will be helpful for further experiments to find more efficient perovskite solar cells.

Conflicts of interest

There are no conflicts to declare.