Next Article in Journal
Hydrophilic and Hydrophobic Mesoporous Silica Derived from Rice Husk Ash as a Potential Drug Carrier
Next Article in Special Issue
Tetranuclear Oxo-Titanium Clusters with Different Carboxylate Aromatic Ligands: Optical Properties, DFT Calculations, and Photoactivity
Previous Article in Journal
Hybrid Top-Down/Bottom-Up Fabrication of a Highly Uniform and Organized Faceted AlN Nanorod Scaffold
Previous Article in Special Issue
New Insights on the Photodegradation of Caffeine in the Presence of Bio-Based Substances-Magnetic Iron Oxide Hybrid Nanomaterials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

First-Principles Investigation on the Electronic and Mechanical Properties of Cs-Doped CH3NH3PbI3

College of Sciences, Northeastern University, Shenyang 110819, China
*
Author to whom correspondence should be addressed.
Materials 2018, 11(7), 1141; https://doi.org/10.3390/ma11071141
Submission received: 2 June 2018 / Revised: 25 June 2018 / Accepted: 25 June 2018 / Published: 5 July 2018
(This article belongs to the Special Issue Novel Photoactive Materials)

Abstract

:
Methylammonium lead iodide, CH3NH3PbI3, is currently a front-runner as light absorber in hybrid solar cells. Despite the high conversion efficiency, the stability of CH3NH3PbI3 is still a major obstacle for commercialization application. In this work, the geometry, electronic structure, thermodynamic, and mechanical property of pure and Cs-doped CH3NH3PbI3 have been systematically studied by first-principles calculations within the framework of the density functional theory (DFT). Our studies suggest that the (CH3NH3)+ organic group takes a random orientation in perovskite lattice due to the minor difference of orientation energy. However, the local ordered arrangement of CH3NH3+ is energetic favorable, which causes the formation of electronic dipole domain. The band edge states of pure and Cs-doped CH3NH3PbI3 are determined by (PbI6) group, while A-site (CH3NH3)+ or Cs+ influences the structural stability and electronic level through Jahn–Teller effect. It has been demonstrated that a suitable concentration of Cs can enhance both thermodynamic and mechanical stability of CH3NH3PbI3 without deteriorating the conversion efficiency. Accordingly, this work clarifies the nature of electronic and mechanical properties of Cs-doped CH3NH3PbI3, and is conducive to the future design of high efficiency and stable hybrid perovskite photovoltaic materials.

1. Introduction

As one of the biggest scientific breakthroughs [1], perovskite solar cell technology has been extensively investigated over the last few years. The conversion efficiency of organic–inorganic lead halide perovskite solar cells has been impressively improved from 3.8% in 2009 to 22.7% recently [2,3,4,5]. In particular, methylammonium lead halide perovskites (ABX3; A = (CH3NH3)+; B = Pb+; X = Cl, Br, I) are regarded as particularly promising light absorbers because of their outstanding photovoltaic properties. With its continuous increase in conversion efficiency, the stability of CH3NH3PbI3 based solar cells has attracted tremendous attentions for the scale-up of industrial applications [6,7,8]. Enhancing the operational stability of CH3NH3PbI3 without weakening the conversion efficiency is still a major challenge for perovskite-type solar cells [9,10,11].
People attempted to control the composition and proportion of dopants by an alloying method to improve the performance of CH3NH3PbI3-based perovskites. For halogen doping, it has been found that mixing Cl/Br in CH3NH3PbI3 can not only realize the continuous tuning of solar absorption, but also improve the carrier mobility and reduce carrier recombination rates [12,13]. However, due to the thermodynamic instability of the solutions, halogen impurities may cause the segregation of CH3NH3PbI3 into iodide-rich minority and bromide-enriched majority domains under light exposure [14].
Besides the substitution of halogens, B-site doping has also been examined. J. Navas et al. performed experimental and theoretical studies on alloying in Pb2+ sites with Sn2+, Sr2+, Cd2+, and Ca2+. They pointed out that the Sn2+, Sr2+, and Cd2+ did not modify the phase structure [15]. Y. Ogom et al. proved that the Sn/Pb halide-based perovskite solar cells can harvest the light in the area up to 1060 nm. Nevertheless, the instability caused by the B-site impurities would lead to the reduction of open circuit voltage for CH3NH3PbI3 solar cell materials [16].
Recently, the substitution for A-site organic group in CH3NH3PbI3 has drawn lots of attentions [17,18,19,20]. It has been reported that [CH(NH2)2]x(CH3NH3)1−xPbI3 has favorable performances in terms of structural stability, and the band gap of the A-site solutions can be tuned between 1.48 eV and 1.57 eV [21,22,23]. Yi et al. demonstrated that a mixed A-site cation Csx[CH(NH2)2]1−xPbI3 exemplifies the potential of high efficiency solar cell material [24]. [CH(NH2)2]0.85Cs0.15PbI3 solution has shown a better performance and device stability than the plain CH(NH2)2PbI3. Compared with the organic (CH3NH3)+ cation, the inorganic Cs+ is far less volatile [25,26]. At present, it seems that A-site doping is an effective scheme to improve stability without degrading the light conversion efficiency of CH3NH3PbI3.
In this work, we first carefully studied the of orientation influence of (CH3NH3)+ group on the total energy and band gap in CH3NH3PbI3. After the determination of ground state geometry of pure and Cs-doped CH3NH3PbI3, a comprehensive investigation on the electronic and mechanical properties was performed by first principles calculations. The present study is conducive to the design of perovskite solar cell material with high stability and efficiency.

