Abstract

Recent experimental results suggest that higher mobility of perovskite-type ligand passivated PbS quantum dots (QDs) could be useful for efficient solar cell applications. However, theoretical understanding of the mechanism through first principal modeling is still lacking. In this study, electronic-, optical-, and temperature-dependent carrier mobility for perovskite ligand passivated PbS QD array is calculated by using the first-principles density functional theory (DFT) combined with the nonequilibrium Green’s function (NEGF) technique and a molecular dynamics (MD)-Landauer approach. It is found that formamidinium (FA)-liganded QDs have higher mobility and enhanced optical absorption comparing to that of Cl-liganded QDs. The difference could be understood through the intermediate band featured electronic structure.

1. Introduction

Harvesting energy directly from solar photons has been a central topic in the quest for sustainable energy. Perovskite-type liganded quantum dots (QDs) solar cells have been investigated extensively for next-generation photovoltaic applications [1, 2]. QD-based intermediate band optoelectronic materials have important effects on the photovoltaic device efficiency [35]. On one hand, introducing the intermediate band through QD nanostructure could introduce extra channel for photon absorption to improve the energy conversion efficiency through the increase of the short-circuit current [4]. On the other hand, the change of mobility gives competing effect at short-circuit condition and open-circuit condition due to charge recombination [6]. Understanding the interplay between optical and transport properties regarding the interface effect and ligand effect provides insight in the trade-off for the design of low-dimensional photoelectronic devices [79]. Recent studies have experimentally reported higher mobility of perovskite-type liganded PbS QDs [7]. However, theoretical principle regarding the enhancement of mobility has not been revealed, making it difficult to further systemically design and tailor the device performance.

Carrier transport property of formamidinium (FA) lead halide perovskite crystal was previously studied for solar cell applications, but the QD array was not considered [10]. Organic ligand for metal halide perovskite solar cells was discussed and potential application to QDs was conjectured [11]. The application of perovskites nanocrystal instead of QDs was also proposed [12]. Optical properties of colloidal PbS QDs were calculated and experimentally measured, but carrier transport property of QD array was not included in these reports [13, 14]. Carrier transport properties of QD arrays were theoretically studied under various idealized situations, which are different from the focus of this study. For example, Green’s functions were calculated for noninteracting systems or small interacting systems, but realistic phonon simulations were not included in these studies [15, 16]. The electronic, optical, and transport properties of PbS QD array are calculated using semiclassical Boltzmann equation for carrier transport [17]. The and drift-diffusion methods were conducted for transportation through QD array junctions [18]. The midgap state effect on the optoelectronic properties of different types of ligand passivated PbS colloidal QDs is studied experimentally with density functional theory (DFT) simulations [19]. Our work is based on molecular dynamics (MD)-Landauer approach, which is different from all the methods mentioned above and will be given in full description in the subsequent section [2022].

In this work, numerical simulation of the electronic properties, optical absorption, as well as temperature-dependent mobility of organic perovskite-type liganded PbS QDs is performed while taking the halogen (here Cl) liganded PbS QDs as a reference. It is observed in simulations that the carrier mobility of FA-liganded QDs is higher than that of Cl-liganded QDs for a wide range of temperature. Further, calculations of electronic structure and optical absorption suggest that the existence or absence of intermediate band plays a role in the change of transport property. Intuitively, the intermediate band provides extra channel of relaxation through scattering, and, hence, could hinder the carrier transport. Our simulation results provide evidence for this intuition and a guide for future development of perovskite-liganded QD solar cells.

This paper is organized as follows. The details of calculation methods and parameters for simulations are listed in Section 2. The main results and discussion are presented in Section 3. The conclusion is given in Section 5.

2. Simulation Methods

2.1. Electronic Structure and Optical Simulation

The geometry of the liganded QDs is optimized through QuantumATK software [23]. In QuantumATK software, global structural optimization is implemented through a genetic algorithm. It works by generating an initial set of random configurations and then evolving them using genetic operators [24]. The electronic structure and optical properties were calculated by using density functional theory-linear combination of atomic orbitals (DFT-LCAO) method embedded in QuantumATK software [25]. The simulation details and initial parameters setting are listed in Table S1.

2.2. MD-Landauer Method for Simulating Mobility

