Next Article in Journal
Remanence Increase in SrFe12O19/Fe Exchange-Decoupled Hard-Soft Composite Magnets Owing to Dipolar Interactions
Previous Article in Journal
Nanoparticulate Perovskites for Photocatalytic Water Reduction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Time–Frequency Signatures of Electronic Coherence of Colloidal CdSe Quantum Dot Dimer Assemblies Probed at Room Temperature by Two-Dimensional Electronic Spectroscopy

by
James R. Hamilton
1,
Edoardo Amarotti
2,†,
Carlo N. Dibenedetto
3,4,
Marinella Striccoli
3,
Raphael D. Levine
5,
Elisabetta Collini
2 and
Francoise Remacle
1,*
1
Department of Theoretical Physical Chemistry, University of Liège, B4000 Liège, Belgium
2
Department of Chemical Sciences, University of Padova, 35131 Padova, Italy
3
CNR-IPCF SS Bari, c/o Chemistry Department, University of Bari Aldo Moro, 70126 Bari, Italy
4
Chemistry Department, University of Bari Aldo Moro, 70126 Bari, Italy
5
The Fritz Haber Research Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
*
Author to whom correspondence should be addressed.
Present address: Division of Chemical Physics and NanoLund, Lund University, P.O. Box 124, 22100 Lund, Sweden.
Nanomaterials 2023, 13(14), 2096; https://doi.org/10.3390/nano13142096
Submission received: 29 June 2023 / Revised: 13 July 2023 / Accepted: 13 July 2023 / Published: 18 July 2023
(This article belongs to the Section Physical Chemistry at Nanoscale)

Abstract

:
Electronic coherence signatures can be directly identified in the time–frequency maps measured in two-dimensional electronic spectroscopy (2DES). Here, we demonstrate the theory and discuss the advantages of this approach via the detailed application to the fast-femtosecond beatings of a wide variety of electronic coherences in ensemble dimers of quantum dots (QDs), assembled from QDs of 3 nm in diameter, with 8% size dispersion in diameter. The observed and computed results can be consistently characterized directly in the time–frequency domain by probing the polarization in the 2DES setup. The experimental and computed time–frequency maps are found in very good agreement, and several electronic coherences are characterized at room temperature in solution, before the extensive dephasing due to the size dispersion begins. As compared to the frequency–frequency maps that are commonly used in 2DES, the time–frequency maps allow exploiting electronic coherences without additional post-processing and with fewer 2DES measurements. Towards quantum technology applications, we also report on the modeling of the time–frequency photocurrent response of these electronic coherences, which paves the way to integrating QD devices with classical architectures, thereby enhancing the quantum advantage of such technologies for parallel information processing at room temperature.

Graphical Abstract

1. Introduction

Semi-conducting nanoparticles, or quantum dots (QDs), are a promising hardware for a wide variety of quantum technologies [1,2,3,4,5,6,7]. Recently, we proposed to exploit the femtosecond fast-beating electronic coherences in small, few-nm QDs and QD dimers for implementing quantum parallel information processing at room temperature [8,9,10,11], using electronic coherences as logic variables. Our scheme offers a significant quantum advantage, as for a set of N coupled quantum states, N2 − 1 coherences can be used to process information in parallel [8,12].
In previous joint theoretical–experimental studies [7,13,14,15], we reported how electronic coherences could be observed and tuned in ensembles of small (mean diameter, D ¯ = 2.5 to 3 nm), size-dispersed ( σ = 5% to 9%) CdSe QDs and QD quasi-homodimers (Figure 1), addressed in a two-dimensional electronic spectroscopy (2DES) BOXCARS setup (see Supplementary Materials Section S1). In this setup, the polarization response is measured as a function of the first two delay times, T1 and T2, and directly in the frequency domain, ω 3 , of the third delay time, T3, using a CCD camera (see Figure 2a). A Fourier transform along the first delay time, T1, brings the maps into the frequency domain, ω 1 , and allows retrieving the conventional 2DES response as a function of ( ω 1 , T 2 , ω 3 ), as shown in Figure 2d [5,16,17,18]. A good agreement was found between the computed and experimental coherence responses of small ≈3 nm QDs and QD dimers along T2 traces of points extracted at different coordinates in the frequency ( ω 1 , ω 3 ) maps [13,14,15].
The ( ω 1 , ω 3 ) frequency maps of the ( ω 1 , T 2 , ω 3 ) response provide ( ω 1 i , ω 3 j ) coordinates for coherences between pairs of excited electronic states (i,j) evolving along T2, where ω 1 i corresponds to the excitation of state i, ω 1 i = ω i ω G S , and ω 3 j is the emission from state j, ω 3 j = ω j ω G S . Bringing T1 to the frequency domain, however, requires a large sampling of this delay time, as well as computationally expensive FFT post-processing [19]. Here, we show, by comparing experimental and computed data for the BOXCARS setup [20], that the electronic coherences can be equally accurately probed along T2 traces in time–frequency (T1, ω 3 ) polarization maps, which require a smaller number of 2DES measurements and less post-processing. In ( T 1 ,   ω 3 ) time–frequency maps, there can be several coherences between excited states that beat along T2 at a specific address ( T 1 , ω 3 ) in the map. The reason is that the coherences along T2 are only partially resolved: they are resolved along ω 3 but not along T1. All the coherences between the excited states i and j, that have a common state j that emits at the frequency ω 3 j and a different state i, will beat at the same address ( T 1 , ω 3 ) in the map. Therefore, by measuring a single point in the ( T 1 , ω 3 ) time–frequency maps, one can characterize a family of coherences that beat along T2. The ability to simultaneously process all these coherences enhances their potential for exploitation in parallel quantum information processing.
The ability to measure electronic coherences in the directly measured data, without the need to compute full frequency maps, greatly enhances the potential for the use of electronic coherences of QDs in room-temperature quantum technologies. Towards applications to quantum information processing, the recently proposed photocurrent action-based setup [21,22,23,24,25,26,27,28] presents several advantages over the BOXCARS one for 2DES: it is a collinear setup, easier to operate, and the output is a photocurrent that can be processed easier and integrated with a classical computer. It is, therefore, of interest to investigate how electronic coherences are probed in the photocurrent response. We show, computationally, that the electronic coherences can also be robustly probed with the action-based photocurrent response.

2. Materials and Methods

2.1. Synthesis of 3 nm CdSe QDs and Assembly into Quasi-Homodimers in Solution

The colloidal CdSe QDs were prepared in solution by mixing cadmium and selenium precursors using a hot-injection technique [29] in the presence of long-alkyl-chain-coordinating agents. Such molecules act as an organic capping layer on the QD surface, helping the controlled and quasi-epitaxial growth of the nanoparticles, allowing their dispersion in organic solvents, and preventing their aggregation. This approach allows controlling the growth parameters, the precursor ratio, and the temperature of the synthesis, for preparing small, 2.5 to 3.5 nm in diameter, QDs with narrow size distributions varying from 5 to 9% in mean diameter, making any further purification procedure unnecessary [30]. More detailed information can be found in the Supplementary Materials Section S2 and in [31].
Quasi-homodimers of CdSe QDs, with mean diameter D ¯ = 3 nm and size distribution σ = 8%, are assembled in solution by bonding pairs of size-dispersed QDs with a short ( 0.5   n m ) propanedithiol ligand (see [31] for details of this procedure and [7] for the details of the synthesis of the QD dimers used in this work). Due to the unavoidable size dispersion, the QD dimers are not identical nano-objects. A schematic representation of a QD quasi-homodimer constructed in this fashion is shown in Figure 1a, and a transmission electronic microscopy (TEM) image of one of the dimers prepared for this work is shown in Figure 1b. A larger TEM micrograph of the dimers is shown in Figure 1c. Not all the QDs in the ensemble will be bonded together with the ligand to form dimers, and therefore, the dimer solution will retain a proportion of single QD monomers. The ratio of dimers to monomers in the sample is estimated at approximately 60:40.

2.2. 2DES Experimental Methods

