Reconstructing Space- and Energy-Dependent Exciton Generation in


Reconstructing Space- and Energy-Dependent Exciton Generation in...

0 downloads 50 Views 1014KB Size

Subscriber access provided by UNIVERSITY OF TOLEDO LIBRARIES

Organic Electronic Devices

Reconstructing Space- and Energy-dependent Exciton Generation in Solution-processed Inverted Organic Solar Cells Yuheng Wang, Yajie Zhang, Guanghao Lu, Xiaoshan Feng, Tong Xiao, Jing Xie, Xiaoyan Liu, Jiahui Ji, Zhixiang Wei, and Laju Bu ACS Appl. Mater. Interfaces, Just Accepted Manuscript • DOI: 10.1021/acsami.7b14698 • Publication Date (Web): 28 Mar 2018 Downloaded from http://pubs.acs.org on March 30, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Reconstructing Space- and Energy-dependent Exciton Generation in Solution-processed Inverted Organic Solar Cells Yuheng Wang1, Yajie Zhang2, Guanghao Lu1,*, Xiaoshan Feng3, Tong Xiao1, Jing Xie4, Xiaoyan Liu5, Jiahui Ji5, Zhixiang Wei2,*, Laju Bu3,*

1 Frontier Institute of Science and Technology, Xi’an Jiaotong University, Xi’an 710054, China

2 CAS key laboratory of nanosystem and hierarchical fabrication, Beijing 100190, China

3 School of Science, Xi’an Jiaotong University, Xi’an 710049, China

4 School of Life Science and Technology, Xi’an Jiaotong University, Xi’an 710049, China

5 School of Electrical Engineering, Xi’an Jiaotong University, Xi’an 710049, China

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 28

KEYWORDS: light absorption, small molecule solar cells, organic photovoltaics, inverted devices, vertical phase separation, optical simulation,

ABSTRACT:

Photon absorption induced exciton generation plays an important role in determining the photovoltaic properties of donor: acceptor organic solar cells with an inverted architecture. However, reconstruction of light-harvesting and thus exciton generation at different locations within organic inverted device are still not well resolved. Here we investigate the film-depth-dependent light absorption spectra in a small molecule donor: acceptor film. Including depth-dependent spectra into an optical transfer matrix method allows us to reconstruct both film-depth- and energy-dependent exciton generation profiles, using which short-circuit current and external quantum efficiency of the inverted device are simulated and compared with the experimental measurements. The film-depth-dependent spectroscopy, from which we are able to simultaneously reconstruct light harvesting profile, depth-dependent composition distribution and vertical energy level variations, provides insights into photovoltaic process. In combination with appropriate materials processing methods and device architecture,

the

method

proposed

in

this

work

will

help

optimize

film-depth-dependent optical/electronic properties for high performance solar cells.

ACS Paragon Plus Environment

Page 3 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

1. INTRODUCTION Organic semiconductors including small molecules and polymeric semiconductors are widely investigated for next generation flexible electronics.1-7 As compared with polymeric semiconductors which are usually multi-dispersed in terms of molecular structure, small organic molecules are of great interests because they have clear chemical structure and could be easily purified for less dispersed performance.8 Currently, organic electronics based on small organic molecules are widely used in organic light emission diodes for display applications. As a comparison, the performance of the state-of-the-art small molecule photovoltaic devices is still too low to be commercially applied. In recent years, the power conversion efficiency of solution-processed small molecule organic photovoltaics has been significantly increased to be higher than 10%, demonstrating a very bright prospect for future energy conversion.8 However, unlike thermally co-evaporated donor: acceptor bulk heterojunction, vertical phase segregation of solution-processed film is commonly observed in donor: acceptor bulk heterojunction,9 which induces a significant variation of light harvesting and thus exciton generation properties along directions of film thickness.10 In order to utilize variable morphology along the vertical direction of the donor: acceptor active layer, two representative device architectures, namely conventional and inverted structures,11-12 are developed to fabricate solar cells. The two types of devices require different materials for cathode, anode, hole/electron transport layers and active layer, which significantly influence the vertical phase evolution of active layer.9 Therefore, for lots of donor: acceptor active layer, inverted

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

