Differentiation of True Nonlinear and Incoherent Mixing of Linear

dynamics during the long measurement time of the action signals can also contribute ... Since its first implementation,1,2 different phase modulation ...
0 downloads 0 Views 900KB Size
Subscriber access provided by UNIV AUTONOMA DE COAHUILA UADEC

A: Kinetics, Dynamics, Photochemistry, and Excited States

Differentiation of True Nonlinear and Incoherent Mixing of Linear Signals in Action-Detected 2D Spectroscopy Alex Arash Sand Kalaee, Fikeraddis Damtie, and Khadga Jung Karki J. Phys. Chem. A, Just Accepted Manuscript • DOI: 10.1021/acs.jpca.9b01129 • Publication Date (Web): 18 Apr 2019 Downloaded from http://pubs.acs.org on April 18, 2019

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 19 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

The Journal of Physical Chemistry

Differentiation of True Nonlinear and Incoherent Mixing of Linear Signals in Action-Detected 2D Spectroscopy Alex Arash Sand Kalaee,∗,† Fikeraddis Damtie,† and Khadga Jung Karki∗,‡ †Mathematical Physics and NanoLund, Lund University, Box 118, 22100 Lund, Sweden ‡Chemical Physics Lund University, P.O. Box 124, 22100 Lund, Sweden E-mail: [email protected]; [email protected]

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

Abstract Phase-modulation and phase-cycling schemes have been commonly used in electronic two-dimensional (2D) spectroscopy where the observables are incoherent signals such as fluorescence or photo-current. Although the methods have distinct advantages compared to the coherent signal detected 2D spectroscopy in sensitivity, possibility to measure spectra from isolated quantum systems and direct visualization of the contributions from the different states to the action signals, ambiguities in interpreting the spectra have emerged. Recent reports have shown that apart from the nonlinear signals from the four pulse interactions, mixing of the linear signals due to nonlinear population dynamics during the long measurement time of the action signals can also contribute to the measured 2D spectra. Exciton-exciton annihilation has been considered to play a major role in the mixing of the linear signals. Thus, it has become important to further characterize the origin of the measured signals. Here, using a non-perturbative simulation of the 2D spectra based on the time evolution of the density matrix in the Lindblad form, we show that the exciton-exciton annihilation contributes to the measured signal only if the quantum yields of the different excited states are not the same. In these cases, the mixed signals can be distinguished from the true nonlinear signals if the phases are measured with respect to the linear signals. In an action-detected 2D spectra, the mixed signals have a π phase shift relative to the true nonlinear signals. A detailed discussion on the experimental implementations of the schemes used in the simulations is also provided.

2

ACS Paragon Plus Environment

Page 2 of 19

Page 3 of 19 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

The Journal of Physical Chemistry

Graphical TOC Entry TOC ENTRY REQUIRED

3

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

Introduction. Since its first implementation, 1,2 different phase modulation schemes 1,3,4 have been used in the Fourier transform spectroscopy of molecules, 5,6 molecular aggregates 7,8 and photoactive devices. 9–12 Electronic 2D spectroscopy 2,7–10 based on four phase modulated beams has recently become popular. Action signals, such as fluorescence and photo-current, have been used in most of the measurements. Similarly, phase cycling schemes 13 have been used in 2D spectroscopy in conjunction with the detection of photo-ions 14 and molecular fragments. 15 One of the far reaching advantages of action-detected 2D spectroscopy is that the method is technically capable of measuring in-situ ultra-fast dynamics in operational devices, 10 isolated systems and even single molecules. However, the research field is still in infancy and relatively few experimental results have been reported. Most of the reported results have shown changes in the 2D spectra at characteristic relaxation time scales of the excited states indicating that the signals originate from nonlinear interactions between the sequence of the four laser pulses and the sample. 5,7,8,10,14,15 Nevertheless, some of the reports have only shown 2D spectra at a single population time, 9 and others have reported no change in the spectra along the population time, 12 which seems to suggest that the applicability of the method to investigate the ultra-fast dynamics may be specific to appropriate systems. In phase modulated 2D spectroscopy, see Fig 1, the phases of four collinear pulsed beams are modulated at four frequencies, φ1 , φ2 , φ3 and φ4 , by using acousto-optic modulators in a dual Mach-Zehnder interferometer. 2,7,10 The interaction of the beams with the sample modulates incoherent signals at frequencies given by ±n1 φ1 ± n2 φ2 ± n3 φ3 ± n4 φ4 , where ni are integers. Out of all the modulated signals, the signals at φ1 − φ2 − φ3 + φ4 and −φ1 + φ2 − φ3 + φ4 , which are also known as the rephasing and the non-rephasing signals, result from the nonlinear interactions with the four beams. Amplitudes and phases of these signals are measured as a function of the inter-pulse time delays to obtain 2D spectra in the time domain. The Fourier transform of the time domain data yields 2D spectra in the frequency domain. 4

