Introduction

In order to improve the photo-voltaic performance of perovskite solar cells (PSCs), we need to further explore the mechanism such as carrier mobility1,2,3, ion migration4,5,6,7,8, density of trap states (DOST) distribution9,10,11,12,13, carrier recombination14,15,16,17,18 and so on. The DOST distribution is a crucial factor that determines the photovoltaic performance of PSCs19,20,21,22,23,24,25,26,27,28,29,30,31,32,33. According to the literatures in recent years19,20,21,22,23, current density–voltage (J–V) hysteresis is caused by the ion migration and trap assisted carrier recombination. The DOST distribution affects the carrier recombination24,25,26,27, influences the open circuit voltage28,29,30, and hinders the enhancement of power conversion efficiency (PCE)31,32,33,34. The DOST distribution cannot be obtained by experimental measurement directly. There are only few methods for the extraction of DOST distribution. The space charge limited current (SCLC) method uses the deconvolution to extract the DOST distribution35,36,37,38,39,40, based on the measured JV data at different temperatures. Walter et al.41,42,43,44,45 put forward the impedance spectroscopy (IS) method to extract the DOST distribution. They extract the DOST distribution based on the plot of equivalent chemical capacitance versus frequency given by IS measurement41,42,43,44,45. They regard the capture and de-capture of carrier by the trap states in the PSCs as charging and discharging of the equivalent chemical capacitance41. They put forward the formula DOST (Eω) = (Vbi/eW)(dC/dω)(ω/kB) to extract the DOST distribution41,42,43,44,45. Here, C is the equivalent chemical capacitance. ω is the angular frequency of the ac signal. Vbi is the built-in electric voltage. W is the depletion width. kB is the Boltzmann constant. e is the elementary charge. The corresponding energy level is calculated by the formula Eω = kBTln(ω0/ω), where T is the ambient temperature and ω0 is the attempt-to-escape frequency41,42,43,44,45. Wang et al.46 proposed a transient photo-voltage (TPV) method for DOST distribution extraction. Based on the hypothesis of exponential type DOST distribution46,47,48,49, multiple-trapping model46,50,51,52, and the zero-temperature approximation46,52, they find that when the DOST distribution is exponential type, the logarithm of carrier lifetime and the photo-voltage satisfy linear relation. They use this relation to extract the DOST distribution based on the TPV result46,47,48. Because of the hypothesis of exponential type distribution46,47,48,49, their method is effective only when the TPV result is linear. The DOST distribution extracted by their method is an exponential type distribution46,47,48.

However, according to the TPV experiments reported in recent years46,47,48, the majority of TPV results are non-linear. In these cases, the method of Wang et al. is not effective to extract the DOST distribution. In this article, we put forward a new technique for extraction of DOST distribution based on the TPV measurement result. We give up the hypothesis of exponential type DOST distribution46,47,48,49 and zero-temperature approximation46,52. The method given in our work is based on the single hypothesis of multiple-trapping model46,50,51,52. Our method is effective for arbitrary TPV results and can be used to extract general type DOST distribution.

Results and discussion

Establishment of theory and method

In this section, we establish the equation for the extraction of DOST distribution based on the TPV result.

Because of the trap states in the perovskite absorber layer, the behavior of trapping and de-trapping of carrier by the trap states determines the recombination rate, which can be described by the multiple-trapping model46,50,51,52. According to the multiple-trapping model, the relation of carrier lifetime τn and the free carrier lifetime τf satisfies46,50,51,52

$$\tau _{{\text{n}}} = (\partial n/\partial n_{\text{c}} )\tau _{{\text{f}}} .$$
(1)

Here n = nt + nc denotes the sum of electron density in trap states nt and electron density in conduct band nc.

Therefore, we have

$$\tau_{\text {n}} = (1 + \partial n_{\text{t}} /\partial n_{\text{c}} )\tau_{\text{f}} ,$$
(2)

which can be rewritten as

$$\tau_{{\text{n}}} = \left( {1 + \frac{{\partial n_{{\text{t}}} /\partial E_{{{\text{Fn}}}} }}{{\partial n_{{\text{c}}} /\partial E_{{{\text{Fn}}}} }}} \right)\tau_{{\text{f}}} .$$
(3)

The density of electron in trap states satisfies46,52