architecture has higher performance as compared with conventional devices. Moreover, inverted organic solar cells are usually more stable and more compatible with industrial processing than conventional devices. In order to correlate efficiency of inverted device with the vertical variation of composition and optical properties, light-harvesting and thus exciton generation at different location within the inverted device are highly required. The bulk heterojunction here is a composite of two-phase blends consisting of a p-type electron donor and an n-type electron acceptor. The p-type material needs to be connected to the anode and the n-type material needs to be connected to the cathode. The p-type material and the n-type material form an interpenetrating network in the film for effective charge transfer and charge transport, thereby achieving high power conversion efficiency. The intermixture morphology of the donor and acceptor and the compositional ratio along the vertical direction of the film directly affect the power conversion efficiency of the device because these parameters are closely related to exciton generation, energy/charge transfer and charge transport. Therefore, in order to optimize the performance of photovoltaic thin films, one needs to analyze how sunlight is attenuated in thin films, as well as the position of exciton generation, and to establish relationship between deep-resolved optical properties and photovoltaic performance.13-18 In fact, for small molecule donor: acceptor film, the variations of light absorption along film-depth direction are related to the facts as follows: 1) Donor and acceptor are usually not fully compatible, and during solution-processing e.g. spincoating, vertical separation between two components is commonly observed. 2) Molecules at

ACS Paragon Plus Environment

Page 4 of 28

Page 5 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

different film-depth might be featured with different molecular conformation and crystallization, which have an important effect on conductance/valance levels and thus light absorption spectra. 3) Donor: acceptor is generally semicrystalline, while amorphous

and

crystalline

materials

have

different

optical

properties.

Film-depth-dependent crystallinity contributes to the film-depth-dependent light absorption. 4) Molecules at the layer-layer interface might have a different molecular or crystal orientation from that inside the active layer, leading to the variable absorption spectra along film-depth. However, correlation the photovoltaic performance of such small molecules organic devices with light harvesting profiles inside the device is still not well resolved. In recent years, film-depth-dependent light absorption spectra have been used to study film-depth-dependent composition, charge transport and light propagation in polymer field-effect transistors and polymer solar cells with conventional architecture. 4, 15-17

UV-Vis absorption spectroscopy is a commonly used analytical technique to

characterize the organic solar cell. This method is easily accessible and can provide insight into the morphology and composition of the film. Oxygen plasma contains a variety of reactive components, including O2 +, O2-, O3, O, O+, O-, ionized ozone, and metastable oxygen excited states, and these highly reactive components can etch the material from the surface of the organic film. Therefore, if the oxygen plasma pressure is high, these active components will diffuse below the surface of the film,15 which may cause chemical reaction with the internal components of the film and change the optical absorption characteristics of the material inside the film, leading to

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 28

many difficulties for resolving the absorption spectrum of sub-layer film. It was found that at low pressure (< 30 Pa), oxygen plasma could incrementally etch the organic film without damaging the materials under the surface.15-16 Therefore, the light absorption spectra of the film could be monitored during the plasma etching. From Beer-Lambert’s law, the light absorption spectra at different film-depth can thus be derived. In this work, we investigate the film-depth-dependent light absorption spectra in a small molecule donor: acceptor film which is subsequently used in inverted solar cells with a high power-conversion efficiency approximately 8%. Using such spectra, composition profile along the vertical direction is obtained. Taking the optical properties of different layers including electrodes and hole/electron transport layers into account, we include the film-depth-dependent absorption spectra of the active layer into an optical interference model to simulate the film-depth-dependent and energy-dependent exciton generation profiles, which were subsequently used to simulate short-circuit current and external quantum efficiency. Our numerical simulations are well consistent with the experimentally measured photovoltaic performance.

2. EXPERIMENTAL SECTION 2-(thiophen-2-yl)thieno[3,2-b]thiophene

as

p-bridges

and

end-capped

with

1H-indene-1, 3(2H)-dione (abbreviated as BTID) was synthesized as reported in our previous work.8 Inverted devices with a structure of glass/ITO/ZnO/active layer/MoOx/Ag is

ACS Paragon Plus Environment

Page 7 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