2. Calculation Method and Model

The general chemical formula for perovskites is ABX3. In the cubic unit cell of a perovskite, the A-atom sits at cube corner positions (0, 0, 0) in 12-fold coordination, while the B-atom sits at body center position (1/2, 1/2, 1/2) and is surrounded by six X-atoms to form an octahedron group. A-site traditionally is occupied by the inorganic cations, while in the present hybrid perovskite A-site is inhabited by the organic (CH3NH3)+ group. The geometry optimization and the electronic structure calculations of (CH3NH3)+ group are investigated by Gaussian 98 code [27], and B3LYP [28,29] method in connection with the 6-311++G basis set [30,31] is chosen as the exchange correlation potential.
All the bulk geometry optimization, electronic structure, and mechanical property of pure and Cs-doped CH3NH3PbI3 are calculated with ab initio total energy and molecular dynamics program VASP (VASP 5.4.1, Faculty of Physics, University of Vienna, Austria) [32]. Perdew–Burke–Ernzerhof (PBE) pseudopotential with vdw and spin–orbit coupling (SOC) correction is adopted during the calculation [33]. The plane wave kinetic energy cutoff is set to 550 eV and Brillouin-zone integration is performed with a 12 × 12 × 12 Monkhorst–Pack k-point mesh. The convergence tolerance for the total energy and Hellmann–Feynman force during the structural relaxation is set to 106 eV and 0.01 eV/Å, respectively.

3. Results and Discussion

3.1. Intrinsic Properties of CH3NH3PbI3: Structure and Band

Different locations of (CH3NH3)+ group in the perovskite lattice would introduce variations in energy and electronic structure due to the intrinsic degree of freedom of the group. Therefore, it is appropriate to start with the study of the nature of (CH3NH3)+. The schematic representation of (CH3NH3)+ group is shown in Figure 1a. N and C atoms are all tetrahedral coordinated, and the length of C-N bond in (CH3NH3)+ group is 1.52 Å, the Van Der Waals volume is 81.67 Å3, the density is 0.65 g/cm3, and the effective radius of the group is 2.69 Å. The electronic property calculations show that (CH3NH3)+ holds a 2.54 Debye intrinsic dipole moment with the direction from N atom pointing to C atom. As can be seen in the following, the direction of spontaneous polarization of (CH3NH3)+ will influence the stability and bandgap of CH3NH3PbI3.
After the full relaxation of (CH3NH3)+ group, we construct (CH3NH3)PbI3 unit cell based on the cubic APbI3 inorganic perovskite frame, in which A-site atom is replaced by the (CH3NH3)+ group. Different from inorganic cation, the organic (CH3NH3)+ group performs a nonspherical shape and can be displayed along different directions in the perovskite lattice [34,35]. In order to determine the energetic favorable configuration, the energy–orientation relationship is studied by fixing the midpoint but rotating the direction of the C-N bond at the center of the cubic cell. Based on the symmetry analysis, there are two independent rotating modes, namely rotating along [110] and [100] directions, for (CH3NH3)+ group and the rotation period of each mode is π/2. The rotating angle is set to zero as the C-N bond lies in the (ab) plane.
As shown in Figure 2a, the lowest energies configuration (at 0 K) appears at the position of ±20° when rotating along [100] direction. However, the maximum rotation barrier caused by the orientation is about 40 meV, which is close to the thermal energy perturbation at room temperature (26 meV). As a result, (CH3NH3)+ group will hold a random orientation in the perovskite lattice at RT. Thus, the calculation results explain the disorder arrangement of the (CH3NH3)+group observed in experiments [36,37].
Although the tiny difference in total energy supports the random orientation, we found that the local order location of (CH3NH3)+ group is energetic favorable in configuration. Since G.R. Berdiyorov et al. reported that 48 atoms contained supercell is large enough to negligible the finite size effects [38]. We built supercells containing four unit cells and the orientation of (CH3NH3)+ group in the adjacent cells is arranged as +−+−, ++−−, and ++++. It is found that the all-aligned ++++ mode is the most stable configuration and ~110 meV lower in total energy than the cross +−+− mode. The local ordered location of (CH3NH3)+ is conducive to the establishment of electric dipole domain in the perovskite lattice, which facilitates the separation of photo-generated carriers and leads to the promotion of photoelectric conversion efficiency. This interesting phenomenon has also been investigated by molecular dynamics simulation and experiment [39,40], and the size of the domains is found to be about 100 nm [39,40].
Furthermore, the variation of bandgap with the group orientation is evaluated (As seen in Figure 2b). According to the results, the fluctuation of bandgap caused by the asymmetry of electronic distribution and the electric dipole moment in (CH3NH3)+ group is 0.1 eV. It is interesting that the lowest total energy and the minimum bandgap appear at the same point, namely at the position of ±20° when rotating along [100] direction. Since that the bandgap shrinks with the intensity of electronic hybridization, it is reasonable to deduce that the system has the lowest energy when the electronic bonding is the strongest.
The point with the lowest energy on the energy–orientation curve corresponds to the ground state of CH3NH3PbI3 (Figure 1b). After geometry relaxations, the band structure of ground state in CH3NH3PbI3, drawn between high symmetry points of the Brillouin zone, has been illustrated in Figure 3. CH3NH3PbI3 is found to have a direct bandgap of 1.68 eV at Γ. This result is well consistent with the reported experimental values [41,42,43]. The partial electronic densities of state (PDOS) indicates that the upper valence bands (VB) mainly consist of I-5p orbital with weak hybridization to Pb-5s state, and the lower conduction bands (CB) are dominated by Pb-6p orbital. While the (CH3NH3)+ group has little contribution to the band edge states. The character of band structure in this work is consistent with recent investigations [44,45,46]. As a result, the band gap of CH3NH3PbI3 will be strongly related to the structure of (PbI6) inorganic framework. We can thus infer that the influence of (CH3NH3)+ group on the electronic structure belongs to Jahn–Teller effect, that is, the change of electronic structure of the system originates from the distortion of (PbI6) octahedral caused by the size and orientation of (CH3NH3)+group, rather than from the direct participation in the electronic structure of (CH3NH3)+. According to our “anion group model” [47], this feature of the electronic structure helps to improve the stability through the element substitution without losing light response and conversion capability.

