Abstract

Femtosecond transient absorption spectroscopy of thin films of two types of morphologies of diketopyrrolopyrrole low band gap polymer/fullerene-adduct blends is presented and indicates triplet state formation by charge recombination, an important loss channel in organic photovoltaic materials. At low laser fluence (approaching solar intensity) charge formation characterized by a 1350 nm band (in ~250 fs) dominates in the two PDPP-PCBM blends with different nanoscale morphologies and these charges recombine to form a local polymer-based triplet state on the sub-ns timescale (in ~300 and ~900 ps) indicated by an 1100 nm absorption band. The rate of triplet state formation is influenced by the morphology. The slower rate of charge recombination to the triplet state (in ~900 ps) belongs to a morphology that results in a higher power conversion efficiency in the corresponding device. Nanoscale morphology not only influences interfacial area and conduction of holes and electrons but also influences the mechanism of intersystem crossing (ISC). We present a model that correlates morphology to the exchange integral and fast and slow mechanisms for ISC (SOCT-ISC and H-HFI-ISC). For the pristine polymer, a flat and unstructured singlet-singlet absorption spectrum (between 900 and 1400 nm) and a very minor triplet state formation (5%) are observed at low laser fluence.

1. Introduction

Since the first observation of photo-induced electron transfer from a conducting polymer to buckminsterfullerene (C60) in 1992 [1], crucial steps towards better materials for plastic solar cells have been made. The introduction of a more soluble fullerene adduct (PCBM) [2], the development of the bulk heterojunction concept [3], and the switch to thiophene-based polymers (P3HT) [4] were paramount. The next steps that have to be taken in order to improve these solar cells are the absorption of a larger part of the available solar radiation and minimizing charge loss [5]. Efficiently harvesting a larger part of the solar spectrum and controlling charge recombination while maintaining the right nanoscale morphology for efficient charge generation and optimizing molecular organization for charge transport are the current challenges in organic photovoltaic (OPV) research. Low band gap materials are therefore anticipated to be a key component in the next generations of organic solar cells [69]. With this type of material combined with [70]PCBM, up to 11.7% power conversion efficiency (PCE) has been obtained [10, 11].

The diketopyrrolopyrrole (DPP) unit has gained a lot of interest as a component for small molecule-based OPV and has also been incorporated into polymers for solar cells [12, 13]. DPP-based polymer materials (see Figure 1) can absorb sunlight up to ~900 nm and it has been shown that these materials can yield OPV devices with a ~5 to ~6% power conversion efficiency [14]. The use of a cosolvent like 1,8-diiodooctane (DIO) or ortho-dichlorobenzene (ODCB) during spin coating is crucial for obtaining the optimized nanostructure and correlated higher efficiency, mainly attributed to a better ordering in the fullerene phase. Optimized [60]PCBM devices were reported with a short circuit current JSC = 9.3 mA cm−2, an open circuit voltage VOC = 0.78 V, a fill factor FF = 0.63, and a PCE = 4.6%. The HOMO and LUMO levels of the polymer were estimated to be −5.35 and −3.53 eV versus vacuum. The absorption onset of the films was reported to be 1.53 eV. The external quantum efficiencies ranged from EQE = 45% (600–750 nm region) to 20% at 530 nm. Furthermore, steady-state photo-induced absorption spectroscopy showed the presence of positive polarons at 1350 nm and triplet states at 1100 nm (see Figure 2), the latter indicating triplet state formation via charge transfer states [14].

Thus, while a larger part of the solar spectrum can be harvested with these DPP-based devices (as compared to, e.g., those based on PCBM/P3HT that have their main P3HT absorption only until 675 nm) [15], the output is not proportionally increased. Clearly, more loss channels are present that reduce the efficiency, and now, triplet state formation plays a role in the charge loss in these materials.

Charge loss can occur through trap states and by Langevin-type recombination at the interface of percolation pathways. It has recently been shown that, through spin statistics, triplet states can be formed through this recombination process on the sub-microsecond to millisecond timescale [16]. Sub-ns timescale charge loss through triplet recombination channels in polymer solar cell materials has also gained a lot of interest in the last years [17]. Furthermore, it is important to note that in a recent report, Howard et al. stated that currently there is no insight into the exact mechanism of charge recombination to the triplet state [18].

