Next Article in Journal
Identification of Chemical Constituents in Blumea balsamifera Using UPLC–Q–Orbitrap HRMS and Evaluation of Their Antioxidant Activities
Previous Article in Journal
Effect of Chemical Refining on the Reduction of β-Carboline Content in Sesame Seed Oil
Previous Article in Special Issue
Time and Angle-Resolved Time-of-Flight Electron Spectroscopy for Functional Materials Science
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretical Description of Attosecond X-ray Absorption Spectroscopy of Frenkel Exciton Dynamics

by
Tim Hansen
1,
Tatiana Bezriadina
1,2 and
Daria Popova-Gorelova
1,2,*
1
I. Institute for Theoretical Physics, Universität Hamburg, Notkestr. 9, 22607 Hamburg, Germany
2
Centre for Ultrafast Imaging, Luruper Chaussee 149, 22671 Hamburg, Germany
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(11), 4502; https://doi.org/10.3390/molecules28114502
Submission received: 26 April 2023 / Revised: 26 May 2023 / Accepted: 29 May 2023 / Published: 1 June 2023
(This article belongs to the Special Issue Ultrafast X-ray Spectroscopy)

Abstract

:
Frenkel excitons are responsible for the transport of light energy in many molecular systems. Coherent electron dynamics govern the initial stage of Frenkel-exciton transfer. Capability to follow coherent exciton dynamics in real time will help to reveal their actual contribution to the efficiency of light-harvesting. Attosecond X-ray pulses are the tool with the necessary temporal resolution to resolve pure electronic processes with atomic sensitivity. We describe how attosecond X-ray pulses can probe coherent electronic processes during Frenkel-exciton transport in molecular aggregates. We analyze time-resolved absorption cross section taking broad spectral bandwidth of an attosecond pulse into account. We demonstrate that attosecond X-ray absorption spectra can reveal delocalization degree of coherent exciton transfer dynamics.

1. Introduction

When a photovoltaic material or a biological complex interacts with light, a bound electron-hole pair, called exciton, can be created. Dynamics of excitons determine processes that are essential for both applications such as solar cells, light-emitting devices and field effect transistors and for the biological significant processes, such as photosynthesis [1,2]. Ultrafast energy transfer after exciton creation governs the functionality of these processes. The role of quantum coherence processes for energy transfer has been heavily debated [3,4,5,6,7,8,9,10,11]. Coherences can occur right after the exciton creation as well as be induced by the environment some time after exciton propagation [12]. Although it is not fully clear, whether coherence processes can last for hundreds of femtoseconds, they certainly play a role on femtosecond time scales before electron-phonon coupling starts to be considerable. Resolving the finest details of pure electronic processes will help to clarify their contribution to the efficiency of light harvesting.
Since recent years, it has become possible to generate attosecond X-ray pulses either at high-harmonic generation sources or at free-electron lasers [13,14,15]. Attosecond X-ray pulses offer new perspectives for revealing mechanisms of excited-state dynamics [16,17,18,19]. Sub-femtosecond temporal resolution provides an access to the ultrashort time scales of electron dynamics [20,21]. In addition, X-ray absorption spectroscopy method encodes information about a probed system with atomic selectivity and chemical sensitivity [22,23,24,25]. Although it is clear that attosecond X-ray absorption spectroscopy carries important information about dynamics in a system, the interpretation of the spectra are challenging. The bandwidths of spectral peaks are inversely proportional to the duration of a probe pulse and become several eV broad for an attosecond probe. Thus, the understanding of the information encoded by attosecond X-ray pulses relies on a thoughtful theoretical analysis. In this article, we study how attosecond X-ray absorption spectroscopy can be applied to probe coherent exciton dynamics in molecular aggregates and demonstrate that X-ray spectra provide information about exciton delocalization in such systems.
Frenkel exciton model [26] is a commonly used and robust model that successfully describes exciton dynamics in organic semiconductors, molecular aggregates or light-harvesting systems [9,27,28,29]. In such systems, interaction with light leads to an excitation of a molecule in an aggregate. Due to electron couplings between neighbouring sites, the excitation becomes delocalized over several molecules, which leads to exciton transfer through the system. Frenkel exciton transport has been studied at time scales of few hundreds of femtoseconds or longer [30,31,32,33,34]. Here, we propose to probe Frenkel exciton transport at much shorter time scales of few femtoseconds, when incoherent processes due to phonons do not start to play a role.
The article is organized as follows. In Section 2, we start with a brief introduction to coherent dynamics of Frenkel excitons. Then, we derive the general expression for time-resolved X-ray absorption cross section for Frenkel exciton system using the first-order time dependent perturbation theory. We proceed with a detailed analysis of the time-resolved X-ray absorption cross section in Section 3. The analysis relies on the study of core-excited final states that can be reached due to X-ray absorption. We connect the time-evolution of the X-ray absorption cross section to the time evolution of Frenkel excitons in Section 3.3.

2. Materials and Methods

2.1. Frenkel Exciton Dynamics

Let us consider a chain of N identical equidistantly aligned molecules. In the Frenkel exciton model, it is assumed that every molecule m can be either in the ground state | ϕ 0 m or in the optically-excited state | ϕ m = c ^ m h ^ m | ϕ 0 m [9,26,27,28,29]. c ^ m ( c ^ m ) are electron creation (annihilation) operator for site m and h ^ m ( h ^ m ) are hole creation (annihilation) operator for site m. It is assumed that the distance between the molecules is large enough to neglect exchange interaction between electrons of different molecules. The ground state of our system is
| Ψ 0 = | ϕ 0 1 | ϕ 0 2 | ϕ 0 N .
We now introduce the basis states of Frenkel Hamiltonian
| Φ m = | ϕ 0 1 | ϕ m | ϕ 0 N .
These are the states, in which one of molecules in the chain is excited and others are in the ground state as shown in Figure 1. The Frenkel Hamiltonian in this basis set is expressed as
H ^ 0 = m U | Φ m Φ m | + V m [ | Φ m Φ m ± 1 | + h . c . ] .
Here, U is the Coulomb energy of an electron-hole pair localized on a single site and V is the coupling between nearest-neighbour molecules due dipole-dipole interaction. The coupling term is responsible for the transfer of the electron-hole pair from site m to the sites m + 1 and m 1 .
The eigenstates of the Frenkel Hamiltonian, Frenkel excitons, are the states, in which the excitation is delocalized over several sites. They can be expressed as a linear combination of the basis states:
| Ψ n = m c n m | Φ m .
with corresponding eigenenergies E n . The Frenkel exciton model assumes that only a single molecule m ˜ gets excited due to interaction with light at time 0. The wave function right after the excitation | Ψ ( t = 0 ) coincides with the basis state of the electronic Hamiltonian and is a superposition of the eigenstates | Ψ ( t = 0 ) = | Φ m ˜ = c ˜ n | Ψ n . Thus, the wave function starts to coherently evolve according to the time-dependent Schrödinger equation Ψ ( t ) = n c ˜ n e i E n t | Ψ n , where this and the following expressions are in atomic units. Representing the eigenstates via the basis states, the time evolution of the wave function can be expressed
| Ψ ( t ) = n , m a n m e i E n t | Φ m ,
where a n m = c n m ˜ c n m . We obtained a n m by connecting the coefficients c ˜ n to the coefficients c n m and considering the initial condition at t = 0 . It holds n c n m c n m ˜ = δ m , m ˜ for the orthonormal basis, which results in c ˜ n = c n m ˜ .
The expression in Equation (5) describes the exciton transport right after the excitation: an exciton is localized at the site m ˜ at time t = 0 and gets delocalized over several molecules due to intersite coupling V. Exciton transport is a coherent process as long as molecular vibrations do not start destroying it. It has been heavily debated, whether coherent processes during exciton transport in molecular systems can survive for several hundred femtoseconds [3,7,8,9,10,11]. The assumption that coherence lasts for at least 20 fs is in any case reasonable [10]. We describe the pump-probe experiment for time delays up to 20 fs and neglect incoherent processes.