The 2DES experiment is implemented using a fully non-collinear setup, in which three fs laser pulses in the visible range are incident upon the dimer solution from different spatial directions, fulfilling the BOXCARS phase-matching conditions [19]. For the 2DES measurements of a QD dimer ensemble assembled into a solid-state multilayer device, see [7]. Additional details about the experimental conditions are reported in the Supplementary Material (SM) Section S1. The experiment probes the optical polarization of the ensemble as a function of the delay times between the pulses, or of their corresponding frequencies. By giving independent control over the delay times, the BOXCARS setup (see Figure S1 of the Supplementary Materials) allows the signal in the different phase-matching directions (PMDs) to be easily extracted. For a description of the experimental setup, see [19], and for details on this implementation, see [7], as well as the Supplementary Materials Section S1.
The three exciting pulses and the final signal observation are separated in time by the delay times, as defined in Figure 2a. The delay time, T0, is the time interval between the arbitrary origin of the time axis, set as the center time of the Local Oscillator (LO), and the first pulse, centered at t1. The LO is the fourth pulse, used as a time reference and employed for heterodyne detection [19].
The first delay time, the excitation or coherence time (T1), is the time separation between the first and second pulses. The second delay time, the population time (T2), is the time between the second and third pulses, and the third delay time, the rephasing or emission time (T3), is the time between the third pulse and the observation. The time at which the measurement is performed is defined as:
t T 0 + T 1 + T 2 + T 3
The polarization response from the ensemble can be measured in specific PMDs as a function of the delay times ( T 1 , T 2 , T 3 ). Repeated measurements with different delays produce a “cube” of data, in which the polarization in a particular PMD is stored as a function of the delay parameters. Figure 2b shows the polarization in the rephasing direction in the time domain, as a function of ( T 1 , T 2 , T 3 ). Each of the delay times ( T 1 , T 2 , T 3 ) can be brought to the frequency domain ( ω 1 , ω 2 , ω 3 ) by Fourier transform, leading to what is sometimes referred to as 3D electron spectroscopy [32]. Figure 2c shows the same response as Figure 2b, only in the time–frequency domain ( T 1 , T 2 , ω 3 ). This is the typical form of the data obtained as raw output of a BOXCARS experimental setup, such as the one used for this work. Indeed, as explained above, in a typical BOXCARS experiment, the signal is measured while scanning the time intervals T1 and T2, whereas the dependence on the third time interval is measured directly in the frequency domain, ω 3 , by the detector [16]. For ease of interpretation, the signal is then typically Fourier-transformed along T1 [5,16,17,18]. This leads to a representation of the same data as a function of the excitation and emission frequencies and the population time ( ω 1 , T 2 , ω 3 ), as shown in Figure 2d. Regardless of the chosen representation, the electronic coherences between the excited electronic states of the QD dimers can be probed by analyzing the data in these cubes as a function of the population time, T 2 , as illustrated in Figure 2b–d. Coherences in the frequency maps ( ω 1 , ω 3 ) (Figure 2d) were studied in our previous work [7], while here, we analyze the ( T 1 , ω 3 ) time–frequency maps.

2.3. Theoretical Methods