prepared.8 A layer of ZnO precursor was spin-coated onto the ITO. The mixture of BTID (Figure 1) and PC71BM with total concentration ca. 18 mg/ml was used for subsequent spin-coating at 40 oC. Finally, 5 nm MoOx and 100 nm Ag layer were evaporated to complete the device fabrication. The J-V curves were measured using Keithley 2400 under AM 1.5G illumination (Newport Thermal Oriel 91159A solar simulator). EQEs were performed with an Oriel Newport system (Model 66902). The film-depth-dependent light absorption spectra were measured in a home-made setup. We integrated the plasma etching machine (PT-15S, SANHOPTT 300 W) with the ultraviolet-visible (UV-Vis) optical fiber spectrometer (PG2000-Pro, Ideaoptics) to ensure that the light is vertically aligned with the film plane and that the in-situ absorption spectra can be acquired during incremental etching. The size of observed region during the absorption measurement is 1 cm (diameter), and we use condenser lens to collect the transmitted light and scattered light. Atomic force microscopy (AFM) with tapping mode (INNOVA, Veeco) is used to characterize the photovoltaic film before and after etching (Figure S1). Raman spectra were acquired at the excitation wavelength 532 nm (LabRAM HR Evolution, HORIBA) (Figure S2). 3. RESULTS AND DISCUSSION In this work we reconstruct the light-induced exciton generation at different positions inside the inverted solar cells where both material and optical properties are very variable along the vertical direction of the donor: accept active layer. The

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 28

continuous Beer-Lambert law and the simplified optical model are used to provide insights into thin-film morphology and light-material interactions, which finally determines the device performance. According to Beer-Lambert's law, the absorbance of the film is defined as

A = log( I / I 0 ) , where I and I 0 are intensity of transmitted and incident light, respectively. We can assume that the whole film is composed of n layers (Figure 1), and each layer has an absorbance A i (i = 1 , 2 , 3 , L , n ) . Therefore, the intensity of incident light after passing through n sublayers is obtained as: n

n

I = I 0 ∏ 10 − Ai = I 0 10



∑ Ai i =1

(1)

i =1

This shows that the total absorbance of the film is equal to the linear superposition of the absorbance of each layer. As is shown in Figure 2(a), through surface-selective etching without underneath damage (Figure S2),15 plasma acts on the film surface and the film thickness decreases gradually. Afterwards, we get in the etching process a series of spectra according to the absorbance of the surface and the reduced amount of each sublayer that is further analysis of each layer component. A series of absorption spectral data obtained before and after etching are subtracted layer by layer, getting each layer’s absorption intensity data as ∆ Ai = Ai +1 − Ai (i=1, 2, 3, L , n-1) , which is shown in Figure 2(b). When the absorbance of each layer is subtracted, the peak dislocation is obviously observed, and the effect of vertical distribution of the component on the active layer is further verified. The variation of the absorption peak along film-depth direction is due to the different composition or crystallinity at different position. The proportion distribution of donor and acceptor in each layer is

ACS Paragon Plus Environment

Page 9 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

then calculated by least square method. Figure 2(c) describes vertical distance distribution of blend composition, obtained from the absorbance of pure samples and non-etched film as shown in Figure 2(d). This vertical phase separation of acceptor and donor reveals the exciton formation site indirectly. Figure 2(e) shows the dependence of absorption on etching time for different wavelengths (330 nm, 480 nm and 600 nm, respectively). Figure 2(f) shows the approach for fitting the sublayer absorption spectrum of the blend by the spectra of pure BTID and pure PC71BM. In principle, the thickness of the film is positively correlated with the intensity of the absorption spectrum. Therefore, we can establish a quantitative relationship between the intensity of the absorption spectrum and the thickness of the film by using a numerical method, according to the change of the absorption spectrum during etching. This means that we do not need a separate film to measure the film thickness after each etching. Since the sample thickness is known, the spectral constants can be used to calculate the optical constants of the samples. For the films, the refractive index and extinction coefficient are satisfied with K-K consistency principle.19 Because of the interference effect, the interference peaks existing in the transmission and reflection spectra (T and R) of the films, the curve T/(1-R) shows a monotonic increasing trend. The absorption coefficient α of the films can be calculated by the monotonicity. Thus, when the wavelength of the incident light is λ, the absorption intensity A and extinction coefficient k of the film are approximately linear with the film thickness d as equations (2)-(5) describe.

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1 d

α =- ln(

α=

T ) 1− R

4π k

λ

I ≈ I 0e−α d

∆Ai (λ ) = − lg(e −4π ki d / λ )

Page 10 of 28

(2) (3) (4) (5)

For organic solar cells, the influence of the refractive index n on the absorption intensity is not large.20 The extinction coefficient k of the sublayer in the film fluctuates in the film thickness direction, which is due to the inhomogeneous distribution of the film in the film-depth direction.6,20 Therefore, we use the transfer matrix method more precisely for the optical simulation of the device. According to the fluctuation of extinction coefficient at each thickness of the film, the active layer film is divided into several sublayers. Each sublayer has a different extinction coefficient k, and we assume that the interior of each sublayer is homogeneous.16 This is the premise of optical analysis and numerical simulation. After the film is divided into multiple sublayers, we assume that each sublayer has a homogeneous medium and is isotropic. Since the equation of the electric field propagation is linear, and the tangential component of the electric field is continuous, the propagation of the photoelectric field in the layered structure of the parallel plane can be described by a 2*2 matrix.6 As is shown in the Figure 1(c), the active layer is divided into a plurality of sublayers, each having a thickness of dj. The optical properties of each layer are determined by the complex refractive index Nj=nj+ikj. Complex refractive index is a complex parameter consisting of refractive index n and