2.2. Interaction with an Ultrashort X-ray Probe Pulse

We assume that the molecular aggregate is probed by an ultrashort X-ray probe pulse with the Gaussian-shaped electric field
E ( r , t ) = ϵ E 0 e 2 ln 2 ( t t p ) 2 τ 2 cos ( k · r ω i n ( t t P ) )
after the time, when Frenkel exciton dynamics has been launched. Here, t p is the delay time between the pump and the probe pulse, ω i n is the central energy, k is the wave vector, ϵ is the polarization of the X-ray pulse and τ is the probe-pulse duration. The choice of the probe-pulse duration depends on the characteristic time scale of coherent electron dynamics. The time evolution of any time-dependent observable O ^ = Ψ ( t ) | O ^ | Ψ ( t ) of the Frenkel exciton system comprises harmonic oscillations with periods 2 π / ( E n E n ) , where n n . The shortest possible beating period is given by the maximum splitting of the eigenvalues of the Frenkel Hamiltonian. It is equal to 4 V for the chain of equivalent molecules. Thus, the duration of the probe pulse τ must be shorter than π / ( 2 V ) to resolve the finest details of exciton dynamics.
The total Hamiltonian in the presence of the X-ray pulse is given
H ^ = H ^ 0 d ^ E ( r , t ) ,
where d ^ is the dipole operator. Time-resolved probability of X-ray absorption from a general coherently evolving electronic system has been derived in Ref. [35] by means of the first-order time-dependent perturbation theory. Applying these results in Appendix A.1, we obtain the absorption cross section as a function of pump-probe time delay and X-ray energy for our system
σ ( t p ) = 4 π 2 ω i n F m Ω F m ( t p ) Φ F | e i k · r ϵ · d ^ | Φ m 2
with the function
Ω F m ( t p ) = π τ 2 8 ln 2 E 0 n a n m e i E n t p e ( E F E n ω i n ) 2 τ 2 8 ln 2 .
Here, the sum runs over all possible final states | Φ F with eigenenergies E F that can be reached after X-ray absorption. | Φ F are eigenstates of a general electronic Hamiltonian, but not of the Frenkel Hamiltonian. The electronic Hamiltonian can be described by the Frenkel Hamiltonian only before the X-ray absorption step. After X-ray absorption, the system is in a core-excited state, which is beyond the Frenkel exciton model.
According to Equation (8), the absorption cross section is determined by the series of peaks with time-dependent amplitudes. The bandwidth and the position of the peaks are encoded in the Gaussian functions e ( ( E F E n ω i n ) 2 τ 2 ) / ( 8 ln 2 ) . The bandwidth of the peaks is inversely proportional to the X-ray probe pulse duration τ . Here, we assumed that the pulse duration is considerably shorter than the shortest beating period of the coherent exciton dynamics. In the spectral domain, this automatically implies that the bandwidth of the probe pulse and, subsequently, the bandwidth of the spectral peaks is considerably larger than the largest energy splitting between the eigenstates of the Frenkel Hamiltonian. The difference in the spectral positions due to different energies E n cannot be distinguished and we can substitute a mean value of the Frenkel Hamiltonian eigenenergies E for any E n in Equation (9).
The time-dependence of the amplitude of the peaks are governed by the terms a n m e i E n t p , which already indicates that time-resolved cross section is sensitive to the coherent Frenkel exciton dynamics. The amplitudes are also determined by the matrix elements of the X-ray transitions Φ F | e i k · r ϵ · d ^ | Φ m . In the following, we will analyze these matrix elements by considering different possible final states after X-ray absorption | Φ F .

3. Results

In X-ray absorption spectroscopy, the energy of the X-ray pulse is tuned to the binding energy of a certain core orbital. The binding energies of core orbitals strongly vary depending on the type of the orbital and atomic species. They can be even considerably different for the same atomic species being in different chemical environments. Atomic specificity and chemical sensitivity makes this technique attractive for the study of a structure and its dynamics.
Let’s consider that an X-ray probe pulse is resonant to an orbital j of a certain atomic type. Core orbitals are strongly localized around the corresponding atom. That is why we can characterize every final many-body state F by the position of a core hole created after X-ray absorption. The electron from a localized orbital is excited to a valence orbital that is delocalized. Thus, a many-body final state can be characterized by a hole in a certain atom j located at a molecule m and a type of an excitation in the valence state. There can be different atoms j of the same type located on the same molecule. For example, the X-ray pulse may be tuned to the binding energy of the 1 s orbital of a Carbon atom and there can be several Carbon atoms in a molecule. Thus, the summation over final states will always imply the summation over orbitals j and molecules m, where the core hole can be located. In the following, we will discuss possible final states and show how X-ray absorption cross section depends on them.

3.1. X-ray Transitions below Fermi Level