3.2. Cs-Doped CH3NH3PbI3: Stability and Electronic Properties

The structure stability of ABX3 perovskites can be estimated by Goldschmidt rule [48], and the tolerance factor t of ABX3 is determined by the expression
t = R A + R X 2 ( R B + R X )
where RA, RB, RX are the respective effective ionic radii of A, B, and X ions. In general, the perovskite structure is considered to be highly stable when the t-factor is between the range of 0.90~1.00 [49]. Although the above rule has been developed for the oxide perovskites, but the criterion is still valid for the structural stability analysis of inorganic–organic hybrid halide perovskite materials [50,51].
To determine the range of doping concentration of Cs, the Goldschmidt’s tolerance factors of Csx(CH3NH3)1−xPbI3 (x from 0 to 1 at an interval of 0.125) are calculated. In our cases, the effective radii of (CH3NH3)+, Pb2+, Cs+, and I in the perovskite-type lattice is 2.69 Å, 1.19 Å, 1.88 Å, and 2.20 Å, respectively. As a result, the t-factor is 1.02 for CH3NH3PbI3 and 0.85 for CsPbI3. According the criterion of Goldschmidt, the stability of CH3NH3PbI3 and CsPbI3 is unsatisfactory. For the Cs doped CH3NH3PbI3, we define the average effective ionic radius of A-site cations as RA-eff = x × RCs + (1 − x) RCNH. When the concentration x increases from 0 to 1 at an interval of 0.125, the corresponding t-factor at each x is 1.0, 0.99, 0.96, 0.92, 0.91, 0.89, and 0.87 for the Cs-doped CH3NH3PbI3. Thus, from the view of Goldschmidt rule, the Cs-concentration of 12.5 at.% to 62.5 at.% is more desirable. In experiments, it has been found that Cs-doped CH3NH3PbI3 can be successfully synthesized with a high solid solubility [52].
To evaluate the function of Cs-dopant on the structural stability, the 12.5 at.% Cs-doped CH3NH3PbI3 has been modeled by a 2 × 2 × 2 CH3NH3PbI3 supercell with one (CH3NH3)+ group replaced by a Cs atom. In this model we do not consider the disordered orientation of (CH3NH3)+ group. The formation energy of CH3NH3PbI3 and Cs0.125(CH3NH3)0.875PbI3 are calculated according to the formula
E f o r m a t i o n = ( E t o t a l j n j E i o n j ) / N t o t a l
where Etotal is the total energy of Cs-doped CH3NH3PbI3, E i o n j is the energy of constituent elements in their respective elemental state, nj is the number of various constituent elements, and Ntotal is the total number of atoms in the supercell. The results show that the formation energy of Cs0.125(CH3NH3)0.875PbI3 is 0.02 eV/atom lower than that of the pure CH3NH3PbI3, indicating that incorporating Cs into CH3NH3PbI3 lattice is exothermic and helps to stabilize CH3NH3PbI3. More interestingly, it is found that the original regular (PbI6) octahedral chain has been distorted in the Cs-doped CH3NH3PbI3, as shown in Figure 4. Thus the symmetry breaking of perfect octahedral chain would also lead to the energy reduction.
As the incorporation of Cs will induce a 2.75 Debye net electric dipole moment, we introduce an electric dipole correction term in the electronic structure calculation of Cs0.125(CH3NH3)0.875PbI3. In Figure 5, energy band diagram clearly shows that the bandgap of 12.5 at.% Cs-doped CH3NH3PbI3 is 1.73 eV, being ~3% wider than that of the pure CH3NH3PbI3. It can be seen that the Cs dopant does not introduce new states at the band edges. This result is consistent with the “anoin group model” in that A-site alkali metal or alkaline earth metal cation in perovskites can hardly influence the electronic state near the Fermi level [47], and this is why the alkali metal doping does not cause the deterioration of photovoltaic performance of CH3NH3PbI3. As a result, A-site participation of alkali metal or alkaline earth metal is a promising way to stabilize the hybrid perovskite photovoltaic materials with reliable photovoltaic properties.