Combining molecular dynamics simulation with nonequilibrium Green’s function (NEGF) for quantum transport gives the so-called MD-Landauer method [2022]. QuantumATK software is used for the calculations [23]. Electron–phonon coupling (EPC) plays a key role in the performance of most electronic devices. However, rigorous simulation of EPC is both challenging for the approaches either based on the first principle or Boltzmann transport equation and for device based on both bulk materials and nano- or atomic-scale materials [2628]. One of the most concerning drawbacks of these conventional approaches, such as the EPC in the Boltzmann equation and the perturbative NEGF methods, is the harmonic approximation of phonons, since introducing anharmonic effects would dramatically increases the computational cost. However, the anharmonic contributions to the phonons are strong at room temperatures and above for many materials, which could not be deal properly with these conventional approaches [22]. In this framework, the electron–phonon interaction is taken into account through the MD calculation, while the carrier transport properties are calculated by DFT-NEGF method and Landauer formula. Figure 1 shows the sketched principle for the proposed methods. For device structure with fixed left and right electrodes, the central transportation region was subjected to MD simulation under different temperature and corresponding NEGF to derive the conductance based on transmittance . Figures 1(a) and 1(b) show the device structure with short and long length of the MD region. Figure 1(c) illustrates the simulation flowchart of MD-Landauer method. Started with a specified length of the MD region, the equilibrium lattice structure is obtained for the central region at the specified temperature . DFT-NEGF method is then applied to calculate the length-dependent transmittance , averaged over the MD samples based on the following Landauer formula to simulate the conductance :

where is the elementary charge, is Planck’s constant, and is the Fermi-level energy. The Landauer formula is used to obtain the conductance and resistance based on the following equations:

As shown in Figure 1(c), this process was repeated for different length of MD region. Here, we have performed four different lengths indexed with and the numbers 0, 3, 4, and 5 represent the number of liganded QDs located between the two electrodes. With the resistance under different length of MD region, the 1 D resistivity is then to be extracted from linear regression over samples of different lengths using the conventional least square regression algorithm.where is defined as random number from normal distribution and is the number of sample data, here . The procedure is repeatedly performed for MD central region with different lengths, so that the resistivity can be obtained by linear regression based on the Formula (4). Using the results of , the mobility can readily be calculated as follows:where is the carrier density. is the bulk resistivity and is calculated with a predetermined device cross-section area . The mobility under specified temperature is obtained and the whole procedure is repeatedly conducted for other temperatures, as shown in Figure 1(c). The parameters used in our MD-Landauer simulations are listed in Table S2 and will be given in more detailed explanation in next section.

3. Results and Discussion

3.1. The Structures of Simulated Systems

The electronic and optical properties of PbS QD arrays without ligand, with FA ligand, and with Cl ligand are studied and compared. Ten ligands are added to one QD for both FA and Cl cases. The structures of the QDs are depicted in Figure 2(a). The geometries of the QDs for both FA ligand and Cl ligand are optimized with DFT-LCAO method from a randomized initial structure, where the detail of the optimization algorithm is presented in the Method section. The QDs are arranged on a square lattice of unit cell with periodic boundary. The distances between QDs are set as 1.7 nm, which has been verified as the maximum distance that shows a coupling effect between QDs in the band structure (see Figure S1 for the results of the QD distance larger than 1.7 nm in Appendix for more details). The superlattice for these two ligands are depicted in Figure 2(b).

3.2. Electronic and Optical Properties of QD with Ligands

Band structure, density of states (DOS), and optical absorption coefficient are calculated using the DFT package embedded in QuantumATK software. Figure 3(a)3(c) shows the band structure of QD arrays with different ligands. For better comparison, we also presented the band structure simulation results for the PbS QDs without ligand, as shown in Figure 3(a). A well-defined Fermi level is found located at the middle of the band gap indicating the intrinsic semiconductor feature of the proposed structures. It is observed, as shown in Figures 3(a) and 3(c), that the intermediate bands located between the energy range of 0.5 and 1.0 eV exist for both pure PbS QD array and and Cl-liganded QD array. However, striking difference was found for the band structure of FA-liganded QD array, as shown in Figure 3(b), where there is no feature indicating the formation of intermediate band and instead band continuum was formed with the lowest band edge at the energy around 0.6 eV. It can also be verified from the DFT band structure simulation that FA-liganded PbS QD array possesses an energy bandgap of 1.28 eV, which is about 0.37 eV narrower than energy bandgap of 1.65 eV for the PbS QD array without ligand. Different from the FA-liganded, the Cl-liganded PbS QD array, which has bandgap energy of 1.48 eV and the intermediate band in the range between 0.5 and 1.0 eV, more closely resembles the electronic band structure to the ones of PbS QD without ligand.