We model the electronic structure of each QD from two-hole, one-electron single-particle states. These single-particle electronic states are calculated using an effective mass-k.p Hamiltonian [33,34,35,36], constructed for CdSe using the size distribution of the ensemble (see [13,14,15]). This approximation defines two-holes, one-electron mono-excitons, 1S and 2S, per QD. When the laser intensity is weak enough, the formation of bi-excitons, two-electron, two-hole states, can be neglected [15,37].
The spin 1 / 2   of the hole of each pair is coupled to the p-type orbitals ( l = 1 ) localized on the Se atoms [38,39]. These spin orbit interactions split each of the S bands into two sub-bands of states, with angular momentum: L = 1 ±   1 / 2 state and L = 1 / 2 , 3 / 2 . The total angular momentum, L ± the spin 1 / 2 of the s-type orbital localized on the Cd atoms, leads to an eight-fold degeneracy of the L = 3 / 2 state and a four-fold degeneracy of the L = 1 / 2 state. These states further undergo a loss of degeneracy due to crystal field and Coulomb interactions [36,38,39,40]. The S 3 / 2 state then forms a band of eight fine-structure (FS), singly excited electronic states, of which five are dark and three are bright, and the S 1 / 2 state forms a band of four FS states, of which all are bright. In this way, 4 bands of 24 FS states are formed per QD with energetic ordering: 1 S 3 / 2 , 1 S 1 / 2 , 2 S 3 / 2 , and 2 S 1 / 2 .  Figure 3d shows the stick spectrum of an ensemble-averaged monomer calculated over a 3 nm/8% ensemble of 4000 QDs [8,13].
When 2 QDs drawn from the size-dispersed ensemble of monomers are covalently bonded to form quasi-homodimers, the 24 FS states of each QD are coupled by interdot electronic Coulomb interactions to create a manifold of 48 FS states per dimer. Since the two QDs in a given dimer slightly differ in size, quasi-homodimers are formed: the quasi-isoenergetic bands of each QD are split by Coulomb interdot interactions into a higher and a lower dimer band. This creates eight bands of FS singly excited states per dimer. These bands are energetically ordered as: 1 S 3 / 2 L ,   1 S 3 / 2 H , 1 S 1 / 2 L ,   1 S 1 / 2 H , 2 S 3 / 2 L ,   2 S 3 / 2 H ,   2 S 1 / 2 L , and 2 S 1 / 2 H , although in dimers made from ensembles with larger D ¯ and/or size dispersion, the FS states of these bands interdigitate and the bands overlap [13,15].
The size differences of the two QDs assembled in each quasi-homodimer mean that they do not obey the optical selection rules of exact homodimers and that all the states will share the oscillator strength. Consequently, all singly excited dimer FS states will be bright, although some FS states will be almost dark. Figure 3b shows the stick spectrum calculated for the 3 nm/8% ensemble averaged over an ensemble of 4000 dimers [8,13].
Figure 3a shows the measured absorption spectra of the dimer solution along with the calculated absorption spectrum for the 3 nm/8% ensemble-averaged dimer and monomer, where the inhomogeneous broadening due to the finite size dispersion is taken into account. This figure also shows the spectral profile of the laser pulse used for the 2DES measurements. Only the 1 S 3 / 2 L and the 1 S 3 / 2 H dimer bands and the 1 S 3 / 2 monomer band fall within the laser pulse energy bandwidth. Tables S1 and S2 in Supplementary Materials, Section S3 present the calculated transition energies, inhomogeneous broadening and corresponding dephasing times, and transition dipole moments from the ground to the 1 S 3 / 2 L and 1 S 3 / 2 H band FS states in the dimer, and from the ground to the 1 S 3 / 2 band FS states in the monomer, averaged over 4000 dimer/monomer ensembles.
Since the laser pulses are short, several electronic states of the 1 S 3 / 2 L and the 1 S 3 / 2 H dimer bands are simultaneously excited, which led to a superposition of several excited FS electronic states. The electronic dimer coherences discussed in the Results Section are coherences between the FS excited states, either within (‘intra’-band) or between (‘inter’-band) the 1 S 3 / 2 L and 1 S 3 / 2 H bands.
We specifically focus on four electronic coherences of the dimer, two intra-band and two inter-band, labeled by the periods of their oscillations, which are governed by the energy difference of the two FS states involved. The first intra-band coherences between the 3rd FS state and the degenerate 4th/5th FS states, all in the 1 S 3 / 2 L band, have a frequency of 680   c m 1 . The second intra-band coherences between the degenerate 1st/2nd FS states and the degenerate 4th/5th FS states, all within the 1 S 3 / 2 L band, have a frequency of 840   c m 1 . The first inter-band coherences between the degenerate 4th/5th FS states in the 1 S 3 / 2 L band and the 11th FS state in the 1 S 3 / 2 H band fall at 450   c m 1 . The second inter-band coherences between the degenerate 1st/2nd FS states in the 1 S 3 / 2 L and the 11th FS state in the 1 S 3 / 2 H band have a frequency of 1300   c m 1 .
All these coherences, both intra- and inter-band, are interdot in character, due to the delocalization of the wavefunctions of the 1 S FS states over the whole dimer [7]. The horizontal double-headed arrows superimposed onto the dimer stick spectrum (Figure 3b) identify these four coherences, and their calculated frequencies, periods, inhomogeneous dephasing times, and emission dipole strengths are presented in Table S3 of Section S3 of the Supplementary Materials. Throughout this work, the 450   c m 1 coherences are labeled in cyan, and the 680   c m 1 , 840   c m 1 , and 1300   c m 1 coherences in purple, green, and orange, respectively.
In the experiments, the pulse excites electronic states in the 1 S 3 / 2 band of the monomers that are also present in the sample, thereby creating monomer electronic coherences. The higher-energy monomer bands do not fall within the laser pulse energy bandwidth. Since five of the eight FS states in the 1 S 3 / 2 band are dark, or practically dark (see Table S2 of the Supplementary Materials, Section S3), only one type of electronic coherence will be produced non-negligibly in the monomers. This is the coherence between the 1st/2nd and 3rd FS states, which has a frequency of 180   c m 1 . This coherence, identified in Figure 3c with a double-headed arrow, will be labeled in yellow throughout this paper. The calculated frequency, period, inhomogeneous dephasing time, and emission dipole strength of this coherence are shown in Table S3.
The size dispersion of the QD ensemble causes an inhomogeneous broadening of the energies of the transitions between electronic states in the dimer and in the monomer, the values of which are presented in Tables S1 and S2 in the Supplementary Materials. In the time domain, the size dispersion leads to the dephasing of the coherences, which have a finite lifetime. In the ( T 1 , ω 3 ) time–frequency domain, the size dispersion, therefore, leads to dephasing along T 1 and to an inhomogeneous broadening of emission bands along ω 3 . The inhomogeneous broadening of the bands corresponding to the four dimer coherences on which we focus is represented as shaded areas on the time–frequency ( T 1 , ω 3 ) map reported in Figure 3d. Traces along T2 for points localized in a given ω 3 band on the ( T 1 , ω 3 ) map will, therefore, exhibit beating periods that are characteristic of coherences (i,j) involving an excited state j, that emits in the range of ω 3 values specified by the inhomogeneously broadened transition frequency ( ω j ω G S ) .
Figure 3e shows the emission dipole strengths of the same four dimer electronic coherences. As was the case in [7], the rather monotonic distribution of these values results from the 8% size dispersion of the ensemble, which breaks the exact homodimer limit of fully dark and bright states. The rather even distribution of dipole strengths of the four coherences means that they will appear in the T 2 traces of points on the time–frequency ( T 1 , ω 3 ) maps with commensurate strength. This is a useful feature of disordered QD quasi-homodimers for applications in quantum technologies because it means that more coherences are available for implementing information processing.
In this work, the partial polarizations in specific PMDs are modeled using a phase-modulated approach [41], which is numerically more straightforward to implement than a full simulation of the BOXCARS setup. In the experimental BOXCARS setup, phase modulation is not needed since the different PMDs are spatially separated by using a non-collinear setup [5,16,17,18]. On the other hand, a phase modulation of the train of pulses is experimentally implemented in 2DES collinear setups that measure action observables such as fluorescence or photocurrent [21,22,23,24,26,27,42].
To simulate the polarization response as measured in the BOXCARS setup, we compute the polarization of the ensemble of QD dimers subject to sequences of three collinear fs phase-modulated laser pulses. The electric field time profile of the pulse sequence is given as:
t = n = 1 3 n ( t )
where n ( t ) is the electric field of each pulse, defined as:
n t = 0 e ( t t n ) 2 2 σ n 2 c o s ( ω n t + ϕ n )
In Equation (3), 0 is the electric field strength, σ n is the width of the Gaussian envelope, ω n is the carrier frequency, and t n is the time at which the nth pulse is centered, as shown in Figure 2a. In all the calculations, we use: 0 = 8.775 × 10 7   W/cm2, ω n = 2.36 eV, and σ n = 3.9   f s , for n = 1 ,   2 ,   3 , in agreement with the experimental values.
We modulate the carrier envelope phase of the pulse, ϕ n , for each set of delay times (T1, T2, T3). The modulation of ϕ n is expressed as ϕ n 2 π k n u , where k n m n / L . By choosing the constants m n of each pulse as integer divisors of L , with m 1 m 2 m 3 , after L repetitions for u varying from 1 to L , each of the carrier envelope phases, ϕ n , will have gone through a different number of complete cycles.
The computations are repeated for ranges of delay times ( T 1 , T 2 , T 3 ), modulating the carrier envelope phases of the pulses for each set of values ( T 1 , T 2 , T 3 ). In this way, the polarization of the ensemble is computed as a function of the delay times and phase modulation, u , P u T 1 , T 2 , T 3 . The wave vector of each PMD is presented as a linear combination of the m n :   k l = l 1 m 1 + l 2 m 2 + l 3 m 3 , where the additional factors of u / L modulate the carrier envelope phase of the pulses. For the rephasing direction: l 1 , l 2 , l 3 = ± 1 , + 1 , + 1 .
The polarization of the ensemble in a particular PMD is extracted by Fourier transforming along the phase modulations, u , and identifying the signal by the value of k l .
P k l T 1 , T 2 , T 3 = u = 1 L P u T 1 , T 2 , T 3 · e i 2 π L k l u
One of these P k l T 1 , T 2 , T 3 values is shown in Figure 2b. P u T 1 , T 2 , T 3 is calculated from the time-dependent ensemble density matrix, ρ e n s , as described in [8,13]:
P u t = T r [ μ   ρ e n s ( t ) ]
where μ is the ensemble dipole matrix. The density matrix, ρ e n s ( t ) , is propagated along time, for each set of delay time and phase modulation parameters, using the ensemble Liouville approach [8,13].
The Liouville equation for the ensemble is given by [8]:
i   t ρ e n s n m t = i j L n m , i j e n s · ρ e n s i j ( t )
where L e n s , the ensemble Liouville matrix, is constructed by averaging the Hamiltonian matrices of the individual size-dispersed dimers over the ensemble and taking into account the size dispersion of the QDs. The ensemble Hamiltonian explicitly includes the interaction of the sequence of three laser pulses in the dipole approximation. Equation (6) is solved numerically using the Cash–Karp Runge–Kutta method.
For a given set of delay times, (T1, T2, T3), the polarization is calculated for given sets, u , of carrier envelope phases, ( ϕ 1 , ϕ 2 , ϕ 3 ) , using Equation (5), and parametrized in terms of these delay times using Equation (1).
P u T 1 , T 2 , T 3 = P u t
The polarization in each PMD is obtained by Fourier transform over u using Equation (4). Here, we focus on the rephasing of PMD. The computed data, P r e p h T 1 , T 2 , T 3 , are the time domain “cube” (Figure 2b) and are converted into the time–frequency domain by Fourier transforming the data along T 3 ω 3 (Figure 2c):
P r e p h ( T 1 , T 2 , ω 3 ) e i ω 3 T 3 P r e p h ( T 1 , T 2 , T 3 ) d T 3
This produces a cube of data for the polarization of the ensemble in the rephasing direction in the time–frequency domain, as shown in Figure 2d. These computed data can be directly compared to the data measured in the BOXCARS setup in the time–frequency domain.
A rotating frame (RF) [43] is applied to the measured and calculated P r e p h ( T 1 , T 2 , ω 3 ) (Equation (8)) [19]. A reference frequency, ω r e f , is subtracted along T 1 , thereby detuning the optical frequency.
P r e p h R F = P r e p h e i ω r e f T 1
In the computations, we take ω r e f ω n , the carrier frequency of the pulse (Equation (3)). Working in the RF makes the electronic coherence between excited FS states clearer by removing the fast-beating coherences between the excited FS states and the GS, which have a much shorter dephasing time. This also allows for a less dense sampling in time along T1 and T3.

3. Results and Discussion

3.1. Comparison of Computed and Experimental Time–Frequency Polarization Maps

