Dominic W.
Ferdani‡
a,
Samuel R.
Pering‡
a,
Dibyajyoti
Ghosh
ab,
Peter
Kubiak
a,
Alison B.
Walker
b,
Simon E.
Lewis
a,
Andrew L.
Johnson
a,
Peter J.
Baker
c,
M. Saiful
Islam
*a and
Petra J.
Cameron
*a
aDepartment of Chemistry, University of Bath, Bath BA2 7AY, UK. E-mail: m.s.islam@bath.ac.uk; p.j.cameron@bath.ac.uk
bDepartment of Physics, University of Bath, Bath BA2 7AY, UK
cISIS Neutron & Muon Source, STFC Rutherford Appleton Laboratory, Didcot OX11 0QX, UK
First published on 15th May 2019
Halide perovskite solar cells containing a mixture of A-site cations are attracting considerable interest due to their improved stability and high power conversion efficiencies. Ionic transport is known to be an important predictor of perovskite behaviour, but the impact of partial A-site substitution on iodide ion diffusion is poorly understood. Here, we combine ab initio modelling, impedance spectroscopy and muon spin relaxation to investigate the effect on iodide ion transport of incorporating a low concentration of each of seven different sized cations (from small rubidium to large guanidinium) into methylammonium lead iodide. Experimental and simulation results are in good agreement, indicating that these cation substitutions increase the activation energy for iodide ion diffusion. We show for the first time that partial guanidinium substitution into methylammonium lead iodide strongly suppresses iodide ion transport. The insights gained from this multi-technique study are important for the future design of mixed-cation perovskite solar cells with enhanced performance.
Broader contextPerovskite solar cells (PSCs) have been generating headlines over the last five years as their solar conversion efficiency has increased rapidly to rival that of multi-crystalline silicon solar cells. PSCs have several advantages over conventional silicon solar cells, including ease of manufacture and projected low energy pay back times. The cells are named for the ‘perovskite’ crystal structure of the light absorbing layer, which is based on an organic–inorganic iodide material. The perovskite layer is a mixed electronic–ionic conductor, containing mobile iodide ions which can strongly influence the measured properties of the solar cell. In this study we show that by moving from a single organic cation perovskite to a mixed-cation perovskite, iodide ion movement is supressed. We have systematically studied seven different cations and shown that the activation energy for iodide ion motion is linked to the size of the cation that is being substituted into the structure. Our results help to understand the improved stability of mixed cation PSCs. |
Several larger A-site cations, including azetidinium (Az) and guanidinium (GA), were first identified by structural factors19 and the mixed materials were subsequently prepared.8,11 MA:Az materials, which retain a perovskite structure up to 5% Az substitution, show enhanced stability and better fill factors compared to MAPbI3.8 Guanidinium (GA) substitution improves the PSC performance by increasing the open cell voltage (VOC) and greatly enhancing stability,20 although the effect on ion transport has not been investigated. Jodlowski et al. showed that GA is phase miscible with MAPbI3 at MA:
GA ratios of up to 75
:
25, and MA0.86GA0.14PbI3 cells can reach efficiencies of over 20%.21
Given the enhanced stability and high efficiencies, what is not fully understood is why mixing A-site cations improves the cell performance. It is widely accepted that iodide ions are mobile in MAPbI3, but the impact of partial A-site cation substitution on iodide ion transport is unclear.22–24
Here, we combine experimental and simulation techniques to investigate, for the first time, the effect on iodide ion diffusion of substituting low concentrations of each of seven different sized cations into MAPbI3 (listed in Table 1). The cation sizes were varied from rubidium and caesium to guanidinium (GA). GA has an ionic radius of 278 pm, indicating significant MA/GA ion size mismatch. All of the cations were included at 5% substitution, giving MA0.95A0.05PbI3. As a result, the compounds retained the same 3D perovskite structure and band gap as MAPbI3, which allowed direct comparison of our results. Ab initio modelling was used to probe the atomic-scale effects of cation substitution on the local structure and ion transport activation energies. Photovoltaic devices were prepared for each mixed-cation composition and investigated by impedance spectroscopy (IS). Finally, powders of cation-substituted perovskites were investigated using muon spin relaxation (μSR) and we show that we can directly detect iodide ion diffusion in MAPbI3 and the doped compounds.
Cation | Structure | Radius/pm19,25,26 |
---|---|---|
Rubidium (Rb) |
![]() |
152 |
Caesium (Cs) |
![]() |
167 |
Methylammonium (MA) |
![]() |
217 |
Azetidinium (Az) |
![]() |
250 |
Formamidinium (FA) |
![]() |
253 |
Dimethylammonium (DM) |
![]() |
272 |
Acetamidinium (Ac) |
![]() |
277 |
Guanidinium (GA) |
![]() |
278 |
Our combined studies show strong agreement between the ion transport results obtained from μSR on perovskite powders, ab initio simulations on the atomic scale and IS on full PSC devices. All the cation substitutions in MAPbI3 increase the activation energy for iodide ion transport. In particular, we show that just 5% GA substitution strongly suppresses iodide ion transport in the halide perovskite, which has implications for solar cell performance.
X-ray diffraction (XRD) was used to confirm that all seven doped materials had the same perovskite crystal structure and that phase separation had not occurred (Fig. 1). UV/vis measurements confirmed that there was minimal difference in the band gap after substitution, with changes of ±0.05 eV relative to MAPbI3 (Fig. S1a, ESI†). These results are supported by our ab initio simulations which also indicate that at 5% cation substitution there is no major effect on the band gap (Table S1, ESI†).
The presence of 5 mol% of the substituted cation was further confirmed by 1H NMR analysis on powders (for Cs, FA and GA substituted MAPbI3; Fig. S2, ESI†). By comparing peak integrals, we show that FA is present at a substitution ratio of 4.35 mol%; GA at 5 mol% and Cs at 4.3 mol%.
The thermodynamic stability of the perovskite phase upon cation substitution was explored using ab initio simulations by calculating the enthalpies of formation of the mixed-cation compound from the component iodide salts and PbI2 (Fig. 1b). As reported previously, pure MAPbI3 shows a very small positive formation enthalpy indicating a degree of metastability.29 More importantly, all the mixed-cation materials exhibit a negative formation enthalpy with GA substitution found to be the most favourable. These results suggest that the thermodynamic stability of MAPbI3 is enhanced by A-site cation substitution, which would suppress the decomposition to the component iodides (for example, MAI, GAI and PbI2). Also, the large cation size of GA (and Ac) promotes hydrogen bonding with iodide anions and these stronger interactions dominate over local strain to help stabilise the structure. Such findings are in good agreement with experimental observations of higher stability against temperature and chemical degradation for mixed-cation perovskites.16
There have been previous ion transport studies on single cation perovskites,22–27 but not the effects of partial substitution of A-site cations as studied here; for example, in terms of comparing single cation perovskites Meloni et al. examined iodide migration in MAPbI3vs. CsPbI3, whereas Oranskaia et al. investigated bromide migration in MAPbBr3vs. FAPbBr3.47 A range of activation energies (Ea) for iodide ion conduction in pure MAPbI3 can be found in the literature (summarised in Table 2);22–24,30–36 most values fall between 0.3 eV and 0.6 eV, with the highest activation energy derived from DFT using the cubic structure.22
Technique | E a (eV) |
---|---|
a Performed on MAPbI3−xClx perovskites. | |
Ab initio/DFT & chronoamperometry22 | 0.58/0.62 |
Impedance and IMVS30![]() |
0.55 |
Thermally stimulated current31 | 0.5 |
Temperature dependent capacitance32![]() |
0.45 |
Ab initio/DFT33 | 0.44 |
Impedance24 | 0.43 |
Temperature dependent JV curves23 | 0.33 |
Temperature dependent current34![]() |
0.31 |
Solid-state NMR35 | 0.17 |
Classical MD, ab initio/DFT36,37 | 0.08 |
However, the precise effects on ion transport of partial A-site substitutions into MAPbI3 have not been characterised and hence remains poorly understood. Here, we extend our related ab initio simulation work39,48,49 to calculate the relative activation energies for the vacancy-mediated migration of iodide ions in three representative systems: pure tetragonal MAPbI3 and Cs- and GA-substituted compounds. Since the simulations are probing local energy barriers, we investigated the MA0.75A0.25PbI3 composition with a higher substitution level allowing us to better examine A-site cation effects.
Energy profiles were derived for iodide vacancy diffusion, and particular care was taken to ensure full structural convergence to avoid unusually low activation energies. Our results are reported in Fig. 2, revealing three key results. First, the activation energy of 0.44 eV for tetragonal MAPbI3 matches well with previous experimental and computational studies (Table 2). Second, there is a small increase in the activation energy for Cs-substituted MAPbI3. Most significantly, we find a relatively steep increase to 0.78 eV for the GA-substituted material (with no major change with different GA orientations; see ESI,† Section C(iv) and Fig. S9). Hence, for long-range iodide ion diffusion through the perovskite material the rate-limiting barrier will be 0.78 eV, which predicts that iodide ion transport would be suppressed.
The simulations also reveal a curved migration path between iodide sites within the Pb/I face for both MAPbI3 and Cs-substituted material, with a small displacement of the neighbouring Pb ion (Fig. 2). In contrast, the presence of the much larger GA cation causes significant steric constraints and distortion of the local inorganic framework leading to a large displacement (0.6 Å) of the neighbouring Pb ion. These structural changes now lead to a curved path that is bowed out of the Pb/I plane (Fig. 2c), resulting in a much higher activation energy in the GA-substituted perovskite. Hence, it is the significant local distortion of the Pb/I cage due to the GA cation size mismatch with MA that leads to the reduction in iodide transport, rather than the exact orientation of the substituted molecular cation. Such atomic-scale mechanistic features are difficult to extract from experiment alone.
To gain insights into relative changes in diffusion rates, we derived iodide ion diffusion coefficients (DI) from the calculated activation energies, assuming a Boltzmann-like barrier-hopping process with a typical attempt frequency of 1012Hz for ionic species at a temperature of 300
K. We find that DI would be approximately one order and five orders of magnitude lower for the Cs- and GA-substituted materials respectively; these results predict the suppression of iodide ion transport with GA substitution, which accords well with the muon spin relaxation and IS results described below.
Temperature-dependent impedance spectroscopy (IS) was carried out on the devices. IS was measured at open circuit under illumination at temperatures between 45 °C and −15 °C. In line with our previous results, three semi-circles were observed in the Nyquist plot for MAPbI3 (Fig. S13, ESI†). We describe these as the high frequency (HF), mid-frequency (MF) and low frequency (LF) responses. The interpretation of impedance spectra for PSCs is still hotly debated in the literature.51–56 Nevertheless, it is generally accepted that the HF response gives the recombination resistance, Rrecomb, and the geometric capacitance, Cgeo. The LF feature is observed between 50 and 0.03 Hz (the time constant, τ, is typically just below 1 s) and the time scale is consistent with slow ionic reorganisation. We have previously attributed the LF response to modulation of Rrecomb as a function of a changing ionic environment. This would occur if, for example, the modulating voltage caused ions to move in such a way to better screen photo-generated charges, which in turn would increase the recombination resistance.
Regardless of the model used to explain the LF response, numerous studies show that it is intrinsically linked to ion diffusion. Bag et al. and Pockett et al. show that an activation energy can be found by plotting the LF time constant as a function of temperature.30,56,57 The activation energy of 0.40 eV for MAPbI3 (Fig. 3b and Table S4, ESI†) is very close to our calculated value and to those measured for iodide ions by a wide range of other techniques (Table 2). It is extremely unlikely that motion of A site cations such as MA+ will be observed on this timescale as their diffusion coefficients are at least three orders of magnitude lower than that for iodide ions.22
Fig. 3 shows representative Nyquist plots for MA0.95Az0.05PbI3 and MA0.95GA0.05PbI3 at 25 °C (see also Fig. S13 and S14, ESI†). Three semi-circles are clearly visible in the Nyquist plot for MA0.95Az0.05PbI3. In contrast, the LF feature disappears completely for MA0.95GA0.05PbI3 (Fig. S14, ESI†) and an activation energy for iodide diffusion can no longer be extracted. These results mirror our ab initio modelling and μSR work and suggest that GA substitution strongly suppresses iodide ion diffusion. Similarly, it was not possible to resolve the LF feature for the second largest cation, Ac. Hence, there are no IS-derived activation energies for either GA or Ac, since ion transport was too slow to be detected.
Activation energies were extracted from the LF impedance response for cells substituted with the other five cations, Rb, Cs, Az, FA, DM (Fig. S15, ESI†), and plotted as a function of ionic radius in Fig. 3(b). The trend with ion size matches well with our calculated activation energies (Fig. 2). The IS results clearly show that partial substitution of MA with smaller or larger cations increases the activation energy for iodide ion transport.
All of the inverted cells prepared for this study showed low hysteresis at the scan rates studied (Fig. S11, ESI†). Despite not being the focus of the present study, iodide ion diffusion in perovskite solar cells has important implications for JV hysteresis, although it is known that hysteresis is also influenced by the nature of the electron- or hole-transport layers. In our study we did not see a reduction in JV-curve hysteresis after cation substitution, but it is important to note that scans were measured at a single scan rate and at room temperature. This is clearly an area that warrants further investigation.
μSR was used to probe ion movement in pure MAPbI3 and 5 mol% guanidinium substituted crystallites. To first investigate whether iodide motion could be detected by μSR, experiments were carried out on MAPbI3 and d6-perdeuterated MAPbI3 powders. Measurements on deuterated and protonated cations helped us to separate cation from anion motion. Data were collected between 50 and 410 K with typical raw data shown in Fig. 4a and Fig. S20–S22 (ESI†). The stability of the powders was checked using three separate techniques (ESI,† Section E(ii)). Fig. S19c (ESI†) shows XRD diffractograms of MAPbI3 before and after muon spectroscopy measurements were carried out. Some very small changes in peak intensity were observed, which is unsurprising as the sample has been subjected to a large change in temperature (40 K to 400 K) over 30 hours. Nevertheless, no evidence of PbI2 formation was observed.
In μSR measurements muons are implanted into the sample where they decay into positrons, which are most likely to be emitted in the muon spin direction at the instant of decay. The effect of local fields on the material on the muon spins is detected by asymmetry in the positron counts measured by detectors around the sample. The asymmetry data at each temperature was fitted to a dynamic Gaussian Kubo–Toyabe function multiplied by an exponential relaxation (λ) (details in ESI,† Section E(i)).60 This allowed calculation of the fluctuation rate (ν, due to nuclear magnetic moments flipping the muon spin as they move nearby) and the static width of the disordered local field distribution (Δ, from nuclear magnetic moments at muon implantation sites).
Fig. 4c shows the temperature dependence of ν for MAPbI3. Δ values are shown in Fig. S23 (ESI†). ν increases from 0.05 MHz at 100 K to 0.6 MHz at 150 K, before falling to 0.2 MHz at 160 K. After 160 K, ν slowly increases until above 350 K, where ν increases more quickly in a thermally activated manner. The two large increases in ν are due to the motion of nuclear moments relative to the muon as different processes occur within the range of timescales to which muon spins are sensitive. The activation energies associated with these two processes in MAPbI3 were derived from an Arrhenius fit of the two activated regions (Fig. S24a and b, ESI†), producing values of 0.172 eV and 0.072 eV for the high-temperature and low-temperature processes respectively.
We have assigned the low temperature process to molecular cation reorientation within the Pb/I cage. As shown previously,61,62 MAPbI3 and d6-MAPbI3 are in a less rigid orthorhombic structure at 100 K with the methylammonium cations fully ordered. The {NH3} groups of MA+ align into the square faces of the unit cell and the {PbI6} octahedra are distorted. MAPbI3 undergoes a phase transition from orthorhombic to tetragonal at about 160 K.62–64 At temperatures close to 100 K, ν starts to increase as the muons are affected by increasing molecular cation motion within the Pb/I cages. Above 150 K the fluctuation rate decreases as the perovskite undergoes a phase transition, above the phase transition cation motion is too fast to be detected by μSR. Our results match neutron scattering studies showing that MA cations reorient within the perovskite lattice at temperatures close to 140 K.63
The high-temperature activation energies of 0.174 eV in MAPI and 0.158 eV in d6-MAPI are within the margin of experimental error (Table 3). Although lower than the DFT and IS derived energies, the muon result agrees well with the Ea value for iodide ion diffusion in MAPbI3 obtained from NMR studies (Table 2). NMR and μSR both probe changes in local fields within the bulk of a material so it would be expected that the two methods produce similar results.35
Perovskite | Low T Ea (eV) | High T Ea (eV) |
---|---|---|
MAPbI3 | 0.072 (±0.005) | 0.174 (±0.010) |
d6-MAPbI3 | 0.082 (±0.003) | 0.158 (±0.019) |
MA0.95GA0.05PbI3 | 0.051 (±0.005) | Not observed |
The hopping rates in all of the samples are indistinguishable between about 150 K and 320 K as the iodide motion is too slow to be measurable by μSR.64 It is only when the diffusion coefficient rises above 1 × 10−13 cm2 s−1 (at about 360 K) that the ionic motion is on the correct timescale for μSR to measure. It is important to stress that this does not mean that the iodide ion is not moving at room temperature, but simply that it is not being detected by μSR.
Our muon study was then extended to MAPbI3 substituted with 5 mol% GA. Fig. 4b shows the raw data collected at 400 K. The retention of the same perovskite crystal structure was confirmed by powder XRD (Fig. S26, ESI†). Our computational predictions show that GA substitution increases the activation energy for iodide ion migration suggesting that iodide ion transport would be supressed. The μSR results support the computational and IS results as after GA substitution iodide motion can no longer be measured in the material. In MA0.95GA0.05PbI3 the increase in fluctuation rate above room temperature disappears completely and no iodide ion diffusion is detected up to 420 K (Fig. 4c). This result is also good evidence that the high-temperature growth in ν for the other samples is due to iodide ion transport, rather than muon motion, which we discuss fully in the ESI,† Section D. The presence of GA also affects the low-temperature fluctuations of the MA cations; the low-temperature Ea are listed in Table 3. Overall, the μSR results agree with our ab initio simulation and impedance spectroscopy results, clearly showing that adding GA to MAPbI3 reduces iodide ion transport.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9ee00476a |
‡ These authors contributed equally to the manuscript. |
This journal is © The Royal Society of Chemistry 2019 |