Open Access Article
This Open Access Article is licensed under a Creative Commons Attribution-Non Commercial 3.0 Unported Licence

Direct triplet sensitization of oligothiophene by quantum dots

Zihao Xu a, Tao Jin a, Yiming Huang b, Karimulla Mulla b, Francesco A. Evangelista *a, Eilaf Egap *b and Tianquan Lian *a
aDepartment of Chemistry, Emory University, 1515 Dickey Dr, Atlanta, GA 30322, USA. E-mail: francesco.evangelista@emory.edu; tlian@emory.edu
bDepartment of Materials Science and NanoEngineering, Department of Chemical and Biomolecular Engineering, Rice University, 6100 Main St, Houston, TX 77005, USA. E-mail: ee30@rice.edu

Received 3rd April 2019 , Accepted 12th May 2019

First published on 13th May 2019


Abstract

Effective sensitization of triplet states is essential to many applications, including triplet–triplet annihilation based photon upconversion schemes. This work demonstrates successful triplet sensitization of a CdSe quantum dot (QD)–bound oligothiophene carboxylic acid (T6). Transient absorption spectroscopy provides direct evidence of Dexter-type triplet energy transfer from the QD to the acceptor without populating the singlet excited state or charge transfer intermediates. Analysis of T6 concentration dependent triplet formation kinetics shows that the intrinsic triplet energy transfer rate in 1[thin space (1/6-em)]:[thin space (1/6-em)]1 QD–T6 complexes is 0.077 ns−1 and the apparent transfer rate and efficiency can be improved by increasing the acceptor binding strength. This work demonstrates a new class of triplet acceptor molecules for QD-based upconversion systems that are more stable and tunable than the extensively studied polyacenes.


Introduction

Triplet excitons in organic materials exhibit long lifetime, extended diffusion length and low-lying energy levels,1–3 and have found promising applications in photon upconversion.4 For example, an optical upconversion layer on a solar cell can capture sub-bandgap photons and emit above bandgap photons,5 increasing the efficiency of the conventional single-junction devices beyond the Shockley–Queisser limit.6 Various strategies for improving the overall efficiency of upconversion systems have been developed.7–10 In recent years, quantum dot (QD) sensitized upconversion systems11–14 have attracted intense interest as a versatile and promising approach because of their large absorption coefficient,15 small singlet-to-triplet energy gap and fast intersystem crossing (ISC) rate,16 and tunable band gap and band alignment.17

In a typical photon-upconversion system, triplet excitons are generated through a sensitizer that undergoes intersystem crossing from an excited singlet state to a triplet state. This process is followed by a sequential triplet energy transfer first to the transmitter then to the emitter. The latter, can undergo triplet–triplet-annihilation (TTA) and emit a higher energy photon.4,11,13,18–22 Thus, the efficiency of a typical sensitizer-emitter-based upconversion system can be represented by the following equation:4

 
ΦUC = ΦISCΦTETΦTTAΦFL(1)

The overall upconversion efficiency (ΦUC) is the product of efficiencies in each step involved, namely the ISC of the sensitizer (ΦISC), the triplet energy transfer (TET)23 from the sensitizer to transmitter to emitter (ΦTET), TTA of the emitter (ΦTTA), and the emitter's fluorescence (ΦFL). For a specific sensitizer and emitter, ΦISC, ΦTTA and ΦFL are determined by the material's properties. A promising area for performance improvement is the design of QD/transmitter complexes to enable efficient TET.24

In most reported QD sensitized upconversion systems, the transmitter/emitter molecules are solely limited to acenes and their derivatives, which have limited structural and energetic tunability and poor stability.11,13,18–22,25 Oligothiophenes have wide range of tunability in energetics and molecular structure,26 making them a very desirable class of triplet acceptors/transmitters in QD-organic hybrid TTA upconversion systems. Although an example of oligothiophene phosphonic acid grafted to cadmium selenide (CdSe) QDs was previously reported,27 the photoluminescence quenching was attributed to charge transfer and TET was not observed. We hypothesize that with appropriate design of the energetics of QDs and oligothiophene acceptors, efficient TET transfer can be facilitated while competing single energy and charge transfer pathways may be suppressed (Fig. 1).