Let us now assume that we tune the energy of the X-ray pulse to match the energy of an electronic transition from a core orbital to a highest occupied molecular orbital (HOMO). If the system were in the ground state, all HOMO orbitals of all molecules would be doubly occupied and such a transition would be forbidden. But since there is an exciton in the system such a transition becomes possible. The electron excited from the core orbital occupies the vacancy in the HOMO orbital. The final core-excited state has all HOMO orbitals doubly occupied, one hole in a core orbital and one electron above the Fermi level, in a LUMO orbital.
We can represent a final state in the basis of many-body core-excited states of the type
| Φ F m , j , m = | ϕ 0 1 | ϕ h j m | ϕ e m | ϕ 0 N ,
where | ϕ h j m denotes the state of a molecule m with a hole in the core orbital located at the atom j as shown in Figure 2. | ϕ e m denotes the state of the molecule m , all HOMO orbitals of which are doubly occupied and one LUMO is singly occupied. The core holes of the final states are strongly localized orbitals. But the analogously to the Frenkel exciton model, LUMO orbitals of the final states maybe delocalized due to coupling between nearest-neighbour molecules. This means that the basis states | Φ F m , j , m may not be the eigenstates of the electronic Hamiltonian.
We assume that the states with a different location of a core hole have the same energy. We also assume that the maximum possible difference between the eigenenergies of the states | Φ F m , j , n and | Φ F m , j , n is similar to the energy difference between the eigenstates of Frenkel Hamiltonian. Thus, the variation of the final-state energies E F would be negligible within the bandwidth of the probe pulse and we substitute E F for any E F . In this case, the substitution of the basis core-excited states for the final states in the expression for the time-resolved cross section in Equation (8) would not change the result (cf. Appendix A.2). We perform this substitution, since it facilitates the analysis.
The derivation of the cross section can now be reduced to the analysis of the matrix elements between the basis states of Frenkel Hamiltonian and the basis core-excited states, Φ F m , j , m | e i k · r ϵ · d ^ | Φ m . The electron in the LUMO orbital in | Φ F m , j , m and | Φ m must be on the same site before and after the X-ray transition:
| ϕ 0 1 | ϕ 0 m | ϕ m | ϕ 0 N
| ϕ 0 1 | ϕ h j m | ϕ e m | ϕ 0 N ,
which means that only matrix elements
Φ F m , j , m | e i k · r ϵ · d ^ | Φ m
are relevant for our further analysis of the cross section.
In the Hartree-Fock approximation, the transition matrix element would reduce to the integral involving the core orbital φ c ( r R m , j ) located at an atom j at the molecule m, where R m , j is the position of the orbital, and the HOMO orbital φ H ( r R m ) of the molecule m centered at R m
Φ F m , j , m | e i k · r ϵ · d ^ | Φ m = e i k · R m , j d 3 r φ c ( r R m , j ) ( ϵ · r ) φ H ( r R m ) .
Here, we use the standard assumption that the wave length of the X-ray pulse is larger than the spatial extend of the core orbital [36], which allowed us to factor out e i k · r from the integral. If the orbitals φ H ( r R m ) and φ c ( r R m , j ) do not have any spatial overlap, the integral above is zero. Thus, the most intuitive assumption to make is that an electron from a core orbital of the molecule m can be excited only to a HOMO orbital of the same molecule
Φ F m , j , m | e i k · r ϵ · d ^ | Φ m = e i k · R m , j δ m , m d j ,
where d j = d 3 r φ c ( r R m , j ) ( ϵ · r ) φ H ( r R m ) . d j does not depend on the position of the molecule. We assume that the orbitals are real functions, which means d j is also real.
This considerably simplifies the expression for the cross section, which becomes a time-independent function as we derive in Appendix A.3
σ ( t p ) = π 3 E 0 2 τ 2 2 ln 2 ω i n e E F E ω i n 2 τ 2 4 ln 2 j d j 2 .
Our assumptions about the transition matrix element led to the appearance of δ n , n in the expression for the cross section, which cancelled out any time-dependent terms. The appearance of δ n , n means that there is no final state, which can be reached from two different eigenstates of Frenkel Hamiltonian. The interference between different transition channels got forbidden and time-dependent information got lost.
Let us revisit the transition matrix element and the integral in Equation (14). We made an assumption that only an orbital on a molecule m can have a spatial overlap with a core orbital located on the same molecule m. Actually, this should not be the case for Frenkel exciton model. The main assumption of this model is that molecules are coupled via the coupling constant V. This implies that a HOMO and the lowest unoccupied molecular orbital (LUMO) of a molecule must extend to its next-neighbour molecules and we can assume
Φ F m , j , m | e i k · r ϵ · d ^ | Φ m = e i k · R m , j ( δ m , m d j + δ m , m ± 1 d ¯ j ) ,
where d ¯ j = d 3 r φ c ( r R m , j ) ( ϵ · r ) φ H ( r R m ± 1 ) , but d ¯ j is considerably smaller than d j .
The time-dependent absorption cross section then becomes (cf. Appendix A.3)
σ ( t p ) = π 3 E 0 2 τ 2 2 ln 2 ω i n e ( E F E ω i n ) 2 τ 2 4 ln 2 × j d j 2 + 2 d j d ¯ j Re m , n , n a n m + 1 a n m + a n m 1 a n m e i ( E n E n ) t p .
The first term in the squared brackets in this expression is time-independent. The second term in the squared brackets contains time-dependent terms with n n and provides time-resolved information about the dynamics. The assumption that a core electron of a molecule can be excited to the delocalized HOMO orbital of a neighbouring molecule with some probability made the interference between different transition channels possible. In the Section 3.3, we will connect properties of exciton dynamics to the time-resolved part of the cross section.

3.2. X-ray Transitions above Fermi Level

Let us now analyze the case, when the X-ray pulse is resonant to the energy of transitions to states above Fermi level, i.e., into LUMO orbitals. A final state after absorption would be a doubly excited state with two holes and two electrons in the LUMO orbitals. In a simple picture, where correlation effects do not play a considerable role, we can consider possible final states in the basis of states
| ϕ 0 1 | ϕ h j m | ϕ e m | ϕ h m | ϕ e m | ϕ 0 N ,
where one hole is in the core orbital located at atom j at molecule m, one hole is in the HOMO orbital of a molecule m , one electron is in the LUMO orbital of molecule m and one electron is in the LUMO orbital of molecule m . If we ignore correlation effects, the HOMO-LUMO excitation at the molecule m is not destroyed after the transition of a core electron into the manifold of unoccupied states (see Figure 3). This means that X-ray absorption transitions are possible only for the final states with m = m or m = m
| ϕ 0 1 | ϕ m | ϕ 0 N | Φ F m , j , m , m = | ϕ 0 1 | ϕ h j m | ϕ e m | ϕ m | ϕ 0 N .
Thus, only matrix elements Φ F m , j , m , m | e i k · r ϵ · d ^ | Φ m matter for the cross section.
With the same arguments as in the previous subsection, the transition matrix element becomes
Φ F m , j , m , m | e i k · r ϵ · d ^ | Φ m = e i k · R m , j ( δ m , m d j + δ m , m ± 1 d ¯ j ) ,
where d j and d ¯ j describe the transition matrix element from a core electron located at the atom j of the molecule m into the LUMO orbital of the same molecule m and into the LUMO orbital of the next-neighbour molecules m ± 1 , respectively.
We derive the time-resolved cross section for resonant transitions above the Fermi level In Appendix A.4 and obtain that it does not provide any time-resolved information about the dynamics
σ ( t p ) = π 3 E 0 2 τ 2 2 ln 2 ω i n e E F E ω i n 2 τ 2 4 ln 2 j N d j 2 + ( 2 N 2 ) d ¯ j 2 .
This happens even despite the fact that we allowed the transitions of core electrons to the nearest-neighbour molecules as in the previous case. The reason for this result is again the absence of interference between different transitions channels. It is not possible to reach the same many-body final state by X-ray transitions from different eigenstates involved in the dynamics and the temporal interference becomes lost.
The cross section for X-ray transitions above Fermi level is ∼N times larger than the cross section for resonant transitions below the Fermi level. This is because there are N unoccupied LUMO orbitals per basis state of the Frenkel exciton model above the Fermi level, whereas there are only one vacant HOMO orbital below the Fermi level. Thus, the time-independent contribution to the cross section would be considerably larger than the time-resolved one. Fortunately, one can exclude the transitions above the Fermi level by tuning the X-ray photon energy to be lower and match only the resonance of the transition to a HOMO orbital. One indeed needs to make sure that the probability of the transitions above the Fermi level does not start to become considerable due to a broad spectral bandwidth of the probe pulse.
We analyzed the transitions above the Fermi level ignoring correlation effects, which may start to play a role due to a doubly excited state of the chain. This may provide some correction to the expression for the cross section in Equation (22) and may even allow for some interference between transition channels. However, such an interference term, if any, would be a small correction to a large time-independent background proportional to N and a resulting time-dependent contribution would be negligible.