3.3. Cs-Doped CH3NH3PbI3: Mechanical Properties

It is believed that the mechanical stability of a material is strongly related to its equilibrium elastic modulus and strength performance [53]. The bulk modulus is a measure of resistance to the volume change due to applied pressure, while the shear strength provides information about the resistance of a material against plastic deformation. Thus, the stability of Cs-doped CH3NH3PbI3 should be basically understood by these parameters. To obtain elastic moduli, elastic constants were first calculated from the stress–strain relation, then the Voigt–Reuss–Hill (VRH) approximation was applied to the CH3NH3PbI3 system, and the effective elastic moduli could be approximated by the arithmetic mean of the Voigt and Reuss limits [54].
In Figure 6, our calculations show that the bulk modulus (B) of pure CH3NH3PbI3 is 11.0 GPa, and The shear modulus (G) is 4.9 GPa, which is comparable to the reported experimental value [55,56]. When the Cs-concentration is below 37.5 at.%, both bulk modulus and shear modulus keep going up and the top value of them is higher than that of pure CH3NH3PbI3 by 3.5 GPa and 0.6 GPa, respectively. When the Cs-concentration increases further, the corresponding elastic modulus decreases with fluctuation. Our calculations suggest that the Cs-concentration should be controlled in the range of 20~35 at.% if we want to achieve the optimal equilibrium mechanical performance. Due to the limitation of solid solubility, the actual Cs-concentration may be lower than this range. However, the elastic moduli still are higher than that of the pure substance, as shown in Figure 6.
According to the work of Pugh, the material is deemed to be ductile if its B/G ratio is greater than the critical value of 1.75 [57]. The B/G ratio of CH3NH3PbI3 is 2.24, while Cs-doped (CH3NH3)PbI3 have higher B/G values indicating that those solutions have a better performance in ductility.
The high value of elastic moduli near the equilibrium does not guarantee the high strength. The fracture feature of Cs-doped (CH3NH3)PbI3 should be basically understood from the ideal shear strength obtained far from the equilibrium. The ideal shear stress–strain curves of (CH3NH3)PbI3 solutions in the (111)<11-2> typical slip system are shown in Figure 7.
The strain of deformation-to-failure ε indicates the maximum deformation that the material can bear and thus is the measure of the brittleness. It is found that CsPbI3 has the highest shear strength (0.78 GPa) but the worst strain capacity, showing a feature of high strength and high brittleness. On the country, (CH3NH3)PbI3 possess the lowest shear strength (0.52 GPa) but the best strain capacity, displaying a feature of low strength and low brittleness. Below 37.5% Cs-concentration, the fracture performance of solutions lies between the boundaries determined by the pure substances CsPbI3 and (CH3NH3)PbI3, that is, the ideal shear strength increases but the maximum bearing deformation decreases. It should be noted that the current strength calculations are strictly related to the defect free lattice, and the obtained ideal might overestimate real shear strengths. Along with the analysis on the equilibrium elastic moduli, we can conclude that below 37.5% Cs-concentration (CH3NH3)PbI3 will possess a desirable performance in mechanical properties including stability, hardness, strength, and ductility.

4. Conclusions

Aiming at developing new efficient and stable perovskite solar cell materials, we report a comprehensive theoretical investigation on the geometry, electronic structure, and mechanical property of pure and A-site Cs-doped CH3NH3PbI3. The main conclusions are drawn as follows:
(1)
The difference in orientation energy of (CH3NH3)+ is comparable to the thermal power at room temperature, which causes a random orientation of (CH3NH3)+ group in the perovskite lattice.
(2)
The local ordered arrangement of (CH3NH3)+ is energetic favorable that facilitates the formation of the electronic dipole domain, which helps to improve the separation and lifetime of photo-generated carriers.
(3)
The band edge states are dominated by (PbI6) anion group in CH3NH3PbI3. A-site (CH3NH3)+ or Cs+ does not directly participate in the construction of the band edge states, but indirectly influences the structural stability and electronic level through Jahn–Teller effect.
(4)
It has been demonstrated that the suitable concentration of Cs can enhance both thermodynamic and mechanical stability of CH3NH3PbI3 without deteriorating the conversion efficiency.
(5)
Goldschmidt’s tolerance factor suggests that the Cs-concentration should be less than 62.5 at.%, while mechanical performance indicates that the optimal Cs-concentration should be less than 37.5%. Below this mark, the mechanical properties including stability, hardness, strength, and ductility will continuously rise with the Cs-concentration.
The adopted research methods, mechanism cognitions and obtained conclusion in the work might help contributing to the future development of efficient and stable hybrid perovskite solar cell materials.

Author Contributions