In our previous work [7], we showed that the electronic coherence response could be consistently characterized in BOXCARS polarization measurements of both solid-state and solution samples. Several electronic coherences between FS states in traces along T2 and their Fourier transforms (FTs) were characterized at specifically addressed points on measured and calculated rephasing frequency maps, P r e p h ω 1 , T 2 , ω 3 (Equation (8)).
Figure 4 compares the real parts of the measured dimer sample (left) and the calculated dimer (right) rephasing ( T 1 , ω 3 ) time–frequency maps at T2 = 20 fs. In this range of T1 values, the main signal in both maps appears at an emission frequency around ω 3 = 18,400 cm−1. This frequency corresponds to the transition energy between the ground state and the strongest dipole FS states in the first 1 S 3 / 2 L dimer band. The experimental and calculated time–frequency maps are in very good agreement. The differences between the measured and calculated maps at emission frequencies below the main signal are attributable to Rayleigh scattering in the measurement.
T2 traces are extracted from the measured and computed maps at points corresponding to the green and pink dots in Figure 4 (see Supplementary Materials Section S4). Figure 4b shows a calculated dimer map, but the same trace was also taken through the calculated monomer maps.
These three traces, i.e., the measured averaged trace and the calculated dimer and monomer traces, are then Fourier-transformed:
e i ω 2 T 2 P r e p h ( T 1 , T 2 , ω 3 ) d T 2 P r e p h ( T 1 , ω 2 , ω 3 )
Three replicates of the 2DES measurements were available. Therefore, we generated three measured FTs (one per replicate), from which a mean and standard deviation are produced. This mean measured FT is plotted with its standard deviation in Figure 5, along with the calculated dimer and monomer FTs.
The mean measured FT shown in Figure 5 has peaks corresponding to all the dimer coherences discussed in the theoretical model section. Their presence is made explicit by the comparison to the calculated dimer FT. The peaks in the mean measured FT at ω 2 = 450 cm−1, 840 cm−1, and 1300 cm−1 are reliably characterized, with a narrow standard deviation. The peak at ω 2 = 680 cm−1 has a larger standard deviation, however it is still distinctive. The region between ω 2 900   a n d   1200   cm 1 has a significantly larger standard deviation and the peaks in this spectral range cannot be reasonably considered. The dimer coherence peaks in the mean measured FT align closely with their corresponding peaks in the calculated dimer FT. The alignment between the experimental and calculated FTs in the frequency domain is very good in the cases of the ω 2 = 450 cm−1, 680 cm−1, and 1300 cm−1 peaks, and within the resolution of the ω 2   points in the measured F T in the case of the 840 cm−1 peak.
The mean measured FT has a broad peak at ω 2 200   cm 1 with a small standard deviation. This signal is primarily caused by the acoustic phonon beating [44]. This phonon signal is not present in the calculated dimer FT because our model does not include the coupling to the two phonon modes. However, in addition to the signal from the acoustic phonon, the 180 cm−1 monomer coherence also contributes to this broad peak, although this contribution is not resolved in the measured data. It should be noted that the coherence between the 1st/2nd and 3rd FS states in the 1 S 3 / 2 is also present in the dimer with a similar frequency. However, as the calculated dimer FT shows, this dimer coherence makes less of a contribution than the monomer coherence.
Figure 5 shows that the four peaks in the measured mean FT at ω 2 = 450 cm−1, 680 cm−1, 840 cm−1, and 1300 cm−1 are caused by coherences in the dimer, as there are no monomer coherences in their ω 2 vicinity.
The dimer coherences discussed above are clearly defined in both the measured and calculated FTs in the time–frequency domain. Producing results in this way, as opposed to using ( ω 1 , ω 3 ) frequency–frequency maps, leads to a drastic reduction in the number of the 2DES measurements needed. The ω 3 coordinate where the measured traces are taken is consistent with the address in the frequency domain and depends upon the transition energies of the excited FS states involved in the coherences to the ground state. The coherences observed in Figure 5 can be best observed within a limited range of T 1 values, because at higher values of T 1 , the fast-beating coherences between the ground state and the mono-excitons generated by the interaction with the first pulse is already de-phased. Hence, intra- and inter-band dimer coherences can be found in the time–frequency maps in a small range of short T 1 values and in a range of ω 3 values which is defined by the transition energies from the FS states involved in the coherences to the GS.
The fact that the same electronic coherences can be characterized in the time–frequency domain as in the frequency–frequency domain of the 2DES experiments, and hence can be exploited in the directly measured data, is advantageous in two respects. The first is that the need for a post-processing step, such as Fourier transforming along T 1 ω 1 , is removed. This yields a computational reduction of O ( N T 1 l o g N T 1 ) , where N T 1 is the number of measurements along T 1 . Following the removal of the FT step, the second advantage gained is the requirement for far fewer measurements along T 1 than are required to sufficiently resolve the FT needed to obtain frequency–frequency maps (see Supplementary Materials Section S5). As it has been discussed, measurements up to T 1 = 10   f s provide an adequate range in which the coherences can be exploited. This range is at least an order of magnitude smaller than the T 1 range required to resolve the FT of T 1 for frequency domain maps, according to the Nyquist limit. As discussed above, in the frequency maps, the ( ω 1 , ω 3 ) address at which a coherence ‘i-j’ will beat along T 2 is given by the values of the ( ω i ω G S ) and ( ω j ω G S ) transition frequencies. In the time–frequency measurements discussed here, since only the third time interval, T3, is Fourier-transformed, only the addresses of the coherences along ω 3 are resolved. All the intra-band and inter-band i-j coherences can be found in traces along T 2 in different ranges of ω 3   but at the same value of T 1 . This partial resolution of the addresses of the coherences significantly reduces the number of time delays that need to be sampled to characterize or exploit them for quantum technology applications.

3.2. Modeling of Action-Based Photocurrent Response

The 2DES action-based fluorescence and photocurrent measurements [21,22,23,24,25,26,27,28,42] are a practical alternative to the polarization BOXCARS ones because they can be implemented with a much simpler collinear setup [5]. In addition, the photocurrent is a more appropriate choice of observable for quantum technology applications, as photocurrent measurements can be directly interfaced with classical electronics. Using the phase modulation approach described above, here, we report on the computed photocurrent response of electronic coherences of an ensemble of monomeric QDs. We set the pulse parameters to access the mono- and bi-exciton states of the QDs and show how electronic coherences involving bi-exciton states can be characterized in the ( T 1 , T 2 , ω 3 ) photocurrent action signal in the double-quantum coherence (DQC) PMD. Note, while the DQC PMD spectra can be easily measured with a phase modulation setup [21], such measurements are harder with the BOXCARS setup because of the lack of fully reliable procedures to correctly phase the signal [45] and the possible contribution of a strong, spurious, non-resonant solvent response [46].
Figure 6 shows the level structure of the mono- and bi-excitonic states of a single CdSe QD of 3 nm. It comprises the ground state 0 , two mono-exciton states, 1 and 2 , the S1 and the S2 states, and the three bi-exciton states, 3 , 4 , and 5 , which correspond to a double-excitation to S1, to an excitation to S1 and to S2, and to a double-excitation to S2, respectively. These five excited states are represented by blue horizontal lines in the figure. The horizontal red dashed lines in Figure 6 show the carrier frequency (and twice the carrier frequency) of the laser pulse (Equation (3)) used in the simulations and the allowed dipole transitions are indicated with green vertical arrows. The calculated transition energies between these states averaged over an ensemble of 4000 monomeric QDs with 8% size dispersion, as well as the energies of the coherences between them, are presented in Table S4, along with the corresponding periods and dipole transition moments. Also presented are the dephasing lifetimes due to the size dispersion.
The energy bandwidth of the laser pulses was selected so that the only transitions energetically allowed are those between the GS and the mono-exciton band and between the mono- and bi-exciton bands, as shown in Figure 6. Intra-band transitions between the states of the mono- or the bi-exciton bands are not resonant with the laser pulse.
Measurements of incoherent actions’ signals, such as fluorescence or, as we show here, photocurrent, require a setup with four pulses [22,23,26,47]. The nonlinear signals in specific PMDs are obtained using the phase modulation approach described above.
The action setup modeled is shown in Figure 7a,b [26]. A sequence of L trains of pulses is incident upon the ensemble, with each train being constituted of four pulses separated in time by the delay times T 1 , T 2 , and T 3 . Each train in the sequence has the same set of delay times, and the trains are spaced apart from one another by the repetition time, t r e p . The total electric field of each train is the same as in Equation (3), only the sum ran over n = 1, 2, 3, and 4. The action signal, from which the photocurrent is calculated, is recorded along t r e p . For a given set of delay times ( T 1 , T 2 , T 3 ) ,   the carrier envelope phases of the four pulses are modulated from u = 1 to L . The sequence is repeated using different values of the delay times, so the total photocurrent response is calculated as a function of the delay times and the phase modulation, P u T 1 , T 2 , T 3 . The total photocurrent is separated into the different PMDs using Equation (4), and these data are Fourier-transformed along T 3 , P P M D T 1 , T 2 , T 3 P P M D ( T 1 , T 2 , ω 3 ) . This post-processing produces a cube of data, in which the photocurrent response in a specific PMD is stored as a function of T 1 , T 2 , ω 3 , as shown in Figure 7c. The analysis of this cube allows identifying and exploiting the coherences contributing to this PMD.
We take here for the pulse parameters: 0 = 5 × 10 6   a.u. (8.775 × 105 W/cm2), ω n = 2.53 eV, and σ n = 3.32   f s , for n = 1 ,   2 ,   3 ,   a n d   4 . The pulse envelope for these parameters in the energy domain is shown in Figure 7d, superimposed onto the stick spectrum of the ground state–mono-exciton and mono–bi-exciton transitions. In the simulations, L = 170 , m 1 = 0 ,   m 2 = 2 , m 3 = 5 ,   and m 4 = 34 .
The observable response being computed is the photocurrent. We took the relaxation times in ranges typical for CdSe QDs [23,47]. The relaxation of the bi-exciton states to the mono-exciton states is in the sub-picosecond range, and the mono-exciton states are assumed to relax to the ground state with a lifetime of a dozen picoseconds. Consequently, the relaxation time from the 3 and 5 bi-exciton states to the 1 and 2 mono-exciton state is fixed to be 318 fs, with corresponding rates of: Γ 3 1 = Γ 5 2 = 1.3 × 10 2   e V . The relaxation time from the 4 bi-exciton state to the 1 and 2 mono-exciton states is taken to be a little longer, at 636 fs, with a corresponding width in energy of: Γ 4 1 = Γ 4 2 = 6.5 × 10 3   eV . The relaxation time from the 1 and 2 mono-exciton states to the 0 ground state is taken to be 15.9 ps, with a corresponding width of: Γ 1 0 = Γ 2 0 = 2.6 × 10 4   eV .
The total photocurrent response is calculated from the density matrix of the ensemble using the approximated Liouville matrix given in Equation (6), to which we added the relaxation rates defined above. Additionally, a decay term was added to the Liouville matrix to account for the dephasing of the coherences caused by phonon coupling with strength γ = 0.005   eV . Both the relaxation of the bi-exciton to the mono-exciton states and of the mono-exciton states to the ground state contribute to the photocurrent signal [23,47]. The action signal response from the relaxation of state m to state n is computed as [47]:
R e s p m n = t 4 t a q u i s i t i o n   d t Γ m n T r n m ρ t m n
The integral of Equation (11) is evaluated along the acquisition time of the experiment, which runs from the time of the fourth pulse in the sequence over the repetition time to the next sequence, t r e p , as shown in Figure 7a. In the simulations, t r e p is approximated as since it is much longer than all the relaxation processes of the mono- and bi-exciton states to the ground state.
The total action signal response is:
R e s p = t 4 t a q u i s i t i o n   d t m , n Γ m n T r n m ρ t m n
The total photocurrent calculated from Equation (12) is separated into its non-linear phase-matching components by Fourier transforming along the phase modulations (Equation (4)). We report here on the DQC PMD: D Q C = m 4 + m 3 m 2 m 1 . As explained before, this is a phase-matching direction, challenging to be reliably measured in the BOXCARS setup, yet very appealing to quantify shifts of the energy correlation between two mono-excitons, in particular to study many-body effects and excited-state landscapes in a wide range of systems, including biomolecules and inorganic materials [19,21,48,49,50].
The calculated cube of data is then converted into the time–frequency domain by Fourier transforming along P D Q C T 1 , T 2 , T 3 P D Q C ( T 1 , T 2 , ω 3 ) .
Figure 8a shows a ( T 1 , T 2 ) time map of the real part of the photocurrent response in the DQC PMD for a value of ω 3 = 20,000   cm 1 . The pathways contributing to the DQC PMD using a third-order perturbative approach are enumerated as double-sided Feynman diagrams by Damtie et al. [47]. We do not use a perturbative approach here, instead we are computing the photocurrent response from the dynamics of the density matrix (Equation (6)). However, the pulse strength used in the simulation is sufficiently weak that the third-order perturbative treatment is a good approximation of the exact time-dependent response and provides good insights into the excitation pathways contributing to the signal. The double-sided Feynman diagrams which contribute to the DQC PMD all have in common that the first pulse excites the ket from the ground state to a mono-exciton state, and the second pulse excites this ket from the mono-exciton state to a bi-exciton state. This means that for the duration of T 1 , the system will be in a coherence between the ground state and a mono-exciton state, and that for the duration of T 2 , the system will be in a coherence between the ground state and a bi-exciton state. These features clearly appear in the time map of Figure 8a. If one compares the number of oscillations along the T 1 axis to those along the T 2 axis, it is clear that there are more beatings in the same time along T 2 than T 1 . This reflects the fact that the ground–bi-exciton coherences have much shorter periods, about two times shorter, than the ground–mono-exciton coherences (see Table S4).
Figure 8b shows the absolute values of the traces along ω 3 of two points on the ( T 1 , T 2 ) time map. The Feynman diagrams [47] show that the third pulse in the train can either relax the ket from the bi-exciton state to a mono-exciton state or excite the bra from the ground state to a mono-exciton state. This means that for the duration of T 3 ,   the system can either be in a coherence between the ground and a mono-exciton state or between a mono-exciton and a bi-exciton state. This is shown in Figure 8b, in which the fast-beating responses corresponding to ground–mono-exciton as well as ground–bi-exciton coherences are identified.
Note that the signal of the two traces plotted in Figure 8b, while being dominated by the same primary frequencies, differs in the specific coherences which can be identified. This small but important variability of the coherences in the ω 3 traces of different ( T 1 , T 2 ) points means that a comprehensive collection of coherence data still requires repetitions for a small number of T 1 and T 2   values, albeit far fewer than would be needed if Fourier transforms along these delay times were required. In addition to the inhomogeneous broadening resulting from the size dispersion of the ensemble, the peaks in Figure 8b are broaden by the coupling to the phonon modes.
The analysis of the computed photocurrent response in an additional PMD, the DQC PMD, shows that coherences involving bi-exciton states can be observed in the action signal photocurrent response of QD ensembles to 2DES. As in the case of the polarization response, directly measured data in the ( T 1 , T 2 , ω 3 ) domain are usable without the need for additional post-processing by Fourier transforms along T 1 and T 2 . A sampling of the delay time, T3, for fixed values of T1 and T2 enables to characterize all the coherences found in the pathways contributing to the DQC PMD. This, in turn, means that fewer measurements along these delay times are needed. Furthermore, these results show that the ( T 1 , T 2 , ω 3 ) is the appropriate domain for looking at coherences in observables in the DQC PMD.