3.3. Connection of the Time-Resolved X-ray Absorption cross Section to Frenkel Exciton Dynamics

The density matrix of the Frenkel exciton system is given by
ρ ^ ( t ) = | Ψ ( t ) Ψ ( t ) |
by the definition. Using the expression in Equation (5), we can represent it in the site basis states
ρ ^ ( t ) = m , m ρ m , m ( t ) | Φ m Φ m |
with
ρ m , m ( t ) = n , n a n m a n m e i ( E n E n ) t .
The terms ρ m , m ( t ) are the relevant property of the Frenkel exciton dynamics, since it determines coherences between different sites. The temporal evolution of the exciton density is determined by the elements of the density matrix (exciton density is a function of space and is not to be confused with the density matrix). The exciton density consists of the spatial distribution localized around the molecules and factored by the populations ρ m , m ( t ) and the spatial distribution in the intersite region and factored by 2 Re ( ρ m , m ) . The spatial distribution in the intersite region exists, if the orbitals of sites m and m have spatial overlap.
Let us consider the observable O ^ m , m = c ^ m h ^ m c ^ m h ^ m + c . c . that describes the degree of exciton delocalization between sites m and m , if m m . Here, c ^ m ( m ) and h ^ m ( m ) are electron and hole annihilation operators at site m ( m ) that we used to define the Frenkel Hamiltonian in Section 2.1. The expectation value of O ^ m , m is given by
O ^ m , m = 2 S Re ρ m , m ( t ) ( δ m , m + 1 + δ m , m 1 ) ,
where S is the overlap integral between site basis states Φ m | c ^ m h ^ m c ^ m h ^ m | Φ m , which is nonzero only for m = m ± 1 . O ^ m , m is the measure of the spatial contribution of the intersite region between molecules m and m to the total Frenkel-exciton density. Using observable O ^ m , m , we can define the degree of exciton intersite delocalization as
h inter m , m O ^ m , m ,
which becomes
h inter ( t ) = 2 S m , n , n Re a n m + 1 a n m + a n m 1 a n m e i ( E n E n ) t .
The time-dependent part of the X-ray absorption cross section in Equation (18) follows exactly the time evolution of the intersite delocalization degree h inter ( t ) .
To illustrate our results, we model coherent Frenkel exciton dynamics and attosecond X-ray absorption spectra of a homoaggregate of eleven equidistantly aligned molecules. We select the intersite coupling V = 0.22 eV and the Coulomb energy of an electron-hole pair localized on a single site U = 2.7 eV, which are typical parameters for Frenkel exciton [37,38,39]. We assume that the molecule in the middle of the chain got excited at time zero, i.e., m ˜ = 6 .
Figure 4 illustrates the time evolution of the site populations, ρ m , m ( t ) , and of the real part of the coherences 2 Re ( ρ m , m + 1 ( t ) ) . We find that Re ( ρ m , m + 1 ( t ) ) are negligible first 5 fs, whereas Im ( ρ m , m + 1 ( t ) ) are considerable during this time. The imaginary part of density matrix elements determine the time-dependent electron current density. This means that the electron current density between sites is considerable first 5 fs and the exciton density starts to fill the intersite regions after 5 fs.
The maximum splitting between eigenenergies of Frenkel Hamiltonian is given by 4 V = 0.88 eV, which means that the shortest beating period of the coherent time evolution of the Frenkel excitons is 4.7 fs. We select the X-ray probe pulse duration τ = 500 as to be short enough to resolve the finest details of the dynamics. We assume that the X-ray probe pulse is resonant with the energy of the excitation of 1 s Carbon orbital to the HOMO orbital, which is approximately 283 eV.
We assume that j d j d ¯ j is ten times smaller as j d j 2 and calculate the cross section as a function of X-ray probe pulse central energy and pump-probe time delay. Figure 5a shows the normalized total cross section at the delay time of 500 as. The cross section is given by a Gaussian peak with FWHM of 3.8 eV centered at ∼284 eV as shown in Figure 5a. The position of the peak is constant, but the amplitude is time-dependent.
The time evolution of the change in the cross section after the delay time of 500 as, σ ˜ ( t p ) = σ ( t p ) σ ( 0.5 fs ) , as a function of energy and time is shown in Figure 5b. The change is negative for our system during the considered pump-probe time delay of 0.5–20 fs. The relative change in the cross section σ ˜ ( t p ) / σ ( t p ) depends only on time and not on energy and is about 1% of the total signal. Figure 5c shows the time evolution of h inter ( t ) obtained with Equation (28). The time evolution of the X-ray signal precisely follows the delocalization degree. Extraction of exciton delocalization is relevant for understanding exciton transport in molecular aggregates [40]. We obtained that attosecond X-ray absorption spectroscopy provides an access to this property in real time.
We can reveal even more details about Frenkel excitons from the X-ray absorption signal. Nonzero signal for Carbon K-edge reveals the p-type structure of the HOMO orbitals, where the hole of Frenkel exciton resides. Another encoded detail is about the bonds connecting the neighbouring molecules. The occurrence of the time-dependent part of the signal is governed by a nonzero transition matrix element d ¯ j . d ¯ j is nonzero, if there is a spatial overlap between orbitals of nearest-neighbour molecules. Thus, it carries information about the bonds responsible for intersite coupling.

4. Conclusions