ACS Paragon Plus Environment

Page 11 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

extinction coefficient k. The photoelectric field at any point in the layer j consists of two components, one component propagating forward along the X axis and the other propagating back along the X axis, denoted respectively as Ej+ and Ej-. The incident matrix for each interface is: 6

I jk =

1 t jk

1 r  jk

r jk  1 

(6)

where tjk and rjk are the final transmission coefficient and reflection coefficient between layers j and k, respectively. The phase shift occurs when the light passes through the layered interface, and the phase shift matrix of the light at layer j is:

 e −ξ j d j Lj =   0

0   ξ d e j j 

(7)

where, ξj = (2π/λ)Nj, ξjdj is the phase change of light passing through layer j, which is positively related to the layer thickness. Thus, the incident matrix and the phase shift matrix of the light are known, and the total transfer matrix from the first layer to the substrate is denoted as S, then the photoelectric field of the whole device can be expressed as:  E 0+   E m+ + 1   − = S  −  E0   E m +1 

(8)

The total transfer matrix S is multiplied by the incident matrix and the phase shift moment of each layer:

S S =  11  S 21

S12   m  =  Π I (ν −1)ν Lν  ⋅ I m ( m +1)  = 1 ν S 22   

(9)

In order to calculate the internal electric field in layer j, the whole layer system can

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 28

be divided into two parts separated by layer j, one from the external environment to the j layer, and one from the substrate to the j layer. Then the total system transfer matrix can be written as a product of both sides:

S = S'j Lj S''j

(10)

S' S 'j =  'j11  S j 21

S 'j12   j −1 L = ΠI S 'j 22   ν =1 (ν −1)ν ν

 S ''j11 S =  ''  S j 21

S ''j12   m  L ⋅ I = Π I S ''j 22   ν = j +1 (ν −1)ν ν  m ( m +1)

'' j

  ⋅ I (j−1) j 

(11)

(12)

Finally, we obtain the intensity of the photoelectric field at any point x in the layer j represented by an electric field formed by the incident wave. −iξ j (d j −x)

Ej (x) =

S''j11 ⋅ e

−iξ j d j

S S ⋅e ' '' j11 j11

iξ j (d j −x)

+ S''j21 ⋅ e

iξ j d j

+ S S ⋅e ' '' j12 j 21

E0+

(13)

Simulation devices were fabricated with a structure of glass/ITO/ZnO/active layer/MoOx/Ag. The active layer is BTID: PC71BM blend (Figure 1). The length parameter and the optical characteristics of other parts are obtained from the optical database. The absorption spectra of the sublayers data measured by etching apply to numerical simulation of photovoltaic devices. We make further efforts to ensure that the result is accurate. Specifically, we put the sublayer spectra obtained by etching into the model to analyze, and then obtain more accurate characteristics of photovoltaic devices by calculation. The experimental design, structural design and construction, and process simulation algorithm provide some references for structural

ACS Paragon Plus Environment

Page 13 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

design of solar cell materials. All in all, the vertical distribution characteristics of the optical constants of materials are fitted by K-K relation,19 and the photoelectric field distribution in the device is calculated by transfer matrix method according to the optical constants. Here, we directly plug in the normal solar test spectrum (AM1.5). Based on the distribution of the photoelectric field, the energy distribution inside the active layer and the density distribution of the exciton are calculated. Then, the short-circuit current density Jsc, the reflection of the device and the absorption of the internal layers and the external quantum efficiency (EQE) of the device are calculated, experiment measured data and simulation result are shown in the Figure 3(b). At around 700 nm, there is some light scattering although we used a series of condenser lens to collect the light, leading to a background value of the simulated EQE. Through the comparison of the experimental values of the device performance parameters Jsc and EQE (Figure 3), the results are basically consistent, so that the accuracy of the simulation is verified by using the multi-layer structure model. Figure 4(a) is TMM simulation of depth-dependent active inner layer light distribution, and Figure 4(b) is energy dissipation contour. Simulation is used to obtain the variation law of the square of the normalized field intensity modulus with the depth distribution of the thin-film device at different wavelengths of light. 6 The calculation of the distribution of light field inside the optical device plays an important role in the study of the internal mechanism of the device. In the film, exciton generation and dissociation occur, resulting in the conversion of solar energy to electrical energy. The absorbed photonic energy is a function of the optical field,

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 28