image file: c9sc01648a-f1.tif
Fig. 1 Cartoon of QD–T6 complexes showing the energetics of the relevant states and associated excited state decay processes.

Results and discussion

In this work, we demonstrate for the first time a successful Dexter-type TET from CdSe QDs to carboxylic acid functionalized oligothiophene (3′′′,4′′-dihexyl-[2,2′:5′,2′′:5′′,2′′′:5′′′,2′′′′:5′′′′,2′′′′′-sexithiophene]-5-carboxylic acid or T6,Fig. 1). In this system, only the triplet energy transfer from the CdSe QD to T6 is energetically favored while the charge and singlet transfer are energetically uphill. We used transient absorption (TA) spectroscopy to first identify the lowest energy singlet (S1) and triplet (T1) state spectral features of free T6 molecules in solution. Then TA study of QD–T6 complexes provide direct evidence for triplet sensitization of T6 by CdSe QD. We also employed Density Function Theory (DFT) computations to characterize the low-lying excited states involved in the triplet–triplet transfer mechanism.

T6 (shown in Fig. 1) was synthesized following a literature procedure28,29 and characterized by 1H-NMR, 13C-NMR and high resolution mass spectrometry (see ESI and Fig. S1–S3 for details). The absorption spectra of four QD–T6 complex samples with increasing T6 concentrations ranging from 12.5 to 50 μM are shown in Fig. 2. These plots show a T6 band at 415 nm and the first excitonic peak of CdSe QDs at 584 nm. In contrast to the absorption spectrum of T6 in toluene (Fig. 2 and S4), the QD-bound T6 shows a red-shift of ∼20 nm, which is attributed to a change in the dielectric environment and has been observed in a similar system.30 The total amounts of T6 in the samples were obtained by subtracting QD only spectra from QD–T6 complex spectra. However, the total T6 amount contains both QD-adsorbed and free T6 molecules, which cannot be easily differentiated by their absorption spectra alone. The amount of QD-bound T6 is determined by fitting the transient kinetics of triplet formation in QD–T6 complexes (see below). QD photoluminescence (PL) intensity decreased with increasing T6 concentration (Fig. 2), suggesting possible energy or charge transfer from the QD to T6.


image file: c9sc01648a-f2.tif
Fig. 2 UV-Vis absorption (solid lines) and emission (dash-dot lines) spectra of QD–T6 samples with increasing T6 concentrations (0 to 50 μM) in toluene. The emission spectra were measured with 520 nm excitation of the QDs. Black lines are QD samples without T6. The adsorbed T6 UV-Vis absorption spectra are shown in dashed lines. The excitation pulse used in the transient spectroscopy study is also shown here. The 520 nm pulse has no overlap with T6 absorbance.

To understand the fluorescence quenching mechanism in QD–T6 complexes, we studied pristine T6 first to characterize its excited states spectral signature via TA spectroscopy. The TA spectra of T6 measured with 400 nm excitation (Fig. 3A and B) show clear evolution from a singlet to triplet excited state on the sub-nanosecond time scale. The TA spectra at <0.5 ns (Fig. 3A) show: (i) a ground (S0) state bleach (GSB) centered at 400 nm caused by the decrease of ground state T6 molecules; (ii) a stimulated emission from the T6 singlet excited state appearing as a negative signal at 500–600 nm; and (iii) a broad singlet (S1) excited state absorption (ESA) from 600 nm to 900 nm. The assignment of GSB and stimulated emission of T6 is based on comparison with the steady state absorption and emission spectra shown in Fig. S4. The (S1) ESA peak is formed at an early delay time (2–5 ps) and its amplitude grows from 2–5 ps (red) to 10–50 ps (yellow) as the stimulated emission band shifts to longer wavelength (by 20 nm) at the same time, which is attributed to a fast relaxation from the initial excited state to the v = 0 vibrational level of the S1 excited state. Similar fast relaxation processes have also been observed in other oligothiophene molecules.26,31–34 Both the S1 ESA and stimulated emission decay on the sub-nanosecond time scale to form the first triplet excited state (T1) with an absorption peak (T1 to Tn transition) centered at 715 nm (Fig. 3B). The lifetime of this species is longer than 10 μs and can be shortened by more than 1000 times in the presence of oxygen (Fig. S7), consistent with the assignment to a T1 state.35