To further understand these phenomena, we investigated the DOS and optical absorption profile. Figure 4(a)4(f) shows the calculated DOS, optical absorption, and orbital-projected PDOS for the FA-liganded and Cl-liganded QDs. It has been reported that PbS QD excitonic peak could smear out due to hybridization of orbitals [29]. In our DOS simulation data, a separated conduction band peak near the Fermi level for Cl-liganded QD is observed (highlighted by red-dashed circle, as shown in Figure 4(d)). This peak does not exist in the FA-liganded QD DOS (Figure 4(a)), and FA-liganded QD has generally higher DOS at conduction band. This difference in the DOS could indeed affect the optical properties of the QD lattice, as confirmed by the optical absorption data (Figures 4(b) and 4(e)). The optical absorption coefficient of Cl-liganded QD also has an isolated peak near 1.5 eV (highlighted by red-dashed circle, as shown in Figure 4(e)), but the peak is not observed in the FA-liganded QD data. The FA-liganded optical coefficient has higher value around 1.75 eV than that of Cl-liganded QD. This difference could be further understood by comparing orbital-projected DOS (Figures 4(c) and 4(f)). The orbital DOS of FA-liganded QD has one extra peak over the energy range of 0.6–2.0 eV in conduction band, while this peak does not exist for Cl ligand. The hybridization with orbital could be the reason of the absence of intermediate band for FA-liganded QD.

3.3. Temperature-Dependent Carrier Mobility for FA-Liganded PbS QD Array

To understand the carrier transport properties of the QD arrays, mobility of QD array junctions is calculated by using the MD-Landauer method for a wide range of temperatures. Figure 5 depicts the structure of the junctions with the element of left and right electrode, as well as the central MD region. The 1D QD array is sandwiched by two CH(NH2)2PbI3 perovskite electrodes. Figure 5(a)5(c) showed single FA-ligand QD, double FA-liganded QDs, and triple FA-liganded QDs. Figure 5(d)5(f) corresponded to the single Cl-liganded QD, double Cl-liganded QDs, and triple Cl-liganded QDs, respectively. This method is described in Section 2 and Formula (1)–(5). The conductance of a junction with MD central region length and temperature is obtained. To adopt the Formula (5), the bulk resistivity has to be calculated in advance using the formula . Here, the area is fixed to the calculated value of in all the simulations. Due to the space limit, we present all the detailed simulation information for deriving the temperature dependent in the Supplementary Materials. Table S2 summarized all the parameters for MD-Landauer calculation. Table S3 provided the carrier density value of used in Formula (5) for both FA- and Cl-liganded PbS QD arrays under different temperatures. The carrier density was calculated by QuantumATK software with parameters: contribution = electron and Fermi shift = 0.13 eV. Table S4 presented the size information characterizing the MD region. Three types of MD region containing one QD, two QDs, and three QDs were prepared and labeled as , respectively, while the PbS QD without ligand was labeled as . Tables S5 and S6 summarized all the resistance data for Cl-liganded PbS QD and FA-liganded PbS QD for performing linear regression and derivation of the resistivity. The length- and temperature-dependent transmittance for bare PbS QD, FA-liganded sample, as well as Cl-liganded sample was summarized in Figure S2. Notice that we have only presented the results simulated at 300 K as an example in Figure S2. In addition, attention must be paid when interpreting the transmittance data, as shown in Figure S2, because the transmittance data alone is not proportional to the carrier mobility and the energy integrated transmittance, namely, the conductance is used to prepare the resistance data for deriving the resistivity , which is in fact inverse proportional to the carrier mobility . Figures 6(a) and 6(b) show the simulation and regression data for FA-liganded QDs and Cl-liganded QDs, respectively, by taking the temperature condition of 300 K as an example. Figure 6(c) showed the log-linear plot of the mobility calculated by the Formula (5) using the parameters listed in Table S2S6. It can be clearly seen from Figure 6(c) that FA-liganded QD (red solid line with circle marker) has higher mobility than that of Cl-liganded QD (blue solid line with square marker) near the room temperature. For FA-liganded QDs at 300 K, the resistivity of was derived by linearly fitting the data, as shown in Figure 6(a), using the Formula (4) and the mobility was simulated as by using Formula (5) where the carrier density was used. For Cl-liganded QDs at 300 K, the resistivity is by fitting the data, as shown in Figure 6(b), using the Formula (4) and the mobility is by using the carrier density . The numerical data of the other temperatures for both FA- and Cl-liganded PbS QDs could be found in Table S3, Table S4, Table S5, and Table S6 in Supplementary Materials.