Here, we report on the primary photophysical events occurring in two types of thin-film blends containing a low band gap polymer (PDPP) and PCBM (see also Figure 1) as well as on the events occurring in the pristine polymer thin film as observed with femtosecond transient absorption spectroscopy (fs-TA). The official name of the PDPP polymer is poly[{2,5-bis(2-hexyldecyl)-2,3,5,6-tetrahydro-3,6-dioxopyrrolo[3,4-c]pyrrole-1,4-diyl}-alt-{[2,2′-(1,4-phenylene)bisthiophene]-5,5′-diyl}]).

At low laser fluence, where singlet-singlet annihilation or processes like “hot singlet fission” [19] do not play a role, a normal excited singlet state is observed for the polymer film.

For the blends of PCBM and the low band gap polymer at low laser fluences (approaching solar intensity) ultra-fast charge generation dominates, but clearly triplet state formation occurs on a sub-ns timescale. Morphology influences charge generation and triplet formation. Fast charge recombination to the triplet state is thus identified as an important charge loss in these blends. The potential gain in efficiency due to the low band gap is counter effected by the lower triplet state of the polymer, which is partly accessible from the charge transfer state. Our results indicate that the efficiency of the plastic solar cells made of these materials will be highly influenced by the solar-light intensity, and is expected to be relatively high at low solar fluence (comparing, e.g., 1 versus 10 sun). Extrapolation to the limit of zero laser fluence projects a (triplet) loss-corrected PCE of 9.6% for these materials.

We present a model for the influence of nanoscale morphology on charge recombination to the triplet state. We use previous studies of molecular electron donor acceptor systems in solution regarding triplet state formation by charge recombination from compact charge transfer states. We combine this with recent results that suggest the occurrence of long range photo-induced electron transfer occurring in optimized photovoltaic blends, implying a large distance between hole and electron.

2. Experimental

2.1. Film Preparation

The synthetic procedures for the preparation of PDPPTPT have been reported [14]. [60]PCBM () was obtained from Solenne B.V. The photoactive layer was deposited on quartz by spin coating a chloroform solution containing 6 mg mL−1 PDPPTPT and 12 mg mL−1 PCBM and 25 mg ortho-dichlorobenzene (or 25 mg 1,8-diiodooctane) per mL solution. Thickness of the active layers (80–90 nm) was measured on a Veeco Dektak 150 profilometer. Optical densities at 530 nm were between 0.098 and 0.11.

2.2. Femtosecond Transient Absorption Spectroscopy

Femtosecond transient absorption experiments were performed with a Spectra-Physics Hurricane Titanium:Sapphire regenerative amplifier system. The full spectrum setup was based on an optical parametric amplifier (Spectra-Physics OPA 800C) as the pump. The residual fundamental light, from the pump OPA, was used for NIR probe-light generation, using a sapphire plate. For femtosecond transient absorption in the NIR region, a Control Development NIR-256L-1.7T1-USB optical spectrometer system and an InGaAs detector with 512 element arrays responding to wavelengths ranging from 900 to 1500 nm were used. The polarization of the pump light was controlled by a Berek Polarization Compensator (New Focus). The Berek Polarizer was always included in the setup to provide the magic angle conditions. The probe light was double passed over a delay line (Physik Instrumente, M-531DD) that provides an experimental time window of 3.6 ns with a maximal resolution of 0.6 fs per step. The OPA was used to generate excitation pulses at 530 nm. The laser output was typically 28 μJ/cm2 per pulse (130 fs FWHM) with a repetition rate of 1 kHz. Laser fluence was adjusted with a neutral density filter wheel. The solution samples were placed into cells of 2 mm path length (Hellma) and were stirred with a downward-projected PTFE shaft using a direct-drive spectro-stir (SPECTRO-CELL).

2.3. Global and Target Analysis Using Glotaran

The analysis of the time-resolved spectroscopic data consisting of a three-dimensional dataset (wavelength, time and absorption difference intensity) was performed with the global and target analysis program Glotaran, http://glotaran.org/.

After a first singular value decomposition of the raw data matrix, the starting number of components was estimated. From the first analyses with a sequential global model, the chirp corrections and the main decay/rise time constants were obtained. The number of components was further optimized. On the basis of the spectral shape (and experience), a target model was set up that contains branching from specific states. Then the ratios of this branching were introduced which influence the amplitude of the species-associated difference spectra (SADS), but not their shape or the decay/rise time constants. The most important spectral assumptions in the target analysis are the equality based on the singlet and triplet features of the polymer and the radicalcation spectrum of the polymer (see also Figure 2).

All spectra (measured at 256 wavelengths) were collated in a matrix, which was globally fitted using a sequential kinetic scheme with increasing lifetimes. From this, the lifetimes and the evolution-associated difference spectra (EADS) were estimated. The quality of the fit was judged by inspection of the singular vectors of the matrix of residuals, which had to be of less structure. The instrument response function was described by a Gaussian shape, and the white-light dispersion over the spectral range was modeled as a second-order polynomial. With increasing lifetimes, and thus decreasing rates, the first of the EADS decays with the first lifetime and corresponds to the difference spectrum at time zero with an ideal, infinitely small instrument response function. The second of the EADS is formed with the first lifetime and decays with the second lifetime and so on. The error in the lifetimes obtained from the fitting procedure does not exceed 10%. EADS may not represent pure species, and they are interpreted as a weighted sum (with only positive contributions) of species-associated difference spectra (SADS).

To resolve the SADS from the EADS, a target analysis was performed on the data. In this target analysis, a kinetic scheme was used to estimate the microscopic rate constants and SADS of the different species.

The matrices of residuals resulting from the target analysis were further analyzed using a singular value decomposition. The first left and right singular vectors show no or little structure, indicating that all kinetics are described satisfactorily.

2.3.1. Quantum Chemical Calculations

We have calculated the triplet level of a DPP tetramer at 0.81 eV with TDDFT/B3LYP with the 6-31G(d) basis set. Full geometry optimization of the DPP4 molecules were performed at the DFT level with Becke’s three-parameter functional and the Lee-Yang-Parr functional (B3LYP) level associated with the 6-31G(d) basis set [20]. Time-dependent (TD) B3LYP calculations are performed for the vertical excitation energies of the lowest 5 singlet excited states, the first triplet state, and their related configuration schemes with the 6-31G(d) basis set. A developmental version of the quantum chemistry package Q-Chem was used for all calculations in the present work [21].

3. Results and Discussion

The steady-state UV-Vis absorption spectra of the pristine polymer and the blends have been reported previously [14] and show an absorption onset at 810 nm. Therefore, fs-TA was performed in the NIR region on a polymer solution (CHCl3), a pristine polymer film, a PCBM blend with a coarse morphology (code COARSE, PCE = 2%, with 200 nm domains), and a PCBM blend with a fine morphology (code FINE; the fine morphology results in a ~2 × higher power conversion efficiency, PCE = 4.6%, with 20 nm domains). A laser excitation fluence of 4 μJ/cm2 (at 530 nm) was used. This corresponds to 1.07 × 1013 photons/cm2 or 1.4 × 1017 photons/cm3 for a 85 nm thick film with an OD at 530 nm of 0.1. Clearly, the initial exciton density approaches charge carrier densities comparable to those in typical organic solar cells under normal operation (approximately 1016 to 1017 cm−3) [22]. In order to simulate normal operating conditions, low photon densities have to be used resulting in conditions where the signal intensity approaches a few milli-OD.

The data were fitted with a global and target analysis method [2325] using the Glotaran software package. Since our time-window of observation is 3.6 ns we approximate the kinetics with exponential functions (using a rate-(K)-matrix) and disregard the second-order kinetic components (see Section 2 for details). Although approaches involving distributed decay kinetics [26] or analytic solutions to rate equations [27] are available, this approximation with exponential functions is often applied for the early timescale [28, 29] and gives a good fit and results. We recently published a target model analysis using Glotaran of a perylene-bis(dicarboximide)/pyrene blend as well as of an annealed PCBM/P3HT film [46].

Figure 3 shows a selection of the spectra of all datasets obtained for the three materials and Tables 1 and 2 display the results of the data analysis. Triplet formation by non-geminate charge recombination of free-charge carriers that are created in ~0.2 ps implies the intermediacy of bound electron hole pairs (also called polaron pairs). In this work, we refer to these states as charge transfer state (CT states).

3.1. Pristine Polymer Films

For the pristine polymer in solution and (at low fluence) for the thin film, we observe a virtually structureless flat excited state absorption spectrum in the 900–1350 nm region, attributed to the singlet excited state (Figure 3(i) and Figure 4(a)). Two singlet excited state lifetimes are observed for the film (2.5 and 41 ps). These are attributed to two excitonic states with different effective conjugation lengths present in the solid film. In the solution, the longest lifetime (253 ps) is attributed to free solvated polymer molecules whereas aggregated species that are also present have a shorter (22 ps) lifetime. Previous studies have indicated that this polymer partly aggregates in solution [14]. The 1 ps component is attributed to solvation processes. A very minor amount of triplet state formation (~5%) is observed for the thin film. The triplet excited state spectrum of the DPPP polymer was previously reported, observed with steady-state PIA, [14] and shows a maximum at ~1100 nm (see also Figure 2). Singlet-singlet, triplet-triplet, CT-CT, and CT-exciton annihilation are known to play a role in fs-TA spectroscopy of thin-film blends [30, 31]. In the simplest approximation, annihilation processes are known to display a quadratic dependence on the laser fluence () [32]. Such a fit (Figure 4(e)) gives a relative efficiency of triplet formation () value of 3%:

The various lifetimes and maximum populations of triplet formation are reported in Table 1. The laser fluence effects on these quantities are visualized in Figure 4(e). The lifetimes of the singlet excited states become longer upon lowering the laser fluence (from 1 to 41 ps). The relative triplet yield shows a similar trend, less triplet is formed at lower laser fluence (from 0.43 to 0.05). The triplet excited state lifetime shows dependence on the fluence, going from 1.1 ns (at 28 μJ/cm2) to more than 3.6 ns (at 4 μJ/cm2) at low fluence. It has to be noted that our observation window is 3.6 ns, and therefore, this latter value is on the limit of our detection. Clearly, the triplet state has a much longer lifetime, as it was observed previously with steady-state PIA for a similar film [14]. The 1.1 ns lifetime at a higher laser fluence is influenced by triplet-triplet annihilation. The singlet excited state lifetime shows a sharper dependence and so does the relative triplet yield. The complex photophysical behavior with respect to the laser fluence of the pristine polymer films is reflected in the properties of the blends.

3.2. Thin-Film Blends

For the blend films at low laser fluence, we observe a short-lived singlet-excitonic state that converts rapidly to a charged state in ~0.2 ps, characterized by the positive polaron absorption. The electrochemically oxidized PDPP has been reported earlier and shows a typical absorption band at 1350 nm (and a shoulder at 900 nm, sloping down from the blue side of the spectral range, see also Figure 2), clearly shifted relative to the triplet state absorption at 1100 nm. Analysis of the data indicates that now, the main part of the charged species recombines to form the triplet excited state of the polymer. As can readily be seen in the spectra shown in Figure 3, the band at 1100 nm develops as the 1350 nm band decays.

The six data sets of the blends (the FINE and COARSE morphology at three light intensities) were analyzed with the global and target analysis method (Glotaran).

The main outcome of the target analysis for the PDPP/PCBM blends is presented in Figure 5 and Table 2. Next to the most important time constants, we obtain the so-called species-associated difference spectra (SADS). These are the spectra belonging to specific species that play a role in the light-induced processes. In order to describe the matrix of time resolved spectra appropriately, we need two species with a “singlet” character, two species with charges, and one triplet state. A characteristic target model is shown in Figure 5(a), in which, for example, the transition from A to E is the singlet fission. The relative contribution of this particular channel will be influenced by laser power and morphology. This is the case for all branching pathways within the target model in Figure 5(a). Figure 5(b) shows the populations versus time of the various excited state species of which the spectra are displayed in Figure 5(c) (species-associated difference spectra, SADS). Within the analysis, the shape of the SADS is a main focus. Separating triplet contributions (1100 nm band) from charge transfer state contributions (1350 nm band) is an important argument when moving from a simple sequential global analysis (e.g., ) to a target model (with branching).

Figure 5 and Table 2 show how the photophysics is influenced by morphology and fluence. The increase in efficiency of charge formation with fluence reduction is stronger for the FINE film, accompanied by a reduction in relative triplet formation. The lifetime of the charges increases more for the FINE film, upon fluence reduction, as compared to the COARSE film. The effects of the laser fluence effects and the morphology can be understood by a kinetic competition of hot singlet fission and CT formation.

If we assume a quadratic dependence on the laser fluence, a fit results in 72% CT and 38% triplet for the FINE morphology, as values extrapolated to zero fluence.

For the COARSE morphology, these values are 63% CT and 48% triplet. Since a negligible amount of triplet is formed in the pristine polymer at low fluence, the ratios of the above values quantify the anticipated charge loss at solar conditions (38/72) × 100 = 52% charge loss for the FINE morphology and 76% charge loss for the COARSE morphology. This would imply loss-corrected PCE values of 9.6 and 8.3%.

Clearly, the ηT and ηCT are strongly influenced by annihilation processes (resulting in triplet formation) that compete with free-charge generation.

A linear fit to the lifetimes gives extrapolated values of 968 and 323 ps for τ(CT) FINE and COARSE. Triplet lifetimes of 6.1 and 3.7 ns were obtained for the FINE and COARSE blends. The linear fit (see Figure 5(e)) indicates that these decays are not governed by annihilation but by, for example, a constant concentration of trap states present within the thin-film blend.

The femtosecond spectroscopy presented here can be correlated to the spectra of the triplet state obtained by ns spectroscopy and spectro-electrochemistry that have been reported before. Figure 2 reports these spectra (from reference [14]) that have been converted to a wavelength scale for clarity.

3.3. Mechanistic Aspects

Non-geminate bimolecular recombination coefficients (β) in polymer/fullerene solar cells have been determined to be 6 × 10−11 cm3 s−1 (for PDPP-5T/PCBM) [33], 3 × 10−12 cm3 s−1 and 2 × 10−13 cm3 s−1 (both for P3HT/PCBM) [34, 35]. This correlates to decay times between 60 ns and 1.6 μs at the lowest laser fluence used in this study (assuming 100% initial charge generation). This is slower than the recombination rates observed here.

Charge recombination to the triplet state has been studied in many molecular electron donor acceptor systems, [3641] and these studies can be correlated to our results. To estimate the energetics of our PDPP/PCBM system, we can use the data Eox = +0.25 V versus Fc/Fc+ in ODCB, Ered = −1.07 V versus Fc/Fc+ in ODCB, and ES1 = 1.53 eV. For the triplet state level, we can use values from a study of molecular model systems [42] giving a triplet level of DPP-T2 at 0.9 eV. We have calculated a triplet level of a DPP tetramer at 0.81 eV with TDDFT/B3LYP with the 6-31G(d) basis set. We thus take ET = 0.85 eV as the average estimated value. The dielectric constant present in the solid can be taken to be 4 [4345].

We can also estimate the energy of the CT density of states by using the open circuit voltage and the “ rule” to give resulting in a value between 1.18 and 1.28 eV. This is just above the anticipated triplet energy level and correlates well with the values obtained for molecule model systems (DPP-T2/C60 Ecss =1.52 to 1.32 eV) [42].

Within the Marcus theory of electron transfer, spin is not taken into account. Clearly, a recombination to the local triplet state must be accompanied by an electron spin flip. The compensating magnetic momentum change in intersystem crossing (ISC) can operate through spin orbital coupling (SOC) or through proton hyperfine interactions (H-HFI).

For molecular diketopyrrolopyrrole-oligothiophene(T)-fullerene triads [42] (DPP-T2-C60 in toluene), we estimated the intrinsic intersystem crossing time constant to be ~15. This was reported to be the rate-limiting step for (triplet) charge recombination in this molecular system studied in solution. For the shorter molecular DPP-T1-C60 (with a larger exchange integral, J, see later) 0.8 ns in ortho-dichlorobenzene was found [42]. Recently, we showed fast and slow rates for charge recombination [46] to the triplet state in perylene bisimide/pyrene thin-film blends ( and ).

For molecular model systems in solution, two mechanisms can be responsible for ISC in charge transfer states: spin-orbit charge-transfer ISC (SOCT-ISC) and radical-pair ISC (RP-ISC) [5257]. In SOCT-ISC, the change of electron magnetic momentum needed for the electron flip is compensated by the change of orbital magnetic momentum that occurs when an electron is moved from one orbital to another. If the symmetry of these orbitals is different or if their main orbital coefficients lie in different regions in space (if they are for instance orthogonal), this is beneficial for SOCT-ISC. Studies by Colvin et al. have given clear evidence for this mechanism [5257]. RP-ISC proceeds via proton hyperfine interactions (H-HFI). The nuclear spin momentum of the hydrogen atoms have to adjust in such a way that magnetic momentum compensation occurs, within the radial pair, converting a singlet CT state into a triplet CT state. This state can then recombine to the local triplet state without spin forbiddenness. So far, we cannot discriminate between these two modes of operation but using a deuterated polymer should be an interesting approach. We suggest a model in which optimized morphology is correlated to a CT state with a smaller exchange integral and a reduced contribution of SOCT-ISC (see Figure 7). This model is in line with recent results obtained by using the photo-induced Stark shift [47], indicating a larger distance of charge separation in optimized blends as studied by Gélinas et al. The correlation between electron and hole distance, the energy difference between singlet and triplet, the relative contribution of spin orbit coupling, and the contribution of the radical pair mechanism was discussed before [48].

The ultrafast (ps) triplet generation we observe here and the strong influence of the laser power on this process indicate that exciton annihilation is responsible [19, 49].

Spin conservation rules dictate that it is spin forbidden to generate one triplet (and a ground state singlet). Thus, hot singlet fission is the only viable option for such ultra-fast triplet state formation. Furthermore it is important to note that, in general, singlet fission occurs in very well ordered crystalline phases. Interestingly, there have also been reports on singlet fission in spin-coated regio-random poly-thiophene films [50].

4. Discussion and Conclusions

At low laser fluence (approaching solar intensity), charge separation dominates in PDPP/PCBM blends. Triplet formation by charge recombination (Figure 6) is an important decay channel of the charge transfer state and morphology influences this process.

We show complex photophysics and fast kinetics of the spin flip dynamics to a triplet excited state occurring from the CT state in a low band gap polymer/PCBM blend at low laser fluence. The charge recombination effect to the polymer triplet state occurs on a ~300 to ~900 ps timescale, depending on film morphology, and plays a role in both PDPP/PCBM blend films. Triplet formation by charge recombination is slower in the FINE blend film of PDPP/PCBM (with 20 nm domains) relative to the COARSE blend film (with 200 nm domains). At the low laser fluences used here, charge generation dominates and triplet formation by charge recombination occurs. If we are able to control the triplet formation by charge recombination, we can anticipate a power conversion efficiency of 9.6% for these types of materials (based on value corrected for the charge loss to the triplet state).

In polymer films (without PCBM), we see that the singlet excited state of the PDPP polymer films shows lifetimes of 2.5 and 41 ps and a virtually flat S1-Sn absorption band in the 900–1400 nm region (at low laser fluence).

It has been estimated that charge recombination to the triplet states with a mode of operation of proton hyperfine interactions occurs with a rate between 1 × 106 and 1.4 × 108 s−1 [51]. However, in molecular D-A systems studied in solution, various examples are available where ns or even sub-ns charge recombination to the triplet occurs. Spin orbit coupling coupled to charge transfer (SOCT) must be an additional important and faster factor for ISC, (e.g., 1 × 109 s−1 to 2.5 × 1010 s−1) [41], especially if the orbitals corresponding to the transition are orthogonal [5257]. Okada et al. observed indications for SOCT-ISC already in 1981 [41], and the basis for these effects correlates to the El-Sayed rules [54]. These types of effects are likely to play a role in the blends studied here.

We present a model that correlates nanoscale morphology with triplet formation by charge recombination based on known mechanisms for intersystem crossing. Fast SOCT-ISC dominates in strongly (coulomb) bound CT states with a large J (and electronic coupling V). Slower H-HFI-ISC becomes important in charged states where V and the exchange integral (J) are small; see Figure 6.

Visualization of the processes occurring in the blends is shown in Figures 6 and 7.

Clearly, we need further studies on ultra-fast charge and triplet generation in low band gap polymer-PCBM blends probing the effects of morphology on triplet charge recombination (also in working devices). To attain a full understanding of the primary events in organic photovoltaics, it is crucial to study and learn to control charge recombination. The use of deuterated polymers and deuterated PCBM (or pristine C60) should modify H-HFI-ISC. This could result in a spin tool to experimentally discriminate the ISC mechanisms occurring in OPV materials. Furthermore, studies of magnetic field effects on these types of materials will be of great importance.

Conflicts of Interest

The authors declare no competing financial interests.

Authors’ Contributions

René A. J. Janssen, Stephan C. J. Meskers, and René M. Williams designed the experiments. Daniele Di Nuzzo prepared the thin films. Hung-Cheng Chen (supervised by René M. Williams) performed the time resolved measurements and analyzed the data. Hung-Cheng Chen performed the DFT calculations. René M. Williams and Hung-Cheng Chen wrote the paper.

Acknowledgments

The authors are grateful to NWO (de Nederlandse Organisatie voor Wetenschappelijk Onderzoek) for the grant for the femtosecond equipment and the UvA (Universiteit van Amsterdam) for the structural support. They thank Caspar Williams and Lisette van Hooijdonk for language corrections.