ACS Paragon Plus Environment

Page 4 of 19

Page 5 of 19 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

The Journal of Physical Chemistry

Figure 1: Schematic of a setup for action-detected 2D spectroscopy. A sequence of four pulses whose phases are modulated at frequencies φ1 , φ2 , φ3 and φ4 interact with the sample imparting modulations in the amplitude of incoherent responses such as fluorescence and photo-current. The modulated response from the sample at the frequencies φ1 − φ2 − φ3 + φ4 and −φ1 + φ2 − φ3 + φ4 are recorded as a function of inter-pulse time delays t21 , t32 and t43 . The Fourier transform of the data along t21 and t43 provides the 2D spectra at different intervals of t32 . The nonlinear response from a system induced by four-wave mixing of light fields can be described by three different interaction pathways commonly called as the ground state bleach (GSB), stimulated emission (SE), and excited state absorption (ESA). In action-detected 2D spectroscopy using phase modulation, 7,10 the ESA consists of two pathways, ESAI and ESAII, with different signs such that the total amplitude of the ESA signal depends on the quantum yields of the signals from the different states. This contrasts with the 2D spectroscopy based on the detection of coherent response, where only one ESA pathway contributes. 16–18 Note that the ESAI pathway corresponds to the bleaching of population in state |1i in model b of Fig 2 (see Methods below) and ESAII corresponds to promoting the population to state |3i in the same model. Particularly in fluorescence detected 2D spectroscopy of excitonic systems with monomolecular recombination kinetics, the quantum yields from both the pathways are often the same such that they cancel each other and the final signal is a sum of the GSB and SE pathways. This is generally not the case in photo-current detected 2D spectroscopy. 10 It is important to consider that the linear signals at frequencies φ21 = φ2 − φ1 and φ43 = φ4 − φ3 , which originate from two-field interactions, are orders of magnitude stronger than the nonlinear signals from four-field interactions. As

5

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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 19

the action-signals are recorded over a significantly longer time compared to the typical ultrafast dynamics in the system, other slow nonlinear population dynamics can incoherently mix the linear signals. The mixed signals appear at the same modulation frequencies as the true nonlinear signals. It has been shown that such effects can dominate the photo-current detected 2D spectra in semiconductors that have high carrier mobility. 19 Similar effects can prevail in any action-detected nonlinear spectroscopy of semiconductors, metal particles and molecular aggregates. Thus, there is an urgent need to develop a methodology to identify if the mixed signals have dominant contributions in the spectra. 19 Here, by explicitly simulating an experiment, we show that the nonlinear signals in the 2D spectra are phase shifted by π with respect to the mixed signals. The phase shifts can be unambiguously identified in an experiment if the linear as well as the nonlinear signals relevant for the 2D spectra are recorded simultaneously without applying any phase correction during the lock-in detection.

a

b 3.10 eV τ21 1.55 eV τ10 0 eV

2 µ21 1 µ10

5 µ54 4 µ43

0 3 Separated system

2.92 eV µ32

τ54

µ31

3 2

µ21

1.46 eV

1

3.01 eV τ32 1.55 meV τ21 1.46 meV τ10

τ43

µ20

0 eV

0 eV 0 Equiv. four lvl. sys.

µ10

Figure 2: Energy diagram for the two model systems, (a) the two-molecule system and (b) an equivalent four level system. The elements µkl indicate the relevant dipole matrix elements in each model and τkl indicate relaxation time scales from level k to l.