In this article, we explored novel capabilities of attosecond X-ray pulses by analyzing how they can probe coherent Frenkel exciton dynamics in molecular aggregates. We focused on Frenkel exciton transport for up to 20 fs, when coherent processes dominate the dynamics. We investigated what attosecond X-ray absorption can reveal about coherent transport. We assumed that the X-ray probe pulse duration is short enough to resolve the shortest beating period of the coherent dynamics. In energy domain, this means that widths of the spectral peaks are considerably broader than energy splittings of electronic eigenstates involved in the dynamics.
We analyzed two different types of X-ray transitions. One case is then the X-ray pulse is resonant to the excitation of an electron from a core orbital to the orbitals below Fermi level, where a hole of the exciton resides. The second case is the transitions to the orbitals above Fermi level, where the electron of the exciton resides. In both cases the spectrum would be formed by a broad spectral peak centered at the central energy of the X-ray pulse. But since the X-ray pulse would be tuned to two different energies, the both peaks would be spectroscopically distinguishable.
Interference effects are necessary for a signal to be time-dependent. The interference of the signal is created, when the same final core-excited state can be reached from several eigenstates involved in the dynamics. We have shown that it is not possible to reach the same final core-excited state from two different eigenstates of Frenkel Hamiltonian for resonant transitions to states above Fermi level. Thus, the part of the X-ray spectrum formed by the peak corresponding to transitions above Fermi level would be simply time-independent and would not carry any information about the coherent dynamics.
The part of the signal corresponding to transitions below Fermi level does carry time-resolved information. The time evolution of the peak amplitude precisely follows the time evolution of the degree of exciton delocalization. The degree of exciton delocalization gives the information about how much of the exciton density is distributed in the intersite region.
Exciton delocalization is at the core of electronic energy transfer of molecular excitons [40]. Electronic energy transfer is the dynamical regime in which pure coherent evolution competes with dissipative evolution in contrast to transfer by incoherent hopping jumps in Förster regime. Electronic energy transfer accompanies photosynthetic processes that result in highly efficient light-energy conversion [41]. Understanding the principles of highly efficient light harvesting in photosynthesis is crucial for design of efficient organic solar cells that also involve molecular excitons for energy transfer. Measuring exciton delocalization will provide the insight into site superposition of excitations present during electronic energy transfer of molecular excitons. Access to this property on attosecond time scale will additionally reveal the role of pure electronic degrees of freedom for electronic energy transfer.
An optical pulse can also be used as an attosecond probe of electron dynamics. But there is one key difference of an optical probe to an X-ray probe. Optical pulses induce electronic transitions from delocalized orbitals to delocalized orbitals, and X-ray pulses induce transitions of electrons from localized orbitals to delocalized orbitals. This specifics of X-ray transitions obviously leads to atomic sensitivity. But other essential advantage for probing exciton transport is that transitions involving localized core states disentangle the degree of delocalization in valence states and map it into the time-resolved part of the X-ray signal.

Author Contributions

Conceptualization, D.P.-G.; methodology, T.H. and D.P.-G.; software, T.H.; validation, T.B. and D.P.-G.; formal analysis, T.H.; investigation, T.H.; writing—original draft preparation, T.H. and D.P.-G.; writing—review and editing, T.B.; visualization, T.H. and T.B.; supervision, T.B. and D.P-G.; project administration, D.P-G. All authors have read and agreed to the published version of the manuscript.

Funding

We acknowledge the funding from the Volkswagen Foundation: 96237 and the funding from the Cluster of Excellence ‘CUI: Advanced Imaging of Matter’ of the Deutsche Forschungsgemeinschaft (DFG)-EXC 2056-project ID 390715994. We acknowledge financial support from the Open Access Publication Fund of Universität Hamburg.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations and notations are used in this manuscript:
HOMOHighest occupied molecular orbital
LUMOLowest unoccupied molecular orbital
E 0 Amplitude of the electric field of a probe pulse
t p Pump-probe time delay
ω i n Central frequency of a probe pulse
k Wave vector of a probe pulse
ϵ Polarization vector of a probe pulse
τ Probe pulse duration
m ˜ Index for the initially excited molecule
jIndex for the atom j
mIndex for the molecule m
UCoulomb energy of an electron-hole pair localized on a single site
VCoupling between nearest-neighbor molecules
h i n t e r Degree of the exciton intersite delocalization
c ˜ n Expansion coefficients of the representation of the time-dependent wave function as eigenstates
c n m Expansion coefficients of the representation of the eigenstates of Frenkel Hamiltonian as basis states
a n m Expansion coefficients of the solution of the time-dependent Schrödinger equation
d j Transition matrix element between core electron orbital of the atom j of molecule m and
HOMO of the same molecule
d ¯ j Transition matrix element between core electron orbital of the atom j of the molecule m and
HOMO of the next-neighbor molecules
d j Transition matrix element between core electron orbital of the atom j of the molecule m and
LUMO of the same molecule
d ¯ j Transition matrix element between core electron of the atom j of the molecule m and
LUMO of the next-neighbor molecules
| ϕ 0 m Ground state of the molecule m
| Ψ 0 Ground state of a chain of molecules
| ϕ m State which corresponds to an electron-hole pair located at the site m
| Φ m Basis state of Frenkel Hamiltonian
| Ψ n Eigenstate of Frenkel Hamiltonian
| ϕ h j m State of the molecule m with core hole at atom j
| ϕ e m State of the molecule m with one electron in LUMO
| ϕ h m State of the molecule m with hole in HOMO
| Φ F Final state after absorption
| Φ F m , j , m Final state with a core hole at the atom j of the molecule m and an electron in LUMO of the molecule m
| Φ F m , j , m , m Final state with a core hole at the atom j of the molecule m, an electron-hole pair at the site m and electron in LUMO at the site m

Appendix A

Appendix A.1. Derivation of Time-Resolved X-ray Absorption cross Section

The absorption cross section within the first-order time-dependent perturbation theory is given by [35]
σ ( t p ) = 4 π 2 ω i n lim t 0 F t 0 + d t Φ F | e i H ^ 0 t E ( r , t t p ) d ^ e i H ^ 0 t | Ψ ( t 0 ) 2
Here, the sum runs over all possible final states | Φ F with eigenenergies E F that can be reached after X-ray absorption. Substituting the time-dependent wave function represented in the basis of the Frenkel basis states in Equation (5), we obtain
σ ( t p ) = 4 π 2 ω i n E 0 2 F + d t Φ F | e i E F t e i k · r ( ϵ · d ^ ) e 2 ln 2 ( t t p ) 2 τ 2 e i ω i n ( t t p ) 2 n , m a n m e i E n t | Φ m 2 .
Performing the integration over time, we obtain the expression for the time-resolved absorption cross section in Equation (8) in Section 2.2.

Appendix A.2. Representation of the X-ray Absorption cross Section via Basis States

Let us represent the final states in Equation (8) via basis states | Φ F m , j , m described in Section 3.1: | Φ F m , j , n = m b n , m | Φ F m , j , m . Since a core orbital is strongly localized, each final eigenstate can be characterized by the position of the core hole in the atom j of a molecule m. The index n refers to a delocalized valence orbital. We assume that the final eigenstates | Φ F m , j , n and | Φ F m , j , n with a different location of a core hole have the same energy. We also assume that the maximum possible difference between the eigenenergies of the states | Φ F m , j , n and | Φ F m , j , n is similar to the energy difference between the eigenstates of Frenkel Hamiltonian. Thus, the variation of the final-state energies E F would be negligible within the bandwidth of the probe pulse and we substitute E F for any E F . We obtain that the expression for the cross section does not depend on the representation of the final states:
σ ( t p ) = 4 π 2 ω i n m , j , n m b n , m Φ F m , j , m | e i k · r ϵ · d ^ | Φ ˜ 2 = 4 π 2 ω i n m , j , n m , m ˜ b n , m b n , m ˜ Φ F m , j , m | e i k · r ϵ · d ^ | Φ ˜ Φ ˜ | e i k · r ϵ · d ^ | Φ F m , j , m ˜ = 4 π 2 ω i n m , j , m Φ F m , j , m | e i k · r ϵ · d ^ | Φ ˜ 2
where we used the notation m Ω F m ( t p ) | Φ m = | Φ ˜ and that n b n , m b n , m ˜ = δ m , m ˜ , since they are coefficients of an orthonormal basis.