image file: c9sc01648a-f3.tif
Fig. 3 TA spectra and kinetics of T6 in toluene. TA spectra of T6 at indicated delay times after 500 μJ cm−2 400 nm excitation: (A) 0–1.6 ns and (B) 0.1–50 μs. Data between 780–820 nm in (A) are cut out due to saturated probe light. 520 nm at 1 mJ cm−2 could not directly excite free T6 as shown in panel A (grey lines, averaged from 2 ps to 1 ns).

We conducted DFT computations to determine the triplet excited state energetics on a model of the T6 molecule by replacing the C6H13 alkyl side chains with methyl groups. All DFT computations were performed with the ORCA36 and Psi4[thin space (1/6-em)]37 programs using the B3LYP38 functional and the def2-SVP basis set.39,40 Equilibrium geometries for the S0 and T1 states are shown in Fig. S10 and S11. From these geometries, we computed vertical excitation energies to higher states by time-dependent density functional theory (TD-DFT) to identify low-lying electronic excited states involved in the triplet–triplet transfer mechanism. At the S0 geometry, T1 lies 1.85 eV above the ground state and is characterized by a spin-flip HOMO → LUMO excitation. Geometric relaxation of the T1 state stabilizes it by 0.47 eV, giving an adiabatic S0 to T1 transition energy of 1.38 eV. The second triplet state corresponds to a HOMO → LUMO + 1 excitation, with vertical transition energy of 2.10 eV. The third and fourth triplet states are predicted to lie at higher energies (2.64 and 2.82 eV) and therefore should not be operative in the observed triplet–triplet energy transfer.

To compute oscillator strengths for the T1 → Tn transitions, we performed a TD-DFT calculation at the optimized geometry of the T1 state using the triplet state as a reference. The vertical excitation spectrum of T1 is dominated by a transition at 1.88 eV (659 nm, fosc = 3.70) with a less intense peak at 2.18 eV (568 nm, fosc = 0.24). These results are in good agreement with the experimental observation of a transient peak at 715 nm in the transient absorption spectra for t > 100 ns.

The mechanism of exciton quenching in CdSe QD/T6 complexes were then studied by transient absorption spectroscopy, which probes directly not only the quenching kinetics, but also the identity of intermediates and products. The TA spectra of the QD–T6 complexes at early delay times (1–10 ps) show that the selective excitation of CdSe QDs at 520 nm (where T6 does not absorb) leads to a QD ground state bleach and excited state absorption41–44 at 420–620 nm (Fig. 4A), similar to those observed in QD only samples (Fig. S5) and indicative of excitons in the QD. The TA spectra at later delay times (Fig. 4A) show the decay of the QD exciton and the formation of a new absorption peak at 640–800 nm (Fig. S8d for all T6 concentration) and a S0 ground state bleach (Fig. 4B). These features are absent in QD only samples and can be assigned to T6 triplet formation by comparing to the T6 triplet TA spectrum (Fig. 4B). Control experiment of T6 only samples at 520 nm excitation (the grey circle in Fig. 4D, grey lines in Fig. 3A) shows negligible T6 triplet signal in the absence of QDs. Here the TA spectra data directly proved a successful triplet sensitization of T6 by QD.


image file: c9sc01648a-f4.tif
Fig. 4 Transient spectra and kinetics of QD/T6 complexes measured at 520 nm excitation. (A) TA spectra of QD–T6 excited at indicated delay times. Inset: expanded view of the triplet spectra at 540–800 nm. (B) Comparison of the averaged TA spectra (from 100 ns to 5 μs) of QD–T6 complexes (excited at 520 nm) and pristine T6 (excited at 400 nm). (C) Comparison of the kinetics of triplet formation of T6 (715 nm, blue circle), T6 GSB (415 nm, grey square, inverted and scaled) and QD exciton bleach (584 nm, blue line, inverted and scaled) in QD–T6 complexes. (D) Kinetics of T6 triplet formation (TF) with increasing T6 concentration (red, yellow, green, blue) in QD–T6 complexes. Solid lines: fits according to a model described in the main text. Grey dots: T6 signal of a control sample with 50 μM T6 and no QD measured under same condition at 520 nm excitation, showing no triplet signal. Excitation pulse energy density is 35 μJ cm−2 for panel A, C and D, 500 μJ cm−2 for panel B.