4. Discussion

In the simulation data, both the mobilities of FA-liganded QD arrays and Cl-liganded QD arrays decrease with the increase of the temperature. Similar temperature-dependent behavior has been reported in the mobilities simulation of the silicon bulk material and silicon nanowire using the same approach [22]. Meanwhile, it is found that the FA-liganded PbS QD arrays have generally higher mobility than that of Cl-liganded QD arrays for a wide range of temperatures. The enhanced mobility appeared in FA-liganded PbS QD could be attributed to the following factors:(1)The difference in transport property could be explained by the corresponding band structure calculation, as shown in Figures 2(a) and 2(b), which indicates the disappearance of intermediate band in FA-liganded QD arrays. This fact also affects the optical absorption data, where FA-liganded QDs have a bulk-like absorption profile with no isolated peak and a higher optical absorption value than that of Cl-liganded QDs. Our DOS data suggests that the disappearance of discrete feature in FA-liganded sample could be related to the hybridization of orbitals. We claim that the FA-liganded QD arrays showed much stronger continuum band featured bulk-like electronic property than the Cl-liganded ones, where the quantum size effect-induced discrete energy features are still dominating. As demonstrated in literatures [22, 30], it has been confirmed both experimentally and theoretically that bulk silicon showed about one order higher mobility than its counterpart of silicon nanowire. In our work, the absolute difference between the mobility of FA- and Cl-liganded sample is quite small, which indicates that although the bulk feature in FA-liganded PbS QD prevails, the effect is still weak.(2)Based on the mobility calculation Formula (5), it is found that the carrier density plays equal role as the resistivity. By closely examining the carrier density simulated for both FA- and Cl-liganded PbS QD array, we found that the carrier density in FA-liganded sample is much less than the Cl-liganded sample. The bandgap value of FA-liganded sample is the smallest when compared to the values of bare PbS QD array and Cl-liganded QD array, as shown in Figures 2(a) and 2(b). One would expect larger carrier density in FA-liganded sample due to the narrower bandgap. However, the simulated carrier density based on rigorous extraction from the DFT simulation showed opposite tendency. We attribute this effect to the multiple energy bands appeared in both the bare PbS QD array sample and Cl-liganded sample. Additional intermediate band provides additional channel for carrier generations, which is exactly the appealing feature for intermediate band solar cell. However, as suggested in the mobility simulation, the enhanced carrier density will inevitably cause the reduction of mobility, which is exactly the primal obstacle, which has been shown both experimentally and theoretically to hinder the realization of the high-efficiency intermediate band solar cell [3133]. In the optimization of power generation of photovoltaic devices, both carrier transport and optical properties should be considered. The calculation in this study provides a fundamental guide for better design of solar cell devices.

5. Conclusions

In this study, we theoretically analyzed the electronic, optical properties and temperature-dependent mobility of QDs with FA ligand by taking the Cl ligand as a reference. FA-liganded PbS QDs showed larger optical absorption coefficients due to the elimination of intermediate levels or band. Mobility simulation results verified that the mobility in FA ligand is larger than the Cl ligand. Our work revealed a novel type of mobility enhancement originated from orbit overlapping, which could provide important information for further device improvement.

Data Availability

All the simulation codes for reproducing the simulation results are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors gratefully acknowledge the funding from New Energy and Industrial Technology Development Organization (NEDO) (Grant Number JPNP20015) and the Ministry of Economy, Trade and Industry (METI), Japan.

Supplementary Materials

Table S1: lists the parameters for DFT-LCAO calculation. Table S2: lists the parameters for MD-Landauer calculation. Table S3: records the carrier density. Table S4: shows the lengths of MD regions. Tables S5 and S6: shows the regression data for Cl ligand and FA ligand, respectively. Figure S1: contains the band structure comparison for different inter-QD distances. Figure S2: shows carrier transmittance for simulating mobility. (Supplementary Materials)