Appendix A.3. Derivation of X-ray Absorption cross Section for Transitions below Fermi Level

Let us consider the time-resolved absorption cross section with the two approximations considered in Section 3.1: X-ray energy is resonant to the transition to a HOMO orbital and the HOMO orbital of a molecule does not have any spatial overlap with a core orbital of an atom located at a different molecule. We substitute the transition matrix element in Equation (15) to the expression for the absorption cross section in Equation (8) and obtain
σ ( t p ) = 4 π 2 ω i n m , j e i k · R m , j m Ω F m ( t p ) d j δ m , m 2 = 4 π 2 ω i n j d j 2 m Ω F m ( t p ) 2 .
We substitute the definition for Ω F m and obtain
σ ( t p ) = π 3 E 0 2 τ 2 2 ln 2 ω i n e E F E ω i n 2 τ 2 4 ln 2 j d j 2 n , n e i ( E n E n ) t p m a n m a n m .
Let us look into the sum m a n m a n m in more details by revisiting the definition of a n m as c n m ˜ c ˜ n m and that c n m are the coefficients of an orthonormal basis set.
m a n m a n m = c n m ˜ c n m ˜ m c n m c n m = c n m ˜ c n m ˜ δ n , n .
Thus, the absorption cross section becomes
σ ( t p ) = π 3 E 0 2 τ 2 2 ln 2 ω i n e E F E ω i n 2 τ 2 4 ln 2 j d j 2 n , n δ n , n | c n m ˜ | 2 e i ( E n E n ) t p .
We note that n | c n m ˜ | 2 = 1 . This considerably simplifies the expression for the cross section, which becomes simply a time-independent function in Equation (16).
We now substitute the transition matrix element in Equation (17) and obtain the new expression for the absorption cross section
σ ( t p ) = 4 π 2 ω i n m , j d j Ω F m ( t p ) + d ¯ j Ω F m ± 1 ( t p ) 2 = 4 π 2 ω i n m , j d j 2 | Ω F m ( t p ) | 2 + 2 d j d ¯ j Re Ω F m ( t p ) Ω F m + 1 ( t p ) + Ω F m ( t p ) Ω F m 1 ( t p ) .
Here, we neglected the term proportional to d ¯ j 2 since it is much smaller than the other terms. From the analysis above, we already obtained that the first term proportional to d j 2 provides a time-independent contribution. Substituting the definition of Ω F m , j , e , we obtain that the time-dependent cross section is given by
σ ( t p ) = π 3 E 0 2 τ 2 2 ln 2 ω i n e ( E F E ω i n ) 2 τ 2 4 ln 2 × j d j 2 + 2 d j d ¯ j Re m , n , n a n m + 1 a n m + a n m 1 a n m e i ( E n E n ) t p .

Appendix A.4. Derivation of X-ray Absorption cross Section for Transitions above Fermi Level

The transition matrix elements in the case, then the core hole and the valence hole are located at the same molecule ( m = m = m ) may deviate from d j due to correlation effects within the same molecule. However, this will not influence the final conclusion and we do not distinguish such transition matrix elements for simplicity.
The matrix elements Φ F m , j , m ˜ , m | e i k · r ϵ · d ^ | Φ m are nonzero for m = m ˜ as discussed in Section 3.2. This means that a Kronecker delta δ m , m ˜ appears in the expression for the time-resolved cross section in Equation (8) and the summation over Frenkel basis states inside the square of absolute values reduces to one term. Substituting the expression for the transition matrix elements in Equation (21) to Equation (8), we obtain
σ ( t p ) = π 3 E 0 2 τ 2 2 ln 2 ω i n m , j , m , m e E F E ω i n 2 τ 2 4 ln 2 δ m , m d j + δ m , m ± 1 d ¯ j n a n m e i E n t p 2 = π 3 E 0 2 τ 2 2 ln 2 ω i n e E F E ω i n 2 τ 2 4 ln 2 m , j , m , m d j 2 δ m , m n a n m e i E n t p 2 + 2 m , j , m , m δ m , m δ m , m ± 1 d j d ¯ j Re n a n m e i E n t p n a n m e i E n t p + m , j , m , m d ¯ j 2 δ m , m ± 1 n a n m e i E n t p 2
The term proportional to d j d ¯ j becomes zero, since δ m , m δ m , m ± 1 = 0 and we obtain
σ ( t p ) = π 3 E 0 2 τ 2 2 ln 2 ω i n e E F E ω i n 2 τ 2 4 ln 2 m , j , m d j 2 δ m , m + m , j , m d j ¯ 2 δ m , m ± 1 m n a n m e i E n t p 2
Applying Equation (A6), we obtain m n a n m e i E n t p 2 = n , n | c ˜ n m ˜ | 2 δ n , n = 1 . The summation over m , j , m δ m , m gives the number of the molecules in the chain and the summation m , j , m ( δ m , m + 1 + δ m , m 1 ) gives twice the number of nearest-neighbour molecules.
It would be more precise to distinguish the case, then m = m = m , i.e., the core excitation and the valence excitation are at the same molecular site in the final state
Φ F m , j , m , m | e i k · r ϵ · d ^ | Φ m = e i k · R m , j ( δ m , m d j + δ m , m ± 1 d ¯ j + δ m , m δ m , m d ˜ j ) .
d ˜ j would be a correction for electronic correlations at the same molecular site. However, this would not change the conclusion that the absorption cross section would remain time-independent and we neglect the correction for simplicity.