X.M. and D.L. built the model. D.L., S.L. and B.F. calculated properties of the compounds. All authors contributed to data analysis. X.M. and D.L. wrote the manuscript with input from all authors. X.M. directed the project.

Funding

This work is financially supported by the National Key Research and Development Program of China 2016YFB0701100, the Fundamental Research Funds for the Central Universities (N150502002, N160208001, and N160206002), and National Natural Science Foundation of China (no. 51525101).

Acknowledgments

Dongyan Liu sincerely thanks Haijun Pan at Northeastern University at Qinhuangdao for helpful discussion.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Science News. Newcomer Juices up the race to harness sunlight. Science 2013, 342, 1438–1439. [Google Scholar]
  2. Kojima, A.; Teshima, K.; Shirai, Y.; Miyasaka, T. Organometal halide perovskites as visible-light sensitizers for photovoltaic cells. J. Am. Chem. Soc. 2009, 131, 6050–6051. [Google Scholar] [CrossRef] [PubMed]
  3. Lee, M.M.; Teuscher, J.; Miyasaka, T.; Murakami, T.N.; Snaith, H.J. Efficient hybrid solar cells based on meso-superstructured organometal halide perovskites. Science 2012, 338, 643–647. [Google Scholar] [CrossRef] [PubMed]
  4. Liu, M.Z.; Johnston, M.B.; Snaith, H.J. Efficient planar heterojunction perovskite solar cells by vapour deposition. Nature 2013, 501, 395–398. [Google Scholar] [CrossRef] [PubMed]
  5. Best Research-Cell Efficiencies (National Renewable Energy Laboratory, 2018). Available online: https://www.nrel.gov/pv/assets/images/efficiency-chart.png (accessed on 26 June 2018).
  6. Conings, B.; Drijkoningen, J.; Gauquelin, N.; Babayigit, A.; D’Haen, J.; D’Olieslaeger, L.; Ethirajan, A.; Verbeeck, J.; Manca, J.; Mosconi, E. Intrinsic thermal instability of methylammonium lead trihalide perovskite. Adv. Energy Mater. 2015, 5. [Google Scholar] [CrossRef]
  7. Christians, J.A.; Herrera, P.A.M.; Kamat, P.V. Transformation of the excited state and photovoltaic efficiency of CH3NH3PbI3 perovskite upon controlled exposure to humidified air. J. Am. Chem. Soc. 2015, 137, 1530–1538. [Google Scholar] [CrossRef] [PubMed]
  8. Graetzel, M.; Janssen, R.A.J.; Mitzi, D.B.; Sargent, E.H. Materials interface engineering for solution-processed photovoltaics. Nature 2012, 488, 304–312. [Google Scholar] [CrossRef] [PubMed]
  9. Pearson, A.J.; Eperon, G.E.; Hopkinson, P.E.; Habisreutinger, S.N.; Wang, J.T.; Snaith, H.J.; Greenham, N.C. Oxygen Degradation in Mesoporous Al2O3/CH3NH3PbI3−xClx Perovskite Solar Cells: Kinetics and Mechanisms. Adv. Energy Mater. 2016, 6. [Google Scholar] [CrossRef]
  10. Alsari, M.; Pearson, A.J.; Wang, J.T.W.; Wang, Z.; Montisci, A.; Greenham, N.C.; Snaith, H.J.; Lilliu, S.; Friend, R.H. Degradation Kinetics of Inverted Perovskite Solar Cells. Sci. Rep. 2018, 8, 5977. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Alsari, M.; Bikondoa, O.; Bishop, J.; Abdi-Jalebi, M.; Ozer, L.Y.; Hampton, M.; Thompson, P.; Hörantner, M.T.; Mahesh, S.; Greenland, C.; et al. In situ simultaneous photovoltaic and structural evolution of perovskite solar cells during film formation. Energy Environ. Sci. 2018, 11, 383–393. [Google Scholar] [CrossRef] [Green Version]
  12. Mosconi, E.; Quarti, C.; Ivanovska, T.; Ruani, G.; Angelis, F.D. Structural and electronic properties of organo-halide lead perovskites: A combined IR-spectroscopy and ab initio molecular dynamics investigation. Phys. Chem. Chem. Phys. 2014, 16, 16137–16144. [Google Scholar] [CrossRef] [PubMed]
  13. Noh, J.H.; Im, S.H.; Heo, J.H.; Mandal, T.N.; Seok, S.I. Chemical management for colorful, efficient, and stable inorganic–organic hybrid nanostructured solar cells. Nano Lett. 2013, 13, 1764–1769. [Google Scholar] [CrossRef] [PubMed]
  14. Hoke, E.T.; Slotcavage, D.J.; Dohner, E.R.; Bowring, A.R.; Karunadas, H.I.; McGehee, M.D. Reversible photo-induced trap formation in mixed-halide hybrid perovskites for photovoltaics. Chem. Sci. 2015, 6, 613–617. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Navas, J.; Sánchez-Coronilla, A.; Gallardo, J.J.; Hernández, N.C.; Piñero, J.C.; Alcántara, R.; Lorenzo, C.F.; De los Santos, D.M.; Aguilar, T.; Calleja, J.M. New insights into organic–inorganic hybrid perovskite CH3NH3PbI3 nanoparticles. An experimental and theoretical study of doping in Pb2+sites with Sn2+, Sr2+, Cd2+ and Ca2+. Nanoscale 2015, 7, 6216–6229. [Google Scholar] [CrossRef] [PubMed]
  16. Ogomi, Y.; Morita, A.; Tsukamoto, S.; Saitho, T.; Fujikawa, N.; Shen, Q.; Toyoda, T.; Yoshino, K.; Pandey, S.; Ma, T.; et al. CH3NH3SnxPb(1−x)I3 perovskite solar cells covering up to 1060 nm. J. Phys. Chem. Lett. 2014, 5, 1004–1011. [Google Scholar] [CrossRef] [PubMed]
  17. Eperon, G.E.; Paternò, G.M.; Sutton, R.J.; Zampetti, A.; Haghighirad, A.A.; Cacialli, F.; Snaith, H.J. Inorganic caesium lead iodide perovskite solar cells. J. Mater. Chem. A 2015, 3, 19688–19695. [Google Scholar] [CrossRef]
  18. Ahmad, M.; Rehman, G.; Ali, L.; Shafiq, M.; Iqbal, R.; Ahmad, R.; Khan, T.; Asadabad, S.J.; Maqbool, M.; Ahmad, I. Structural, electronic and optical properties of CsPbX3 (X = Cl, Br, I) for energy storage and hybrid solar cell applications. J. Alloy Compd. 2017, 705, 828–839. [Google Scholar] [CrossRef]
  19. Protesescu, L.; Yakunin, S.; Bodnarchuk, M.I.; Krieg, F.; Caputo, R.; Hendon, C.H.; Yang, R.X.; Walsh, A.; Kovalenko, M.V. Nanocrystals of cesium lead halide perovskites (CsPbX3, X = Cl, Br, and I): Novel optoelectronic materials showing bright emission with wide color gamut. Nano Lett. 2015, 15, 3692–3696. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Trots, D.M.; Myagkota, S.V. High-temperature structural evolution of caesium and rubidium triiodoplumbates. J. Phys. Chem. Solids 2008, 69, 2520–2526. [Google Scholar] [CrossRef] [Green Version]
  21. Eperon, G.E.; Beck, C.E.; Snaith, H.J. Cation exchange for thin film lead iodide perovskite interconversion. Mater. Horiz. 2016, 3, 63–71. [Google Scholar] [CrossRef]
  22. Jeon, N.J.; Noh, J.H.; Yang, W.S.; Kim, Y.C.; Ryu, S.; Seo, J.; Seok, S., II. Compositional engineering of perovskite materials for high-performance solar cells. Nature 2015, 517, 476–480. [Google Scholar] [CrossRef] [PubMed]
  23. Binek, A.; Hanusch, F.C.; Docampo, P.; Bein, T. Stabilization of the trigonal high-temperature phase of formamidinium lead iodide. J. Phys. Chem. Lett. 2015, 6, 1249–1253. [Google Scholar] [CrossRef] [PubMed]
  24. Yi, C.; Luo, J.; Meloni, S.; Boziki, A.; Astani, N.A.; Gratzel, C.; Zakeeruddin, S.M.; Rothlisberger, U.; Gratzel, M. Entropic stabilization of mixed A-cation ABX3 metal halide perovskites for high performance perovskite solar cells. Energy Environ. Sci. 2016, 9, 656–662. [Google Scholar] [CrossRef]
  25. Li, Z.; Yang, M.; Park, J.S.; Wei, S.H.; Berry, J.J.; Zhu, K. Stabilizing perovskite structures by tuning tolerance factor: Form ation of formamidinium and cesium lead iodide solid-state alloys. Chem. Mater. 2016, 28, 284–292. [Google Scholar] [CrossRef]
  26. Kulbak, M.; Cahen, D.; Hodes, G. How important is the organic part of lead halide perovskite photovoltaic cells? Efficient CsPbBr3 cells. J. Phys. Chem. Lett. 2015, 6, 2452–2456. [Google Scholar] [CrossRef] [PubMed]
  27. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Zakrzewski, V.G.; Montgomery, J.A.; Stratmann, R.E.; Burant, J.C.; et al. Gaussian 98; Gaussian Inc.: Pittsburgh, PA, USA, 1998. [Google Scholar]
  28. Becke, A.D. Density functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648–5652. [Google Scholar] [CrossRef]
  29. Lee, C.; Yang, W.; Parr, R.G. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 1988, 37, 785. [Google Scholar] [CrossRef]
  30. McLean, A.D.; Chandler, G.S. Contracted Gaussian basis sets for molecular calculations. I. Second row atoms, Z = 11–18. J. Chem. Phys. 1980, 72, 5639–5648. [Google Scholar] [CrossRef]
  31. Perdew, J.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865. [Google Scholar] [CrossRef] [PubMed]
  32. Kresse, G.; Furthmuller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169. [Google Scholar] [CrossRef]
  33. Krishnan, R.; Binkley, J.S.; Seeger, R.; Pople, J.A. Self-consistent molecular orbital methods. XX. A basis set for correlated wave functions. J. Chem. Phys. 1980, 72, 650–654. [Google Scholar] [CrossRef]
  34. Bakulin, A.A.; Selig, O.; Bakker, H.J.; Rezus, Y.L.; Muller, C.; Glaser, T.; Lovrincic, R.; Sun, Z.; Chen, Z.; Walsh, A.; et al. Real-time observation of organic cation reorientation in methylammonium lead iodide perovskites. J. Phys. Chem. Lett. 2015, 6, 3663–3669. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Frost, J.M.; Butler, K.T.; Brivio, F.; Hendon, C.H.; Schilfgaarde, M.V.; Walsh, A. Atomistic origins of high-performance in hybrid halide perovskite solar cells. Nano Lett. 2014, 14, 2584–2590. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. La-o-vorakiat, C.; Salim, T.; Kadro, J.; Khuc, M.; Haselsberger, R.; Cheng, L.; Xia, H.; Gurzadyan, G.G.; Su, H.; Lam, Y.M.; et al. Elucidating the role of disorder and free-carrier recombination kinetics in CH3NH3PbI3 perovskite films. Nat. Commun. 2015, 6, 7903. [Google Scholar] [CrossRef] [PubMed]
  37. Stroppa, A.; Quarti, C.; De Angelis, F.; Picozzi, S. Ferroelectric polarization of CH3NH3PbI3: A detailed study based on density functional theory and symmetry mode analysis. J. Phys. Chem. Lett. 2015, 6, 2223–2231. [Google Scholar] [CrossRef] [PubMed]
  38. Berdiyorov, G.R.; Madjet, M.E.; El-Mellouhi, F.; Peeters, F.M. Effect of crystal structure on the electronic transport properties of the organometallic perovskite CH3NH3PbI3. Sol. Energy Mater. Sol. Cells 2016, 148, 60–66. [Google Scholar] [CrossRef]
  39. Frost, J.M.; Walsh, A. What is moving in hybrid halide Perovskite solar cells? Acc. Chem. Res. 2016, 49, 528–535. [Google Scholar] [CrossRef] [PubMed]
  40. Kutes, Y.; Ye, L.; Zhou, Y.; Pang, S.; Huey, B.D.; Padture, N.P. Direct observation of ferroelectric domains in solution-processed CH3NH3PbI3 Perovskite Thin Films. J. Phys. Chem. Lett. 2014, 5, 3335–3339. [Google Scholar] [CrossRef] [PubMed]
  41. Jeng, J.Y.; Chiang, Y.F.; Lee, M.H.; Peng, S.R.; Guo, T.F.; Chen, P.; Wen, T.C. CH3NH3PbI3 perovskite/fullerene planar heterojunction hybrid solar cells. Adv. Mater. 2013, 25, 3727–3732. [Google Scholar] [CrossRef] [PubMed]
  42. Qiu, J.; Qiu, Y.; Yan, K.; Zhong, M.; Mu, C.; Yan, H.; Yang, S. All-solid-state hybrid solar cells based on a new organometal halideperovskite sensitizer and one-dimensional TiO2 nanowire arrays. Nanoscale 2013, 5, 3245–3248. [Google Scholar] [CrossRef] [PubMed]
  43. Umari, P.; Mosconi, E.; Filippo, D. Relativistic GW calculations on CH3NH3PbI3 and CH3NH3SnI3 perovskites for solar cell applications. Sci. Rep. 2014, 4, 4467. [Google Scholar] [CrossRef] [PubMed]
  44. Mosconi, E.; Amat, A.; Nazeeruddin, M.; Gratzel, M.; Angelis, F.D. First-principles modeling of mixed halide organometal perovskites for photovoltaic applications. J. Phys. Chem. C 2013, 117, 13902–13913. [Google Scholar] [CrossRef]
  45. Umebayashi, T.; Asai, K.; Kondo, T.; Nakao, A. Electronic structures of lead iodide based low-dimensional crystals. Phys. Rev. B 2003, 67, 155405. [Google Scholar] [CrossRef]
  46. Brivio, F.; Walker, A.B.; Walsh, A. Structural and electronic properties of hybrid perovskites for high-efficiency thin-film photovoltaics from first-principles. APL Mater. 2013, 1, 042111. [Google Scholar] [CrossRef] [Green Version]
  47. Meng, X.Y.; Liu, D.Y.; Qin, G.W. Band engineering of multicomponent semiconductors: A general theoretical model on the anion group. Energy Environ. Sci. 2018, 11, 692–701. [Google Scholar] [CrossRef]
  48. Goldschmidt, V.M. Crystal structure and chemical correlation. Ber. Dtsch. Chem. Ges. 1927, 60, 1263–1268. [Google Scholar] [CrossRef]
  49. Goldschmidt, V.M. Die Gesetze der Krystallochemie. Naturwissenschaften 1926, 14, 477–485. [Google Scholar] [CrossRef]
  50. Kieslich, G.; Sun, S.; Cheetham, A.K. Solid-state principles applied to organic–inorganic perovskites: New tricks for an old dog. Chem. Sci. 2014, 5, 4712–4715. [Google Scholar] [CrossRef]
  51. Stoumpos, C.C.; Kanatzidis, M.G. The renaissance of halide perovskites and their evolution as emerging semiconductors. Acc. Chem. Res. 2015, 48, 2791–2802. [Google Scholar] [CrossRef] [PubMed]
  52. Niemann, R.G.; Gouda, L.; Hu, J.; Tirosh, S.; Gottesman, R.; Cameron, P.J.; Zaban, A. Cs+ incorporation into CH3NH3PbI3perovskite: Substitution limit and stability enhancement. J. Mater. Chem. A 2016, 4, 17819–17827. [Google Scholar] [CrossRef]
  53. Born, M.; Huang, K. Dynamical Theory of Crystal Lattices; Clarendon Press: Oxford, UK, 1954. [Google Scholar]
  54. Hill, R. The elastic behaviour of a crystalline aggregate. Proc. Phys. Soc. A 1952, 65, 349–354. [Google Scholar] [CrossRef]
  55. Beecher, A.N.; Semonin, O.E.; Skelton, J.M.; Frost, J.M.; Terban, M.W.; Zhai, H.W.; Alatas, A.; Owen, J.S.; Walsh, A.; Billinge, S.J.L. Direct observation of dynamic symmetry breaking above room temperature in methylammonium lead iodide perovskite. ACS Energy Lett. 2016, 1, 880–887. [Google Scholar] [CrossRef]
  56. Rakita, Y.; Cohen, S.R.; Kedem, N.K.; Hodes, G.; Cahen, D. Mechanical properties of APbX3 (A = Cs or CH3NH3; X = I or Br) perovskite single crystals. MRS Commun. 2015, 5, 623–629. [Google Scholar] [CrossRef]
  57. Pugh, S.F. Relations between the elastic moduli and the plastic properties of polycrystalline pure metals. Philos. Mag. Ser. 1954, 45, 823–843. [Google Scholar] [CrossRef]