4. Conclusions

We have shown that the electronic coherences resulting from the excitation of ensembles of size-dispersed QDs and QD dimers by sequences of fs broad-bandwidth laser pulses, as in 2DES, can be observed and characterized in the directly measured time–frequency domain. This is the case for both polarization and action-based (here, photocurrent) measurements. For the case of the cube of data that depends on three delay times, the electronic coherences are characterized by robust and distinct beating patterns in the traces of the signal, as a function of one delay time at a single point in the time–frequency domain of the other two delay times. In the case of polarization measurements, we obtained very good agreement between the modeled and experimental ( T 1   ,   ω 3 ) time–frequency maps in the rephasing PMD, and in the T 2 / ω 2   traces along these maps.
Compared to the conventional ( ω 1   ,   ω 3 ) frequency maps, the advantage of directly exploiting time–frequency data leads to a considerable reduction of the number of time intervals that need to be sampled. This is because, in a time–frequency map, the addresses of the coherences are only partially resolved. All the coherences between excited states i and j, which have a common state j, beat in the traces of points with the same T1 value. By fixing ω 3 to be within the inhomogeneous broadening of the transition j to the GS, and recording a trace at a point (T1, ω 3 ) along T2, one can characterize, at the same time, all the beating frequencies of a family of coherences that involve the same emitting state j in the chosen PMD, because the absorbing states i are not resolved along T1.
That several coherences are simultaneously accessible in the directly measured time–frequency domain is a huge advantage for the exploitation of QD electronic coherences in quantum technologies for parallel information processing. The savings in fewer measurements and less computation time will greatly enhance the advantage of logical operations encoded onto coherences over classical logic operations. Action-based photocurrent measurements of electronic coherences are a step further in coherence exploitation since they pave the way for QD devices to be integrated into classical architectures.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nano13142096/s1, Section S1: Experimental details on the 2DES set-up. Section S2: Materials and method for the synthesis of CdSe QDs and dimers preparation. Section S3: Tables of the computed energy levels of the QD dimers. Section S4: Extraction of T2 traces and their corresponding ω3 FTs from the measured and calculated time-frequency maps. Section S5: Scaling of data processing for time-frequency responses (T1, T2, ω3) and frequency-frequency responses (ω1, T2, ω3). References [7,19,20,30,31] are cited in the supplementary materials.

Author Contributions

Conceptualization, F.R., E.C. and R.D.L.; sample synthesis and preparation, C.N.D. and M.S.; 2DES experiments, E.C. and E.A.; theory and modeling, F.R., R.D.L. and J.R.H.; software, J.R.H.; validation, F.R., J.R.H. and E.C.; formal analysis, F.R., J.R.H. and E.C.; writing—review and editing, F.R., R.D.L., J.R.H., E.C., M.S. and C.N.D.; funding acquisition, F.R., R.D.L., E.C. and M.S. All authors have read and agreed to the published version of the manuscript.

Funding