and the equation of which is: 21-22

Q ( x) =

2πε 0ckn E ( x)

2

λ

(14)

where |E(x)|2 is the square of the optical field modulus, n is the refractive index, k is the extinction coefficient, λ is the wavelength of incident light, c is the speed of light in vacuum, and ɛ0 is the dielectric constant. The equation for the rate of exciton production is: 21-22

G (x) =



800 nm 300nm

λ hc

Q ( x, λ )d λ (15)

According to the equation, the distribution of exciton production rate is calculated, as shown in the Figure 5(a). Subsequently, the short-circuit current density (Jsc) of the solar cell is calculated. For clarity, here we assume that one photon generates one exciton and is immediately translated into one free electron and one free hole, both of which are finally collected by both side electrodes. Internal quantum efficiency (IQE) is set as 1 for each wavelength.23 The device is assumed to be working under AM 1.5 condition and is consistent with the actual situation. The short-circuit current calculation equation is:

Jsc = ηIQE q ∫

100nm

0

G( x)dx

(16)

Where, q is elementary charge. The as-calculated Jsc for the device with 100 nm active layer is 16.0 mA/cm-2, which is close to the measured value 14.0 mA/cm-2.7 We explain that this tiny difference between simulated and measured Jsc attributed to the solar spectrum does not match the solar simulator. Meanwhile, the EQE of the device is emulated. Since IQE is set to be 1 for each wavelength, EQE is then the number of

ACS Paragon Plus Environment

Page 15 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

absorbed photons in active layer that occupies the total incident photons number proportion, as shown in the Figure 3(b) for comparison with experimentally measured data (Table S1). On the other hand, because light absorption spectra are correlated to electronic excitation from HOMO to LUMO levels, the depth-dependent light absorption profiles as shown in Figure 2b mean that the distributions of such energy levels are also film-depth dependent. From the depth-dependent light absorption spectra shown in Figure 2b, we found that the BTID component at different film-depth is featured with different absorption peak varying from ca. 500 nm at top part of the film to ca. 600 nm in the middle region and ca. 550 nm at the bottom surface. This variation is correlated to the variation of energetic distribution of localized states and thus definitely influences the exciton diffusion and free charge transport after charge transfer from donor to acceptor.24,25 Moreover, the light absorption of BTID at the top part of the film is much broader, as compared with the well-evolved absorption peaks in the bottom part of the film, implying a significant distribution of energetic disorder of localized states at each film-depth (lateral) in the top part of the film. The energetic disorder is found to induce low open-circuit voltage and fill factor.24,25 Optimizing the vertical and lateral distribution of localized states could be helpful to improve the photovoltaic performance. The exciton distribution shown in this manuscript is not the optimal case. Optimizing the sample preparation method (solvent type and its evaporation rate, processing temperature, thermal/solvent annealing or other post-treatment26) and

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

device architecture could be used to tune the exciton distribution, which could not only tune the light harvesting behavior, but also help optimize energy level distributions to facilitate charge transport. It can be considered that the plasma in-situ etching and absorption spectra can be used to analyze the sublayer of the film. Compared with other existing methods, it has the advantages of simple operation and accurate results. Sublayer absorption spectra and the optical model for active layer, compared to the traditional method based on the heterogeneity within the active layer, can make the results more accurate; meanwhile, the design of photovoltaic devices, which is expected to analyze optical multilayer structure organic solar cells, helps to further enhance the performance of the device to provide a theoretical basis and reference value. Our innovation lies in the continuous etching of the films to obtain continuous spectroscopic data of the films, and further analysis of the thin films is based on this continuous detection process. Compared to the traditional model, the innovation point of the new optical model lies in considering the heterogeneity within the active layer, thus improving the calculation accuracy. Therefore, the internal mechanism can be qualitatively and quantitatively analyzed, and the internal structure of the organic film is diagnosed. Thus, the optical analysis of multilayer structures of photovoltaic devices and the improvement of device performance of organic solar cells are solved.

4. CONCLUSION In summary, we investigated the film-depth-dependent light absorption spectra in a

ACS Paragon Plus Environment

Page 16 of 28