To determine the mechanism of triplet formation, we compare the kinetics of T6 triplet state growth and exciton decay in QDs. The kinetics of triplet formation (Fig. 4C) can be monitored by both the triplet ESA (TF, blue circle) and S0 ground state bleach (GSB, black square). Here, triplet ESA kinetics was obtained by averaging TA signals from 650 nm to 780 nm to increase the signal/noise ratio. The T6 GSB is obtained by subtracting the overlapping QD contribution at 415 nm from the total signal. As reported previously, the CdSe XB signal is dominated by state filling of electrons in the conduction band edge.41–45 As shown in Fig. 4C, the triplet formation kinetics agree well with the exciton bleach (XB, blue line) up to ∼20 ns, consistent with a direct Dexter type TET from the QD to T6. These kinetics deviate after ∼20 ns, which suggests the presence of other exciton decay pathways. As discussed below, we attribute this to the presence of QDs without adsorbed T6, in which excitons decay by radiative and nonradiative pathways within the QD. It has been previously reported that QD sensitized triplet formation can also occur via sequential charge transfer pathways: electron (hole) transfer followed by hole (electron) transfer.46–48 These pathways would produce T6+ or T6 intermediates, which are not observed in our study. This observation is also supported by the energetics alignment that prohibits charge transfer or singlet energy transfer (Fig. S12). Therefore, we attribute this sensitization process to be a direct Dexter-type triplet energy transfer from a triplet exciton state in the QD.

To examine whether TET occurs via static or dynamic quenching mechanisms, we measure TET kinetics as a function of total T6 concentration in the sample. As shown in Fig. 4D, the final amplitude and apparent formation rate of the T6 triplet state increase with the total T6 concentration, until it reaches saturation. Similar acceptor concentration dependent kinetics has been observed in previous studies of QD/electron acceptor complexes.44,49–52 This kinetics can be well described by a model that assumes that the number of adsorbed acceptor molecules in QD/acceptor complexes follows a Poisson distribution and the rate of transfer is proportional to the number of adsorbed acceptors. Assuming the same model for TET in QD/acceptor complexes, we can derive an expression of the kinetics of triplet formation (eqn (2)) and QD exciton bleach recovery (see ESI).

 
image file: c9sc01648a-t1.tif(2)

In eqn (2), P(m,i) is the Poisson distribution probability of finding QD/T6 complexes with i adsorbed T6 molecules in a sample with m average number of acceptors. kTET is the intrinsic energy transfer rate in a 1[thin space (1/6-em)]:[thin space (1/6-em)]1 QD/T6 complex. The decay of QD exciton bleach through radiative and nonradiative pathways within the QD is described by multiple exponential decay with amplitude aj and rate constant k1j. [QD(0)]* is the initial concentration of excited QDs. Further details of the model and derivation can be found in the ESI. Because the decay rate of T6 triplet excited state (k3) is three orders of magnitude slower than triplet formation (Fig. S7), it is not considered in the fitting of the triplet formation kinetics. The T6 triplet formation (Fig. 4D) and QD exciton bleach recovery (Fig. S8a) can be well fit by this model with an intrinsic TET of 0.077 ± 0.002 ns−1. The fit also yields the average number of adsorbed T6 molecules as a function of total T6 concentration, ranging from 0.16 to 0.54. This dependence can be fit to a Langmuir isotherm, as shown in Fig. S8b, to obtain a maximum number of adsorption sites and binding constant of T6 on CdSe QDs of θ = 1.4 ± 0.8, and K = 0.014 ± 0.012 μM−1, respectively. It is important to note that in the range of m (from 0.16 to 0.54) studied in this work, over 90% of the QD population in the samples has zero or one adsorbed T6 (Table S2) and the triplet formation kinetics is mostly contributed by QD with one bounded T6. This is confirmed by the comparison of normalized triplet formation kinetics for different samples (Fig. S8c), which are identical within the S/N ratio of this measurement. The increased T6 signal amplitude in Fig. 4D reflects a larger percentage of T6 bounded QDs in these samples.