$$n_{{\text{t}}} = \int \rho_{{\text{t}}} \left( E \right)f\left( E \right)dE.$$
(4)

Here \(\rho_{{\text{t}}} \left( E \right)\) denotes the DOST distribution. \(f\left( E \right) = \frac{1}{{exp\left( {\frac{{E - E_{{{\text{Fn}}}} }}{{k_{B} T}}} \right) + 1}}\) denotes the Fermi–Dirac distribution53.

Therefore, we have

$$\frac{{\partial n_{{\text{t}}} }}{{\partial \left( {\frac{{E_{{{\text{Fn}}}} }}{{k_{{\text{B}}} T}}} \right)}} = \frac{\partial }{{\partial \left( {\frac{{E_{{{\text{Fn}}}} }}{{k_{{\text{B}}} T}}} \right)}}\left( {\int \rho_{{\text{t}}} \left( E \right)f\left( E \right)dE} \right).$$
(5)

According to the Fermi–Dirac distribution53, we have

$$\frac{{\partial n_{{\text{t}}} }}{{\partial \left( {\frac{{E_{{{\text{Fn}}}} }}{{k_{{\text{B}}} T}}} \right)}} = \int \rho_{{\text{t}}} \left( E \right)f\left( E \right)\left( {1 - f\left( E \right)} \right)dE.$$
(6)

We rewrite Eq. (6) as

$$\frac{{\partial n_{t} }}{{\partial E_{Fn} }} = \frac{1}{{k_{B} T}}\int \rho_{{\text{t}}} \left( E \right)f\left( E \right)\left( {1 - f\left( E \right)} \right)dE.$$
(7)

The carrier density in conductor band satisfies46,52,53

$$n_{{\text{c}}} = N_{{\text{c}}} exp\left( {\frac{{E_{{{\text{Fn}}}} - E_{{\text{c}}} }}{{k_{{\text{B}}} T}}} \right).$$
(8)

Here Nc is the density of effective states in conduction band. Ec is the conduction band energy level position53.

Therefore, we have

$$\partial n_{\text{c}} /\partial E_{\text{Fn}} = n_{\text{c}} /k_{\text{B}} T.$$
(9)

According to Eqs. (3), (7), and (9), we have

$$\tau_{{\text{n}}} = \left( {1 + \frac{{\int \rho_{{\text{t}}} \left( E \right)f\left( E \right)\left( {1 - f\left( E \right)} \right)dE}}{{n_{{\text{c}}} }}} \right)\tau_{{\text{f}}} .$$
(10)

We rewrite Eq. (10) as

$$\int \rho_{{\text{t}}} \left( E \right)f\left( E \right)\left( {1 - f\left( E \right)} \right)dE = \left( {\frac{{\tau_{{\text{n}}} }}{{\tau_{{\text{f}}} }} - 1} \right)n_{{\text{c}}} .$$
(11)

Substituting the Fermi–Dirac distribution into Eq. (11), we have

$$\int \rho_{{\text{t}}} \left( E \right)\frac{1}{{exp\left( {\frac{{E - E_{{{\text{Fn}}}} }}{{k_{B} T}}} \right) + 1}}\left( {1 - \frac{1}{{exp\left( {\frac{{E - E_{{{\text{Fn}}}} }}{{k_{B} T}}} \right) + 1}}} \right)dE = \left( {\frac{{\tau_{{\text{n}}} }}{{\tau_{{\text{f}}} }} - 1} \right)n_{{\text{c}}} .$$
(12)

We rewrite Eq. (12) as

$$\int \rho_{{\text{t}}} \left( E \right)\frac{1}{{exp\left( {\frac{{ - \left( {E_{{{\text{Fn}}}} - E} \right)}}{{k_{B} T}}} \right) + 1}}\left( {1 - \frac{1}{{exp\left( {\frac{{ - \left( {E_{{{\text{Fn}}}} - E} \right)}}{{k_{B} T}}} \right) + 1}}} \right)dE = \left( {\frac{{\tau_{{\text{n}}} }}{{\tau_{{\text{f}}} }} - 1} \right)n_{{\text{c}}} .$$
(13)

We define a derivation factor

$$g\left( {E_{{{\text{Fn}}}} - E} \right) = \frac{1}{{exp\left( {\frac{{ - \left( {E_{{{\text{Fn}}}} - E} \right)}}{{k_{B} T}}} \right) + 1}}\left( {1 - \frac{1}{{exp\left( {\frac{{ - \left( {E_{{{\text{Fn}}}} - E} \right)}}{{k_{B} T}}} \right) + 1}}} \right).$$
(14)

and rewrite Eq. (13) as

$$\int \rho_{{\text{t}}} \left( E \right)g\left( {E_{{{\text{Fn}}}} - E} \right)dE = \left( {\frac{{\tau_{{\text{n}}} }}{{\tau_{{\text{f}}} }} - 1} \right)n_{{\text{c}}} .$$
(15)

According to the TPV result, the carrier lifetime is a function of photo-voltage46,47,48,49. Therefore, we rewrite Eq. (15) as

$$\int \rho_{{\text{t}}} \left( E \right)g\left( {E_{{{\text{Fn}}}} - E} \right)dE = \left( {\frac{{\tau_{{\text{n}}} \left( {V_{{{\text{ph}}}} } \right)}}{{\tau_{{\text{f}}} }} - 1} \right)n_{{\text{c}}} .$$
(16)

Equation (16) is the fundamental equation of our method. The right-hand side of Eq. (16) can be obtained from experimental measurement. τn(Vph) can be measured from TPV experiment. nc can be obtained from differential charging method5 or from SCLC under different intensity of illumination39. We can use absorbance spectrum, Kelvin probe (KP), ultraviolet photoemission spectroscopy (UPS), and X-ray photoelectron spectroscopy (XPS) to get conduction band energy level position Ec, valence band energy level position Ev, Fermi energy level position EF0 and band gap Eg35, respectively. After getting the conductor band electron density in dark n0 and in different intensity of illumination nc, according to the relation EFn = EF0 + kBTln(nc/n0)53, we obtain the electron quasi Fermi energy level EFn corresponding to the photo-voltage Vph in the different intensity of illumination. The free carrier lifetime τf can be measured or can be calculated according to the equationτf = 1/CnNt54. Here, Nt is the electron trap concentration. Cn denotes the capture coefficients for electrons. Since the left-hand side of Eq. (16) is the convolution integral of the DOST distribution and the derivative factor, we can get the DOST distribution by deconvolution. We use the numerical deconvolution function ([q,r] = deconv(u,v)) of MATLAB to solve the equation for DOST distribution.

For intrinsic perovskite, the electron density in conductor band satisfies53

$$n_{{\text{c}}} = N_{{\text{c}}} exp\left( {\frac{{E_{{{\text{Fn}}}} - E_{{{\text{cb}}}} }}{{k_{{\text{B}}} T}}} \right) = n_{0} exp\left( {\frac{{V_{{{\text{ph}}}} e}}{{2k_{{\text{B}}} T}}} \right).$$
(17)

Here \(n_{0} = N_{{\text{c}}} exp\left( {\frac{{E_{{{\text{F}}0}} - E_{{{\text{cb}}}} }}{{k_{{\text{B}}} T}}} \right)\) is the electron density in conductor band at dark state53.

Substituting Eq. (17) into Eq. (16), we have

$$\int \rho_{{\text{t}}} \left( E \right)g\left( {E_{{{\text{Fn}}}} - E} \right)dE = \left( {\frac{{\tau_{{\text{n}}} \left( {V_{{{\text{ph}}}} } \right)}}{{\tau_{{\text{f}}} }} - 1} \right)n_{0} exp\left( {\frac{{V_{{{\text{ph}}}} e}}{{2k_{{\text{B}}} T}}} \right).$$
(18)

Equation (18) is the fundamental equation for DOST distribution extraction of intrinsic perovskite. Similarly, we obtain the DOST distribution using numerical deconvolution function in MATLAB. The parameters used for calculation are list in Table 127,35.

Table 1 Parameters used for DOST distribution calculation27,35.

Exponential type DOST distribution

In this section, we explore the exponential type DOST distribution. Wang et al.46 find that for the exponential type DOST distribution, the logarithm of carrier lifetime is proportional to the photo-voltage (TPV result). The proof of this relation is shown as follows.

The exponential type DOST distribution satisfies46,47,48,49

$$\rho_{{\text{t}}} \left( E \right) = \frac{{N_{T} }}{{E_{{\text{B}}} }}exp\left( {\frac{{E - E_{{\text{c}}} }}{{E_{{\text{B}}} }}} \right).$$
(19)

Here, EB is the characteristic energy and NT is the total density of the trapped state. Substituting formula (19) into formula (4), and taking zero-temperature approximation46,52, we have

$$n_{{\text{T}}} = N_{{\text{T}}} \left( {exp\left( {\frac{{E_{{{\text{Fn}}}} - E_{{\text{c}}} }}{{E_{{\text{B}}} }}} \right) - exp\left( {\frac{{E_{{\text{v}}} - E_{{\text{c}}} }}{{E_{{\text{B}}} }}} \right)} \right).$$
(20)

Making the approximation of n ≈ nT50, we rewrite the multiple-trapping model46,50,51,52 as

$$\tau_{\text{n}} = (\partial n_{\text{T}} /\partial n_{\text{c}} )\tau_{\text{f}} ,$$
(21)

which is equivalent to

$$\tau_{{\text{n}}} = \frac{{\partial n_{{\text{T}}} /\partial V_{{{\text{ph}}}} }}{{\partial n_{{\text{c}}} /\partial V_{{{\text{ph}}}} }}\tau_{{\text{f}}} .$$
(22)

The electron density in conductor band satisfies53

$$n_{{\text{c}}} = N_{{\text{c}}} exp\left( {\frac{{E_{{{\text{Fn}}}} - E_{{\text{c}}} }}{{k_{{\text{B}}} T}}} \right).$$
(23)

According to Eqs. (20), (22), and (23) and from the relation of EFn = EFp + eVph, we have

$$ln \tau_{\text{n}} = \left( {e/E_{\text{B}} - e/k_{\text{B}} T} \right)V_{\text{ph}} + lnA.$$
(24)

Here \(A = \frac{{N_{{\text{T}}} k_{{\text{B}}} T}}{{N_{{\text{c}}} E_{{\text{B}}} }}exp\left( {\frac{{E_{{{\text{Fp}}}} - E_{{\text{c}}} }}{{E_{{\text{B}}} }} - \frac{{E_{{{\text{Fp}}}} - E_{{\text{c}}} }}{{k_{{\text{B}}} T}}} \right)\tau_{{\text{f}}}\).

We write Eq. (24) as linear mathematical form lnτn = aVph + b.

Here, a = e/EB − e/kBT, b = lnA.

Therefore, we finish the proof of this relation. Note that for the intrinsic perovskite, Eq. (24) can be written as

$$ln \tau_{\text{n}} = \left( {e/2E_{\text{B}} - e/2k_{\text{B}} T} \right)V_{\text{ph}} + lnA.$$
(25)

Here \(A = \frac{{N_{{\text{T}}} k_{{\text{B}}} T}}{{N_{{\text{c}}} E_{{\text{B}}} }}exp\left( {\frac{{E_{{{\text{F}}0}} - E_{{\text{c}}} }}{{E_{{\text{B}}} }} - \frac{{E_{{{\text{F}}0}} - E_{{\text{c}}} }}{{k_{{\text{B}}} T}}} \right)\tau_{{\text{f}}}\). We write Eq. (25) as linear form lnτn = aVph + b. Here, a = e/2 EB− e/2kBT, b = lnA.

Similarly, we can also prove that when TPV result is linear, the extracted DOST distribution is exponential type. Details of the proof are given in the supporting information. We can use this relation to extract the DOST distribution when the TPV result is linear. Below, this method is called analytic method.

We can use the analytic method to verify the validity of our numerical method. We use both our numerical method and the analytic method to extract the DOST distribution and make a comparison. In Fig. 1a–c, we set the b = − 5 and the slope parameter a as − 1.5, − 2 and − 2.5, respectively. Using the analytic method, we derive EB and NT (see Table 2). Figure 1d–f shows the extracted DOST distributions from Fig. 1a–c using our numerical method. As expected by the analytic method, the DOST distribution is exponential type. In order to compare with the DOST distribution extracted by analytic method, we use the exponential fitting (f(E) = cexp(dE)) to calculate the EB and NT (see Table 2). As shown in Table 2, the EB and NT calculated by our numerical method are consistent with the EB and NT calculated by the analytic method, indicating that the numerical algorithm to do the deconvolution in our calculation is reliable. We explain the slight deviation of EB and NT extracted by the two methods as follows. The analytic method is established based on the three hypothesis of multiple-trapping model46,50,51,52, zero-temperature approximation46,52, and exponential type DOST distribution46,47,48,49. Our numerical method is established based on the unique hypothesis of multiple-trapping model46,50,51,52. Therefore, the slight deviation of EB and NT extracted by the two methods is attributed to the zero-temperature approximation and the fitting error.

Figure 1
figure 1

(a)–(c) Plots of carrier lifetime versus photo-voltage. We set b =  − 5 and a as − 1.5, − 2 and − 2.5, respectively. (d)–(f) DOST distributions extracted from (a)–(c), respectively. The blue dots represent the extracted DOST distributions, the red lines represent the exponential fitting of extracted DOST distributions.

Table 2 Exponential fitting coefficients, distribution coefficients calculated from Fig. 2.

Non-exponential type DOST distribution

In this section, we investigate the non-exponential type DOST distribution. We take the TPV data in Ref.47 as an example. As shown in Fig. 2a,b, the TPV result is non-linear. Hence, we cannot use the analytic method to extract the DOST distribution.

Figure 2
figure 2

Plots of carrier lifetime versus photo-voltage given by TPV measurement47. The black lines represent the result of TPV measurement. The red lines represent linear fitting in the subintervals. (a) shows the linear fittings of TPV in two subintervals, respectively. (b) shows the linear fittings of TPV in three subintervals, respectively.

To overcome this difficulty, Wang et al.47 made a linear fitting in the subintervals of 0.05–0.6 V and 0.6–0.91 V (see Fig. 2a). They used formula (24) to get two values of EB in these two subintervals, respectively47. They explained these two EB as two types of exponential type DOST distribution (They called them as deep trap type and shallow trap type)47. However, there are some difficulties in their method. (1) Equation (24) is derived from one exponential type DOST distribution, not the two types of DOST distribution (deep trap state type DOST distribution and shallow trap state type DOST distribution). We cannot derive Eq. (24) based on these two types of distribution. (2) There is no clear boundary for deep and shallow trap. So the definition of deep and shallow trap is ambiguous. (3) It is more accurate to take linear fittings in three subintervals, respectively (see Fig. 2b). According to the differential theory, we can divide the photo-voltage interval into infinite differential subintervals and use Eq. (24) to calculate the EB in these differential subintervals, respectively. These EB cannot be explained by the concept of deep and shallow trap.

Our method is effective for arbitrary TPV results, which can be used to extract general type DOST distribution. Using our method, we extract the DOST distribution (as illustrated in Fig. 3). It can be seen that the extracted DOST distribution is not an exponential type distribution. Compared to the method given by Wang et al., our method is more accurate and it can give the fine structure of DOST distribution.

Figure 3
figure 3

DOST distribution calculated by our method given in this article. The blue dots represent the extracted DOST distribution. The red, yellow and purple lines represent the exponential, double exponential and Gauss fittings for the extracted DOST distribution, respectively.

We make the exponential, double exponential and Gauss fittings for the extracted DOST distribution, respectively (see Fig. 3). The R-square (Coefficient of determination) of these fittings are given in Table 3. It can be see that the R-square of double exponential fitting is larger than the exponential fitting and Gauss fitting, indicating that the DOST distribution is more consistent with the double exponential type distribution than with the exponential type distribution. Therefore, when the TPV result is non-linear, the extracted DOST distribution is non-exponential type.

Table 3 Different type fitting for extracted DOST distribution.

Calculation example

In this section, we give some examples of DOST distribution extraction using our method. We take the TPV data in Ref.47 as an example. Figure 4a–c shows the TPV results for meso-structured perovskite solar cells with large, middle and small size of perovskite grain, respectively. Using our method, we extract the DOST distributions from Fig. 4a–c, respectively, (as shown in Fig. 4d–f).

Figure 4
figure 4

(a)–(c) TPV results for meso-structured perovskite solar cells with large, middle and small size of perovskite grain47. (d)–(f) DOST distributions extracted by our method from (a)–(c), respectively.

For comparing, we plot the extracted DOST distributions in Fig. 5. Using the formula \({T}_{\mathrm{total}}=\int {\rho }_{\mathrm{t}}\left(E\right)dE\), we calculate the total amount of trap states for three solar cells. We derived Ttotal = 5.4431 × 1020 cm−3 (large size), Ttotal = 5.9096 × 1020 cm−3 (middle size) and Ttotal = 4.5628 × 1020 cm−3 (small size). It can be seen that the solar cell with small size of perovskite grain has the least amount of trap states. This result could give a guidance for preparing perovskite solar cells with less trap states.

Figure 5
figure 5

DOST distributions for meso-structured perovskite solar cells with large, middle and small size of perovskite grain extracted by our method given in this paper.

Comparison to SCLC method

SCLC measurement gives the relation of current density and voltage at different temperatures (j = j (U, T)). For SCLC method, we extract the DOST distribution adopting the equations listed as follows35,36,37,38,39,40

$$E_{\text{a}} = - dlnj/d\left( {k_{\text{B}} T} \right)^{ - 1} ,$$
(26)
$$m = d\left( {lnj} \right)/d\left( {lnU} \right),$$
(27)
$$B = - \left[ {dm/d\left( {lnU} \right)} \right]/\left[ {m\left( {m - 1} \right)\left( {2m - 1} \right)} \right],$$
(28)
$$C = \left( {B\left( {2m - 1} \right) + B^{2} \left( {3m - 2} \right) + d\left( {ln\left( {1 + B} \right)} \right)/d\left( {lnU} \right)} \right)/\left( {1 + B\left( {m - 1} \right)} \right),$$
(29)
$$\int \rho_{{\text{t}}} \left( E \right)f\left( E \right)\left( {1 - f\left( E \right)} \right)dE = \frac{1}{{k_{{\text{B}}} T}}\frac{\varepsilon U}{{eL^{2} }}\frac{2m - 1}{{m^{2} }}\left( {1 + C} \right).$$
(30)

Here j denotes the current density. U denotes the voltage. Ea is the activation energy. m, B, C are the parameters for calculation. L is the thickness of perovskite absorber layer. e is the elementary charge. ε is the dielectric constant of perovskite absorber layer. Based on the SCLC measurement result, we need to calculate Ea, m, B, C and finally use Eq. (30) to extract the DOST by deconvolution35,36,37,38,39,40. The calculations of Ea, m, B, C are complicated. For our method, we only need to use formula (18) to extract the DOST distribution by deconvolution. Therefore, our method takes less computation.

Conclusion

In conclusion, this article presents a new technique for DOST distribution extraction based on the TPV measurement result. The approach given in this paper is effective for extraction of general type DOST distribution. We prove that when the TPV result is linear, the DOST distribution is exponential type and vice versa. Our method needs less computation than the SCLC method. The obtained results provide a guidance for preparing perovskite solar cells with less trap states.

Methods

The transient photo-voltage measurement is an effective technique for the study of carrier recombination3,46,47,48,49. Figure 6a shows the typical TPV measurement result. Figure 6b shows the mechanism of TPV experiment. Photo-voltaic device is held under open circuit. At the start of the TPV experiment, we use a steady-state bias light to illuminate the device until the equilibrium between generation and recombination is established. The steady-state bias light produces a bias photo-voltage Vph. As a result, the Fermi energy level EF0 changes to electron quasi Fermi energy level EFn and hole quasi Fermi energy level EFp. Thereafter, we apply an additional small light pulse to the device. With the perturbation of small light pulse, the photo-voltage increases to Vph + ΔV. The electron quasi Fermi energy level EFn shifts to \({E}_{\mathrm{Fn}}^{*}\), and the hole quasi Fermi energy level EFp to \({E}_{\mathrm{Fp}}^{*}\). After switching off the small light pulse source, due to carrier recombination, the photo-voltage decays exponentially until it reaches the bias value Vph. The electron quasi Fermi energy level comes back to EFn from \({E}_{\mathrm{Fn}}^{*}\) and the hole quasi Fermi energy level comes back to EFp from \({E}_{\mathrm{Fp}}^{*}\).The carrier lifetime τn is defined as the time taken for the photo-voltage to decay from Vph + ΔV to Vph + ΔV/e. Here, e is the natural constant. By tuning the intensity of steady-state bias light I, we get the relation between the photo-voltage Vph and carrier lifetime τn, which is written as τn = τn(Vph) (TPV result)46,47,48,49.

Figure 6
figure 6

(a) Result of TPV measurement. (b) Mechanism of TPV experiment.