Page 17 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

small molecule donor: acceptor film with a high photovoltaic efficiency (~8%). Depth-dependent absorption spectra are subsequently utilized to study vertical phase segregation behavior. Including film-depth-dependent spectra of active layer into an optical interference method, with the inputs of optical properties of electrodes and hole/electron transport layers, allows us to reconstruct the film-depth-dependent and energy-dependent exciton generation profiles in an inverted solar cell, from which short-circuit current and external quantum efficiency of the device are simultaneously simulated. Our numerical simulations are well consistent with the experimentally measured photovoltaic performances. The film-depth-dependent spectroscopy, from which we are able to simultaneously reconstruct light harvesting profile, depth-dependent composition distribution and vertical variations of energy levels, provides insights into photovoltaic process. In combination with optimizing materials processing method, the approach proposed in this work helps improve the power conversion efficiency upon finely optimizing the film-depth-resolved photon harvesting profiles and simultaneously constructing appropriate vertical-variation of energy levels.

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(a)

Page 18 of 28

(b)

BTID

PC71BM

(c)

Figure 1. (a, b) Chemical structures of thiophene-substituted benzodithiophene (BTID) and [6,6]-phenyl-C71-butyric acid methyl ester (PC71BM). (c) Device configuration of the small molecule solar cell with inverted configuration. The donor: acceptor active layer is assumed to be composed of n sublayers.

ACS Paragon Plus Environment

Page 19 of 28

(a)

(b)

0.6

A

A

0.4

0.2

0.0 200

400

600

800

Wavelength (nm)

(c)

400

BTID PC71BM

800

0.6 0.5

0.6 A

0.4

0.4

0.3 Blend BTID PC71BM

0.2 0.2 0.0

600

Wavelength (nm) 0.7

0.8 content (vol %)

top to bottom

nm 0-8 8-16 16-30 30-38 38-57 57-74 74-86 86-92 92-100 0 200

(d)

1.0

0.1 0

20

40 60 Depth (nm)

80

0.0

100

(e)

400

600 Wavelength (nm)

800

(f)

0.7

Blend BTID PC71BM

0.6 0.5

Fit

A

0.4 A

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

0.3 0.2 0.1 0.0

330 nm 480 nm 600 nm 0 20 40 60 80 100 120 140 160 180 Time (s)

400

600 Wavelength (nm)

800

Figure 2. (a) Real time absorption spectra of blend film (100 nm) during etching by oxygen plasma, where the etching step is 20 s. (b) The calculated absorption spectrum of each sublayer, and the film-depth range for each sublayer is shown. (c) Distribution of BTID and PC71BM along the vertical direction as obtained from the absorption spectra. (d) Absorption spectra of BTID, PC71BM and BTID: PC71BM blends. (e) Dependence of absorption on etching time at different wavelength (330 nm, 480 nm and 600 nm) as obtained form (a). (f) Fitting the sublayer absorption spectrum of the blend by the spectra of pure BTID and pure PC71BM films. The black line is from the blend film after etching 60 s, and the line is shifted vertically for comparison.

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces

(a)

(b) 80

10

70 60

5 EQE (%)

2

Current density (mA/cm )

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 28

0 -5

50 40 30 experiment simulation

20

-10

10

-15 -0.5

0

0.0 0.5 Voltage (V)

1.0

400

500 600 Wavelength (nm)

700

Figure 3. Measured photovoltaic characteristics of BTID: PC71BM solar cell with inverted structure. (a) Measured J-V curves. (b) Measured and simulation EQE profile.

ACS Paragon Plus Environment

Page 21 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

(a)

(b)

Figure. 4 (a) Distribution of |E|2 inside the active layer. The incident solar light is from the left side. (b) Energy dissipation Q contour. (unit, mWcm-3nm-1 at its position and wavelength).

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

MoOX ZnO

Ag

Figure. 5 Simulated exciton generation rate (s-1cm-3) via reconstruction of film-depth-dependent and energy-dependent exciton generation profile.

ACS Paragon Plus Environment

Page 22 of 28

Page 23 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

ASSOCIATED CONTENT

Supporting Information.

Atomic Force Microscopy height images and Raman spectroscopy of the active layer before and after etching; Comparison of simulated and experimented results at specific wavelengths.

AUTHOR INFORMATION Corresponding Author *[email protected]

*[email protected]

*[email protected]