References

  1. Cheng, Y.C.; Fleming, G.R. Dynamics of Light Harvesting in Photosynthesis. Annu. Rev. Phys. Chem. 2009, 60, 241–262. [Google Scholar] [CrossRef] [PubMed]
  2. Dimitriev, O.P. Dynamics of Excitons in Conjugated Molecules and Organic Semiconductor Systems. Chem. Rev. 2022, 122, 8487–8593. [Google Scholar] [CrossRef] [PubMed]
  3. Duan, H.G.; Prokhorenko, V.I.; Cogdell, R.J.; Ashraf, K.; Stevens, A.L.; Thorwart, M.; Miller, D. Nature does not rely on long-lived electronic quantum coherence for photosynthetic energy transfer. Proc. Natl. Acad. Sci. USA 2017, 114, 8493–8498. [Google Scholar] [CrossRef] [PubMed]
  4. Cui, X.Y.; Yan, Y.J.; Wei, J.H. Theoretical study on the exciton dynamics of coherent excitation energy transfer in the phycoerythrin 545 light-harvesting complex. Chin. Phys. B 2021, 31, 018201. [Google Scholar] [CrossRef]
  5. Wong, C.Y.; Alvey, R.M.; Turner, D.B.; Wilk, K.E.; Bryant, D.A.; Curmi, P.M.; Silbey, R.J.; Scholes, G.D. Electronic coherence lineshapes reveal hidden excitonic correlations in photosynthetic light harvesting. Nat. Chem. 2012, 4, 396–404. [Google Scholar] [CrossRef]
  6. Wang, L.; Allodi, M.A.; Engel, G.S. Quantum coherences reveal excited-state dynamics in biophysical systems. Nat. Rev. Chem. 2019, 3, 477–490. [Google Scholar] [CrossRef]
  7. Engel, G.S.; Calhoun, T.R.; Read, E.L.; Ahn, T.K.; Mančal, T.; Cheng, Y.C.; Blankenship, R.E.; Fleming, G.R. Evidence for wavelike energy transfer through quantum coherence in photosynthetic systems. Nature 2007, 446, 782–786. [Google Scholar] [CrossRef]
  8. Lambert, N.; Chen, Y.N.; Cheng, Y.C.; Li, C.M.; Chen, G.Y.; Nori, F. Quantum biology. Nat. Phys. 2013, 9, 10–18. [Google Scholar] [CrossRef]
  9. Fassioli, F.; Dinshaw, R.; Arpin, P.C.; Scholes, G.D. Photosynthetic light harvesting: Excitons and coherence. J. R. Soc. Interface 2014, 11, 20130901. [Google Scholar] [CrossRef]
  10. Cassette, E.; Pensack, R.D.; Mahler, B.; Scholes, G.D. Room-temperature exciton coherence and dephasing in two-dimensional nanostructures. Nat. Commun. 2015, 6, 6086. [Google Scholar] [CrossRef]
  11. Cao, J.; Cogdell, R.J.; Coker, D.F.; Duan, H.G.; Hauer, J.; Kleinekathöfer, U.; Jansen, T.L.C.; Mancal, T.; Miller, R.J.D.; Ogilvie, J.P.; et al. Quantum biology revisited. Sci. Adv. 2020, 6, eaaz4888. [Google Scholar] [CrossRef]
  12. Fassioli, F.; Olaya-Castro, A.; Scholes, G.D. Coherent Energy Transfer under Incoherent Light Conditions. J. Phys. Chem. Lett. 2012, 3, 3136–3142. [Google Scholar] [CrossRef] [PubMed]
  13. Rossi, G.M.; Mainz, R.E.; Yang, Y.; Scheiba, F.; Silva-Toledo, M.A.; Chia, S.H.; Keathley, P.D.; Fang, S.; Mücke, O.D.; Manzoni, C.; et al. Sub-cycle millijoule-level parametric waveform synthesizer for attosecond science. Nat. Photonics 2020, 14, 629–635. [Google Scholar] [CrossRef]
  14. Duris, J.; Li, S.; Driver, T.; Champenois, E.G.; MacArthur, J.P.; Lutman, A.A.; Zhang, Z.; Rosenberger, P.; Aldrich, J.W.; Coffee, R.; et al. Tunable isolated attosecond X-ray pulses with gigawatt peak power from a free-electron laser. Nat. Photonics 2020, 14, 30–36. [Google Scholar] [CrossRef]
  15. Buades, B.; Picón, A.; Berger, E.; León, I.; Di Palo, N.; Cousin, S.L.; Cocchi, C.; Pellegrin, E.; Martin, J.H.; Manas-Valero, S.; et al. Attosecond state-resolved carrier motion in quantum materials probed by soft X-ray XANES. Appl. Phys. Rev. 2021, 8, 011408. [Google Scholar] [CrossRef]
  16. Garratt, D.; Misiekis, L.; Wood, D.; Larsen, E.W.; Matthews, M.; Alexander, O.; Ye, P.; Jarosch, S.; Ferchaud, C.; Strüber, C.; et al. Direct observation of ultrafast exciton localization in an organic semiconductor with soft X-ray transient absorption spectroscopy. Nat. Commun. 2022, 13, 3414. [Google Scholar] [CrossRef] [PubMed]
  17. Zinchenko, K.S.; Ardana-Lamas, F.; Seidu, I.; Neville, S.P.; van der Veen, J.; Lanfaloni, V.U.; Schuurman, M.S.; Wörner, H.J. Sub-7-femtosecond conical-intersection dynamics probed at the carbon K-edge. Science 2021, 371, 489–494. [Google Scholar] [CrossRef] [PubMed]
  18. Géneaux, R.; Kaplan, C.J.; Yue, L.; Ross, A.D.; Bækhøj, J.E.; Kraus, P.M.; Chang, H.T.; Guggenmos, A.; Huang, M.Y.; Zürch, M.; et al. Attosecond Time-Domain Measurement of Core-Level-Exciton Decay in Magnesium Oxide. Phys. Rev. Lett. 2020, 124, 207401. [Google Scholar] [CrossRef]
  19. Geneaux, R.; Marroux, H.J.B.; Guggenmos, A.; Neumark, D.M.; Leone, S.R. Transient absorption spectroscopy using high harmonic generation: A review of ultrafast X-ray dynamics in molecules and solids. Philos. Trans. A Math. Phys. Eng. Sci. 2019, 377, 20170463. [Google Scholar] [CrossRef]
  20. Krausz, F.; Ivanov, M. Attosecond physics. Rev. Mod. Phys. 2009, 81, 163–234. [Google Scholar] [CrossRef]
  21. Popova-Gorelova, D. Imaging Electron Dynamics with Ultrashort Light Pulses: A Theory Perspective. Appl. Sci. 2018, 8, 318. [Google Scholar] [CrossRef]
  22. Saes, M.; Bressler, C.; Abela, R.; Grolimund, D.; Johnson, S.L.; Heimann, P.A.; Chergui, M. Observing Photochemical Transients by Ultrafast X-ray Absorption Spectroscopy. Phys. Rev. Lett. 2003, 90, 047403. [Google Scholar] [CrossRef] [PubMed]
  23. Bhattacherjee, A.; Leone, S.R. Ultrafast X-ray Transient Absorption Spectroscopy of Gas-Phase Photochemical Reactions: A New Universal Probe of Photoinduced Molecular Dynamics. Accounts Chem. Res. 2018, 51, 3203–3211. [Google Scholar] [CrossRef] [PubMed]
  24. Khakhulin, D.; Otte, F.; Biednov, M.; Bömer, C.; Choi, T.K.; Diez, M.; Galler, A.; Jiang, Y.; Kubicek, K.; Lima, F.A.; et al. Ultrafast X-ray Photochemistry at European XFEL: Capabilities of the Femtosecond X-ray Experiments (FXE) Instrument. Appl. Sci. 2020, 10, 995. [Google Scholar] [CrossRef]
  25. Shakya, Y.; Inhester, L.; Arnold, C.; Welsch, R.; Santra, R. Ultrafast time-resolved X-ray absorption spectroscopy of ionized urea and its dimer through ab initio nonadiabatic dynamics. Struct. Dyn. 2021, 8, 034102. [Google Scholar] [CrossRef]
  26. Frenkel, J. On the Transformation of light into Heat in Solids. I. Phys. Rev. 1931, 37, 17–44. [Google Scholar] [CrossRef]
  27. Saikin, S.K.; Eisfeld, A.; Valleau, S.; Aspuru-Guzik, A. Photonics meets excitonics: Natural and artificial molecular aggregates. Nanophotonics 2013, 2, 21–38. [Google Scholar] [CrossRef]
  28. Bardeen, C.J. The Structure and Dynamics of Molecular Excitons. Annu. Rev. Phys. Chem. 2014, 65, 127–148. [Google Scholar] [CrossRef]
  29. Agranovich, V. Excitations in Organic Solids; Oxford University Press: Oxford, UK, 2008. [Google Scholar]
  30. Akselrod, G.M.; Deotare, P.B.; Thompson, N.J.; Lee, J.; Tisdale, W.A.; Baldo, M.A.; Menon, V.M.; Bulović, V. Visualization of exciton transport in ordered and disordered molecular solids. Nat. Commun. 2014, 5, 3646. [Google Scholar] [CrossRef]
  31. Wan, Y.; Guo, Z.; Zhu, T.; Yan, S.; Johnson, J.; Huang, L. Cooperative singlet and triplet exciton transport in tetracene crystals visualized by ultrafast microscopy. Nat. Chem. 2015, 7, 785–792. [Google Scholar] [CrossRef]
  32. Zhu, T.; Wan, Y.; Huang, L. Direct Imaging of Frenkel Exciton Transport by Ultrafast Microscopy. Accounts Chem. Res. 2017, 50, 1725–1733. [Google Scholar] [CrossRef] [PubMed]
  33. Stadtmüller, B.; Emmerich, S.; Jungkenn, D.; Haag, N.; Rollinger, M.; Eich, S.; Maniraj, M.; Aeschlimann, M.; Cinchetti, M.; Mathias, S. Strong modification of the transport level alignment in organic materials after optical excitation. Nat. Commun. 2019, 10, 1470. [Google Scholar] [CrossRef] [PubMed]
  34. Emmerich, S.; Hedwig, S.; Arnoldi, B.; Stöckl, J.; Haag, F.; Hemm, R.; Cinchetti, M.; Mathias, S.; Stadtmüller, B.; Aeschlimann, M. Ultrafast Charge-Transfer Exciton Dynamics in C60 Thin Films. J. Phys. Chem. C 2020, 124, 23579–23587. [Google Scholar] [CrossRef] [PubMed]
  35. Santra, R.; Yakovlev, V.S.; Pfeifer, T.; Loh, Z.H. Theory of attosecond transient absorption spectroscopy of strong-field-generated ions. Phys. Rev. A 2011, 83, 033405. [Google Scholar] [CrossRef]
  36. Santra, R. Concepts in X-ray physics. J. Phys. B At. Mol. Opt. Phys. 2008, 42, 023001. [Google Scholar] [CrossRef]
  37. Albert, V.V.; Badaeva, E.; Kilina, S.; Sykora, M.; Tretiak, S. The Frenkel exciton Hamiltonian for functionalized Ru(II)–bpy complexes. J. Lumin. 2011, 131, 1739–1746. [Google Scholar] [CrossRef]
  38. Farahvash, A.; Lee, C.K.; Sun, Q.; Shi, L.; Willard, A.P. Machine learning Frenkel Hamiltonian parameters to accelerate simulations of exciton dynamics. J. Chem. Phys. 2020, 153, 074111. [Google Scholar] [CrossRef]
  39. Nematiaram, T.; Padula, D.; Troisi, A. Bright Frenkel Excitons in Molecular Crystals: A Survey. Chem. Mater. 2021, 33, 3368–3378. [Google Scholar] [CrossRef]
  40. Scholes, G.D.; Smyth, C. Perspective: Detecting and measuring exciton delocalization in photosynthetic light harvesting. J. Chem. Phys. 2014, 140, 110901. [Google Scholar] [CrossRef]
  41. Chenu, A.; Scholes, G.D. Coherence in Energy Transfer and Photosynthesis. Annu. Rev. Phys. Chem. 2015, 66, 69–96. [Google Scholar] [CrossRef]