Methods. We consider two different models of an excitonic system that consist of two closely spaced excited states, see Fig. 2. The energy levels at 1.55 eV and 1.46 eV are similar to the levels 6

ACS Paragon Plus Environment

Page 7 of 19 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

The Journal of Physical Chemistry

corresponding to the outer and inner rings, respectively, of bacteriochlorophylls in the lightharvesting complex of purple bacteria. 20 Recent measurements of fluorescence-detected 2D spectra have shown peaks across the diagonal, 5,21 the origin of which has been debated. It has been suggested that either a correlated excitation of the two rings 21 or exciton-exciton annihilation after the energy transfer from the outer ring to the inner ring could result in the off-diagonal peaks. 22 We have used model a (Fig. 2 a) to simulate the 2D spectra that includes contributions from the exciton-exciton annihilation as a representative process that can mix the linear signals. The model has two independent systems, each with three eigenstates that are denoted by |0i, |1i and |2i, and |3i, |4i and |5i, respectively. The transition dipole matrix elements are adjusted such that µkl E0 = 8 meV for all the included elements, E0 is the electric field strength, and the corresponding relaxation rates are τ10 = τ43 = 1 ns and τ21 = τ54 = 100 fs. Similarly, we have used model b (Fig. 2 b) with four eigenstates, |0i, |1i, |2i and |3i, to simulate the spectra of correlated excitation of the two rings. The model assumes that all the states can be excited from a common ground state. The transition matrix elements have the same size as in a and the relaxation rates are τ10 = 160 ps, τ21 = 160 fs and τ32 = 100 fs. The density matrix, ρˆ(t), is propagated in time by using the Lindblad Master Equation. 23 The simulations are done following the method developed by Damtie et al. 24 and the details are given in the Supplementary Information (SI).

Results and Discussion. In previous simulations, the fluorescence signal has been calculated by summing the recombination events. 24 This method, however, cannot be used to simulate the mixing of linear signals that is proportional to the product of the population densities in the low lying excited states of model a, Sx ∝ ρ11 × ρ44 . 25 Thus, in order to obtain the contributions of such processes in the 2D spectra, we have used the modulations in the density of populations imparted by the train of phase modulated pulses. In order to rigorously show that the two

7

ACS Paragon Plus Environment

The Journal of Physical Chemistry

signals are equivalent, we compare the 2D spectra obtained from them. Fig. 3 shows the spectra obtained for model a from the recombination signals Rel10 (a) and Rel43 (b) with the spectra obtained from the corresponding population density signals ρ11 (c) and ρ44 (d). Only the analysis of real part of the total correlation spectra, i.e. the sum of the rephasing and non-rephasing spectra, are presented throughout this work. Apart from a difference in scale, the shape and sign of the spectral surface plot for the two pairs of spectra share same topological characteristics, which indicates similar dynamics. Hence, we consider the modulations in the density of populations as a viable proxy for the externally measured signals including the fluorescence.

667 0 -667 -1334 -1334 -667 0 667 1334 Wavenumber k1 (cm−1) 1334

×10−6 0.0 1334 Wavenumber k3 (cm−1)

(a)

×10−6 0.0 −0.5 −1.0 −1.5 −2.0 −2.5 −3.0 −3.5 −4.0 −4.5

667 0 -667 -1334

−0.6

667

−1.2

0

−2.4

−1.8 −3.0

-667

−3.6 −4.2 -1334 -667 0 667 1334 Wavenumber k1 (cm−1)

0.00 −0.25 −0.50 −0.75 −1.00 −1.25 −1.50 −1.75 −2.00 −2.25

(c)

(b)

-1334

1334 Wavenumber k3 (cm−1)

Wavenumber k3 (cm−1)

1334

Wavenumber k3 (cm−1)

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 19

(d)

667 0 -667 -1334

−4.8

0.0 −0.3 −0.6 −0.9 −1.2 −1.5 −1.8 −2.1 −2.4

-1334 -667 0 667 1334 Wavenumber k1 (cm−1)

-1334 -667 0 667 1334 Wavenumber k1 (cm−1)

Figure 3: Comparison of the signals for (a) relaxations Rel10 and (b) Rel43 with the corresponding population densities (c) ρ11 and (d) ρ44 . The mixing of linear signals due to exciton-exciton annihilation in a molecular aggregate