Notes The authors declare no competing financial interest.

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Funding Sources This work is financially supported by Natural Science Foundation of China (Grant No. 21574103) and China Postdoctoral Science Foundation (2015M580841 and 2016T90910). L.B. thanks the Fundamental Research Funds for the Central Universities. G. L. thanks Cyrus Tang Foundation.

ACKNOWLEDGMENT The authors are grateful to Dongfan Li, Zihao Wang, Jinde Yu, Tianzhu Yao, Prof. Dingxin Liu and Prof. Yongquan Qu for fruitful discussions. The authors thank Yu Wang in Instrument Analysis Center of Xi'an Jiaotong University for assistance with the Raman experiments, and Nan Zhu and Kun Gao in State Key Laboratory for Manufacturing Systems Engineering for Infrared Spectroscopy and Atomic Force Microscopy measurements. We also sincerely thank Zexin Deng for the Scanning Electron Microscopy measurement.

ACS Paragon Plus Environment

Page 24 of 28

Page 25 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

REFERENCES (1) Trost, S.; Becker, T.; Polywka, A.; Görrn, P.; Oszajca, M. F.; Luechinger, N. A.; Rogalla, D.; Weidner, M.; Reckers, P.; Mayer, T.; Riedl, T. Avoiding Photoinduced Shunts in Organic Solar Cells by the Use of Tin Oxide (SnOx) as Electron Extraction Material Instead of ZnO. Adv. Energy Mater. 2016, 6 (15), 1600347. (2) Kondo, M.; Uemura, T.; Matsumoto, T.; Araki, T.; Yoshimoto, S.; Sekitani, T. Ultraflexible and ultrathin polymeric gate insulator for 2 V organic transistor circuits.

Appl. Phys. Express 2016, 9 (6), 061602. (3) Luo, G.; Ren, X.; Zhang, S.; Wu, H.; Choy, W. C.; He, Z.; Cao, Y. Recent Advances in Organic Photovoltaics: Device Structure and Optical Engineering Optimization on the Nanoscale. Small 2016, 12 (12), 1547-1571. (4) Bu, L.; Qiu, Y.; Wei, P.; Zhou, L.; Lu, W.; Li, S.; Lu, G. Manipulating Transistor Operation via Nonuniformly Distributed Charges in a Polymer Insulating Electret Layer. Phys. Rev. Applied 2016, 6 (5). 054022. (5) Lai, L. H.; Speirs, M. J.; Chang, F. K.; Piveteau, L.; Kovalenko, M. V.; Chen, J. S.; Wu, J. J.; Loi, M. A. Increasing photon absorption and stability of PbS quantum dot solar cells using a ZnO interlayer. Appl. Phys. Lett. 2015, 107 (18), 183901. (6) Pettersson, L. A. A.; Roman, L. S.; Inganäs, O. Modeling photocurrent action spectra of photovoltaic devices based on organic thin films. J. Appl. Phys. 1999, 86 (1), 487-496. (7) Gsänger, M.; Bialas, D.; Huang, L.; Stolte, M.; Würthner, F. Organic Semiconductors based on Dyes and Color Pigments. Adv. Mater. 2016, 28 (19), 3615-3645. (8) Deng, D.; Zhang, Y.; Zhang, J.; Wang, Z.; Zhu, L.; Fang, J.; Xia, B.; Wang, Z.; Lu, K.; Ma, W. Fluorination-enabled optimal morphology leads to over 11% efficiency for inverted small-molecule organic solar cells. Nature Commun. 2016, 7, 13740. (9) Finck, B. Y.; Schwartz, B. J. Drift-Diffusion Modeling of the Effects of Structural Disorder and Carrier Mobility on the Performance of Organic Photovoltaic Devices.