Figure 1. Structure model of (a) (CH3NH3)+ group (the arrow represents the direction of electric dipole moment in the group), and (b) CH3NH3PbI3 ground state (the PbI6 octahedral is rendered).
Figure 1. Structure model of (a) (CH3NH3)+ group (the arrow represents the direction of electric dipole moment in the group), and (b) CH3NH3PbI3 ground state (the PbI6 octahedral is rendered).
Materials 11 01141 g001
Figure 2. The fluctuation of total energy and bandgap with the orientation of (CH3NH3)+ group in the perovskite lattice.
Figure 2. The fluctuation of total energy and bandgap with the orientation of (CH3NH3)+ group in the perovskite lattice.
Materials 11 01141 g002
Figure 3. Band level and PDOS of CH3NH3PbI3.
Figure 3. Band level and PDOS of CH3NH3PbI3.
Materials 11 01141 g003
Figure 4. Regular (PbI6) octahedral chain in CH3NH3PbI3 (a) and distorted chain in 12.5 at.% Cs-doped CH3NH3PbI3 (b).
Figure 4. Regular (PbI6) octahedral chain in CH3NH3PbI3 (a) and distorted chain in 12.5 at.% Cs-doped CH3NH3PbI3 (b).
Materials 11 01141 g004
Figure 5. Band level and PDOS of Cs0.125(CH3NH3)0.875PbI3.
Figure 5. Band level and PDOS of Cs0.125(CH3NH3)0.875PbI3.
Materials 11 01141 g005
Figure 6. The variation of bulk modulus and shear modulus with Cs-concentration in CH3NH3PbI3.
Figure 6. The variation of bulk modulus and shear modulus with Cs-concentration in CH3NH3PbI3.
Materials 11 01141 g006
Figure 7. DFT calculated ideal shear stress–strain curves for the (CH3NH3)PbI3 solutions.
Figure 7. DFT calculated ideal shear stress–strain curves for the (CH3NH3)PbI3 solutions.
Materials 11 01141 g007

Share and Cite

MDPI and ACS Style

Liu, D.; Li, S.; Bian, F.; Meng, X. First-Principles Investigation on the Electronic and Mechanical Properties of Cs-Doped CH3NH3PbI3. Materials 2018, 11, 1141. https://doi.org/10.3390/ma11071141

AMA Style

Liu D, Li S, Bian F, Meng X. First-Principles Investigation on the Electronic and Mechanical Properties of Cs-Doped CH3NH3PbI3. Materials. 2018; 11(7):1141. https://doi.org/10.3390/ma11071141

Chicago/Turabian Style

Liu, Dongyan, Shanshan Li, Fang Bian, and Xiangying Meng. 2018. "First-Principles Investigation on the Electronic and Mechanical Properties of Cs-Doped CH3NH3PbI3" Materials 11, no. 7: 1141. https://doi.org/10.3390/ma11071141

APA Style

Liu, D., Li, S., Bian, F., & Meng, X. (2018). First-Principles Investigation on the Electronic and Mechanical Properties of Cs-Doped CH3NH3PbI3. Materials, 11(7), 1141. https://doi.org/10.3390/ma11071141

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop