• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    In situ studies on the positive and negative effects of 1,8-diiodoctane on the device performance and morphology evolution of organic solar cells

    2021-07-02 09:29:04DaHuangChunXiaHongJianHuaHanNieZhaoXiuHongLiXiChangBaoFengGangBianYuYingHuangChunMingYang
    Nuclear Science and Techniques 2021年6期

    Da Huang · Chun-Xia Hong · Jian-Hua Han· Nie Zhao ·Xiu-Hong Li · Xi-Chang Bao · Feng-Gang Bian · Yu-Ying Huang ·Chun-Ming Yang

    Abstract The introduction of solvent additives is one of the most common approaches for enhancing the power conversion efficiency of organic solar cells (OSCs). However,the use of solvent additives has some negative effects,and an understanding of how solvent additives affect OSCs is currently limited. In this study, we developed an in situ grazing incidence wide-angle X-ray scattering (GIWAXS)technique in the SAXS beamline(BL16B1)at the Shanghai Synchrotron Radiation Facility, and the additive effects of 1,8-diiodoctane(DIO)on the performance and morphology evolution of the PTB7-Th/PC71BM device was investigated in depth. The results revealed that the crystal size increased with the volume ratio of DIO, and a drastic evolution of lattice space and crystal coherence length was observed during thermal annealing for the first time,to our knowledge. The discrete PC71BM molecules dissolved by DIO have an effect similar to that of the nucleating agent for PTB7-Th, boosting the crystallization of PTB7-Th,reducing phase separation, and inducing more drastic morphological evolution during thermal annealing. Our results provide a deep perspective for the mechanism of solvent additives, while also showing the significance and feasibility of the in situ GIWAXS technique we developed at BL16B1.

    Keywords Grazing incidence wide-angle X-ray scattering · Grazing incidence small-angle X-ray scattering·Organic solar cells·Solvent additive·Thermal stability

    1 Introduction

    Organic solar cells (OSCs), which are lightweight,flexible, and environmentally friendly [1], have power conversion efficiencies (PCEs) of over 17.5% [2]. The morphology of the active layer plays a crucial role in the performance of OSCs [3], which can be manipulated by several approaches, such as the selection of donor and acceptor materials [4], selection of the used solvent [5],thermal annealing [6], solvent vapor annealing [7], upside down solvent annealing treatment [8], and the addition of solvent additives [9]. For the solvent additive approach,1,8-diiodoctane (DIO) is one of the most commonly used solvent additives to enhance OSC device performance[10].However, a series of negative effects, such as light instability, thermal instability, and long-term instability, can be induced by the use of DIO [11]. However, the mechanism of both the positive and negative effects of DIO usage on morphology is not yet completely understood. Grazing incidence wide-angle X-ray scattering (GIWAXS), which provides statistical information in a large region, is sensitive to the surface and inner structure information of the active layer [12]. More importantly, in situ GIWAXS can monitor the dynamic evolution of active layer morphology,which is of great importance for related mechanism research[13].Most of the current research was based on ex situ GIWAXS methods to characterize the active layer morphology [14], and therefore, the understanding of the mechanism remains limited.In this study,we developed an in situ GIWAXS technique in the small-angle X-ray scattering beamline (BL16B1) at the Shanghai Synchrotron Radiation Facility (SSRF) and probed the dynamic morphology evolution of the active layer of a typical fullerene system of OSCs with donor poly[(ethylhexyl-thiophenyl)-benzodithiophene-(ethylhexyl)-thienothiophene] (PTB7-Th) and acceptor [6,6]-phenyl C71-butyric acid methyl ester (PC71BM). A phenomenon in which the evolution was enlarged with an increase in the proportion of DIO being observed for the first time.The mechanism by which the discrete PC71BM molecules dissolved by DIO had a similar effect to that of the nucleating agent for PTB7-Th,leading to a larger crystal size of PTB7-Th, thereby, promoting the morphology evolution during thermal annealing, was revealed, and is of great significance for the further development of OSCs.

    2 Experimental section

    2.1 Development of the in situ GIWAXS technique

    The small-angle X-ray scattering (SAXS) beamline(BL16B1) at the SSRF [15], which mainly engages in the application of X-ray scattering [16], provides perfect conditions for the development of the GIWAXS method.Figure 1 shows the developed GIWAXS setup for BL16B1. PILATUS3 X 2 M was employed as a 2D area detector, and the pixel size was 172 μm × 172 μm, the number of pixels was 1475 × 1679, and the maximum frame rate was 250 Hz.The high frame rate guaranteed the ability to collect data during time-resolved measurements.The incident angle could be adjusted by the sample stage,which comprised four Kohzu stages, from 0° to 2°. The four Kohzu stages control the height and shift in the YZ plane,rotation on the y-axis,and rotation on the x-axis.The glove box, as shown in Fig. 1, provided a gaseous atmosphere for the samples and reduced background scattering.The sample temperature was controlled with a homemade temperature controller using a thermocouple that was thermally connected to the top of the sample stage. The temperature range was 25-200 °C.

    Fig. 1 (Color online) Photograph of grazing incidence wide-angle X-ray scattering setup at BL16B1

    2.2 GIWAXS and GISAXS measurements

    For GIWAXS and GISAXS measurements, samples were fabricated by spin-coating the respective solutions of active layers onto 1.5 × 1.5 cm2silicon substrates with pre-coated PEDOT: PSS (30 nm). Silicon substrates were cleaned in sequence with detergent; deionized water; acetone; and isopropyl alcohol in an ultrasonic bath for 15 min, dried by ultra-pure N2, and then treated with O2plasma for 8 min. The GIWAXS and GISAXS measurements were carried out using beamline BL16B1 at the SSRF. The beam size was 0.4 × 0.5 mm2, and the photon flux at the sample was 1 × 1011phs/s. The X-ray wavelength was 0.124 nm (E = 10 keV),and the incident angle was set at 0.15°, which was close to the critical angle, to maximize the intensity. The sample-to-detector distance was 266.2 and 1997.6 mm for GIWAXS and GISAXS measurements,respectively, which was calibrated with the standard sample of silver behenate. For the static measurements, the exposure time of individual measurements was 60 s. For in situ measurements, samples were heated with a hotplate under an air atmosphere. Single images were collected every 60 s to obtain a high signal-to-noise ratio and high sampling rate. For each set of experiments,the pixel coordinates of the images were converted into q values after correction of the beam center by standard samples. One-dimensional profiles in the out-of-plane direction were obtained by integrating at the azimuthal angle around 90°, while the profiles in the in-plane direction were extracted with an azimuthal angle of approximately 0°.

    2.3 Materials

    All reagents and solvents were obtained from commercial sources and used directly without further purification.N, N-Dioxide of bis(N,N′-dimethylaminopropyl)perylene diimide (PDINO) was purchased from Solarmer Materials Inc. (Beijing, China). Poly(3,4 ethylenedioxythiophene):poly(styrenesulfonate) (PEDOT:PSS) (PH1000)was obtained from Heraeus Inc. PTB7-Th was purchased from One-Material Corp. The PC71BM, chlorobenzene,and 1,8-diiodooctane (DIO) used in the device fabrication processes were obtained from Sigma-Aldrich.

    2.4 Device fabrication

    The OSC devices were fabricated with the configuration ITO/PEDOT:PSS/active layer/PDINO/Al. ITO substrates were cleaned with detergent,deionized water,acetone,and isopropyl alcohol in an ultrasonic bath sequentially for 15 min, then dried by ultra-pure N2, and finally treated with O2plasma for 8 min. Then, a 30 nm PEDOT:PSS layer was spin-coated onto the ITO-coated glass at 4,000 rpm and baked at 150 °C for 20 min in an oven.The PTB7-Th/PC71BM solution was prepared by dissolving blends with a weight ratio of 1:1.5 in chlorobenzene (total concentration 20 mg/mL, with 0, 1, 3, and 8 vol% DIO additive).When the substrates were transferred into a glove box, the solutions were stirred overnight and spin-coated on the PEDOT: PSS layer to form active layers (the optimal thickness was approximately 100 nm).Subsequently,a PDINO solution(in CH3OH,1 mg/mL)was spin-coated at 3,000 rpm for 15 s to form the electron transfer layer.Finally, the Al (100 nm) metal electrode was thermally evaporated under approximately 4 × 10-4Pa, and the device area was 0.1 cm2, defined by a shadow mask.

    2.5 Annealing process

    To perform in situ investigations on device PCE and active layer morphology, an annealing strategy that involved placing devices on the hot plate at room temperature (25 °C) and then heating them to the annealing temperature (100 °C) at a heating rate of 20 °C/min was employed to simulate the conventional thermal annealing strategy that involved directly placing devices on a hot plate that had been heated to a specified temperature.Hence, it would remain for 3.75 min in the heating stage process from 25 °C to 100 °C, 10 min in the temperature holding process at 100 °C, and 3.75 min in the cooling process from 100 °C to 25 °C.

    3 Results and discussion

    3.1 DIO addition effects on device PCE and thermal stability

    Fig. 2 (Color online) Device configuration and the molecular structures of PTB7-Th and PC71BM

    Fig.3 (Color online)a J-V curves and b photovoltaic parameters of organic solar cell devices based on PTB7-Th/PC71BM devices with various ratios of 1,8-diiodoctane ranging from 0 to 8 vol%. The average parameters were calculated from 20 independent devices.DIO, 1,8-diiodoctane; VOC, open circuit voltage; JSC, low shortcircuit current density; FF, fill factor; PCE, power conversion efficiency

    Table 1 Photovoltaic parameters of PTB7-Th/PC71BM devices with various ratios of 1,8-diiodoctane

    To demonstrate the effects of using DIO on the PTB7-Th/PC71BM device, devices with various ratios of DIO ranging from 0 to 8 vol% were fabricated with a conventional device structure of ITO/PEDOT:PSS/active layer/PDINO/Al, as shown in Fig. 2. The J-V curves are shown in Fig. 3a, and the detailed photovoltaic parameters for devices with various ratios of added DIO are shown in Fig. 3b and Table 1.As shown in Fig. 3a,a drastic change in the J-V curve of the PTB7-Th/PC71BM device was induced by the addition of DIO. The J-V curve for the device without DIO showed an unfavorable shape with a low short-circuit current density (JSC) and fill factor (FF)compared with the values for devices with added DIO.While the values of open circuit voltage (VOC) for the device without the addition of DIO were basically the same as the value for devices with different ratios of DIO,a trend of decreasing with the increase in DIO was observed. As shown in Fig. 3b, the PCE was only 5.49% when no DIO was added.As the ratio of DIO increased,the value of PCE first increased to 8.48% for the device with 1 vol% DIO,and a maximum value of 9.17% was achieved for devices with 3 vol% DIO. However, the PCE value decreased to 8.36% when the DIO ratio increased to 8 vol%. These results demonstrate that the appropriate addition of DIO can improve the PCE of PTB7-Th/PC71BM devices, while the excessive addition of DIO could harm the device performance. A similar phenomenon in PTB7/PC71BM was reported by Kim et al. [17].

    As shown in Fig. 3b,for devices with no added DIO,the values of VOC, JSC, and FF were 0.83 V, 15.37 mA/cm2,and 43.37%, respectively. After adding DIO, the value of VOCshowed a trend of decreasing with an increase in DIO proportion,where the median value of 0.82 V was obtained from the device with 3 vol% DIO added. For JSCand FF,the values significantly increased after the addition of DIO,and maximum values of 17.77 mA/cm2and 63.33% were achieved,respectively, when the ratio of DIO was 3 vol%,leading to the optimum PCE. However, when the ratio of DIO was increased to 8 vol%, both the FF and JSCvalues decreased, resulting in a decrease in PCE. The increase in JSCand FF appeared to be the main reason for the improvement in PCE induced by the addition of DIO. For OSCs, VOCis generally influenced by the morphology of the active layer, especially large-scale phase separation[18]. JSCand FF are usually dominated by the generation,diffusion, and dissociation of photogenerated excitons and the transport and collection of charge, which are strongly dependent on the donor/acceptor interface and charge transport channel in the active layer [19]. Therefore, the photovoltaic parameters changed with various ratios of DIO, suggesting that there was a morphological evolution in the active layer.

    Fig.4 (Color online) a J-V curves and b normalized photovoltaic parameters of organic solar cell devices based on PTB7-Th/PC71BM devices with various ratios of 1,8-diiodoctane ranging from 0 to 8 vol% after 10 min of thermal annealing at 100 °C. The average parameters were calculated from 20 independent devices. DIO, 1,8-diiodoctane; VOC, open circuit voltage; JSC, low short-circuit current density; FF, fill factor; PCE, power conversion efficiency

    Typical PTB7-Th/PC71BM devices suffer severe PCE degradation even after short-term thermal annealing at temperatures ranging from 70 to 150 °C[20].The addition of DIO can further harm the thermal stability of OSC devices [21]. However, the mechanism by which DIO can harm the thermal stability of PTB7-Th/PC71BM devices is not yet completely understood. Figure 4a, b show the JV curves and normalized photovoltaic parameters obtained after 10 min of thermal annealing at 100 °C for PTB7-Th/PC71BM blends with various ratios of added DIO. As shown in Fig. 4a,the J-V curve after thermal annealing for the device without DIO showed an insignificant change compared to that of the J-V curve before annealing(Fig. 3a). However, drastic changes were observed for devices with various proportions of DIO, where the changes were primarily induced by the decrease in JSCand FF.As shown in Fig. 4b, after thermal annealing, the residual JSC, FF, and PCE values were 99.3%, 92.4%, and 90.4%,respectively,of their initial values for devices without DIO.However, only 95.1%, 82.9%, and 77.6% of JSC, FF, and PCE, respectively, were retained for the devices with 1 vol% DIO. For devices with 3 vol% DIO, the JSC, FF, and PCE values were 91.4%, 83.2%, and 75.8%, respectively.When the proportion of DIO was further increased to 8 vol%, the residual JSC, FF, and PCE values decreased to 88.7%,82.1%,and 71.4%,respectively.The degradation in JSC, FF, and PCE increased after the addition of DIO. In addition,these three parameters showed a decreasing trend as the proportion of DIO increased. However, the value of residual VOCwas approximately unchanged with the increase in DIO,which demonstrated that the effect of DIO on VOCcould be negligible and no large-scale phase separation was induced during thermal annealing. The photovoltaic parameters after thermal annealing are summarized in Table S1. These results demonstrate the negative effect of DIO addition on the thermal stability of the PTB7-Th/PC71BM devices.

    3.2 Effects of addition of DIO on the morphology of the active layer to improve device performance

    To further explore the mechanism by which addition of various ratios of DIO improves device performance,GIWAXS measurements were carried out to monitor the evolution of the active layer morphology after adding various ratios of DIO. The 2D GIWAXS patterns of the blends with various ratios of DIO are shown in Fig. S1.Figure 5a, b show the 1D profiles integrated from the 2D patterns in the out-of-plane and in-plane directions,respectively.For the typical PTB7-Th/PC71BM blends with 3 vol% DIO, the polymer alkyl stacking peaks showed an almost ring-like scattering pattern at q ≈2.90 nm-1with an associated lattice spacing of d ≈2.17 nm, which corresponds to the (100) Bragg reflection and is consistent with that of previous results [22]. The almost ring-like scattering pattern of PC71BM was observed at q ≈13.60 nm-1, which demonstrated that the PC71BM molecules were in an amorphous state. A weak π-π stacking peak,which could hardly be distinguished from the ring of PC71BM, occurred in the out-of-plane direction, implying that the PTB7-Th lamellae had a preference for face-on orientation. The lamellar structure was likely to form nanofibrils under common processing conditions, which could be inclined relative to the substrate, causing the peaks to have an angle spread in the scattering patterns.

    For the (100) peak of PTB7-Th shown in Fig. 5a, b, a weak peak occurred in the blend film without DIO owing to the poor crystallinity of PTB7-Th [22]. The peak intensity increased with decreasing peak width in both the out-ofplane and in-plane directions as the ratio of DIO increased.These results suggested that the crystallization of PTB7-Th could be enhanced by increasing the proportion of DIO.The changes in the peak position and full width at half maximum (FWHM) were also observed after adding various proportions of DIO, indicating that DIO could influence the lattice space (d-spacing) and crystal coherence length(CCL)of PTB7-Th.For the peak of PC71BM,slight changes in FWHM were also observed for blends with various ratios of DIO, while the peak position was approximately constant.

    Fig. 5 (Color online) 1D profiles integrated from 2D patterns in the a out-of-plane and b in-plane directions.Crystal coherence length and d-spacing values of c(100)peaks and d[6,6]-phenyl C71-butyric acid methyl ester rings of blends with various ratios of 1,8-diiodoctane.CCL, crystal coherence length; DIO, 1,8-diiodoctane

    Table 2 Crystal coherence length and lattice space values of the(100)peaks and[6,6]-phenyl C71-butyric acid methyl ester rings in the out-ofplane direction for blends with various ratios of 1,8-diiodoctane

    The d-spacing and CCL values of the (100) peaks are shown in Fig. 5c. The detailed values integrated from the out-of-plane and in-plane directions are summarized in Table 2 and S2, respectively. The CCL values were determined using the Scherre equation,as shown in Eq. (1)[23]:where B is the FWHM of the peaks,K is a constant with a value of 0.89,γ is the X-ray wavelength,and θ is the Bragg diffraction angle. For the data processing process, the software Fit2D [24] and SG-Tool [25] were used to subtract background scattering and convert 2D patterns into 1D profiles,and the software Peakfit[26]was employed to fit the peak position and FWHM from 1D profiles. The detailed data processing process is shown in the Supporting Information (Sect. 2). The original experimental profiles with background and(100)profiles from Fig. 5a are shown in Fig. S2 and S3, respectively.

    As shown in Fig. 5c, the d-spacing increased as the proportion of DIO increased, which suggested that the packing distance between the lamellae was enlarged. The CCL also increased with increasing ratio of DIO,indicating that a larger crystal size was obtained when adding more DIO.When the ratio of DIO was 0 vol%,the CCL value of PTB7-Th was approximately 3.98 nm. The values of CCL for blends with various ratios of DIO were significantly larger (5.10, 5.71, and 6.15 nm when adding 1, 3, and 8 vol%DIO,respectively).The values of CCL obtained from the out-of-plane and in-plane directions were approximately the same regardless of the amount of DIO added,indicating that the orientation of the PTB7-Th crystal was not affected by DIO.

    The d-spacing and CCL values for PC71BM are shown in Fig. 5d. The value of d-spacing for PC71BM remained approximately unchanged at 0.46 nm with various ratios of DIO. The CCL value first decreased and then increased as the proportion of DIO increased, where the CCL values were 2.89 nm, 2.84 nm, 2.77 nm, and 2.83 nm. This phenomenon indicated that adding an appropriate proportion of DIO decreases the size of the PC71BM phase, while excessive addition of DIO increases it.

    To explore the effect of adding DIO on the aggregation of the amorphous phase in the active layer, GISAXS analyses were performed for blends with 0 and 3 vol%DIO to monitor the large-scale morphological information in the active layer. Figure 6a, b show the 2D GISAXS patterns for blends with 0 and 3 vol% DIO. The 1D GISAXS profiles are shown in Fig. 6c, which were integrated from the 2D patterns in the in-plane direction at the value of qznear the reflected position to obtain enhanced signals.The fitting parameters obtained from GISAXS data were fitted with the Debye-Anderson-Brumberger term [27] in the low q region(0.025-0.5 nm-1)as shown in Table 3.The fitting equation is shown in Eq. (2):DIO, 1,8-diiodoctane; ξ,domain size where the acceptor is dispersed in polymer; D, fractal dimension of the structure; η, the correlation length of the fractal-like network

    Table 3 GISAXS fitting parameters of blends with 0 and 3 vol%1,8-diiodoctane

    In polymer fullerene systems, the first term is used to model the scattering of dispersive PC71BM in the amorphous domain, and the second term represents the contribution of fractal-like aggregations of PC71BM, where q is the scattering wave vector, A1 and A2 are independent fitting parameters, and ξ is traditionally used to estimate the size of the domain where the acceptor is dispersed in polymer, which can be easily identified from the ‘‘shoulder’’ occurring in the 1D GISAXS profile. In the second term, P(q, R) is the form factor modeled as spheres, and S(q, R, η, D) is the structure factor of the cluster domain modeled as a fractal-like network expressed in Eq. (3):

    where R is the mean spherical radius of the primary particles,η is the correlation length of the fractal-like network,D is the fractal dimension of the structure, and η is the correlation length of the fractal-like network. 2Rg-polymerand 2Rg-PCBMare the Guinier radii employed to approximately characterize the average domain sizes of the aggregated cluster phases of PTB7-Th and PC71BM,which are calculated according to Eqs. (4) and (5):

    Fig. 6 (Color online) 2D Grazing incidence small angle X-ray scattering patterns for blends with a 0 and b 3 vol% 1,8-diiodoctane. c The 1D profiles integrated from 2D patterns and the fitting curves. DIO, 1,8-diiodoctane

    As shown in Table 3, the values of ξ were approximately constant for blends with 0 and 3 vol% DIO, while the CCL values of blends without DIO were apparently smaller than those of blends with 3 vol%DIO,as shown in Fig. 5c. This indicated that the addition of DIO had no effect on the size of the PTB7-Th cluster, but both d spacing and crystal size were enhanced with the addition of DIO. The stronger aggregation of PC71BM molecules was obtained in blends without DIO, which was demonstrated by the significantly larger η and 2Rg-PCBMvalues shown in Table 3. However, there were more discrete PC71BM molecules in the active layer with 3 vol%DIO,in which a smaller 2Rg-polymervalue of the aggregation size of PTB7-Th was obtained (Table 3). Therefore, the addition of DIO weakened the phase separation in the active layer by weakening the aggregation of PC71BM molecules. A homogeneous phase separation scale was obtained,accompanied by smaller clusters of both PC71BM and PTB7-Th.Schematic diagrams of the 3D nanostructures for blends with 0 and 3 vol%DIO were drawn according to the results of GIWAXS and GISAXS measurements,as shown in Fig. 7.

    For the blend without DIO shown in Fig. 7a, the CCL value of PTB7-Th was significantly lower than that of blends with 3 vol%DIO,according to the GIWAXS results(Fig. 5c), resulting in poor charge transport efficiency due to the lack of charge transport channels.The larger D value(2.97) of the blends without DIO suggested that molecules tended to form spherical clusters, which was also harmful to the charge transport.In addition,the maximum values of CCL and 2Rg-PCBMfor PC71BM were obtained, indicating that stronger PC71BM aggregation was formed in this case,which resulted in a relatively large phase separation in the active layer. Moreover, a stronger phase separation was also demonstrated by the larger 2Rg-polymervalue. Thus,exciton dissociation was harmed because of the lack of a donor/acceptor interface caused by the strong phase separation. Owing to these two factors, devices without the addition of DIO showed significantly lower JSCand FF values than devices with added DIO, leading to an apparently lower PCE of 5.49% (Fig. 3). After the addition of DIO(Fig. 7b),the value of CCL for PTB7-Th increased to 5.70 nm. The value of D decreased to 2.04,which implied that molecules preferred to form planar aggregations accompanied by a large crystal size, which was beneficial to the improvement of charge transport efficiency, hence,increasing JSCand FF. Simultaneously, the values of CCL for PC71BM, 2Rg-PCBM, and 2Rg-polymerdecreased with the addition of DIO, implying the formation of weaker phase separation with smaller aggregation clusters, which was also favorable for the enhancement of exciton dissociation and the improvement of JSCand FF by forming more donor/acceptor interfaces.Therefore,the PCE of the device can be improved (Fig. 3).

    Fig. 7 (Color online).Schematic diagrams of the 3D nanostructures drawn according to the results of grazing incidence wide-angle X-ray scattering and GISAXS measurements for blend films with a 0 and b 3 vol% 1,8-diiodoctane. DIO, 1,8-diiodoctane; PTB7-Th,poly[(ethylhexyl-thiophenyl)-benzodithiophene-(ethylhexyl)-thienothiophene]; PC71BM, 6,6-Phenyl C71 butyric acid methyl ester

    3.3 In situ investigation of morphology evolution and correlations with the degradation of thermal stability

    In situ GISAXS measurements were performed for PTB7-Th/PC71BM blend films with 0 and 3 vol% DIO to explore the large-scale morphological evolution during the thermal annealing process,as shown in Figs.S4 and S5.No significant change in the 1D curves was observed during thermal annealing of the blend films, regardless of the proportion of DIO added. This indicates that the effect of short-term thermal annealing on the aggregation of amorphous molecules in the active layer was insignificant. The size of the aggregation, size of the fractal network, and fractal dimension were almost unchanged. No large-scale phase separation was induced after thermal annealing,which is consistent with the results in Fig. 4b, where the values of VOCremained approximately constant for devices after thermal annealing, regardless of the addition proportion of DIO. Therefore, the degradation of the device PCE was attributed to the small-scale morphological evolution,which was closely related to the JSCand FF of the devices.

    To explore the mechanism of increasing DIO harming device thermal stability, in situ GIWAXS measurements were performed during the thermal annealing process for PTB7-Th/PC71BM blends with various ratios of DIO. To simulate the conventional thermal annealing strategy that directly places devices on a hot plate that has been heated to a specified temperature, an annealing strategy that placed devices on the hot plate at room temperature(25 °C), which was then heated to the annealing temperature (100 °C) at a heating rate of 20 °C/min, was employed. Figure 8 shows the 1D profiles of GIWAXS during the thermal annealing process obtained in the outof-plane direction for blend films with different ratios of DIO. A similar transition in the (100) peak was observed for all blend films, regardless of the proportion of DIO.Taking the blend film with 3 vol% DIO as an example, as the annealing time increased, the FWHM of the original(100)peak in the out-of-plane direction gradually increased from 0.99 to 1.23 nm-1, and the peak intensity gradually decreased, which led to the gradual disappearance of the original(100)peak.Simultaneously,a newly formed(100)peak occurred at the position with a higher q value(q = 3.04 nm-1)accompanied by a gradual decrease in the FWHM to 0.78 nm-1, and peak intensity gradually increased. The detailed values of the peak position and FWHM are summarized in Table S3-S6. This indicates that PTB7-Th crystallization could be refactored during the temperature holding process, where the original crystallization was destroyed and new crystallization was formed,accompanied by molecular nucleation and crystal growth processes.It can be concluded from Fig. 8 that the transition of the(100)peak became more drastic as the proportion of DIO increased. In addition, as shown in Table S7-S10, it remained approximately unchanged for both the peak width and peak position of PC71BM during thermal annealing,which implies that the effect of short-term thermal annealing on the PC71BM crystal could be ignored.

    Fig. 8 (Color online) The 1D profiles of grazing incidence wideangle X-ray scattering obtained in the out-of-plane direction for blend films with a 0,b 1,c 3,and d 8 vol%1,8-diiodoctane during thermal annealing at 100 °C for 10 min, where 0 represents the temperature before annealing (25 °C). DIO, 1,8-diiodoctane

    Fig. 9 (Color online) Normalized a lattice space and b crystal coherence length obtained from the (100) peak in the out-of-plane direction for blends with various ratios of 1,8-diiodoctane during thermal annealing at 100 °C. ‘‘0’’ represents the value obtained at 25 °C.The heating and cooling rate was 20 °C/min.d-spacing,lattice space; CCL, crystal coherence length; DIO, 1,8-diiodoctane

    Table 4 Values of crystal coherence length and lattice space of(100)peaks integrated from the out-of-plane direction after thermal annealing at 100 °C for 10 min

    The evolution of the normalized d-spacing integrated from the out-of-plane (100) peak for blend films with various ratios of DIO during thermal annealing at 100 °C is shown in Fig. 9a. The detailed values of d-spacing and CCL after thermal annealing are summarized in Table 4.As shown in Fig. 9a,the same trends occurred,wherein the values of d-spacing increased initially and then decreased for all blend films with various ratios of DIO. The difference was that both the increasing and decreasing rates of dspace increased as the proportion of DIO increased. The maximum values of d-spacing were 102.8%, 102.9%,103.2%, and 103.4% of their initial values, and the final values were 99.3%, 98.0%, 97.3%, and 96.5% of their initial values for blends with 0, 1, 3, and 8 vol% DIO,respectively.It was indicated that the blend film with more DIO added would suffer more drastic evolution in the dspacing. Moreover, the time at which d-spacing started to decrease and CCL started to increase was advanced as the proportion of DIO increased, where the time was approximately 2 min earlier for blends with 8 vol% DIO (about 8 min after thermal annealing started) than for the blends without DIO (approximately 10 min after thermal annealing started), which also suggested that the addition of DIO could facilitate the morphological evolution during the annealing process.

    The evolution of the normalized value of CCL integrated from the out-of-plane (100) peak for blend films with various ratios of DIO during thermal annealing at 100 °C is shown in Fig. 9b.The same trends wherein CCL first decreased and then increased were obtained from blend films, regardless of the amount of DIO. Similarly,larger increasing and decreasing rates of CCL were observed as the proportion of DIO increased. The minimum values of CCL were 82.1%, 81.4%, 81.0%, and 79.5%,and the final values were 111.4%,118.9%,126.6%,and 137.8% of their initial values for blends with 0, 1, 3,and 8 vol% DIO, respectively. These results indicate that the morphological stability of the blends decreased with an increase in the proportion of DIO, which behaved as a faster morphological evolution by a large margin,as shown in Figs. 8 and 9, resulting in more severe degradation of device thermal stability as more DIO was added (Fig. 3).

    In the heating stage, as shown in Fig. 9, the CCL first decreased, which resulted in the reduction of charge transport channels.Simultaneously,a looser packing with a larger d-spacing value was formed due to the destruction of the original crystal, which allowed the PC71BM molecules to diffuse into the polymer domains to form more interfaces between the donor and acceptor, thereby enhancing exciton dissociation. In the temperature holding stage, the CCL continued to decrease at a relatively slower rate compared to the decreasing rate in the heating stage.Meanwhile,the d-spacing increased at a slower rate.These results demonstrated that the destruction of polymer crystalline phase was still in progress, while the slower evolution slowed the decay rate of charge transport, which may have balanced the effects of weakened charge transport and improved the donor/acceptor interface. Subsequently,as the annealing time increased,the value of CCL increased, accompanied by newly formed charge transport channels, which facilitated the electron mobility and hole mobility in the active layer, thereby enhancing the device PCE and decreasing the recombination loss. However, the interface between the donor and acceptor was weakened owing to the enlarged phase separation due to crystal growth. In the cooling stage, the CCL continued to increase,while the growth rate gradually slowed.However,the final crystal size became too large, resulting in unfavorable phase separation that harmed the photogenerated exciton dissociation, thereby decreasing JSCand FF.

    When the devices were thermally annealed, the varying rates of CCL and d-spacing in all stages increased as the proportion of DIO increased, which suggested that the addition of DIO causes the active layer to become unstable and facilitate the rate of morphological evolution. The blends with more DIO suffered a more drastic morphological evolution,which was apparently unfavorable for the thermal stability of the device.Therefore,the residual PCE gradually decreased as the proportion of DIO increased(Fig. 4).

    3.4 Mechanism of analysis of the effect of DIO on the morphology-performance relationship

    Figure 10 shows a schematic diagram of the morphological evolution of blends with 0 and 3 vol% DIO during the thermal annealing processes.DIO is a selective solvent additive that can selectively dissolve PC71BM [28]. In addition, the boiling point of DIO (332 °C) [29] is significantly higher than that of the solvent chlorobenzene(130.9 °C) [30]. As shown in Fig. 10a, when no DIO was added, large-scale sphere PC71BM clusters were formed,leading to a strong phase separation between PTB7-Th and PC71BM, which was demonstrated by the larger values of 2Rg-PCBMand 2Rg-polymerlisted in Table 3.When DIO was added (Fig. 10b), part of the clustered PC71BM was dissolved into discrete molecules,and the shape of the cluster became planar, accompanied by decreased values of CCL for PC71BM and 2Rg-PCBM(Fig. 5d and Table 3). Discrete PC71BM molecules could diffuse into the polymer-rich domains and even the crystal domains of PTB7-Th during the film-forming process,owing to the slower volatilization of DIO,leading to an increase in the d-spacing of PTB7-Th as the ratio of DIO increased (Fig. 5c). Subsequently, a heterogeneous nucleation and crystal growth process was observed, where the discrete PC71BM molecules had an effect similar to that of the nucleating agent for PTB7-Th.The d-spacing, and crystal size of PTB7-Th increased as the proportion of DIO increased (Fig. 5c). When the DIO ratio was increased to 8 vol%, many PC71BM molecules were dissolved in DIO. Owing to the higher boiling point of DIO, chlorobenzene volatilized before DIO during the film-forming process. Afterward, the excessive PC71BM molecules dissolved in DIO formed a new large-scale PC71BM cluster during the volatilization of DIO,leading to an enlarged CCL of PC71BM, as shown in Fig. 5d.

    Fig.10 (Color online)Schematic diagram of the morphological evolution of blends with 0 and 3 vol%1,8-diiodoctane during thermal annealing.DIO, 1,8-diiodoctane; PC71BM, [6,6]-phenyl C71-butyric acid methyl ester

    For the morphological evolution during thermal annealing, there were two stages; the crystal destruction process and recrystallization process, which corresponded to the CCL decreasing and increasing stages, respectively,as shown in Fig. 9b. Organic photovoltaic molecules have long flexible side chains, allowing molecular chains to relax at temperatures below the glass transition temperature(Tg) [31], which is approximately 125 °C for PTB7-Th[32]. When the annealing process started, the original crystal was destroyed first because of the strong relaxation of the molecules at this temperature. The addition of DIO can decrease the Tgof PTBT-Th[32];therefore,a stronger molecular relaxation occurred at the annealing temperature for blends with added DIO,which led to a more significant decrease in CCL in the first stage of thermal annealing(Fig. 9). In the second stage, molecular interaction and movement caused amorphous molecules to nucleate and crystallize again to form new crystals with gradually decreased d-spacing and increased CCL. For blends with DIO, there were many discrete PC71BM molecules distributed in the polymer domain, which could facilitate the nucleation and crystallization of PTB7-Th (Fig. 10b).Hence, the crystal growth rate and growth starting time both increased as the proportion of DIO increased(Fig. 9).Therefore, more drastic morphological evolutions induced when the proportion of DIO was increased resulted in morphological instability and subsequently, led to a relatively poor thermal stability compared to that of devices without DIO.

    4 Conclusion

    In this study, the effect of the DIO additive on the PCE and thermal stability of the PTB7-Th/PC71BM device was investigated in depth using an in situ GIWAXS technique that we developed in the beamline BL16B1 at SSRF. The PCE first increased and then decreased as the proportion of DIO increased,achieving the highest value of 9.17%with 3 vol% DIO. In addition, gradual degradation in thermal stability with an increasing DIO proportion occurred,where the residual PCEs after thermal annealing were 90.4%, 77.6%, 75.8%, and 71.4% for devices with 0, 1, 3,and 8 vol% DIO, respectively. GIWAXS and GISAXS measurements showed that the device with 3 vol% DIO could form a favorable active layer morphology with appropriate donor crystal size, aggregation size of amorphous molecules, and phase separation, which offered plenty of donor/acceptor interfaces and charge transport channels to enhance photogenerated exciton dissociation and electron and hole transport. Hence, the optimum PCE was obtained for the device with 3 vol% DIO. However,when the proportion of DIO was 8 vol%, excessive PC71-BM was dissolved in DIO, which would form larger aggregates during the film-forming process to harm the device PCE. In situ GIWAXS measurements revealed that the rate of morphological evolution increased as the proportion of DIO increased. The acceptor of PC71BM was partly dissolved by DIO into discrete molecules, which diffused into polymer-rich domains and had an effect similar to that of the PTB7-Th nucleating agent,leading to enhanced crystallization. However, a fast morphological evolution was also induced by the addition of DIO during thermal annealing, thereby resulting in poorer thermal stability. These results also showed the significance,importance, and feasibility of the in situ GIWAXS technique that we developed at BL16B1.

    AcknowledgementsThe authors thank the colleagues at the SAXS beamlines (BL10U1, BL19U2, and BL14B1) of the Shanghai Synchrotron Radiation Facility for their support and discussions.

    Author contributionsAll authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by Da Huang, Chun-Xia Hong, Jian-Hua Han, and Chun-Ming Yang.The first draft of the manuscript was written by Da Huang and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.

    97人妻精品一区二区三区麻豆| 国产成人a∨麻豆精品| 最近的中文字幕免费完整| 欧美高清性xxxxhd video| 男女无遮挡免费网站观看| a级一级毛片免费在线观看| 欧美老熟妇乱子伦牲交| 免费观看在线日韩| 黄片wwwwww| 成人特级av手机在线观看| 2021少妇久久久久久久久久久| 亚洲精品乱码久久久v下载方式| 在线精品无人区一区二区三 | 国产av不卡久久| 亚洲aⅴ乱码一区二区在线播放| 欧美3d第一页| 欧美高清成人免费视频www| 交换朋友夫妻互换小说| 国产欧美另类精品又又久久亚洲欧美| 国产国拍精品亚洲av在线观看| 欧美xxxx黑人xx丫x性爽| 王馨瑶露胸无遮挡在线观看| 亚洲美女搞黄在线观看| 久久鲁丝午夜福利片| 久久99精品国语久久久| 日本av手机在线免费观看| 三级国产精品片| 亚洲在线观看片| 国产有黄有色有爽视频| 国产欧美另类精品又又久久亚洲欧美| 下体分泌物呈黄色| 97超碰精品成人国产| 天天躁夜夜躁狠狠久久av| 国产亚洲91精品色在线| 亚洲成人久久爱视频| 国产综合懂色| 国产黄色视频一区二区在线观看| 久久久久久久久久久丰满| 黄色配什么色好看| 最新中文字幕久久久久| 亚洲综合精品二区| 亚洲丝袜综合中文字幕| 国产女主播在线喷水免费视频网站| 18禁裸乳无遮挡动漫免费视频 | 成人毛片60女人毛片免费| 亚洲一级一片aⅴ在线观看| 最后的刺客免费高清国语| av线在线观看网站| 大陆偷拍与自拍| 亚洲国产av新网站| 日本欧美国产在线视频| 精品久久久久久久人妻蜜臀av| 亚洲精品aⅴ在线观看| 亚洲av福利一区| 亚洲最大成人手机在线| 成人美女网站在线观看视频| 午夜福利网站1000一区二区三区| 夫妻午夜视频| 在线免费观看不下载黄p国产| 精品视频人人做人人爽| 亚洲一区二区三区欧美精品 | 色综合色国产| 七月丁香在线播放| 国产亚洲精品久久久com| 久久久a久久爽久久v久久| 精品国产一区二区三区久久久樱花 | 伦精品一区二区三区| 国产精品国产三级专区第一集| 国产精品爽爽va在线观看网站| 久久精品夜色国产| 亚洲最大成人中文| av免费在线看不卡| 久久久色成人| 国产成人午夜福利电影在线观看| 欧美日韩视频精品一区| 一区二区三区四区激情视频| 国产乱人偷精品视频| 欧美日韩综合久久久久久| 亚洲国产精品专区欧美| 国产人妻一区二区三区在| 亚洲丝袜综合中文字幕| 亚洲熟女精品中文字幕| 午夜福利视频精品| 国产色婷婷99| 一区二区三区精品91| 舔av片在线| 一级黄片播放器| 欧美老熟妇乱子伦牲交| 看非洲黑人一级黄片| 91久久精品国产一区二区成人| 国产久久久一区二区三区| 久久女婷五月综合色啪小说 | 亚洲精品国产av蜜桃| h日本视频在线播放| 九九久久精品国产亚洲av麻豆| 亚洲va在线va天堂va国产| 又黄又爽又刺激的免费视频.| 制服丝袜香蕉在线| 欧美国产精品一级二级三级 | av女优亚洲男人天堂| 秋霞伦理黄片| 97人妻精品一区二区三区麻豆| 欧美zozozo另类| av线在线观看网站| 欧美精品一区二区大全| 欧美性猛交╳xxx乱大交人| 成人二区视频| 在线a可以看的网站| 国产精品一及| 国产老妇伦熟女老妇高清| 国产久久久一区二区三区| 内射极品少妇av片p| 免费播放大片免费观看视频在线观看| 国语对白做爰xxxⅹ性视频网站| 国产精品麻豆人妻色哟哟久久| 三级国产精品欧美在线观看| 男女国产视频网站| 亚洲精品乱码久久久v下载方式| 99视频精品全部免费 在线| 97超视频在线观看视频| 秋霞在线观看毛片| 国产黄片美女视频| 日韩大片免费观看网站| 人妻少妇偷人精品九色| 搞女人的毛片| 国产欧美另类精品又又久久亚洲欧美| 国产亚洲一区二区精品| 国产高清不卡午夜福利| 久久ye,这里只有精品| 成人鲁丝片一二三区免费| 国产精品一区二区三区四区免费观看| 国产成年人精品一区二区| 简卡轻食公司| 久久精品国产亚洲av涩爱| 免费播放大片免费观看视频在线观看| 简卡轻食公司| 成人免费观看视频高清| 午夜激情久久久久久久| 成年免费大片在线观看| 人体艺术视频欧美日本| av一本久久久久| 免费看日本二区| av天堂中文字幕网| 啦啦啦在线观看免费高清www| 国产成人aa在线观看| 人人妻人人爽人人添夜夜欢视频 | 免费av不卡在线播放| 精品午夜福利在线看| 99热6这里只有精品| 十八禁网站网址无遮挡 | 精品国产乱码久久久久久小说| 国产免费一区二区三区四区乱码| 国产淫语在线视频| 日韩中字成人| 在现免费观看毛片| 看免费成人av毛片| 国产国拍精品亚洲av在线观看| 国产69精品久久久久777片| 色综合色国产| 一边亲一边摸免费视频| 久久久久久久大尺度免费视频| 大陆偷拍与自拍| 欧美另类一区| 极品少妇高潮喷水抽搐| 亚洲av欧美aⅴ国产| 色视频在线一区二区三区| 国产v大片淫在线免费观看| 高清午夜精品一区二区三区| 国产亚洲午夜精品一区二区久久 | 色吧在线观看| 午夜激情福利司机影院| 久久久久久久久久成人| 亚洲国产av新网站| 99热全是精品| 国产午夜福利久久久久久| 亚洲最大成人av| 中文字幕制服av| 精品熟女少妇av免费看| 观看美女的网站| 亚洲人成网站在线播| 亚洲第一区二区三区不卡| av在线天堂中文字幕| 日本三级黄在线观看| 成人毛片60女人毛片免费| 午夜亚洲福利在线播放| 联通29元200g的流量卡| 亚洲熟女精品中文字幕| 2021少妇久久久久久久久久久| 日日摸夜夜添夜夜添av毛片| 亚州av有码| 国产视频内射| 免费av不卡在线播放| 国产一区二区亚洲精品在线观看| 欧美成人精品欧美一级黄| av.在线天堂| av在线天堂中文字幕| 极品少妇高潮喷水抽搐| 婷婷色麻豆天堂久久| 天天躁日日操中文字幕| 97超碰精品成人国产| 精品久久久久久电影网| 51国产日韩欧美| 国产成人精品福利久久| 国产高清国产精品国产三级 | 九九在线视频观看精品| 久久综合国产亚洲精品| 男人添女人高潮全过程视频| 欧美bdsm另类| 少妇熟女欧美另类| 亚洲成人av在线免费| 久久久久久久午夜电影| 国产久久久一区二区三区| 亚洲精品亚洲一区二区| 99久久精品国产国产毛片| 久久久久久久久久成人| 久久久久久国产a免费观看| 蜜桃久久精品国产亚洲av| av网站免费在线观看视频| 91久久精品国产一区二区成人| 日韩一区二区视频免费看| 日韩精品有码人妻一区| 建设人人有责人人尽责人人享有的 | 人人妻人人看人人澡| 成人国产麻豆网| 黄色视频在线播放观看不卡| av.在线天堂| 亚洲精品久久午夜乱码| 秋霞伦理黄片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产男女超爽视频在线观看| 午夜免费鲁丝| 在线播放无遮挡| 一级毛片 在线播放| 99热这里只有是精品在线观看| 国产高清不卡午夜福利| 欧美一区二区亚洲| 人妻 亚洲 视频| 日韩不卡一区二区三区视频在线| 一级av片app| 两个人的视频大全免费| 亚洲欧美清纯卡通| 黄色配什么色好看| 少妇裸体淫交视频免费看高清| 亚洲av成人精品一区久久| 自拍偷自拍亚洲精品老妇| 国产亚洲av嫩草精品影院| 搞女人的毛片| 欧美性猛交╳xxx乱大交人| 在线观看一区二区三区激情| 亚洲欧美成人综合另类久久久| 久久久久性生活片| 国产男女超爽视频在线观看| 免费观看无遮挡的男女| 在线精品无人区一区二区三 | 最近的中文字幕免费完整| 美女脱内裤让男人舔精品视频| 超碰av人人做人人爽久久| 亚洲人成网站高清观看| 高清日韩中文字幕在线| 日本三级黄在线观看| 国产黄频视频在线观看| 国产片特级美女逼逼视频| 久久久久久久久久久免费av| 超碰av人人做人人爽久久| 日韩电影二区| 菩萨蛮人人尽说江南好唐韦庄| 天堂俺去俺来也www色官网| 永久网站在线| 18禁裸乳无遮挡免费网站照片| 欧美激情国产日韩精品一区| 国产精品精品国产色婷婷| 青春草视频在线免费观看| 搡女人真爽免费视频火全软件| 99热这里只有是精品在线观看| 国产精品一区二区性色av| 欧美成人a在线观看| 亚洲国产精品999| 777米奇影视久久| 久久精品国产鲁丝片午夜精品| 久久99蜜桃精品久久| 国产人妻一区二区三区在| 精品国产露脸久久av麻豆| 观看美女的网站| 波多野结衣巨乳人妻| 国产成年人精品一区二区| 欧美成人午夜免费资源| 精品国产乱码久久久久久小说| 少妇熟女欧美另类| 亚洲一区二区三区欧美精品 | 深夜a级毛片| 狂野欧美激情性bbbbbb| 亚洲av欧美aⅴ国产| 美女cb高潮喷水在线观看| 精品酒店卫生间| 听说在线观看完整版免费高清| 日韩电影二区| 亚洲内射少妇av| 真实男女啪啪啪动态图| 国产精品一及| 国产亚洲av片在线观看秒播厂| 美女主播在线视频| 婷婷色av中文字幕| 日本与韩国留学比较| 国产精品99久久99久久久不卡 | 韩国高清视频一区二区三区| 久久精品人妻少妇| 51国产日韩欧美| 麻豆久久精品国产亚洲av| 一区二区三区乱码不卡18| 国产欧美亚洲国产| 国产精品精品国产色婷婷| a级一级毛片免费在线观看| 大码成人一级视频| 精品久久久噜噜| 精品久久久久久电影网| 久久精品国产鲁丝片午夜精品| eeuss影院久久| 日本黄色片子视频| 亚洲成人精品中文字幕电影| 亚洲欧美日韩另类电影网站 | 亚洲熟女精品中文字幕| 免费看日本二区| 只有这里有精品99| 99热这里只有是精品在线观看| 男女边摸边吃奶| 亚洲欧美中文字幕日韩二区| 国产爱豆传媒在线观看| 国产黄色视频一区二区在线观看| 91午夜精品亚洲一区二区三区| 久久久久网色| 欧美xxⅹ黑人| 精品一区二区三区视频在线| 日韩欧美 国产精品| 最近手机中文字幕大全| 亚洲电影在线观看av| 丝袜脚勾引网站| 国产成人精品久久久久久| videossex国产| 国产毛片a区久久久久| 视频中文字幕在线观看| 尾随美女入室| 综合色丁香网| 人妻制服诱惑在线中文字幕| 亚洲,一卡二卡三卡| 麻豆成人av视频| 舔av片在线| 午夜免费男女啪啪视频观看| 激情五月婷婷亚洲| 国产成人午夜福利电影在线观看| 青春草亚洲视频在线观看| 成人午夜精彩视频在线观看| 啦啦啦在线观看免费高清www| 亚洲人成网站高清观看| 又粗又硬又长又爽又黄的视频| 国产成人a区在线观看| 日日撸夜夜添| 新久久久久国产一级毛片| 亚洲怡红院男人天堂| 日本三级黄在线观看| 亚洲怡红院男人天堂| 成人二区视频| 少妇丰满av| 精品99又大又爽又粗少妇毛片| 日本色播在线视频| 熟妇人妻不卡中文字幕| 国产伦精品一区二区三区视频9| 成年av动漫网址| 五月伊人婷婷丁香| 欧美成人午夜免费资源| 亚洲av男天堂| 国产免费又黄又爽又色| 久久久久久国产a免费观看| 久久久精品免费免费高清| 一边亲一边摸免费视频| 久久久精品免费免费高清| 亚洲欧洲日产国产| 国产精品不卡视频一区二区| 日本av手机在线免费观看| 国产成人freesex在线| 日韩成人av中文字幕在线观看| 十八禁网站网址无遮挡 | 久久久久久久午夜电影| 日日摸夜夜添夜夜添av毛片| 91在线精品国自产拍蜜月| 看免费成人av毛片| 久久午夜福利片| 成人亚洲欧美一区二区av| 18禁裸乳无遮挡免费网站照片| 一级毛片黄色毛片免费观看视频| 2021少妇久久久久久久久久久| 日本黄大片高清| 亚洲国产精品成人久久小说| 国产黄色免费在线视频| 嘟嘟电影网在线观看| 成人黄色视频免费在线看| 国产精品伦人一区二区| 免费播放大片免费观看视频在线观看| 亚洲天堂国产精品一区在线| 亚洲国产欧美在线一区| 小蜜桃在线观看免费完整版高清| 久久人人爽av亚洲精品天堂 | 在线观看一区二区三区| 新久久久久国产一级毛片| 丰满人妻一区二区三区视频av| 国产精品麻豆人妻色哟哟久久| 国产男人的电影天堂91| 亚洲av中文字字幕乱码综合| 国产精品女同一区二区软件| 成人免费观看视频高清| 高清欧美精品videossex| 身体一侧抽搐| 在线免费观看不下载黄p国产| 麻豆久久精品国产亚洲av| av卡一久久| 久久久久国产网址| 天天躁日日操中文字幕| 日本熟妇午夜| 久久影院123| kizo精华| 69av精品久久久久久| 亚洲经典国产精华液单| 大香蕉97超碰在线| 在线看a的网站| 国产黄色视频一区二区在线观看| 91aial.com中文字幕在线观看| 香蕉精品网在线| 亚洲美女视频黄频| 午夜精品一区二区三区免费看| 国国产精品蜜臀av免费| 欧美日本视频| 女人十人毛片免费观看3o分钟| 成年版毛片免费区| 99久国产av精品国产电影| 国内少妇人妻偷人精品xxx网站| 韩国av在线不卡| 免费观看a级毛片全部| 美女视频免费永久观看网站| 美女国产视频在线观看| 美女被艹到高潮喷水动态| 国产精品精品国产色婷婷| 国产一区二区三区综合在线观看 | 新久久久久国产一级毛片| 在线观看三级黄色| 亚洲精华国产精华液的使用体验| 女人十人毛片免费观看3o分钟| 三级男女做爰猛烈吃奶摸视频| 最近的中文字幕免费完整| 午夜福利在线在线| 亚洲人成网站在线播| 一级毛片 在线播放| 免费黄频网站在线观看国产| 丰满乱子伦码专区| 欧美激情久久久久久爽电影| av国产久精品久网站免费入址| 身体一侧抽搐| 久久精品久久久久久久性| 最新中文字幕久久久久| 国产成人精品一,二区| 欧美少妇被猛烈插入视频| 在线精品无人区一区二区三 | 国产亚洲av片在线观看秒播厂| 亚洲内射少妇av| 91狼人影院| 一级二级三级毛片免费看| 国产在线男女| 欧美+日韩+精品| 久久精品人妻少妇| 男人和女人高潮做爰伦理| 插逼视频在线观看| 亚洲欧美日韩卡通动漫| 91狼人影院| 麻豆国产97在线/欧美| 亚洲图色成人| 国产伦精品一区二区三区视频9| 欧美丝袜亚洲另类| 久久这里有精品视频免费| 国产亚洲91精品色在线| 国产亚洲av片在线观看秒播厂| 久久99热这里只有精品18| 欧美极品一区二区三区四区| 搡女人真爽免费视频火全软件| 欧美成人午夜免费资源| 亚洲综合色惰| 新久久久久国产一级毛片| 国产一区二区三区av在线| 久久国内精品自在自线图片| 久久久久久久久久久免费av| 视频中文字幕在线观看| 久久久a久久爽久久v久久| av国产免费在线观看| 欧美丝袜亚洲另类| 国产免费一级a男人的天堂| 黄色视频在线播放观看不卡| 亚洲精品aⅴ在线观看| 在线 av 中文字幕| 日本黄大片高清| 国内精品宾馆在线| 纵有疾风起免费观看全集完整版| 一本久久精品| 日本-黄色视频高清免费观看| 国产淫片久久久久久久久| 在线观看av片永久免费下载| 亚洲精品国产av成人精品| 高清午夜精品一区二区三区| 亚洲四区av| 欧美三级亚洲精品| 国内精品宾馆在线| 中文字幕av成人在线电影| 人妻系列 视频| 国产黄频视频在线观看| 秋霞在线观看毛片| 久久午夜福利片| 99久久九九国产精品国产免费| 精品少妇久久久久久888优播| 国产精品99久久久久久久久| av国产精品久久久久影院| 久久久久久久久久人人人人人人| 婷婷色综合www| 日韩人妻高清精品专区| 三级国产精品片| av在线观看视频网站免费| 美女视频免费永久观看网站| 国产精品成人在线| 欧美成人精品欧美一级黄| 日本wwww免费看| 新久久久久国产一级毛片| 久久久久久久久大av| 日本一本二区三区精品| 一级毛片电影观看| 国产成人一区二区在线| 视频中文字幕在线观看| 色5月婷婷丁香| 欧美zozozo另类| 老司机影院毛片| 禁无遮挡网站| 91午夜精品亚洲一区二区三区| 成人二区视频| 亚洲国产精品专区欧美| 在线观看人妻少妇| 久久人人爽人人爽人人片va| 亚洲在线观看片| 亚洲最大成人av| 少妇高潮的动态图| 日本欧美国产在线视频| 亚洲精品视频女| 亚洲国产成人一精品久久久| 亚洲一级一片aⅴ在线观看| 男女无遮挡免费网站观看| 亚洲一级一片aⅴ在线观看| 人妻 亚洲 视频| 亚洲av电影在线观看一区二区三区 | 亚洲欧美一区二区三区国产| 欧美精品人与动牲交sv欧美| 久久久久网色| 亚洲激情五月婷婷啪啪| 性色av一级| 搡女人真爽免费视频火全软件| 又黄又爽又刺激的免费视频.| 亚洲精品国产色婷婷电影| 国产欧美亚洲国产| 2021少妇久久久久久久久久久| 神马国产精品三级电影在线观看| 中文字幕免费在线视频6| 久久久久久伊人网av| 老师上课跳d突然被开到最大视频| 成人国产av品久久久| 久久久久性生活片| 蜜臀久久99精品久久宅男| 成年av动漫网址| 久久人人爽人人爽人人片va| 色5月婷婷丁香| 精品久久久久久久末码| 成人毛片a级毛片在线播放| a级毛色黄片| 国产一区二区在线观看日韩| 亚洲av在线观看美女高潮| 国产男女内射视频| 久久久欧美国产精品| 国产精品精品国产色婷婷| 熟女av电影| 久久精品久久精品一区二区三区| 人人妻人人澡人人爽人人夜夜| 久久久久精品性色| av免费观看日本| 免费少妇av软件| 久久人人爽av亚洲精品天堂 | 亚洲国产高清在线一区二区三| 99热6这里只有精品| 日本午夜av视频| 亚洲精品中文字幕在线视频 | 国产色爽女视频免费观看| 亚洲一级一片aⅴ在线观看| 男人添女人高潮全过程视频| 一级a做视频免费观看| 午夜免费男女啪啪视频观看| 久久久久久久精品精品| 老司机影院毛片| 麻豆成人午夜福利视频| 国产精品偷伦视频观看了| 内地一区二区视频在线| 免费不卡的大黄色大毛片视频在线观看| 久久精品熟女亚洲av麻豆精品| 青春草亚洲视频在线观看| 伦理电影大哥的女人| 国产亚洲av嫩草精品影院| 亚洲精品亚洲一区二区| 特大巨黑吊av在线直播| 日韩三级伦理在线观看| 日韩欧美一区视频在线观看 | 91久久精品电影网| 一级av片app| 国产成人免费观看mmmm| 80岁老熟妇乱子伦牲交| 久久6这里有精品| 在线看a的网站|