The work of J.R.H. and F.R. is supported by the Fonds National de la Recherche, F.R.S.—FNRS (Belgium), #T.0205.20. Computational resources have been provided by the Consortium des Equipements de Calcul Intensif (CECI), funded by the F.R.S.—FNRS under Grant #2.5020.11. E.C. acknowledges partial support by “CQ-TECH” Supporting Talent in Research at the University of Padova (STARS), Grant 2019 (2019-UNPD0Z9-0166571). E.C. and M.S. acknowledge the Project NEST—NETWORK 4 ENERGY SUSTAINABLE TRANSITION SPOKE 9: ENERGY-SUSTAINABLE ADVANCED MATERIALS, code PE0000021, Concession Decree No. 1561 of 11 October 2022 by Ministero dell’Università e della Ricerca (MUR) (CUP C93C22005230007—University of Padua, CUP B53C22004060006—CNR).

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank Rosaria Brescia of the Italian Institute of Technology for the high-resolution TEM image in Figure 1b, and Hugo Gattuso for discussions on the modeling of the electronic structure of the QDs.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rebentrost, P.; Stopa, M.; Aspuru-Guzik, A. Förster Coupling in Nanoparticle Excitonic Circuits. Nano Lett. 2010, 10, 2849–2856. [Google Scholar] [CrossRef] [Green Version]
  2. Spittel, D.; Poppe, J.; Meerbach, C.; Ziegler, C.; Hickey, S.G.; Eychmüller, A. Absolute Energy Level Positions in Cdse Nanostructures from Potential-Modulated Absorption Spectroscopy (Emas). ACS Nano 2017, 11, 12174–12184. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Kagan, C.R.; Bassett, L.C.; Murray, C.B.; Thompson, S.M. Colloidal Quantum Dots as Platforms for Quantum Information Science. Chem. Rev. 2021, 121, 3186–3233. [Google Scholar] [CrossRef]
  4. Koley, S.; Cui, J.; Panfil, Y.E.; Banin, U. Coupled Colloidal Quantum Dot Molecules. Acc. Chem. Res. 2021, 54, 1178–1188. [Google Scholar] [CrossRef] [PubMed]
  5. Collini, E. 2D Electronic Spectroscopic Techniques for Quantum Technology Applications. J. Phys. Chem. C 2021, 125, 13096–13108. [Google Scholar] [CrossRef]
  6. Dibenedetto, C.N.; Fanizza, E.; De Caro, L.; Brescia, R.; Panniello, A.; Tommasi, R.; Ingrosso, C.; Giannini, C.; Agostiano, A.; Curri, M.L. Coupling in Quantum Dot Molecular Hetero-Assemblies. Mater. Res. Bull. 2022, 146, 111578. [Google Scholar] [CrossRef]
  7. Hamilton, J.R.; Amarotti, E.; Dibenedetto, C.N.; Striccoli, M.; Levine, R.D.; Collini, E.; Remacle, F. Harvesting a Wide Spectral Range of Electronic Coherences with Disordered Quasi-Homo Dimeric Assemblies at Room Temperature. Adv. Quantum Technol. 2022, 5, 2200060. [Google Scholar] [CrossRef]
  8. Gattuso, H.; Levine, R.D.; Remacle, F. Massively Parallel Classical Logic Via Coherent Dynamics of an Ensemble of Quantum Systems with Dispersion in Size. Proc. Natl. Acad. Sci. USA 2020, 117, 21022. [Google Scholar] [CrossRef]
  9. Komarova, K.; Gattuso, H.; Levine, R.D.; Remacle, F. Quantum Device Emulates the Dynamics of Two Coupled Oscillators. J. Phys. Chem. Lett. 2020, 11, 6990–6995. [Google Scholar] [CrossRef] [PubMed]
  10. Komarova, K.; Gattuso, H.; Levine, R.D.; Remacle, F. Parallel Quantum Computation of Vibrational Dynamics. Front. Phys. 2020, 8, 486. [Google Scholar] [CrossRef]
  11. Remacle, F.; Levine, R.D. A Quantum Information Processing Machine for Computing by Observables. Proc. Natl. Acad. Sci. USA 2023, 120, e2220069120. [Google Scholar] [CrossRef] [PubMed]
  12. Fresch, B.; Hiluf, D.; Collini, E.; Levine, R.D.; Remacle, F. Molecular Decision Trees Realized by Ultrafast Electronic Spectroscopy. Proc. Natl. Acad. Sci. USA 2013, 110, 17183–17188. [Google Scholar] [CrossRef] [PubMed]
  13. Collini, E.; Gattuso, H.; Levine, R.D.; Remacle, F. Ultrafast Fs Coherent Excitonic Dynamics in Cdse Quantum Dots Assemblies Addressed and Probed by 2d Electronic Spectroscopy. J. Chem. Phys. 2021, 154, 014301. [Google Scholar] [CrossRef] [PubMed]
  14. Collini, E.; Gattuso, H.; Bolzonello, L.; Casotto, A.; Volpato, A.; Dibenedetto, C.N.; Fanizza, E.; Striccoli, M.; Remacle, F. Quantum Phenomena in Nanomaterials: Coherent Superpositions of Fine Structure States in Cdse Nanocrystals at Room Temperature. J. Phys. Chem. C 2019, 123, 31286–31293. [Google Scholar] [CrossRef] [Green Version]
  15. Collini, E.; Gattuso, H.; Kolodny, Y.; Bolzonello, L.; Volpato, A.; Fridman, H.T.; Yochelis, S.; Mor, M.; Dehnel, J.; Lifshitz, E. Room-Temperature Inter-Dot Coherent Dynamics in Multilayer Quantum Dot Materials. J. Phys. Chem. C 2020, 124, 1622–16231. [Google Scholar] [CrossRef]
  16. Gelzinis, A.; Augulis, R.; Butkus, V.; Robert, B.; Valkunas, L. Two-Dimensional Spectroscopy for Non-Specialists. Biochim. Et Biophys. Acta (BBA)—Bioenerg. 2019, 1860, 271–285. [Google Scholar] [CrossRef]
  17. Brańczyk, A.M.; Turner, D.B.; Scholes, G.D. Crossing Disciplines—A View on Two-Dimensional Optical Spectroscopy. Ann. Der Phys. 2014, 526, 31–49. [Google Scholar] [CrossRef] [Green Version]
  18. Hybl, J.D.; Albrecht, A.W.; Faeder, S.M.G.; Jonas, D.M. Two-Dimensional Electronic Spectroscopy. Chem. Phys. Lett. 1998, 297, 307–313. [Google Scholar] [CrossRef] [Green Version]
  19. Bolzonello, L.; Volpato, A.; Meneghin, E.; Collini, E. Versatile Setup for High-Quality Rephasing, Non-Rephasing, and Double Quantum 2d Electronic Spectroscopy. J. Opt. Soc. Am. B 2017, 34, 1223–1233. [Google Scholar] [CrossRef]
  20. Eckbreth, A.C. Boxcars: Crossed-Beam Phase-Matched Cars Generation in Gases. Appl. Phys. Lett. 2008, 32, 421–423. [Google Scholar] [CrossRef]
  21. Mueller, S.; Draeger, S.; Ma, X.; Hensen, M.; Kenneweg, T.; Pfeiffer, W.; Brixner, T. Fluorescence-Detected Two-Quantum and One-Quantum–Two-Quantum 2d Electronic Spectroscopy. J. Phys. Chem. Lett. 2018, 9, 1964–1969. [Google Scholar] [CrossRef] [Green Version]
  22. 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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Karki, K.J.; Widom, J.R.; Seibt, J.; Moody, I.; Lonergan, M.C.; Pullerits, T.; Marcus, A.H. Coherent Two-Dimensional Photocurrent Spectroscopy in a Pbs Quantum Dot Photocell. Nat. Commun. 2014, 5, 5869. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. 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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Bian, Q.; Ma, F.; Chen, S.; Wei, Q.; Su, X.; Buyanova, I.A.; Chen, W.M.; Ponseca, C.S., Jr.; Linares, M.; Karki, K.J. Vibronic Coherence Contributes to Photocurrent Generation in Organic Semiconductor Heterojunction Diodes. Nat. Commun. 2020, 11, 617. [Google Scholar] [CrossRef] [Green Version]
  26. Bolzonello, L.; Bernal-Texca, F.; Gerling, L.G.; Ockova, J.; Collini, E.; Martorell, J.; Van Hulst, N.F. Photocurrent-Detected 2d Electronic Spectroscopy Reveals Ultrafast Hole Transfer in Operating Pm6/Y6 Organic Solar Cells. J. Phys. Chem. Lett. 2021, 12, 3983–3988. [Google Scholar] [CrossRef] [PubMed]
  27. Bruschi, M.; Gallina, F.; Fresch, B. Simulating Action-2d Electronic Spectroscopy of Quantum Dots: Insights on the Exciton and Biexciton Interplay from Detection-Mode and Time-Gating. Phys. Chem. Chem. Phys. 2022, 24, 27645–27659. [Google Scholar] [CrossRef] [PubMed]
  28. Karki, K.J.; Ciappina, M.F. Advances in Nonlinear Spectroscopy Using Phase Modulated Light Fields: Prospective Applications in Perturbative and Non-Perturbative Regimes. Adv. Phys. X 2022, 7, 2090856. [Google Scholar] [CrossRef]
  29. Murray, C.B.; Norris, D.J.; Bawendi, M.G. Synthesis and Characterization of Nearly Monodisperse Cde (E = Sulfur, Selenium, Tellurium) Semiconductor Nanocrystallites. J. Am. Chem. Soc. 1993, 115, 8706–8715. [Google Scholar] [CrossRef]
  30. Roberge, A.; Dunlap, J.H.; Ahmed, F.; Greytak, A.B. SizeDependent Pbs Quantum Dot Surface Chemistry Investigated Via Gel Permeation Chromatography. Chem. Mater. 2020, 32, 6588–6594. [Google Scholar] [CrossRef]
  31. Dibenedetto, C.N.; Fanizza, E.; Brescia, R.; Kolodny, Y.; Remennik, S.; Panniello, A.; Depalo, N.; Yochelis, S.; Comparelli, R.; Agostiano, A.; et al. Coupling Effects in Qd Dimers at Sub-Nanometer Interparticle Distance. Nano Res. 2020, 13, 1071–1080. [Google Scholar] [CrossRef]
  32. Irgen-Gioro, S.; Spencer, A.P.; Hutson, W.O.; Harel, E. Coherences of Bacteriochlorophyll a Uncovered Using 3d-Electronic Spectroscopy. J. Phys. Chem. Lett. 2018, 9, 6077–6081. [Google Scholar] [CrossRef]
  33. Luttinger, J.M.; Kohn, W. Motion of Electrons and Holes in Perturbed Periodic Fields. Phys. Rev. 1955, 97, 869–883. [Google Scholar] [CrossRef] [Green Version]
  34. Brus, L. Electronic Wave Functions in Semiconductor Clusters: Experiment and Theory. J. Phys. Chem. 1986, 90, 2555–2560. [Google Scholar] [CrossRef]
  35. Zunger, A. Electronic Structure Theory of Semiconductor Quantum Dots. MRS Bull. 1998, 23, 35–42. [Google Scholar] [CrossRef]
  36. Efros Al, L.; Rosen, M. The Electronic Structure of Semi-Conducting Nanocrystal. Ann. Rev. Mater. Sci. 2000, 30, 475–521. [Google Scholar] [CrossRef]
  37. Righetto, M.; Bolzonello, L.; Volpato, A.; Amoruso, G.; Panniello, A.; Fanizza, E.; Striccoli, M.; Collini, E. Deciphering Hot-and Multi-Exciton Dynamics in Core–Shell Qds by 2d Electronic Spectroscopies. Phys. Chem. Chem. Phys. 2018, 20, 18176–18183. [Google Scholar] [CrossRef] [Green Version]
  38. Efros, A.L.; Rosen, M.; Kuno, M.; Nirmal, M.; Norris, D.J.; Bawendi, M. Band-Edge Exciton in Quantum Dots of Semiconductors with a Degenerate Valence Band: Dark and Bright Exciton States. Phys. Rev. B 1996, 54, 4843–4856. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Norris, D.J.; Bawendi, M.G. Measurement and Assignment of the Size-Dependent Optical Spectrum in Cdse Quantum Dots. Phys. Rev. B 1996, 53, 16338–16346. [Google Scholar] [CrossRef] [Green Version]
  40. Wong, C.Y.; Scholes, G.D. Using Two-Dimensional Photon Echo Spectroscopy to Probe the Fine Structure of the Ground State Biexciton of Cdse Nanocrystals. J. Lumin. 2011, 131, 366–374. [Google Scholar] [CrossRef]
  41. Kjellberg, P.; Brüggemann, B.; Pullerits, T. Two-Dimensional Electronic Spectroscopy of an Excitonically Coupled Dimer. Phys. Rev. B 2006, 74, 024303. [Google Scholar] [CrossRef]
  42. Draeger, S.; Roeding, S.; Brixner, T. Rapid-Scan Coherent 2d Fluorescence Spectroscopy. Opt. Express 2017, 25, 3259–3267. [Google Scholar] [CrossRef]
  43. Calhoun, T.R.; Ginsberg, N.S.; Schlau-Cohen, G.S.; Cheng, Y.C.; Ballottari, M.; Bassi, R.; Fleming, G.R. Quantum Coherence Enabled Determination of the Energy Landscape in Light-Harvesting Complex Ii. J. Phys. Chem. B 2009, 113, 16291–16295. [Google Scholar] [CrossRef]
  44. Mork, A.J.; Lee, E.M.; Tisdale, W.A. Temperature Dependence of Acoustic Vibrations of Cdse and Cdse–Cds Core–Shell Nanocrystals Measured by Low-Frequency Raman Spectroscopy. Phys. Chem. Chem. Phys. 2016, 18, 28797–28801. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Nemeth, A.; Milota, F.; Mančal, T.; Pullerits, T.; Sperling, J.; Hauer, J.; Kauffmann, H.F.; Christensson, N. Double-Quantum Two-Dimensional Electronic Spectroscopy of a Three-Level System: Experiments and Simulations. J. Chem. Phys. 2010, 133, 094505. [Google Scholar] [CrossRef] [PubMed]
  46. Gellen, T.A.; Bizimana, L.A.; Carbery, W.P.; Breen, I.; Turner, D.B. Ultrabroadband Two-Quantum Two-Dimensional Electronic Spectroscopy. J. Chem. Phys. 2016, 145, 064201. [Google Scholar] [CrossRef]
  47. 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. [Google Scholar] [CrossRef] [Green Version]
  48. Kim, J.; Mukamel, S.; Scholes, G.D. Two-Dimensional Electronic Double-Quantum Coherence Spectroscopy. Acc. Chem. Res. 2009, 42, 1375–1384. [Google Scholar] [CrossRef] [Green Version]
  49. Stone, K.W.; Gundogdu, K.; Turner, D.B.; Li, X.; Cundiff, S.T.; Nelson, K.A. Two-Quantum 2d Ft Electronic Spectroscopy of Biexcitons in Gaas Quantum Wells. Science 2009, 324, 1169. [Google Scholar] [CrossRef]
  50. Christensson, N.; Milota, F.; Nemeth, A.; Pugliesi, I.; Riedle, E.; Sperling, J.; Pullerits, T.; Kauffmann, H.F.; Hauer, J. Electronic Double-Quantum Coherences and Their Impact on Ultrafast Spectroscopy: The Example of Β-Carotene. J. Phys. Chem. Lett. 2010, 1, 3366–3370. [Google Scholar] [CrossRef]