Figure 1. A many-body basis state | Φ m of Frenkel Hamiltonian.
Figure 1. A many-body basis state | Φ m of Frenkel Hamiltonian.
Molecules 28 04502 g001
Figure 2. A many-body basis state | Φ F m , j , m , in which a core electron occupied a vacant HOMO orbital after absorption of an X-ray photon.
Figure 2. A many-body basis state | Φ F m , j , m , in which a core electron occupied a vacant HOMO orbital after absorption of an X-ray photon.
Molecules 28 04502 g002
Figure 3. A many-body basis state | Φ F m , j , m , m , in which a core electron occupied a vacant LUMO orbital after absorption of an X-ray photon.
Figure 3. A many-body basis state | Φ F m , j , m , m , in which a core electron occupied a vacant LUMO orbital after absorption of an X-ray photon.
Molecules 28 04502 g003
Figure 4. Populations ρ m , m for all molecules in the chain and the real part of the coherences 2 Re ( ρ m , m + 1 ) for all nearest neighbours depending on time.
Figure 4. Populations ρ m , m for all molecules in the chain and the real part of the coherences 2 Re ( ρ m , m + 1 ) for all nearest neighbours depending on time.
Molecules 28 04502 g004
Figure 5. (a) Normalized cross section at the time delay of 500 as as a function of energy. X-ray probe pulse is assumed to be resonant with the transition energy of an electron from the 1 s orbital of Carbon into outermost orbitals below Fermi level. (b) Normalized time evolution of the change in the cross section after the delay time of 500 as, σ ˜ ( t p ) = σ ( t p ) σ ( 0.5 fs ) , as a function of energy and time. (c) Time evolution evolution of the normalized delocalization degree h inter ( t ) .
Figure 5. (a) Normalized cross section at the time delay of 500 as as a function of energy. X-ray probe pulse is assumed to be resonant with the transition energy of an electron from the 1 s orbital of Carbon into outermost orbitals below Fermi level. (b) Normalized time evolution of the change in the cross section after the delay time of 500 as, σ ˜ ( t p ) = σ ( t p ) σ ( 0.5 fs ) , as a function of energy and time. (c) Time evolution evolution of the normalized delocalization degree h inter ( t ) .
Molecules 28 04502 g005
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

Hansen, T.; Bezriadina, T.; Popova-Gorelova, D. Theoretical Description of Attosecond X-ray Absorption Spectroscopy of Frenkel Exciton Dynamics. Molecules 2023, 28, 4502. https://doi.org/10.3390/molecules28114502

AMA Style

Hansen T, Bezriadina T, Popova-Gorelova D. Theoretical Description of Attosecond X-ray Absorption Spectroscopy of Frenkel Exciton Dynamics. Molecules. 2023; 28(11):4502. https://doi.org/10.3390/molecules28114502

Chicago/Turabian Style

Hansen, Tim, Tatiana Bezriadina, and Daria Popova-Gorelova. 2023. "Theoretical Description of Attosecond X-ray Absorption Spectroscopy of Frenkel Exciton Dynamics" Molecules 28, no. 11: 4502. https://doi.org/10.3390/molecules28114502

APA Style

Hansen, T., Bezriadina, T., & Popova-Gorelova, D. (2023). Theoretical Description of Attosecond X-ray Absorption Spectroscopy of Frenkel Exciton Dynamics. Molecules, 28(11), 4502. https://doi.org/10.3390/molecules28114502

Article Metrics

Back to TopTop