8

ACS Paragon Plus Environment

Page 9 of 19 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

The Journal of Physical Chemistry

shown in model a can be understood as a two-step process where the first step is an energy transfer from state |1i to state |4i, which at the same time de-excites the first molecule and excites the second molecule to the high lying state |5i as shown by the solid red arrows in Fig. 4. The excitation transferred to |5i is given by S5x ∝ ρ11 × ρ44 . As the populations in the states |1i and |4i are modulated at frequencies φ21 and φ43 , S5x exhibits modulations at the frequencies φ43 ± φ21 , which are identical to the modulations in the non-rephasing and rephasing nonlinear signals. Note that the population in |5i is transferred from |4i, a process which also imparts modulations in the population density of |4i, S4x ∝ ρ11 × ρ44 , at φ43 ±φ21 . However, the modulations in |4i and |5i are in anti-phase (see Fig.4). Similar to the two ESA pathways, whether the exciton-exciton annihilation contributes to the measured signal depends on the quantum yields, η5 and η4 , of the states |5i and |4i, respectively. Some authors have wrongly assumed a universal contribution from the annihilation process in fluorescence detected 2D spectra of molecular dimers by analyzing only the population of |5i, 22 which is not based on an experimental observable. For example, typical quantum yields of fluorescence from both the states in a molecular system is the same such that η5 − η4 = 0; consequently the annihilation does not contribute to the measured nonlinear signals. On the other hand, quantum yields of photo-current or photo-ion from the two states are generally not equal and in most of the cases η5 > η4 . Therefore, the contributions of the annihilation in such signals are given by Sx ∝ ηS5x , where 0 ≤ η. Fig. 5 (a) shows the 2D spectra of the mixed signal Sx from model a obtained by setting η = 1 and (b) shows the corresponding spectra obtained from the nonlinear population signal ρ11 from model b. Considering Fig. 5 (a), we note strong off-diagonal peaks from exciton-exciton annihilation. As model b only allows a single exciton to be present, excitonexciton annihilation is not possible. Hence, the obtained signal is the true non-linear signal, which provides a reference for the mixed signal and its spectral line shape. Comparison with Fig. 5 (b) reveals that the inclusion of exciton-exciton annihilation shifts the phase of the signal by π such that the mixed 2D spectra has an opposite sign. Similarly, the amplitudes

9

ACS Paragon Plus Environment

The Journal of Physical Chemistry

Figure 4: Schematic depiction of mixing of linear transitions in model a. In exciton-exciton annihilation energy is transferred from |1i to |4i, which de-excites the first molecule and excites the second molecule to the high lying state |5i. The second molecule subsequently relaxes back to |4i via fast non-radiative vibrational relaxation. The modulations in the population densities at |4i and |5i at the mixed frequencies are in anti-phase. of the peaks in (a) are significantly smaller than in (b). From the relative sizes of the offdiagonal peaks in the two spectra we observe that if the quantum yield of Sx is about an order of magnitude larger than that of the true non-linear signal then the annihilation signals dominate in the off-diagonal peaks of the total spectrum. Thus, in a typical measurement of a system where Sx and true non-linear are of similar magnitude, the resulting 2D spectra has dominant diagonal peaks and possibly weak off-diagonal peaks with an opposite sign. 0

0.64

(a)

0.56

667

0.48 0.40

0

0.32 0.24

-667

0.16 0.08

-1334 -1334 -667 0 667 1334 Wavenumber k1 (cm−1)

1334 Wavenumber k3 (cm−1)

1334 Wavenumber k3 (cm−1)

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 10 of 19

0.00

(b)

−2

667

−4

0

−8

−6 −10 −12

-667

−14

-1334 -1334 -667 0 667 1334 Wavenumber k1 (cm−1)

−16

Figure 5: Comparison of (a) combined signals Sx and (b) population density ρ11 for model b. Finally we outline a procedure that can be used in experiments to distinguish between the mixed and the nonlinear signals. This distinction is not possible in the current im-

10

ACS Paragon Plus Environment

Page 11 of 19 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

The Journal of Physical Chemistry