Figure 1. (a) Schematic representation of a CdSe quasi-homodimer. The dimer is assembled by covalently bonding two QDs with a 0.55 nm S C H 3 2 S ligand. The QDs are drawn from an ensemble with a mean diameter, D ¯ , of 3 nm and size dispersion, σ . (b) A high-resolution transmission electronic microscopy (HR-TEM) image of one of the quasi-homodimers prepared for this work. (c) A larger TEM micrograph of the CdSe QDs dimers.
Figure 1. (a) Schematic representation of a CdSe quasi-homodimer. The dimer is assembled by covalently bonding two QDs with a 0.55 nm S C H 3 2 S ligand. The QDs are drawn from an ensemble with a mean diameter, D ¯ , of 3 nm and size dispersion, σ . (b) A high-resolution transmission electronic microscopy (HR-TEM) image of one of the quasi-homodimers prepared for this work. (c) A larger TEM micrograph of the CdSe QDs dimers.
Nanomaterials 13 02096 g001
Figure 2. (a) Time ordering of the three fs laser pulses as used in the 2DES experiment and in the modeling. (bd) Different possible representations of the final three-dimensional signal obtained after a 2DES measurement (or of modeled data). The final cube of polarization response data can be cast as a function of: (b) the three time delays between pulses (T1, T2, T3) (Equation (4)), and (c) the first two time delays between pulses and the third delay in the frequency domain (T1, T2, ω 3 ) (Equation (8)). These are the data studied in this work. (d) The first and third delays in the frequency domain and the second delay in the time domain ( ω 1 , T2,   ω 3 ). This is the representation typically used when 2DES measurements are published, as was the case in [7]. In panels (bd), the green, purple and orange lines schematically represent the beatings of three different electronic coherences along T2.
Figure 2. (a) Time ordering of the three fs laser pulses as used in the 2DES experiment and in the modeling. (bd) Different possible representations of the final three-dimensional signal obtained after a 2DES measurement (or of modeled data). The final cube of polarization response data can be cast as a function of: (b) the three time delays between pulses (T1, T2, T3) (Equation (4)), and (c) the first two time delays between pulses and the third delay in the frequency domain (T1, T2, ω 3 ) (Equation (8)). These are the data studied in this work. (d) The first and third delays in the frequency domain and the second delay in the time domain ( ω 1 , T2,   ω 3 ). This is the representation typically used when 2DES measurements are published, as was the case in [7]. In panels (bd), the green, purple and orange lines schematically represent the beatings of three different electronic coherences along T2.
Nanomaterials 13 02096 g002
Figure 3. (a) Measured (violet) absorption spectra of the QD dimer sample and the calculated absorption spectra of the QD dimers (azure) and monomers (green), with the spectral profile of the laser pulse (orange) used for the 2DES measurement. (b) Calculated stick spectrum for the averaged dimer computed from an ensemble of 4000 3 nm/8% dimers. (c) Calculated stick spectrum for the averaged monomer computed from an ensemble of 4000 3 nm/8% monomers. The five electronic coherences between specific electronic FS states discussed in this work (four in the dimer, one in the monomer) are identified with horizontal, double-arrowed lines between their constituent FS states. (d) Inhomogeneous broadening of the four dimer electronic coherences in frequency, ω3. Patterns are superimposed onto a (T1, ω3) map, indicating the range of ω3 in which the four dimer electronic coherences can be found. The vertical grey pattern indicates the range corresponding to the inhomogeneous broadening of the 450 cm−1 and 680 cm−1 coherences, and the horizontal grey pattern indicates the range of the 840 cm−1 and 1300 cm−1 coherences. All four dimer coherences will be found in the overlapping square pattern region. (e) Emission dipole strengths of the four dimer electronic coherences. The color labeling of the FS bands and coherences in this figure is used throughout.
Figure 3. (a) Measured (violet) absorption spectra of the QD dimer sample and the calculated absorption spectra of the QD dimers (azure) and monomers (green), with the spectral profile of the laser pulse (orange) used for the 2DES measurement. (b) Calculated stick spectrum for the averaged dimer computed from an ensemble of 4000 3 nm/8% dimers. (c) Calculated stick spectrum for the averaged monomer computed from an ensemble of 4000 3 nm/8% monomers. The five electronic coherences between specific electronic FS states discussed in this work (four in the dimer, one in the monomer) are identified with horizontal, double-arrowed lines between their constituent FS states. (d) Inhomogeneous broadening of the four dimer electronic coherences in frequency, ω3. Patterns are superimposed onto a (T1, ω3) map, indicating the range of ω3 in which the four dimer electronic coherences can be found. The vertical grey pattern indicates the range corresponding to the inhomogeneous broadening of the 450 cm−1 and 680 cm−1 coherences, and the horizontal grey pattern indicates the range of the 840 cm−1 and 1300 cm−1 coherences. All four dimer coherences will be found in the overlapping square pattern region. (e) Emission dipole strengths of the four dimer electronic coherences. The color labeling of the FS bands and coherences in this figure is used throughout.
Nanomaterials 13 02096 g003
Figure 4. (a) Measured dimer sample and (b) calculated dimer ( T 1 , ω 3 ) maps at T 2 = 20 fs after the rotating frame has been applied. The green and pink dots correspond to the locations where traces along T2 were taken, from which the FTs in Figure 5 were produced.
Figure 4. (a) Measured dimer sample and (b) calculated dimer ( T 1 , ω 3 ) maps at T 2 = 20 fs after the rotating frame has been applied. The green and pink dots correspond to the locations where traces along T2 were taken, from which the FTs in Figure 5 were produced.
Nanomaterials 13 02096 g004
Figure 5. FTs of the time traces through the measured (Equation (10)) and calculated maps. The blue and green lines show the calculated dimer and monomer FTs, respectively, transformed along T 2 at the coordinate ( T 1 = 7.8   f s , ω 3 = 17,186   cm 1 ) on the time–frequency maps. The pink line is the mean of FTs from three measurements of the dimer solution sample, and the shaded area is their standard deviation, see Supplementary Materials Section S4 for more details.
Figure 5. FTs of the time traces through the measured (Equation (10)) and calculated maps. The blue and green lines show the calculated dimer and monomer FTs, respectively, transformed along T 2 at the coordinate ( T 1 = 7.8   f s , ω 3 = 17,186   cm 1 ) on the time–frequency maps. The pink line is the mean of FTs from three measurements of the dimer solution sample, and the shaded area is their standard deviation, see Supplementary Materials Section S4 for more details.
Nanomaterials 13 02096 g005
Figure 6. Energy level diagram for a monomeric QD excitonic system. The blue lines show the mono- and bi-exciton electronic states of the QD. The red dashed line shows the carrier frequency (and twice the carrier frequency) of the pulses (Equation (3)) used to compute the 2DES response.
Figure 6. Energy level diagram for a monomeric QD excitonic system. The blue lines show the mono- and bi-exciton electronic states of the QD. The red dashed line shows the carrier frequency (and twice the carrier frequency) of the pulses (Equation (3)) used to compute the 2DES response.
Nanomaterials 13 02096 g006
Figure 7. (a) Sequence of L phase modulation of 4 pulse trains used in the simulation. Each train has a different value of u, from 0 to L, which modulated the carrier envelope phase, and the trains are spaced by the laser repetition time, trep. (b) The 4 pulses in each train are separated by the delay times, T1, T2, and T3, which are fixed in each sequence of L trains. The measurement is repeated using sequences with different sets of delay times. (c) The observable response in a specific PMD stored in a cube of data as a function of T 1 , T 2   ,   a n d   ω 3 . (d) The calculated stick spectrum with ground state–mono-exciton transitions in blue and mono-exciton–bi-exciton transitions in pink. The pulse envelope in the energy domain is superimposed onto the stick spectrum in green.
Figure 7. (a) Sequence of L phase modulation of 4 pulse trains used in the simulation. Each train has a different value of u, from 0 to L, which modulated the carrier envelope phase, and the trains are spaced by the laser repetition time, trep. (b) The 4 pulses in each train are separated by the delay times, T1, T2, and T3, which are fixed in each sequence of L trains. The measurement is repeated using sequences with different sets of delay times. (c) The observable response in a specific PMD stored in a cube of data as a function of T 1 , T 2   ,   a n d   ω 3 . (d) The calculated stick spectrum with ground state–mono-exciton transitions in blue and mono-exciton–bi-exciton transitions in pink. The pulse envelope in the energy domain is superimposed onto the stick spectrum in green.
Nanomaterials 13 02096 g007
Figure 8. (a) ( T 1 , T 2 ) Time map of the real part of the photocurrent response in the DQC PMD for a value of ω 3 = 20,000   cm 1 . (b) Traces along ω 3 at the points indicated in (a) with a green dot, ( T 1 , T 2 ) = ( 6.0   f s , 6.0   f s ), and a burgundy dot, ( T 1 , T 2 ) = ( 8.0   f s , 8.0   f s ). The signal corresponding to ground–mono-exciton state coherences is identified with blue arrows, and the signal corresponding to mono–bi-exciton state coherences is identified with pink arrows.
Figure 8. (a) ( T 1 , T 2 ) Time map of the real part of the photocurrent response in the DQC PMD for a value of ω 3 = 20,000   cm 1 . (b) Traces along ω 3 at the points indicated in (a) with a green dot, ( T 1 , T 2 ) = ( 6.0   f s , 6.0   f s ), and a burgundy dot, ( T 1 , T 2 ) = ( 8.0   f s , 8.0   f s ). The signal corresponding to ground–mono-exciton state coherences is identified with blue arrows, and the signal corresponding to mono–bi-exciton state coherences is identified with pink arrows.
Nanomaterials 13 02096 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hamilton, J.R.; Amarotti, E.; Dibenedetto, C.N.; Striccoli, M.; Levine, R.D.; Collini, E.; Remacle, F. Time–Frequency Signatures of Electronic Coherence of Colloidal CdSe Quantum Dot Dimer Assemblies Probed at Room Temperature by Two-Dimensional Electronic Spectroscopy. Nanomaterials 2023, 13, 2096. https://doi.org/10.3390/nano13142096