From the relative amplitude of the triplet state absorption and QD exciton bleach and their extinction coefficients, we can estimate a TET efficiency of 15.4 ± 0.6% in the sample with 50 μM T6 concentration (see ESI for details). The low efficiency can be attributed to both the lower number of adsorbed T6 molecules and slow intrinsic TET rate. The average number of adsorbed T6 per QD is only around 0.5 at the highest coverage level, presumably limited by the competitive adsorption between T6 and the stabilization oleic acid ligand on the QD surface. In comparison, the Castellano group reported CdSe QD to 9-anthracene carboxylic acid (ACA) Dexter TET with an overall rate of 2 × 109 s−1 with an average number of ACA per QD of 12.12 This yields a similar intrinsic TET rate constant of 0.17 ns−1. To verify this point, we vigorously sonicated the sample with excess amount of T6 (375 μM) and measured the TA spectrum within an hour. This temporarily created a high loading which leads to more efficient TET at 31.8 ± 1.2% as shown in Fig. S9. Thus, a promising approach to increase the overall TET rate and efficiency in the T6/QD complex is to improve the binding strength of T6, through perhaps a different anchoring group.

Conclusions

We have demonstrated the first example of triplet energy sensitization of an oligothiophene (T6) by CdSe QDs. By TA spectroscopy we confirmed this process occur via direct Dexter-type triplet energy transfer. The efficiency of TET increases with the concentration of T6 in solution, reaching a value of 15.4 ± 0.6% in the sample with 50 μM T6 and 31.8 ± 1.2% in the sample with 375 μM of T6, although the latter exhibits poor long-term stability. The T6 concentration-dependent TET kinetics can be understood by a model that assumes a Poisson distribution of adsorbed T6 molecule on the QD. This analysis suggests that TET rates are similar in the 1[thin space (1/6-em)]:[thin space (1/6-em)]1 QD–T6 and previously reported QD-anthracene complexes, but T6 acceptors have relatively weaker binding constant on the surface, which leads to a smaller average number of acceptors and slower apparent TET rate. Our finding suggests that TET efficiency QD–T6 complexes can likely be further improved by increasing the acceptor binding strength (through modification of the acceptor and QD surface ligands) and fine-tuning energetic levels. Finally, this work identifies a new class of stable and highly tunable/functionalizable organic TET acceptors.

Conflicts of interest

There are no conflicts to declare.

Acknowledgements

T. L. gratefully acknowledges the financial support from the National Science Foundation (CHE-1709182). E. A. acknowledges support from the National Science Foundation (CHE-1821863). F. A. E. acknowledges support from the U.S. Department of Energy under Award No. DE-SC0016004, the Camille and Henry Dreyfus Foundation under award TC-18-045, and the Alfred P. Sloan Foundation under award FG_20166748.

Notes and references

  1. C. J. Bardeen, Annu. Rev. Phys. Chem., 2014, 65, 127–148 CrossRef CAS PubMed.
  2. A. Köhler and H. Bässler, J. Mater. Chem., 2011, 21, 4003–4011 RSC.
  3. A. Köhler and H. Bässler, Mater. Sci. Eng., R, 2009, 66, 71–109 CrossRef.
  4. T. N. Singh-Rachford and F. N. Castellano, Coord. Chem. Rev., 2010, 254, 2560–2573 CrossRef CAS.
  5. J. de Wild, A. Meijerink, J. K. Rath, W. G. J. H. M. van Sark and R. E. I. Schropp, Energy Environ. Sci., 2011, 4, 4835–4848 RSC.
  6. W. Shockley and H. J. Queisser, J. Appl. Phys., 1961, 32, 510–519 CrossRef CAS.
  7. J. Peng, X. Guo, X. Jiang, D. Zhao and Y. Ma, Chem. Sci., 2016, 7, 1233–1237 RSC.
  8. T. C. Wu, D. N. Congreve and M. A. Baldo, Appl. Phys. Lett., 2015, 107, 031103 CrossRef.
  9. N. Yanai, M. Kozue, S. Amemori, R. Kabe, C. Adachi and N. Kimizuka, J. Mater. Chem. C, 2016, 4, 6447–6451 RSC.
  10. N. Yanai and N. Kimizuka, Acc. Chem. Res., 2017, 50, 2487–2495 CrossRef CAS PubMed.
  11. Z. Huang, X. Li, M. Mahboub, K. M. Hanson, V. M. Nichols, H. Le, M. L. Tang and C. J. Bardeen, Nano Lett., 2015, 15, 5552–5557 CrossRef CAS PubMed.
  12. C. Mongin, S. Garakyaraghi, N. Razgoniaeva, M. Zamkov and F. N. Castellano, Science, 2016, 351, 369–372 CrossRef CAS PubMed.
  13. M. Wu, D. N. Congreve, M. W. B. Wilson, J. Jean, N. Geva, M. Welborn, T. Van Voorhis, V. Bulović, M. G. Bawendi and M. A. Baldo, Nat. Photonics, 2016, 10, 31–34 CrossRef CAS.
  14. K. Okumura, K. Mase, N. Yanai and N. Kimizuka, Chem.–Eur. J., 2016, 22, 7721–7726 CrossRef CAS PubMed.
  15. W. W. Yu, L. Qu, W. Guo and X. Peng, Chem. Mater., 2003, 15, 2854–2860 CrossRef CAS.
  16. G. D. Scholes and G. Rumbles, Nat. Mater., 2006, 5, 683–696 CrossRef CAS PubMed.
  17. J. Jasieniak, M. Califano and S. E. Watkins, ACS Nano, 2011, 5, 5888–5902 CrossRef CAS PubMed.
  18. Z. Huang, X. Li, B. D. Yip, J. M. Rubalcava, C. J. Bardeen and M. L. Tang, Chem. Mater., 2015, 27, 7503–7507 CrossRef CAS.
  19. M. Mahboub, Z. Huang and M. L. Tang, Nano Lett., 2016, 16, 7169–7175 CrossRef CAS PubMed.
  20. V. Gray, P. Xia, Z. Huang, E. Moses, A. Fast, D. A. Fishman, V. I. Vullev, M. Abrahamsson, K. Moth-Poulsen and M. Lee Tang, Chem. Sci., 2017, 8, 5488–5496 RSC.
  21. M. Mahboub, H. Maghsoudiganjeh, A. M. Pham, Z. Huang and M. L. Tang, Adv. Funct. Mater., 2016, 26, 6091–6097 CrossRef CAS.
  22. Z. Huang, D. E. Simpson, M. Mahboub, X. Li and M. L. Tang, Chem. Sci., 2016, 7, 4101–4104 RSC.
  23. D. L. J. Dexter, A Theory of Sensitized Luminescence in Solid, 1953 Search PubMed.
  24. Z. Huang and M. L. Tang, J. Am. Chem. Soc., 2017, 139, 9412–9418 CrossRef CAS PubMed.
  25. A. R. Reddy and M. Bendikov, Chem. Commun., 2006, 1179–1181,  10.1039/b513597d.
  26. R. S. Becker, J. S. deMelo, A. L. Macanita and F. Elisei, J. Phys. Chem., 1996, 100, 18683–18695 CrossRef CAS.
  27. R. Stalder, D. Xie, R. Zhou, J. Xue, J. R. Reynolds and K. S. Schanze, Chem. Mater., 2012, 24, 3143–3152 CrossRef CAS.
  28. S. Dufresne, G. S. Hanan and W. G. Skene, J. Phys. Chem. B, 2007, 111, 11407–11418 CrossRef CAS PubMed.
  29. Y. Huang, Z. Xu, S. Jin, C. Li, K. Warncke, F. A. Evangelista, T. Lian and E. Egap, Chem. Mater., 2018, 30, 7840–7851 CrossRef CAS.
  30. Z. Huang, Z. Xu, M. Mahboub, X. Li, J. W. Taylor, W. H. Harman, T. Lian and M. L. Tang, Angew. Chem., 2017, 129, 16810–16814 CrossRef.
  31. R. S. Becker, J. S. Demelo, A. L. Macanita and F. Elisei, Pure Appl. Chem., 1995, 67, 9–16 CAS.
  32. J. P. Yang, W. Paa and S. Rentsch, Chem. Phys. Lett., 2000, 320, 665–672 CrossRef CAS.
  33. J. P. Yang, W. Paa and S. Rentsch, Synth. Met., 1999, 101, 624–625 CrossRef CAS.
  34. D. V. Lap, D. Grebner and S. Rentsch, J. Phys. Chem. A, 1997, 101, 107–112 CrossRef CAS.
  35. C. Grewer and H.-D. Brauer, J. Phys. Chem., 1994, 98, 4230–4235 CrossRef CAS.
  36. F. Neese, Wiley Interdiscip. Rev.: Comput. Mol. Sci., 2012, 2, 73–78 CAS.
  37. R. M. Parrish, L. A. Burns, D. G. A. Smith, A. C. Simmonett, A. E. DePrince, E. G. Hohenstein, U. Bozkaya, A. Y. Sokolov, R. Di Remigio, R. M. Richard, J. F. Gonthier, A. M. James, H. R. McAlexander, A. Kumar, M. Saitow, X. Wang, B. P. Pritchard, V. Prakash, H. F. Schaefer, K. Patkowski, R. A. King, E. F. Valeev, F. A. Evangelista, J. M. Turney, T. D. Crawford and D. Sherrill, J. Chem. Theory Comput., 2017, 13, 3185–3197 CrossRef CAS PubMed.
  38. A. D. Becke, J. Chem. Phys., 1993, 98, 5648–5652 CrossRef CAS.
  39. F. Weigend and R. Ahlrichs, Phys. Chem. Chem. Phys., 2005, 7, 3297–3305 RSC.
  40. F. Weigend, J. Comput. Chem., 2008, 29, 167–175 CrossRef CAS PubMed.
  41. J. Q. Chen, K. F. Wu, B. Rudshteyn, Y. Y. Jia, W. D. Ding, Z. X. Xie, V. S. Batista and T. Q. Lian, J. Am. Chem. Soc., 2016, 138, 884–892 CrossRef CAS PubMed.
  42. H. Zhu, N. Song and T. Lian, J. Am. Chem. Soc., 2010, 132, 15038–15045 CrossRef CAS PubMed.
  43. J. Huang, Z. Huang, Y. Yang, H. Zhu and T. Lian, J. Am. Chem. Soc., 2010, 132, 4858–4864 CrossRef CAS PubMed.
  44. J. Huang, D. Stockwell, Z. Huang, D. L. Mohler and T. Lian, J. Am. Chem. Soc., 2008, 130, 5632–5633 CrossRef CAS PubMed.
  45. V. I. Klimov, Annu. Rev. Phys. Chem., 2007, 58, 635–673 CrossRef CAS PubMed.
  46. S. Garakyaraghi, C. Mongin, D. B. Granger, J. E. Anthony and F. N. Castellano, J. Phys. Chem. Lett., 2017, 8, 1458–1463 CrossRef CAS PubMed.
  47. S. S. Skourtis, C. Liu, P. Antoniou, A. M. Virshup and D. N. Beratan, Proc. Natl. Acad. Sci. U. S. A., 2016, 113, 8115–8120 CrossRef CAS PubMed.
  48. D. J. Weinberg, S. M. Dyar, Z. Khademi, M. Malicki, S. R. Marder, M. R. Wasielewski and E. A. Weiss, J. Am. Chem. Soc., 2014, 136, 14513–14518 CrossRef CAS PubMed.
  49. A. Boulesbaa, A. Issac, D. Stockwell, Z. Huang, J. Huang, J. Guo and T. Lian, J. Am. Chem. Soc., 2007, 129, 15132–15133 CrossRef CAS PubMed.
  50. M. Tagliazucchi, D. B. Tice, C. M. Sweeney, A. J. Morris-Cohen and E. A. Weiss, ACS Nano, 2011, 5, 9907–9917 CrossRef CAS PubMed.
  51. A. Issac, S. Jin and T. Lian, J. Am. Chem. Soc., 2008, 130, 11280–11281 CrossRef CAS PubMed.
  52. S. Jin, J. C. Hsiang, H. Zhu, N. Song, R. M. Dickson and T. Lian, Chem. Sci., 2010, 1, 519–526 RSC.

Footnote

Electronic supplementary information (ESI) available. See DOI: 10.1039/c9sc01648a

This journal is © The Royal Society of Chemistry 2019