plementations of the experimental setups mainly because of the adjustment of the phase during data processing. The phases of the signals are adjusted such that the real parts are maximized during the overlap of the pulses. This renders the main features of the total correlation spectra to be positive irrespective of whether it originates from the mixed or the true nonlinear signal. As we have mentioned previously, the linear signals at φ21 and φ43 are orders of magnitude stronger than the nonlinear signals at φ43 ± φ21 . Typically in a measurement, the linear signals are filtered out and only the nonlinear signals are recorded by using lock-in amplifiers. These signals, however, are important in assigning the absolute phase of the nonlinear signals. In order to do so, a different method of data acquisition has to be implemented where the signal from the detector is digitized and the data is post-processed using algorithms of fast-Fourier transforms or the generalized lock-in amplifiers. 26–28 Time domain 2D spectra are then constructed from the amplitudes and phases of the linear signals at φ21 and φ43 , and the nonlinear signals at φ43 ± φ21 . Multiplication of the 2D spectra from the linear signals give the reference spectra at φ43 ± φ21 . We have tested the method using the data from the simulations. Fig. 6 compares the 2D spectra of model a (see Fig. 2) obtained from the annihilation signal Sx (a), the true nonlinear signal (b), ESAII pathway (c) and the reference (d). The spectra from the true nonlinear signal in (b) is obtained from the nonlinear modulations of the populations in the lowest excited states ρ11 and ρ44 . These signals arise from the GSB, SE and ESAI pathways. Similarly, the spectra in (c) are obtained from the nonlinear modulations of the populations in the high-lying excited states ρ22 and ρ55 . These signals arise from the ESAII pathways. Although the two spectra have similar features, they have different signs. Note that the spectra measured in the experiments is the sum of (b) and (c), and generally if the quantum yields of the observable from high lying states are the same as that from the lowest excited state, the two ESA pathways cancel resulting in the spectra that has the same sign as in (b). When we compare the reference spectra in (d) with the nonlinear spectra in (b) and (c), we find that the sign of the spectra from the ESAII pathways is the same as the reference spectra. Thus, from the comparison

11

ACS Paragon Plus Environment

The Journal of Physical Chemistry

we can identify if the measured spectra is dominated by the ESAII pathways or the usual GSB, SE and ESAI pathways. Similarly the sign of the spectra in (a) is similar to that of the reference (d), which also allows us to distinguish the annihilation contributions from the GSB, SE and ESAI pathways. However, the sign alone is not sufficient to distinguish between exciton-exciton annihilation and ESAII. In this case, the prevalence of annihilation can be inferred from the prominent off-diagonal peaks in the spectra. 0.64 0.56

667

0.48 0.40

0

0.32 0.24

-667

0.16 0.08

-1334 -1334 -667 0 667 1334 Wavenumber k1 (cm−1) 1334

1334 Wavenumber k3 (cm−1)

(a)

0.56 0.48

667

0.40 0

0.32 0.24

-667

0.16 0.08

-1334 -1334 -667 0 667 1334 Wavenumber k1 (cm−1)

0 -667 -1334

0.0 −0.3 −0.6 −0.9 −1.2 −1.5 −1.8 −2.1 −2.4

-1334 -667 0 667 1334 Wavenumber k1 (cm−1)

0.64

(c)

(b)

667

0.00

1334 Wavenumber k3 (cm−1)

Wavenumber k3 (cm−1)

1334

Wavenumber k3 (cm−1)

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 19

0.00

(d)

667 0 -667 -1334 -1334 -667 0 667 1334 Wavenumber k1 (cm−1)

2250 2000 1750 1500 1250 1000 750 500 250 0

Figure 6: 2D spectra from the mixed signals due to (a) exciton-exciton annihilation and (b) the true nonlinear signals based on the populations ρ11 + ρ44 , (c) population signal ρ22 + ρ55 that corresponds to ESAII and (d) the reference spectra obtained from the modulations at φ21 =500 kHz and φ43 =800 kHz of the populations ρ11 + ρ44 . See the text for the details.

12

ACS Paragon Plus Environment

Page 13 of 19 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

The Journal of Physical Chemistry

Conclusion To conclude, we have shown that exciton-exciton annihilation does not contribute significantly in the 2D spectra of molecular systems where quantum yields of the observable are same for all the excited states. In systems with different quantum yields, the signatures of annihilation can be identified if the signs (or phase) of the signals are compared with the reference spectra constructed from the linear signals. The results address an important and unresolved ambiguity in identifying the origin of the signals measured in action-detected 2D spectroscopy. The methodology we have presented can be applied in any experimental setup of action-detected 2D spectroscopy to distinguish the contributions from the mixing of linear signals, ESAII, and GSB and SE pathways.

Supporting Information Available The time evolution of the system under the Lindblad master equation approach is given in the Supplementary Information. This material is available free of charge via the Internet at http://pubs.acs.org/.

Acknowledgement Financial support from the Swedish Research Council (VR), the Crafoord Foundation, the Knut and Wallenberg Foundation and NanoLund is gratefully acknowledged. Authors also thank Prof. Andreas Wacker for guidance in implementing the Lindblad Master Equations.

References (1) Tekavec, P. F.; Dyke, T. R.; Marcus, A. H. Wave packet interferometry and quantum state reconstruction by acousto-optic modulation. J. Chem. Phys. 2006, 125, 194303.

13

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

(2) Tekavec, P. F.; Lott, G. A.; Marcus, A. H. Fluorescence-detected two-dimensional electronic coherence spectroscopy by acousto-optic phase modulation. J. Chem. Phys. 2007, 127, 214307. (3) Karki, K. J.; Kringle, L.; Marcus, A. H.; Pullerits, T. Phase-synchronous detection of coherent and incoherent nonlinear signals. J. Opt. 2016, 18, 015504. (4) Osipov, V. A.; Shang, X.; Hansen, T.; Pullerits, T.; Karki, K. J. Nature of relaxation processes revealed by the action signals of intensity-modulated light fields. Phys. Rev. A 2016, 94, 1–9. (5) Tiwari, V.; Matutes, Y. A.; Gardiner, A. T.; Jansen, T. L. C.; Cogdell, R. J.; Ogilvie, J. P. Spatially-resolved fluorescence-detected two-dimensional electronic spectroscopy probes varying excitonic structure in photosynthetic bacteria. Nat. Commun. 2018, 9, 4219. (6) Kumar, P.; Karki, K. J. Two-photon excitation spectroscopy of 1, 5–diphenyl-1, 3, 5-hexatriene using phase modulation. J. Phys. Commun. 2019, 3, 035008 (7) Perdomo-Ortiz, A.; Widom, J. R.; Lott, G. A.; Aspuru-Guzik, A.; Marcus, A. H. Conformation and electronic population transfer in membrane supported self-assembled porphyrin dimers by 2D fluorescence spectroscopy. J. Phys. Chem. B 2012, 116, 10757– 10770. (8) Widom, J. R.; Johnson, N. P.; von Hippel, P. H.; Marcus, A. H. Solution conformation of 2-aminopurine dinucleotide determined by ultraviolet two-dimensional fluorescence spectroscopy. New J. Phys. 2013, 15, 1–16. (9) Nardin, G.; Autry, T. M.; Silverman, K. L.; Cundiff, S. T. Multidimensional coherent photocurrent spectroscopy of a semiconductor nanostructure. Opt. Express 2013, 21, 28617–28627.

14

ACS Paragon Plus Environment

Page 14 of 19

Page 15 of 19 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

The Journal of Physical Chemistry