AMA Style

Hamilton JR, Amarotti E, Dibenedetto CN, Striccoli M, Levine RD, Collini E, Remacle F. Time–Frequency Signatures of Electronic Coherence of Colloidal CdSe Quantum Dot Dimer Assemblies Probed at Room Temperature by Two-Dimensional Electronic Spectroscopy. Nanomaterials. 2023; 13(14):2096. https://doi.org/10.3390/nano13142096

Chicago/Turabian Style

Hamilton, James R., Edoardo Amarotti, Carlo N. Dibenedetto, Marinella Striccoli, Raphael D. Levine, Elisabetta Collini, and Francoise Remacle. 2023. "Time–Frequency Signatures of Electronic Coherence of Colloidal CdSe Quantum Dot Dimer Assemblies Probed at Room Temperature by Two-Dimensional Electronic Spectroscopy" Nanomaterials 13, no. 14: 2096. https://doi.org/10.3390/nano13142096

APA Style

Hamilton, J. R., Amarotti, E., Dibenedetto, C. N., Striccoli, M., Levine, R. D., Collini, E., & Remacle, F. (2023). Time–Frequency Signatures of Electronic Coherence of Colloidal CdSe Quantum Dot Dimer Assemblies Probed at Room Temperature by Two-Dimensional Electronic Spectroscopy. Nanomaterials, 13(14), 2096. https://doi.org/10.3390/nano13142096

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop