Reducing Energy Loss by Developing Luminescent Triphenylamine Functionalized Electron Acceptor For High-Performance Organic Solar Cells
Reducing Energy Loss by Developing Luminescent Triphenylamine Functionalized Electron Acceptor For High-Performance Organic Solar Cells
Supplementary Information
Zhang*c, Jiaxin Gaod, Xuelin Wangd, Xunchang Wange, Zheng Tangd, Renqiang Yange, Feng Gaoc, and
Yanming Sun*
Sweden
E-mail: [email protected]
dState Key Laboratory for Modification of Chemical Fibers and Polymer Materials, Center for
1
1.Experiment section
1.1 Materials
PM6 was purchased from Solarmer Materials (Beijing, China). Compound 1(2-(5-bromo-3-oxo-2,3-
[1,2,5]thiadiazolo[3,4-e]thieno[2'',3'':4',5']thieno[2',3':4,5]pyrrolo
from Energy Chemical. [2-(9H-Carbazol-9-yl)ethyl] phosphonic acid (2PACz) was purchased from
Aladdin. All purchased materials were used without further purification. The detailed synthetic
procedures and characterizations of the Z-Tri were described in the following sections.
Compound 3: Compound 1 (273 mg, 1.00 mmol) and compound 2 (1446 mg, 5.00 mmol) were
dissolved into tetrahydrofuran (45 ml) in a three-neck flask. K2CO3 (484 mg, 3.50 mmol) dissolved in
water (6 ml) was added to the three-neck flask. The solution was flushed with nitrogen for 10 min.
Then, add Pd (Pph3)4 (110 mg) to the solution, then let the solution temperature to 95°C for 5 h.
Washed with saturated salt water and ethyl acetate. The solvent was removed under reduced pressure.
The crude product was subsequently purified by column chromatography on silica gel to afford
compound 3 as an orange solid (417 mg, 93% yield). 1H NMR (500 MHz, CD2Cl2) δ 8.63 (d, J = 8.3
Hz, 1H), 8.19 – 8.00 (m, 2H), 7.60 (d, J = 8.8 Hz, 2H), 7.33 (dd, J = 8.5, 7.3 Hz, 4H), 7.22 – 7.07 (m,
Z-Tri: Compound 4 (L8-BO-CHO, 120 mg, 0.113 mmol) and compound 3 (396 mg, 0.904 mmol)
were dissolved into dry chloroform (40 mL) in a three-neck flask. The solution was flushed with
nitrogen for 15 min. After 0.6 mL pyridine was added, the mixture was stirred at 65°C overnight. After
cooling to room temperature, the reaction mixture was poured into water and extracted several times
with chloroform. Then the solvent was removed under reduced pressure, and the crude product was
2
purified by column chromatography on silica gel to yield Z-Tri as a black solid (187 mg, 87% yield)
.1H NMR (400 MHz, Chloroform-d) δ 9.15 (s, 2H), 8.73 (d, J = 8.4 Hz, 2H), 8.14 (d, J = 1.8 Hz, 2H),
7.97 (dd, J = 8.4, 1.9 Hz, 2H), 7.61 (d, J = 8.8 Hz, 4H), 7.34 (t, J = 8.5 Hz, 8H), 7.19 (m, 10H), 7.16 –
7.10 (m, 4H), 4.78 (t, J = 7.8 Hz, 2H), 3.19 (d, J = 7.5 Hz, 2H), 2.11 (m, 4H), 1.39 (m, 4H), 1.24 (m,
32H), 1.03 (m, 10H), 0.85 (dt, J = 13.8, 6.9 Hz, 12H), 0.76 (td, J = 7.5, 3.8 Hz, 6H), 0.64 (td, J = 7.2,
4.0 Hz, 6H); 13C NMR (126 MHz, CDCl3, δ): 188.49, 160.60, 152.31, 149.18, 147.60, 147.09, 146.45,
144.94, 137.99, 137.76, 137.74, 134.96, 134.09, 133.54, 132.23, 131.02, 129.70, 129.57, 127.92,
125.64, 125.33, 123.97, 122.48, 121.58, 120.26, 115.90, 115.39, 113.35, 67.22, 55.73, 40.35, 39.85,
34.52, 33.56, 33.30, 31.84, 29.73, 29.66, 28.87, 28.85, 27.67, 26.61, 26.59, 23.36, 23.33, 23.01, 22.88,
22.86, 22.67, 14.11, 14.07, 13.79, 13.77, 10.40; MS (MALDI-TOF) m/z: [M + H]+ calcd for
S
HO OH N N
CN
B C 4H 9 C 4H 9
NC CN NC S S
90℃ C6H13 C6H13
N
+ N Pd(Pph3)4,THF, O + OHC S
C 4H 9
N N
S
C 4H 9
CHO
O
C 2H 5 C 2H 5
Br Na2CO3
1 2 3 4
S
N N
C 4H 9 C 4H 9
S S
C6H13 C6H13
CN N N NC
S S
60℃ NC
C 4H 9 C 4H 9 CN
O C 2H 5 C 2H 5 O
CHCl3, pyridine
N N
Z-Tri
All the OSCs devices were fabricated by using a conventional structure of ITO/2PACz/active
water, acetone, and isopropyl alcohol for 20 min in an ultrasonic cleaning machine. The cleaned ITO-
3
coated substrates were dried in an oven at 100 °C overnight. Before use, the ITO-coated substrates
were pretreated by a plasma cleaner for 1 min under a vacuum condition below 10 Pa. Then a
monolayer of 2PACz (0.3 mg ml-1 in ethanol) was firstly deposited on the top of the ITO substrates at
a rate of 3000 rpm for 30 s, followed by thermal annealing at 100 oC for 10 min under ambient
conditions. The PM6:L8-BO (1:1.2 w/w) was dissolved in chloroform at the donor concentration of
7.3mg/mL with 165 wt% DBCl as the solid additive. The PM6:Z-Tri (1:1.2 w/w) was dissolved in
chloroform at the donor concentration of 7.3mg/mL with 165 wt% DBCl as the solid additive. The
PM6:L8-BO: Z-Tri (1:1:0.2 w/w/w) was dissolved in chloroform at the donor concentration of
7.3mg/mL with 165 wt% DBCl as the solid additive. Then the active layers with different amounts of
Z-Tri were spin-cast on the 2PACz layer surface in an N2-filled glove box at 2800 rpm for 40 s with a
total concentration of 16 mg mL−1 in chloroform (CF). The active layer solutions were stirred at 80 oC
for 4h before used. Then, the active layers were thermally annealed at 80 °C for 8 min. Subsequently,
a thin layer of PNDIT-F3N (1.2 mg ml-1 in methanol with 0.5% acetic acid, v/v) was spin-coated onto
the active layer at a rate of 4200 rpm for 30 s. Finally, a 120 nm-thick silver layer was thermally
deposited under a vacuum condition of 2×10-4 Pa. The active area of devices is 5.12 mm2.
The devices were tested through a mask with an area of 3.15 mm2. The device performance was
measured by using a 510 Air Mass 1.5 Global (AM1.5G) solar simulator (SS-F5-3A, Enlitech) with
an irradiation intensity of 100 mW cm-2, which was determined by using a calibrated silicon solar cell
(SRC2020, Enlitech). The J-V characteristics were measured by using a Keithley 2400 Source
Measure Unit. EQE spectra were performed by using a QE-R3011 Solar Cell EQE measurement
system (Enlitech).
All calculations were optimized using the ORCA 5.0.4 program package. The R2SCAN-3c composite
method proposed by S.Grimme was used to optimize the structure and ωB97X-D3 hybrid functional
combined with Def2-TZVP was used for single point calculations. The solvent effect of chloroform
4
was considered, using the continuum polarizable continuum model (CPCM). All post-calculation
analysis was carried out using the Multiwfn 3.8.4 Dev software package and the results were visualized
UV-vis absorption spectra of the pristine and blend films were acquired with a UV-vis
Cyclic voltammogram (CV) measurement was carried out on a CHI660E electrochemical workstation
with three electrodes, employing glassy carbon discs as the working electrode, Pt wire as the counter
electrode, and Ag/AgCl electrode as the reference electrode at a scanning rate of 100 mV s-1 under a
acetonitrile (HPLC grade) was adopted as the supporting electrolyte. PBDF-NS was drop-cast the
electrode from chloroform solutions to form thin films. Potentials were referenced to the
ferrocene/ferrocenium (Fc/Fc+) redox couple by using ferrocene as the external standard in acetonitrile
solutions (EFc/Fc+=0.42 eV). EHOMO =−[(4.8− EFc/Fc+) +Eox] (eV); ELUMO =−[(4.8− EFc/Fc+) +Ered] (eV).
The SCLC method was used to study the charge transport properties of neat film and blend films. The
charge carrier mobility was determined by fitting the corresponding J-V characteristics by using the
equation of J=9ε0εrµV2/8L3, where J is the current density, ε0 is the permittivity of free space, εr is the
permittivity of the active layer, µ is the charge carrier mobility (μh and μe), V is the effective voltage
(V=Vappl-Vbi-VR, in which Vappl is the applied voltage, Vbi is the built-in potential and VR is the voltage
loss from series resistance of the device) and L is the film thickness of the neat film or blend film.
5
sEQE measurements.
The sEQE measurements were conducted using a halogen lamp light source, a monochromator
(Newport CS260), a Stanford SR830 lock-in amplifier, and a Stanford SR570 current amplifier. The
EL measurements
The EL spectra were measured by using a source meter (Keithley 2400) to inject electric current, and
EQEEL measurements
EQEEL was determined by a homebuilt setup using a Keithley 2400 to inject current into the solar
cells. The Si detector (Hamamatsu s1337-1010BQ) and a Keithley 6482 picoammete were used to
The water and glycerol contact angle images of neat films were recorded by using a KRÜSS DSA 100
GIWAXS measurements were conducted at beamline 7.3.3 of the Advanced Light Source. Samples
were prepared on Si substrates utilizing blend solutions identical to those employed in the fabrication
of devices. A 10 keV X-ray beam was incident at a grazing angle ranging from 0.12° to 0.16°, chosen
to optimize the scattering intensity from the samples. The scattered X-rays were captured by a Dectris
AFM measurement was performed on a Dimension Icon AFM (Bruker) in a tapping mode under
ambient conditions. (GIWAXS) patterns were acquired by detector Pilatus3R 1M, Dectris (X-ray
6
PiFM is a multimodal Atomic Force Microscopy (AFM) technique that integrates the high resolution
of AFM with infrared spectroscopy to acquire topographical and molecular information with sub-10
nanometer resolution. All measurements were conducted using Molecular Vista's Visascope
microscope and Block Engineering's Laser Tune QCL, with an adjustable range from 775 to 1885 cm-1
and a spectral linewidth of 2 cm-1. PiFM images were created with a resolution of 256×256 pixels, line
speed of 0.89 lines/s, and a size of 1 µm2, using wavenumbers 1648 cm-1, 1278 cm-1, and 1539 cm-1
pump-probe system in collaboration with a regenerative amplified laser system from Coherent. An
800 nm pulse with a repetition rate of 1kHz, a length of 100 fs, and an energy of 7 mJ/pluse, was
generated by an Ti:sapphire amplifier (Astrella, Coherent). The time delay between pump and probe
was controlled by a motorized optical delay line with a maximum delay time of 8 ns. The samples
films were spin-coated onto the 1 mm-thick quartz plates and are encapsulated by epoxy resin in
nitrogen-filled glove box to resist water and oxygen in the air. The pump pulse is chopped by a
mechanical chopper with 500 Hz and then focused on to the mounted sample with probe beams. The
probe beam was collimated and focused into a fiber-coupled multichannel spectrometer with CCD
sensor. The energy of pump pulse was measured and calibrated by a power meter (PM400, Thorlabs)
The glass transition temperature (Tg) was determined using ultraviolet-visible spectroscopy.
Absorption spectra measurements were conducted on the receptor thin film, with the annealing
temperature ranging from 20°C to 240°C. Then, the deviation metric (DMT) of the each absorption
spectra was calculated, following the method from Samuel E. Root et al. DMT= ∑λmax λmin [IRT (λ) -
IT (λ)]2where λ is the wavelength, λmax and λmin are the upper and lower bounds of the optical sweep,
respectively, IRT(λ) and IT(λ) are the normalized absorption intensities of the as-cast (room
temperature) and annealed films, respectively. Then, Tgs were determined to points where the two
7
8
2.Figures and Tables
CF CB
9
(b)
L8-BO
Z-Tri
Deviation Metric
Tg= 106 ℃
Tg= 139 ℃
10
Fig. S3. UV-vis absorption spectra of blend films (a) in the solutions and (b) in the films.
11
(a) (b)
——PM6 ——L8-BO
LUMO=-3.48 eV LUMO=-3.88 eV
HUMO=-5.44 eV HUMO=-5.69 eV
Current
Current
-0.50 1.31
-0.90 1.06
-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5
(c)
——Z-Tri
LUMO=-3.52 eV
HUMO=-5.45 eV
Current
-0.86 1.07
Potential (V)
Fig. S4. Cyclic voltammetry curves of (a) PM6, (b) L8-BO, and (c) Z-Tri.
12
(a) (b)
4
10 Hole-only device Electron-only device
104
Current Density (mA cm-2)
102
102
101
101
100
100
Fig. S5. (a) SCLC curves of hole-only devices based on binary and ternary devices. (b) SCLC curves of
13
v(a)
1.0
Normalized PCE
0.8
PM6:L8-BO
0.6 PM6:Z-Tri
PM6:L8-BO:Z-Tri
0.4
0.2
0.0
0 150 300 450 600 750 900
Time (h)
Figure S6. The photostability of the solar cells under continuous 100 mW/cm² illumination in air.
14
v(a)
1.0
0.8
Normalized PCE
0.6
PM6:L8-BO
PM6:Z-Tri
0.4 PM6:L8-BO:Z-Tri
0.2
0.0
0 100 200 300 400 500
Time (h)
Figure S7. The thermal stability of the solar cells under continuous heating at 80°C.
15
L8-BO Z-Tri
LUMO=-3.45 eV LUMO=-3.24 eV
HUMO=-5.31eV HUMO=-5.25 eV
Fig. S8. The HOMO/LUMO structures of L8-BO and Z-Tri are based on DFT (B3LYP/6-31G level)
calculation. (Alkyl side chains were replaced with methyl groups for simplifying calculation.)
16
Heat Maps Calculation
In the donor-acceptor (D-A) interface of OSCs, heat maps are crucial for demonstrating the charge
transfer process. They vividly display the electron transfer from the donor to the acceptor and the hole
accumulation in the donor. On the heat map, the electron concentration gradually increases on the acceptor
side, indicating successful electron migration from the donor material, which is a key step in the charge-
separation process of OSCs. Meanwhile, the hole concentration rises on the donor side as electrons leave,
leaving behind holes that accumulate in this region. Such visualization by heat maps directly reflects the
charge transfer direction. The movement of electrons and the corresponding hole accumulation define the
charge-transfer pathway. Additionally, the degree of change in electron and hole concentrations on the
heat map can be used to infer the charge transfer efficiency. A steeper increase in electron concentration
on the acceptor side and a more significant rise in hole concentration on the donor side may imply a more
efficient process, whereas a less pronounced change may suggest barriers or inefficiencies in the
interface's charge-transfer mechanism. Overall, heat maps are powerful tools for understanding the
17
(a)
PM6:L8-BO
LE
CT
(b)
PM6:Z-Tri
LE
CT
Fig. S9. Optimized geometry of the PM6:L8-BO and PM6:Z-Tri complexes, electron-hole distribution
diagram.
18
Luminescence Lifetime Calculation
1.5
𝜏=
𝑓·∆𝐸2
(1)
This formula plays a crucial role in the field of luminescence research, enabling us to quantitatively
determine the duration for which a luminescent material emits light. In this formula, the variable 𝜏
represents the luminescence lifetime, and its unit is seconds. The luminescence lifetime is a key parameter
that characterizes the decay behavior of the excited - state population in a luminescent system. It provides
valuable information about the stability and efficiency of the luminescence process. On the other hand,
ΔE is the energy gap. It is measured in the unit of cm-1, which is a commonly used unit in spectroscopy
to describe energy differences between energy levels. This energy gap is related to the electronic or
vibrational transitions within the luminescent material. Based on the mathematical form of the formula 1,
we can clearly observe an inverse - square relationship between the luminescence lifetime 𝜏 and the energy
2
gap ΔE. As the value of the energy gap ΔE becomes smaller, the denominator 𝑓·∆𝐸 of the formula
decreases. Since the numerator is a constant value of 1.5, a decrease in the denominator leads to an
increase in the value of the entire fraction. Consequently, the luminescence lifetime 𝜏 becomes longer.
19
(a)
PM6:L8-BO
ΔE=0.519 eV
E=2.399 eV PM6:L8-BO
(b)
PM6-Z-Tri
ΔE=0.362 eV
PM6:Z-Tri
E=2.283 eV
Fig. S10. Optimized geometry of the (a) PM6:L8-BO and (b) PM6:Z-Tri excitation process diagram, and
20
Fig. S11. Contact angles of PM6, L8-BO, and Z-Tri neat films based on water and glycerol drops.
21
Energy Loss Calculation
As shown in equation 2, Eloss can be segmented into three distinct components: the radiative
recombination loss from absorption above the bandgap (ΔE1), the additional radiative recombination loss
from the absorption below the bandgap (ΔE2), and the non-radiative recombination loss (ΔEnr). ΔE1
SQ SQ
defined as Eg - qV oc (qV oc is maximum Voc predicted by the Shockley–Queisser limit) strongly
depends on the Eg, which is considered as unavoidable loss in any types of solar cells. ΔE2 can be obtained
SQ rad rad
from the equation of qV oc - qV oc , where V oc is V when only considering radiative recombination.
oc
rad
By fitting sEQE and EL spectra, we determined the V oc of the device.
SQ SQ rad rad
Eloss = Eg - qVoc = ( Eg - qV oc )+(qV oc - qV oc )+(qV oc - qVoc)
= ∆𝐸1+∆𝐸2+ΔEnr
(2)
ΔEnr has been recognized as the key factor to limit the Voc of OSCs, which can be calculated from the
quation 3:
ΔEnr = - KT ln(EQEEL)
(3)
where K is the Boltzmann constant, T is the Kelvin temperature. The Eg is determined by identifying the
22
(a) 1.2 (b) 1.2
PM6:L8-BO PL PM6:L8-BO:Z-Tri PL
Absorption Absorption
1.0 1.0
Normalized PL/Abs
Normalized PL/Abs
0.8 0.8
0.6 0.6
0.2 0.2
0.0 0.0
1.0 1.2 1.4 1.6 1.8 2.0 1.0 1.2 1.4 1.6 1.8 2.0
(c) 1.2
PM6:Z-Tri PL
Absorption
1.0
Normalized PL/Abs
0.8
0.6
0.4 Eg =1.546 eV
0.2
0.0
Energy (eV)
Fig. S12. Normalized reduced PL and absorption spectra of (a) PM6:L8-BO (b)PM6:L8-BO: Z-Tri, and
23
(i)
102
101
100
10-1
10-2
EQE
10-3
10-4
10-8
1.0 1.5 2.0 2.5 3.0 3.5
Energy (eV)
24
(a) (b)
1010 1010
Blank Blank
109 L8-BO on 109 Z-Tri on
L8-BO off Z-Tri off
108 108
107
Count
107
Count
106 106
105 105
104 104
103 103
700 800 900 1000 1100 700 800 900 1000 1100
(c)
1010
Blank
109 L8-BO:Z-Tri on
L8-BO:Z-Tri off
108
107
Count
106
105
104
103
700 800 900 1000 1100
Wavelength (nm)
Fig. S14. Schematic for a PLQY measurement of (a) L8-BO, (b) Z-Tri, and (c) L8-BO:Z-Tri. (The red
line is the light on the sample, directly exciting the sample, blank substrate beside; The blue line is the
light hitting the blank substrate, scattering the light to excite the sample; The black lines are only two
25
(e)
1.0
Linear summation
0.8 PM6:L8-BO:Z-Tri PM6:Z-Tri
Normalized EL
+
PM6:L8-BO
0.6
0.4
0.2
Fig. S15. The EL spectra of PM6:L8-BO:Z-Tri and binary blends (16.7% PM6:L8-BO+83.3% PM6:Z-
Tri).
26
(a) (b)
1000 1000
ΔA (mOD) ΔA (mOD)
L8-BO:Z-Tri L8-BO
Pump @ 650 nm -2.0 -0.50 1.0
Pump @ 650 nm
-2.0 -0.50 1.0
100 100
10 10
1 1
0.1 0.1
600 700 800 900 600 700 800 900
Wavelength (nm) Wavelength (nm)
(c)
1000
ΔA (mOD)
Z-Tri
-2.0 -0.50 1.0
Pump @ 650 nm
100
Delay time (ps)
10
0.1
600 700 800 900
Wavelength (nm)
Fig. S16. 2D TAS images of (a) L8-BO:Z-Tri, (b) L8-BO, and (c) Z-Tri films with 650 nm excitation.
27
Fig. S17. TA spectra of the Z-Tri pumped at 900 nm.
28
Fig. S18. TA spectra of the L8-BO pumped at 900 nm.
29
Fig. S19. TA spectra of the L8-BO:Z-Tri blends pumped at 900 nm.
30
Fig. S20. TA spectra of the PM6:L8-BO:Z-Tri blends pumped at 900 nm.
31
Fig. S21. TA spectra of the PM6:L8-BO blends pumped at 900 nm.
32
Fig. S22. TA spectra of the PM6:Z-Tri blends pumped at 900 nm.
33
(c)
1.0 L8-BO
Z-Tri
Normalized Absorption (a.u.)
L8-BO:Z-Tri
0.8
0.6
0.4
0.2
0.0
400 500 600 700 800 900 1000
Wavelength(nm)
34
(a) (b)
PM6:L8-BO PM6:Z-Tri
(c)
PM6:L8-BO:Z-Tri
Fig. S24. AFM height images of (a) PM6:L8-BO, (b) PM6:Z-Tri, and (c) PM6:L8-BO:Z-Tri blend films.
35
(a) (b) (c)
PM6 (d)
L8-BO (e)
Z-Tri
OOP
IP
(d) (e)
10-2 PM6 PM6
IP OOP
L8-BO L8-BO
Z-Tri Z-Tri
10-3
Intensity
10-3
Intensity
10-4
10-4
Fig. S25. The GIWAXS patterns of neat films of donor and acceptors. (a-c) The 2D GIWAXS patterns
of the films; (d) In-plane and (e) out-of-plane cuts of the films with GIWAXS patterns.
36
Fig. S26. 1H-NMR spectrum of compound 3 in CD2Cl2.
37
1892.978
Z-Tri
38
Z-Tri
39
Z-Tri
40
Table S1. Optical and electrochemical properties of the photoactive materials.
a) The maximum absorption in solutions; b) The maximum absorption in thin films; c) Absorption onsets
in thin films; d) Optical band gap (Eg) was calculated through the onset of the thin-film absorption edge;
e) Obtained by UV/Vis DMT results; f) Calculated by EHOMO = -4.80 eV - Eox; g) ELUMO = -Ered - 4.80
eV.
41
Table S2. Summary of photovoltaic parameters of the ternary devices with different Z-Tri contents.
Active Layer Voc (V) Jsc (mA cm-2) FF (%) PCE (%)a
42
Table S3. The hole mobilities and electron mobilities of the binary and ternary devices.
μ ha μe a
Active Layer μh/μe
(10−4 cm2 V−1 s −1) (10−4 cm2 V−1 s −1)
43
Table S4. Contact angle and surface tension data of PM6, L8-BO, and Z-Tri neat films.
44
Table S5. Crystal coherence length and the d-spacings of (010) peaks in OOP direction and (100) peaks
a) Obtained by the equation of d = 2π/q, in which q corresponds to the diffraction peak positions; b)
Calculated using the equation: CCL = 2 πK/FWHM, in which w is the full width at half maxima and K is
45
Table S6. The photovoltaic performance for PM6:L8-BO:Z-Tri (1:1:0.2) based OSCs with different
Thickness (nm) Voc (V) Jsc (mA cm-2) FF (%) PCE (%)
46
Table S7. Summary of photovoltaic parameters for PM6:L8-BO:Z-Tri (1:1:0.2) devices with different
annealing temperatures.
Temperature
Voc (V) Jsc (mA cm-2) FF (%) PCE (%)
(℃)
47
Table S8. Comparison of efficiency and energy loss between this work and earlier references.
48
PM6:Y6:ZY-4Cl 0.198 0.863 0.520 18.69 22
3.References
1. C. Li, J. Zhou, J. Song, J. Xu, H. Zhang, X. Zhang, J. Guo, L. Zhu, D. Wei, G. Han, J. Min, Y.
Zhang, Z. Xie, Y. Yi, H. Yan, F. Gao, F. Liu and Y. Sun, Nat. Energy, 2021, 6, 605-613.
2. H. Chen, Y. Zou, H. Liang, T. He, X. Xu, Y. Zhang, Z. Ma, J. Wang, M. Zhang, Q. Li, C. Li, G.
Long, X. Wan, Z. Yao and Y. Chen, Sci. China Chem., 2022, 65, 1362-1373.
3. Y. Zhang, J. Deng, S. You, X. Huang, J. Liu, Y. Cheng, B. Huang, X. Chen, Z. Sun, C. Yang, Q.
4. F. Meng, Y. Qin, Y. Zheng, Z. Zhao, Y. Sun, Y. Yang, K. Gao and D. Zhao, Angew. Chem. Int.
5. X. Si, Y. Huang, W. Shi, R. Wang, K. Ma, Y. Zhang, S. Wu, Z. Yao, C. Li, X. Wan and Y.
6. C. He, Q. Shen, B. Wu, Y. Gao, S. Li, J. Min, W. Ma, L. Zuo and H. Chen, Adv. Energy Mater.,
8. J. Wan, Y. Wu, R. Sun, J. Qiao, X. Hao and J. Min, nergy Environ. Sci., 2022, 15, 5192-5201.
9. K. Liu, Y. Jiang, F. Liu, G. Ran, F. Huang, W. Wang, W. Zhang, C. Zhang, J. Hou and X. Zhu,
10. R. Sun, Y. Wu, X. Yang, Y. Gao, Z. Chen, K. Li, J. Qiao, T. Wang, J. Guo, C. Liu, X. Hao, H.
11. C. He, Y. Pan, Y. Ouyang, Q. Shen, Y. Gao, K. Yan, J. Fang, Y. Chen, C.-Q. Ma, J. Min, C.
Zhang, L. Zuo and H. Chen, Energy Environ. Sci., 2022, 15, 2537-2544.
12. J. Gao, N. Yu, Z. Chen, Y. Wei, C. Li, T. Liu, X. Gu, J. Zhang, Z. Wei, Z. Tang, X. Hao, F.
13. J. Song, C. Zhang, C. Li, J. Qiao, J. Yu, J. Gao, X. Wang, X. Hao, Z. Tang, G. Lu, R. Yang, H.
Yan and Y. Sun, Angew. Chem. Int. Ed., 2024, 63, e202404297.
14. J. Fu, Q. Yang, P. Huang, S. Chung, K. Cho, Z. Kan, H. Liu, X. Lu, Y. Lang, H. Lai, F. He, P.
W. K. Fong, S. Lu, Y. Yang, Z. Xiao and G. Li, Nat. Commun., 2024, 15, 1830.
15. C. He, Z. Chen, T. Wang, Z. Shen, Y. Li, J. Zhou, J. Yu, H. Fang, Y. Li, S. Li, X. Lu, W. Ma, F.
Gao, Z. Xie, V. Coropceanu, H. Zhu, J.-L. Bredas, L. Zuo and H. Chen, Nat. Commun., 2022, 13,
2598.
16. G. Zhang, H. Ning, H. Chen, Q. Jiang, J. Jiang, P. Han, L. Dang, M. Xu, M. Shao, F. He and Q.
17. X. Ma, J. Wang, J. Gao, Z. Hu, C. Xu, X. Zhang and F. Zhang, Adv. Energy Mater., 2020, 10,
2001404.
18. Q. An, J. Wang, X. Ma, J. Gao, Z. Hu, B. Liu, H. Sun, X. Guo, X. Zhang and F. Zhang, Energy
19. Q. Ma, Z. Jia, L. Meng, J. Zhang, H. Zhang, W. Huang, J. Yuan, F. Gao, Y. Wan, Z. Zhang and
50
20. J. Song, C. Li, L. Zhu, J. Guo, J. Xu, X. Zhang, K. Weng, K. Zhang, J. Min, X. Hao, Y. Zhang,
21. Y. Cai, Y. Li, R. Wang, H. Wu, Z. Chen, J. Zhang, Z. Ma, X. Hao, Y. Zhao, C. Zhang, F. Huang
22. X. Duan, W. Song, J. Qiao, X. Li, Y. Cai, H. Wu, J. Zhang, X. Hao, Z. Tang, Z. Ge, F. Huang
23. C. Xiao, X. Wang, T. Zhong, R. Zhou, X. Zheng, Y. Liu, T. Hu, Y. Luo, F. Sun, B. Xiao, Z. Liu,
24. D. Sun, Z. Chen, J. Zhang, W. Song, J. Shi, J. Zhu, Y. Meng, F. Jin, S. Yang and Z. Ge, Sci.
25. H. Xia, Y. Zhang, K. Liu, W. Deng, M. Zhu, H. Tan, P. W. K. Fong, H. Liu, X. Xia, M. Zhang,
T. A. Dela Peña, R. Ma, M. Li, J. Wu, Y. Lang, J. Fu, W.-Y. Wong, X. Lu, W. Zhu and G. Li,
26. C. Li, G. Yao, X. Gu, J. Lv, Y. Hou, Q. Lin, N. Yu, M. S. Abbasi, X. Zhang, J. Zhang, Z. Tang,
Q. Peng, C. Zhang, Y. Cai and H. Huang, Nat. Commun., 2024, 15, 8872.
27. R. Sun, T. Wang, Q. Fan, M. Wu, X. Yang, X. Wu, Y. Yu, X. Xia, F. Cui, J. Wan, X. Lu, X.
28. C. Zhang, X. Zhong, X. Sun, Y. Wei, H. Wang, F. Wang, Y. Ji, X. Ding, J. Lv, S. Pang, H. Hu,
29. H. Xia, Y. Zhang, W. Deng, K. Liu, X. Xia, C. J. Su, U. S. Jeng, M. Zhang, J. Huang, J. Huang,
C. Yan, W. Y. Wong, X. Lu, W. Zhu and G. Li, Adv. Mater., 2022, 34, 2107659.
30. X. Yan, J. Wu, J. Lv, L. Zhang, R. Zhang, X. Guo and M. Zhang, J. Mater. Chem. A, 2022, 10,
15605-15613.
31. Y. Jiang, S. Sun, R. Xu, F. Liu, X. Miao, G. Ran, K. Liu, Y. Yi, W. Zhang and X. Zhu, Nat.
32. H. Lu, D. Li, W. Liu, G. Ran, H. Wu, N. Wei, Z. Tang, Y. Liu, W. Zhang and Z. Bo, Angew.
51
Chem. Int. Ed., 2024, 63, e202407007.
52