(10) Karki, K. J.; Widom, J. R.; Seibt, J.; Moody, I.; Lonergan, M. C.; onu Pullerits, T.; Marcus, A. H. Coherent two-dimensional photocurrent spectroscopy in a PbS quantum dot photocell. Nat. Commun. 2014, 5 . (11) Bakulin, A. A.; Silva, C.; Vella, E. Ultrafast spectroscopy with photocurrent detection: Watching excitonic optoelectronic systems at work. J. Phys. Chem. Lett. 2016, 7, 250– 258. (12) Vella, E.; Li, H.; Grégoire, P.; Tuladhar, S. M.; Vezie, M. S.; Few, S.; Bazán, C. M.; Nelson, J.; Silva-Acuña, C.; Bittner, E. R. Ultrafast decoherence dynamics govern photocarrier generation efficiencies in polymer solar cells. Sci. Rep. 2016, 6, 29437. (13) Tian, P.; Keusters, D.; Suzaki, Y.; Warren, W. S. Femtosecond phase-coherent two dimensional spectroscopy. Science 2003, 300, 1553–1555. (14) Aeschlimann, M.; Brixner, T.; Fischer, A.; Kramer, C.; Melchior, P.; Pfeiffer, W.; Schneider, C.; Strüber, C.; Tuchscherer, P.; Voronine, D. V. Coherent Two-Dimensional Nanoscopy. Science 2011, 333, 1723–1726. (15) Roeding, S.; Brixner, T. Coherent two-dimensional electronic mass spectrometry. Nat. Commun. 2018, 9, 2519. (16) Jonas, D. M. Two-dimensional femtosecond spectroscopy. Annu. Rev. Phys. Chem. 2003, 54, 425–463. (17) Cho, M. Coherent Two-dimensional optical spectroscopy. Chem. Rev. 2008, 108, 1331– 1418. (18) Hochstrasser, R. M. Two-dimensional spectroscopy at infrared and optical frequencies. Proc. Natl. Acad. Sci. U.S.A. 2007, 104, 14190–14196.

15

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

(19) Grégoire, P.; Srimath Kandada, A. R.; Vella, E.; Tao, C.; Leonelli, R.; Silva, C. Incoherent population mixing contributions to phase-modulation two-dimensional coherent excitation spectra. The Journal of Chemical Physics 2017, 147, 114201. (20) McDermott, G.; Prince, S. M.; Freer, A. A.; Hawthornthwaite-Lawless, A. M.; Papiz, M. Z.; Cogdell, R. J.; Isaacs, N. W. Crystal structure of an integral membrane light-harvesting complex from photosynthetic bacteria. Nature 1995, 374, 517–521. (21) Karki, K. J.; Chen, J.; Sakurai, A.; Shi, Q.; Gardiner, A. T.; Kühn, O.; Cogdell, R. J.; Pullerits, T. Unexpectedly large delocalization of the initial excitation in photosynthetic light harvesting. ArXiv e-prints 2018, (22) Maly, P.; Mancal, T. Signatures of Exciton Delocalization and Exciton-Exciton Annihilation in Fluorescence-Detected Two-Dimensional Coherent Spectroscopy. The Journal of Physical Chemistry Letters 2018, 9, 5654–5659. (23) Lindblad, G. On the generators of quantum dynamical semigroups. Commun. Math. Phys. 1976, 48, 119–130. (24) Damtie, F. A.; Wacker, A.; Pullerits, T.; Karki, K. J. Two-dimensional action spectroscopy of excitonic systems: Explicit simulation using a phase-modulation technique. Phys. Rev. A 2017, 96, 053830. (25) van Amerongen, H.; van Grondelle, R.; Valkunas, L. Photosynthetic excitons; World Scientific: Singapore, 2000. (26) Karki, K. J.; Torbjörnsson, M.; Widom, J. R.; Marcus, A. H.; Pullerits, T. Digital cavities and their potential applications. J. Instrum. 2013, 8, 1–12. (27) Fu, S.; Sakurai, A.; Liu, L.; Edman, F.; Pullerits, T.; Öwall, V.; Karki, K. J. Generalized lock-in amplifier for precision measurement of high frequency signals. Rev. Sci. Instrum. 2013, 84 . 16

ACS Paragon Plus Environment

Page 16 of 19

Page 17 of 19 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

The Journal of Physical Chemistry

(28) Jin, A.; Fu, S.; Sakurai, A.; Liu, L.; Edman, F.; Pullerits, T.; Öwall, V.; Karki, K. J. Note: High precision measurements using high frequency gigahertz signals. Rev. Sci. Instrum. 2014, 85 .

17

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

Figure 7

Action-detected 2D spectra depicting exciton-exciton annihilation (bottom-left) and nonlinear ground-state bleach and stimulated emission (bottom-right) from isolated molecular systems with two energy levels (top). The spectra have a relative phase shift of π.

18

ACS Paragon Plus Environment

Page 18 of 19

Page 19 ofThe 19 Journal of Physical Chemistry 1 2 3 4 5 6 7

ACS Paragon Plus Environment