Phys. Rev. Applied 2015, 4 (3), 034006.

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(10) Fabiano, S.; Himmelberger, S.; Drees, M.; Chen, Z.; Altamimi, R. M.; Salleo, A.; Loi, M. A.; Facchetti, A. Charge Transport Orthogonality in All‐Polymer Blend Transistors, Diodes, and Solar Cells. Adv. Energy Mater. 2014, 4 (6), 1301409. (11) Albrecht, S.; Schäfer, S.; Lange, I.; Yilmaz, S.; Dumsch, I.; Allard, S.; Scherf, U.; Hertwig, A.; Neher, D. Light management in PCPDTBT:PC70BM solar cells: A comparison of standard and inverted device structures. Organ. Electron. 2012, 13 (4), 615-622. (12) Jia, X.; Jiang, Z.; Chen, X.; Zhou, J.; Pan, L.; Zhu, F.; Sun, Z.; Huang, S. Highly efficient and air stable inverted polymer solar cells using LiF-modified ITO cathode and MoO3/AgAl alloy anode. ACS Appl. Mater. Interfaces 2016, 8 (6), 3792–3799. (13) Sun, Y.; Sheng, P.; Di, C.; Jiao, F.; Xu, W.; Qiu, D.; Zhu, D. Organic thermoelectric materials and devices based on p- and n-type poly(metal 1,1,2,2-ethenetetrathiolate)s. Adv. Mater. 2012, 24 (7), 932-937. (14) Guo, Y.; Di, C. A.; Ye, S.; Sun, X.; Zheng, J.; Wen, Y.; Wu, W.; Yu, G.; Liu, Y. Multibit Storage of Organic Thin-Film Field-Effect Transistors. Adv. Mater. 2010, 21 (19), 1954-1959. (15) Lu, G.; Di Pietro, R.; Kölln, L. S.; Nasrallah, I.; Zhou, L.; Mollinger, S.; Himmelberger, S.; Koch, N.; Salleo, A.; Neher, D. Dual‐Characteristic Transistors Based on Semiconducting Polymer Blends. Adv. Electron. Mater. 2016, 2 (10), 1600267. (16) Lu, G.; Blakesley, J.; Himmelberger, S.; Pingel, P.; Frisch, J.; Lieberwirth, I.; Salzmann, I.; Oehzelt, M.; Di Pietro, R.; Salleo, A. Moderate doping leads to high performance of semiconductor/insulator polymer blend transistors. Nature Commun. 2013, 4 (3), 1588. (17) Gao, S.; Bu, L.; Zheng, Z.; Wang, X.; Wang, W.; Zhou, L.; Hou, J.; Lu, G. Probing film-depth-related light harvesting in polymer solar cells via plasma etching.

AIP. Adv. 2017, 7 (4), 045312. (18) Finck, B. Y.; Schwartz, B. J. Drift-Diffusion Studies of Compositional Morphology in Bulk Heterojunctions: The Role of the Mixed Phase in Photovoltaic Performance. Phys. Rev. Applied 2016, 6 (5), 054008. ACS Paragon Plus Environment

Page 26 of 28

Page 27 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

(19) Cardona, M. Infrared Dielectric Constant and Ultraviolet Optical Properties of Solids with Diamond, Zinc Blende, Wurtzite, and Rocksalt Structure. J. Appl. Phys. 1965, 36 (7), 2181-2186. (20) Burkhard, G. F.; Hoke, E. T.; Mcgehee, M. D. Accounting for Interference, Scattering, and Electrode Absorption to Make Accurate Internal Quantum Efficiency Measurements in Organic and Other Thin Solar Cells. Adv. Mater. 2010, 22 (30), 3293-3297. (21) Persson, N. K.; Schubert, M.; Inganäs, O. Optical modelling of a layered photovoltaic device with a polyfluorene derivative/fullerene as the active layer. Sol.

Energy Mater. Sol. Cells. 2004, 83 (2–3), 169-186. (22) Moulé, A. J.; Meerholz, K. Minimizing optical losses in bulk heterojunction polymer solar cells. Appl. Phys. B. 2007, 86 (4), 721-727. (23) Park, H.; An, J.; Song, J.; Lee, M.; Ahn, H.; Jahnel, M.; Im, C. Thickness-dependent internal quantum efficiency of narrow band-gap polymer-based solar cells. Sol. Energy Mater. Sol. Cells. 2015, 143, 242-249. (24) Blakesley, J. C.; Neher, D. Relationship between energetic disorder and open-circuit voltage in bulk heterojunction organic solar cells. Phys. Rev. B. 2011, 84 (7),075210. (25) Zhang, Y.; Deng, D.; Wang, Z.; Wang, Y.; Zhang, J.; Fang, J.; Yang, Y.; Lu, G.; Ma, W.; Wei, Z. Enhancing the Photovoltaic Performance via Vertical Phase Distribution Optimization in Small Molecule:PC71BM Blends. Adv. Energy Mater. 2017, 7 (22), 1701548. (26) Yang, X.; Zheng, F.; Xu, W.; Bi, P.; Feng, L.; Liu, J.; Hao, X. Improving the Compatibility of Donor Polymers in Efficient Ternary Organic Solar Cells via Post-Additive Soaking Treatment. ACS Appl. Mater. Interfaces 2017, 9 (1), 618–627

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

TOC

ACS Paragon Plus Environment

Page 28 of 28