Abstract
We theoretically propose a biphoton entanglement-enhanced multidimensional spectroscopic technique as a probe for the dissipative polariton dynamics in the ultrafast regime. It is applied to the cavity-confined monomeric photosynthetic complex that represents a prototypical multi-site excitonic quantum aggregate. The proposed technique is shown to be particularly sensitive to inter-manifold polariton coherence between the two and one-excitation subspaces. It is demonstrated to be able to monitor the dynamical role of cavity-mediated excitonic correlations, and dephasing in the presence of phonon-induced dissipation. The non-classicality of the entangled biphoton sources is shown to enhance the ultra-fast and broadband correlation features of the signal, giving an indication about the underlying state correlations responsible for long-range cavity-assisted exciton migration.
1 Introduction
The quantum aggregates consisting of multiple centers of electronic excitations, e.g., the naturally occurring light-harvesting photosynthetic complexes, artificially designed molecular light-harvesters offer a uniquely favorable testing ground for the entangled photon-induced dynamics. These systems intrinsically host collective excitations, frequently termed as molecular excitons, extending over several excitation centers/sites and often within the coherence domain of the spatial (temporal) length (time) scales. Thus, it provides opportunities for the external modulation of delocalized excitons and testing the limits of coherent dynamics. Further, the vibrational motions associated with these structures often give rise to a collective dephasing mechanism for the delocalized excitons. Additionally, the vibrational motions modulate the energy gradient of the kinetics, leading to situations where coherent excitons undergo inter-site transport and eventually localize. The spectroscopic investigation and control of such exciton kinetics offers insights into the microscopic nature of the coherent energy transfer mechanism and prescribes guidelines for bio-mimetic engineering which builds on the operational equivalence [1–6]. In a separate line of development, there have been a series of studies that have demonstrated the effective role of external electromagnetic cavities in manipulating material excitations. These studies range from the control of electronic excitations [7–16], vibrational modes [17, 18], collective mode responses [19, 20] to the cavity mode-assisted modulation of dynamical resonances [21–25]. Cavity mode interacting resonantly with narrow-band of excitons and off-resonantly with the rest may modulate the excitation dispersal by tuning the delocalization properties and influence the dephasing properties via the spectral weight modulation [26–28]. Hence, combining these two developments to the case of light-harvesting quantum aggregates offers an opportunity to investigate the role of the cavity in controlling a prototypical, extended yet aperiodic system that hosts collective excitons.
Previous spectroscopic studies of the cavity-modulated dynamics [29–35] have increasingly focused their attention on the ultrafast, nonlinear techniques analogs e.g. pump-probe here. In the majority of these studies, participation of the higher-order excitations remain relatively unexplored. The dynamics involving the higher-order nonlinear excitations introduce possibilities for cavity-modulated exciton-exciton annihilation, cavity-assisted exciton fusion, and many more correlated mechanisms of practical interest. On a fundamental level, the external dielectric modification of the Coulomb interaction between different excitons underlies all of these coherent mechanisms. However, the nonlinear spectroscopic signatures studied so far depend dominantly on the associated vibrational (or vibronic) processes [37, 38]. Therefore, in order to explore the complexity of cavity-modulated nonlinear exciton dynamics, it is desirable to combine ultrafast spectroscopic tools that are sensitive to high-order exciton correlations in the presence of collective vibrational dephasing. Among the multi-pulse nonlinear spectroscopic techniques available for mapping out the correlation between two excitons and probing the two-exciton state-specific dephasing in the ultrafast regime, double-quantum multidimensional correlation spectra (DQC) have been proven useful [39–42]. The ultrafast nature of the exciton dynamics in aggregates results in difficulty in measuring the role of energetically-distant states even without the presence of the cavity mode. These delocalized multi-exciton states often participate in dynamics within the same timescale due to dynamical dephasing properties. Associated with the normal exciton-number conserving dephasing process, one encounters a dephasing mechanism that couples the energy-manifolds with the different number of excitons. The presence of cavity adds additional cavity-exciton hybrid states, introduces multiple dynamical timescales, and gives rise to novel cavity-exciton-phonon coupling mechanisms. These excitations are of principal interest in this communication. The investigation of them requires the deployment of probes that specifically excite spectrally narrow-band states (creation of non-linear polarization), allow the associated dynamics to evolve in time in the presence of dephasing (evolution of the polarization), and project them to desired frequency components of interest (projection of polarization to signal components). Robust, correlated state excitation can be controllably achieved by deploying the entangled photon pairs, namely entangled biphotons, via a scheme being referred to as entangled two-photon absorption [43–47]. In comparison to the shaped laser pulses, the entangled biphotons have been shown to improve the spectral resolution while scaling favorably with the intensity of the sources. The non-classical correlation properties can also be utilized to obtain favorable spectral-temporal resolution in the probing via selective state projections of the nonlinear polarization [29, 33, 48–55]. The latter requires the selection of a few states whose correlation properties are of particular interest from a manifold. In this regard, a combination of the aforementioned DQC signal measurement scheme with the entangled photon sources may provide a technique that studies cavity-modulated correlated exciton kinetics in the ultrafast regime with higher spectral resolution.
In what follows, in Section 2, we introduce the Frenkel exciton Hamiltonian and describe the underlying model. It will be used to describe the dissipative exciton-polariton phenomenology and obtain the relevant Green’s functions. Subsequently, in Section 2 we introduce the DQC signal in a modular manner, discuss the nature of the signal and present the simulation results in the relevant parameter regime. Section 4 discusses limitations, the scope of the presented signal within the broader scope of entangled photon-enhanced spectroscopies, and the outlook.
2 Dissipative Exciton-Polariton Phenomenology
Here we present the exciton Hamiltonian concurrently interacting with a cavity mode and the phonon reservoir, which will be used to build up the phenomenology using a quasi-particle picture. It is given by,where we set the ℏ to unity. The components of the Hamiltonian are explained below.
2.1 Exciton, Cavity and Exciton-Cavity Interactions
The first three terms constitute the bare exciton Hamiltonian where are the m-th site exciton creation (annihilation) operators with the respective commutation relation, and ) (where κm = μ21/μ10, is given in terms of ratio of transition dipoles of ij-th level of the m-th site) [56–58]. is the on-site excitation energies and Jmn is the inter-site Coulomb-mediated hopping. The presence of the interacting terms in the exciton Hamiltonian indicates that the single and two exciton states are delocalized over sites. The two-excitation manifold which is composed of permutable composition of pure local and non-local two-exciton, pure two cavity excitations and joint one exciton-one photon excitation is given by, . The terms further dictates the interaction among higher-order excitonic states which gives rise to static, local energy shifts, exciton-exciton scattering. The values of is taken as , with κ = μ21/μ10 which accounts for the local energy shifts for two-exciton states and . Bare two-exciton energy is given by (where the Δm is the two-exciton state anharmonicity from the Harmonic limit). The parameter values for the Hamiltonian are obtained from the semi-empirical simulation in [59–62] where these data-set have been shown to simultaneously fit the spectral data obtained from the linear absorption, fluorescence, and circular dichroism spectra. The combined cavity and the cavity-exciton interaction Hamiltonian is be given by the fourth and fifth term where the gc;m,α is the cavity-exciton coupling strengths describing the resonant dipolar interaction. The effects arising from the non-resonant cavity-exciton interaction could be significant in the static limit if the exciton-cavity coupling strengths are comparable to the bare Rabi frequencies of the multi-exciton system [63–66]. The cavity-exciton parameter considered in this communication corresponds to moderately weak and within the close to resonant regime which allows this term being neglected. The cavity photonic modes are described in the oscillator basis with the corresponding creation (annihilation) operators denoted as , for the α-th mode. Within the scope of this study, a single-mode limit (i.e. α = 1) and two-excitation per-mode (with he fundamental frequency denoted by ωc) has been considered. Furthermore, we assume an uniform cavity-exciton coupling strengths for all sites, i.e. gc;m,α = gc = 100cm−1 is independent of index m, although the analysis presented is not limited by that choice. We note that a detailed estimation of the cavity-exciton couplings requires investigation regarding the mode volume of the cavity and the quality factor which is beyond the scope of the article. The Coulombic origin inter-site hopping in the Hamiltonian is of different magnitude, indicating a different propensity for delocalization in the absence of cavity interaction. The cavity coupling modulates site-delocalization at different extents for one-exciton, local and non-local two-excitons. This gives rise to the possibility of creating a different admixture of cavity-matter excitations. In addition, it offers cavity mediates coupling between these configurations, thereby coupling diverse sets of excitations and inter-site processes. Since the cavity accommodates two-photon excitations, the analysis is capable of describing resonant two-polariton processes.
2.2 Phonon and Exciton-Phonon Interactions
The exciton-phonon interactions originate from the inter and intra-molecular vibrational motions associated with the relative nuclear motions of the aggregate. Normal modes of the low-energy vibrational degrees of freedom related to the collective vibrational coordinates are assigned as phonon modes and are mapped onto an infinite set of Harmonic oscillators. It is given by the free phonon Hamiltonian appearing in the sixth term where υk is the mode frequency associated with the k-th normal mode whose creation (annihilation) operators are denoted via . The phonon operators follow the free-boson commutation relations, and . The seventh term presenting the exciton-phonon interactions is taken in the site-uncorrelated, local form and is characterized by the distribution of the corresponding coupling functions, . These phonon modes are responsible for exciton dephasing, relaxation and exciton-transport phenomena. The extent of these phenomena are governed by the site-dependent exciton-phonon coupling strengths gm,j which are taken as 1 (1.4) for the sites identified to represent Chl-A (Chl-b). While constructing the two-exciton phonon interaction Hamiltonian these coupling strengths are taken as 0.6 (gm,j + gn,j). Due to the differential coupling strengths, the excitations undergo dephasing at different rates, the one and two-exciton transport occur to a different extent over the sites leading to time-windowed interferences at a broader energy window. The inter-band dephasing, a quantity of particular interest in this article has contributions from the elements of the relaxation tensor that also participates in the cavity mediated one-exciton and two-exciton transport. Within this model, the single and the double excitons are coupled to a common set of phonon modes which can be characterized by the discrete distribution function, from which the spectral density function is obtained in the continuum frequency limit. The spectral density function, presented as, describes Nb = 48 multi-mode Brownian oscillator modes and one over-damped oscillator mode. The number of multi-mode Brownian oscillators (Nb), the respective spectral shift parameters (λ0, λj = υjϒj where υj and ϒj are the j-th oscillator frequencies and Huang-Rhys parameters respectively), and the damping parameter (γ0, γj) which are required to optimally describe the equilibrium spectral density function have been obtained from the work of [35, 59–61]. The numerical values corresponding to those parameters are enlisted in Appendix A2.
2.3 Polariton Manifolds
The polariton states are obtained as number-conserving manifolds via exact-diagonalization of the field-free Hamiltonian subspaces, by using , where we define operators as projectors onto the eigenstates of the polariton Hamiltonian (denoted by index n). The three distinct manifolds thus obtained corresponds to ground n = 0, single n = 1 and double n = 2 polariton manifolds with , , polariton states respectively (Figure 1). The manifolds are composed of the cluster of states which depend on an interplay of the hopping and cavity coupling parameters. The proposed spectroscopic technique utilizing entangled biphotons, in the parameter regime of interest, interact only via two-quantum interactions, which allows the joint-excitation manifold to be truncated at the level of the double polaritons. For the simulation, we have chosen the cavity mode frequency and the coupling strength parameters as ωc (cm−1) = 15.4 × 103 and gc (cm−1) = 0.1 × 103. The cavity frequency is resonant with the narrow frequency band around the frequency and off-resonant with the rest of the cluster of states in the bare absorption spectra of the excitonic system. The cavity coupling results in the polariton states to encompass a wider range than the cavity-free counterpart with the close-to-resonance states becoming affected to a greater degree. It can also be noted that the two-polariton states accommodate more states within the comparable bandwidth which is determined by the cavity coupling parameter. These two observations, combined, affirm that the exciton-cavity hybridization dominantly takes around the resonances, as expected. The spectral weights of the resultant delocalized polariton states have respective contributions from both the exciton and cavity modes, as determined by the exciton-cavity coupling matrix elements.
FIGURE 1
The progression of dephasing-broadened bare polariton band with the cavity frequencies and a fixed value of cavity coupling strength are presented. For the dephasing a representative set of values (generated from a random Gaussian distribution of mean and variance of μ1 = 25, σ1 = 10 and μ2 = 35 respectively) are chosen for clarity. In the displayed one-polariton (left panel) and the two-polariton band (right panel), the increase in the density of states in the latter within a comparable energy-window (e.g., ) is noticeable. The fundamental cavity frequencies along the x-axis are given in terms of detuning from the central region of the one-polariton band.
2.4 Polariton-Laser Interactions
The interaction between the external field i.e. biphoton sources and polaritons are treated within the optical dipolar interaction limit and within the rotating wave approximation. Corresponding Hamiltonian is given by,where external photon mode creation operators aj and dipole-weighted inter-manifold polariton transition operators, were defined. Also, the mode quantization volume V and the Fourier expansion frequency ωj for the photon modes have been introduced. These operators will be used to derive the signal expressions in the next sections. It is notable that in the spirit of the weak laser driving limit suitable for spectroscopy only the one-photon transition operators had been employed.
2.5 Phonon Induced State Broadening and Polariton Green’s Functions
In this section, we introduce the framework to obtain the state-dependent dephasing timescales and obtain the polariton Green’s functions required for the signal expressions. It is carried out by seeking an integral solution of the generalized master equation (written in the multi-polariton basis) obtained in the Markovian and the secular limit (often termed as the Redfield equation). The secular approximation limits the polariton-phonon mode interactions to be describable within the resonant cases only. The latter suffices our treatment of acoustic phonons mediate interactions are mediated by displacive perturbations of exciton states simultaneously dressed by the cavity interactions that are comparatively much stronger. The kinetic equation is given by, where we denote σ(n) as number resolved reduced multi-polariton density operator. We defined the super-operator as, . The polariton-phonon memory kernel capable of describing the phonon-induced polariton relaxation and dephasing (neglecting the effects of the driving field on the relaxation) is given by, . The kernel contains the dissipation operator D(t) = X (t − t′)Cb(t − t′) where we have the convolution of polariton-phonon coupling operator and time-domain correlation function given by,
In the above we have , , and, matsubara frequencies, νn = n (2π/β) (with nm = 20). Furthermore, we defined β = 1/κT where the Boltzmann constant κ and temperature T is denoted. The energy domain version of the relaxation kernel was used to extract the dephasing parameters. We define the line-broadening functions as, where we have Cb(Ω) = ∫dt exp iΩtC±(t) and ) which can be evaluated using the previous expressions. With the help of the line-broadening functions we estimate the inter-manifold dephasing parameters as, . Finally, we obtain the inter-manifold Greens functions relevant for dephasing as, . The advanced Green’s functions are defined likewise.
3 Multidimensional Double-Quantum Coherence Signal With Entangled Biphoton Sources
The theoretical development mentioned in the preceding sections provides us the ingredients to introduce the Double Quantum Coherence (DQC) signal. The DQC signal, as hinted earlier, involves excitation of state-selective (or narrow-band around the selected state) two-polariton coherence that is followed by projection of the oscillating polarization components onto two plausible sets of inter-manifold coherence. Provided that the participant two-polariton states are sufficiently correlated with the states onto which they are being projected, the desired nonlinear polarization have dominant component in the phase-matched direction ks = k1 + k2 − k3. Below we present a modular derivation of the signal starting from the time-dependent dual-perturbation scheme, corroborate the derivation to the Keldysh-Schwinger loop diagrams Figure 2 and introduce the entangled biphoton sources that serves as probe.
FIGURE 2
(A) Schematic description of the dynamical phenomena (in WMEL diagrammatic convention [36]) that can be addressed by the proposed DQC signal measurement technique. The first two interactions (in red) are common to both the pathways. The last two interactions (in black and gray) are distinguished by two alternative possibilities indicated by , pair. Entangled photons affect both these pair interactions {k1, k2} and { − k3, ks}. (B) Associated Keldysh-Schwinger diagrams that describe the corresponding dynamical pathways in the multi-polariton basis contributing to the signal and which are distinct in the last two interactions in time.
3.1 Double Quantum Coherence Signal
The DQC signal is typically generated by inducing four external field-matter interactions. In the case of time-domain (frequency-domain) classical field sources, three external fields with controllable delays (relative phases) are allowed to interact with the matter. The radiation field emitted by the time-dependent nonlinear polarization is registered, typically, via suitable heterodyning after another delay (spectrally dispersing) the signal. Typically, the deployment of quantum fields e.g., the biphoton sources to the measurement of DQC signals requires considerable care. These field sources are often parametrically scanned via schemes akin to multi-pulse phase-cycling. An involved discussion regarding the details of the deployment and measurement scheme is beyond the scope of the article [67, 68]. However, we assume that the biphoton generation scheme is capable of producing two sets of entangled photon pairs that lend themselves to external manipulation via central frequencies and delays. Before proceeding further, we introduce the two sets of pathways involved in the DQC signal generation as,
where the first and the second one have been written in the Heisenberg representation and correspond to diagrams in Figure 2). These pathways, notably, differing in the last two components signify dynamical spectral weights. Interference features between these pathway components depend on the nature of the polariton correlation and dephasing properties. We also introduce the time-domain, four-point external field correlation function as D (τ4, τ3, τ2, τ1) which is capable of incorporating the generalized nature of the external field. With the help of these definitions, we present the signal expression in the time-domain as, where Cs represents the coefficients that arises from the perturbative expansion and have been taken as a scaling factor. The Heaviside functions are defined in reference to the interaction times denoted along the loop diagrams. In going forward, the interaction times i.e. the loop time-instance variables were transformed to the loop-delay variables and mapped onto the real-time parameters that are externally tunable using, θ(s3)θ(s2)θ(s1) → θ(t3)θ(t2)θ(t1) for the diagram I and θ(s3)θ(s2)θ(s1) → θ(s3)θ(s2 − s3)θ(s1) = θ(t3)θ(t2)θ(t1) for the diagram II. These parameters are experimentally realizable. We also use the time-domain phonon-averaged Green’s functions expanded in the multi-polariton basis and write the field correlation function in the frequency domain. Following these exercises, we obtain the generalized signal as,
This expression is valid for a general class of DQC signal measurement which may use different kinds of external field sources beyond biphotons and simple Gaussian classical fields. The field correlation function acts as a convolutional probing function for the bare signal. Additional possibilities for the external manipulation of the field correlation function extend the applicability of DQC signals to a wide range of scenarios. Further, we aim to introduce a two-dimensional frequency-domain representation of the signal in order to facilitate a visualization of the correlation features contained in the matter correlation functions. In order to allow such representations to be generated via real-time delay-scanning protocols, we introduce integral transform, . We also define the following variable mapping of the parameters, , for the diagram I and , for the diagram II (the is the delays between the centering times). Scanning of these set of parameters and obtaining joint Fourier transforms w. r.t the delays generate the desired two-dimensional correlation plots. The final expression can be presented as the following expression.where the functions are specified as,
These functions encode the polariton dynamical resonances. These resonances show up, as predicted, during the scan of the Fourier transformed parameter. The field correlation function encodes the information about the ability to manipulate the spectral weights of the matter excitations and reveal desired dynamical resonances.
3.2 Entangled Biphoton Properties
The principal aim of using the entangled biphoton sources is to avail the non-classical relation between the joint time of arrival and frequency pairs of the biphotons. This in turn allows one to excite relatively short-lived two-polariton states (i.e., within an ultra-short time window) that are outside the excitation energy window of the classical two-photon laser pulses. These constraints remain difficult to surpass via independent variable manipulation, even in the case of multiple classical pulses. Below we present some basic features of the entangled photon source properties that were used in the simulation and describe their correlation features. The entangled biphoton field is traditionally generated via the spontaneous parametric down-conversion (SPDC) process (in the weak down-conversion limit) by pumping the source material with an ultra-short classical laser. The pump pulse bandwidth and the central frequency determine the correlation properties and time-frequency regime of the generated pairs. An effective Hamiltonian procedure which has been used to derive the correlation properties as outlined previously [47, 69–71] is avoided here for succinctness. Following a similar derivation, the entangled biphoton field correlation function can be obtained as, where we have, The function and the temporal entanglement parameter of the pairs , via , quantifies the spectral-temporal properties of the entangled biphotons [71–73]. The term A0 denotes the amplitude of the pump. This form of factorization underlies the fact that the signal scales linearly with the intensity. The temporal entanglement parameter can be viewed as an estimator of upper bound of delay between the time of generation of the entangled photon pairs inside the SPDC source material. The entanglement time parameter arises from the phase matching function Δk (ω1, ω2) and parametrically depends on the group velocity of propagation inside the SPDC material. The phase mismatching function under the linearization approximation around the central frequencies of the beams, for the collinear case leads to the identification of parameters, Tj = 1/vp − 1/vj where j ∈ {1, 2} and vp/j denotes the group velocity of the pump, and two biphotons inside the material. These values can be estimated from the inverse of the marginal distribution function of the joint spectral amplitude function in the frequency domain [51, 73]. In contrast, the classical field, in a similar weak-field limit factorizes into the product of amplitudes as, and scales quadratically with the intensity. The biphoton field correlation properties for different typical parameter regimes can be examined by plotting the joint-spectral amplitude which describes the frequency-dependent correlation of the same as shown in Figure 3. The bottom-right plot (i.e., (d)) in Figure 3 corresponds to a classical field scenario. The simulation has the freedom of selecting frequency pairs from the plot region where the function has finite support.
FIGURE 3
Field correlation functions for different entangled photon states corresponding to parameter variation of temporal entanglement parameter Tent (fs) and pump width τp (fs) for uniform pump frequency 31.0, ×, 103cm−1 and biphoton frequencies ω1 = ω2 = 15.5 × 103cm−1. The bottom right corresponds to the typical case of a classical field two-photon pulse (A) corresponds to the Tent = 10; τp = 20 (B) corresponds to Tent = 10; τp = 50 while (C) assumes values Tent = 50; τp = 20. The (D) corresponds to a Gaussian pair-pulse with τg (fs) = 10 with ω1,g = ω2,g = 15.5 × 103cm−1.
3.3 Simulation
The correlated two-polariton excitations via entangled biphoton sources may focus on several experimental configurations which will be of particular interest to the condensed phase spectroscopies. The excitation of specific two-polariton states via higher-energy sectors of the one-polariton manifold and contrasting them with those via the lower-energy sectors may give information about polariton scattering, delocalization, and dephasing. In other words, the specific two-polariton states may have dominant contributions from certain one-polariton states which are distant on the site basis but energetically closer. Alternatively for the same two-polariton excitation, projecting to the higher and lower-energy one-polariton sector offers insights into the state compositions. Combining two strategies may provide important insight into the state resolved polariton correlations. These features can be probed as shown in the upper and lower panel of Figure 4. Particularly the polariton states (e.g., ) that are specifically prone to phonon-induced dephasing (higher dephasing induced broadened) have been excited while the temporal entanglement parameter of the probe, projected one-polariton sector has been varied.
FIGURE 4
Multidimensional correlation plots (A–F) for signal with parametric variation of the frequencies and the temporal entanglement parameters of the entangled biphotons. The short-lived, intermediate one-polariton resonances along the Ω3 axis can be explored and associated with the two-polariton resonances along the Ω2 axis. For parameters corresponding to each plot see text.
In the upper panel of Figure 4, we present, along the rows, three sets of results for the variation of the temporal entanglement parameter . Each of them correspond to fixed and (therefore exciting the target at ) with pump width τp (fs) = 20.0. Therefore, the two polariton excitation occurs via the middle sector of the one-polariton band. It also projects the two-polariton coherence to the mixed-energy region of the one-polariton band by choosing the corresponding frequencies as and . It is noticed that the strong entanglement between the photon pairs, as we move from (A) to (C) corresponding to values 60.0, 50.0, 40.0 respectively, allows correlated signal features to develop. These emergent features may not be visible in the signal obtained by using classical pulse pairs of comparable spectral-temporal properties. Time-frequency correlated excitation gives the freedom of simultaneously choosing the narrow-band target while ensuring that system remains less affected by the high dephasing components of the intermediate states. The latter has the possibility of making the short-time kinetics in the one-polariton manifold more accessible. Even an individually controlled classical field two-pulse scenario may offer less advantage because the temporal and spectral components are bound.
The bottom row ((d) to (f)) accomplishes the aforementioned goal of exploring excitation via different energy sectors of the one-polariton manifold. Here the (d) and (e) allow excitation via middle-sector and (f) lower-sector while projecting all of them to the same mixed-energy sectors as the above panel. With the increase in temporal entanglement parameter in going from (d) to (e) ( values 40.0, 10.0, respectively) we find features shows distinguishable increase. The last panel (f) whose excitation via lower energy sector (excitation via and ) which reveals many more correlation features. It reveals the higher participation of the particular set of one-polariton states to chosen target in the two-polariton manifold throughout the simulation. The short temporal entanglement parameter also certifies the capability of the biphoton sources to map out correlation involving energetically distant states. The success of the parameter regime and overall strategy of the simulation can be traced back to the ability to choose the frequency pairs from a broader distribution. In addition, the fact that they are also bound by the temporal constraints allow excitation via fast dephasing components in the one-polariton band (during the first time delay i.e., T1) while simultaneously projecting the resultant two-polariton coherence to short-lived coherences (during the last time delay i.e., T3).
4 Conclusion and Outlook
In this article, we have proposed a theoretical protocol that is suitable for the investigation of the inter-manifold coherence properties associated with the two-polariton manifold. We observe that the proposal exclusively focused on the cavity control of exciton correlation and modulation of exciton-phonon dephasing via the former. In other words, the cavity affects the exciton transitions directly and redistributes the excitonic spectral weights. The extent to which this redistribution occurs is encoded in the action of the polariton transformation matrices. As a result, the novel hybridized polariton states interact with the phonons rather differently than that of the cavity-free case. We demonstrated that biphoton sources are capable of studying the ultrafast signatures of the related dynamics without losing the state specificity. In the process, it is capable of mapping out the cavity modulated exciton correlation. In this direction, two further extensions namely, the detailed study involving the parametric variation of the cavity-free case and a comparison employing controlled classical fields are worth looking at. They will be part of future communication. However, we note that the signal expressions presented in this communication will be sufficient for such extended analysis.
The role of cavity coupling has been included non-perturbatively with the same quasi-particle excitation picture. Here the real-space coupling variations were neglected for convenience, in the spirit of the first simulation. Also, the role of the mean number of photons in the cavity has not been investigated and the role of the cavity has been confined to the coherence created between the states within the proposed regime of operation.
In comparing and contrasting the present technique to the transmission mode pump-probe measurements several features distinguish the present technique. DQC measures one specific component of the nonlinear polarization of the cavity polariton, unlike the pump-probe analog. The pump-probe technique also includes the pathways analogous to polaritonic Raman scattering pathways. Thus the DQC signal is more specific to the purpose of this article. The two-polariton coherence is explicitly monitored by isolating the signal components as suggested by the plausible implementation via phase-cycling schemes. It can be also highlighted that the present technique works by projecting the two-polariton coherence in two competing coherence components of lower order. The relevant processes leading to the signal occur during the last two time intervals. The degree of discrimination of the pathways, reflected in the associated dynamical spectral weights decides the magnitude of the signal. In contrast to the pump-probe studies, the off-diagonal spectral signatures in the correlation plots required to be interpreted differently. In the latter, the cooperative features appearing in the above-mentioned sector carry less specific information regarding the physical origin of the polaritonic correlation due to non-discrimination between pure two-one polariton coherence and two-polariton-one polariton coherence.
The study can be extended to accommodate the explicit two-polariton transport phenomenology by studying the fluorescence-detected phase cycling protocols [74, 75]. The latter is a four-wave mixing analog but aimed at investigating the longer-time state correlation properties in the presence of phonon-induced dynamical population redistribution. The longer-time dynamical information thereby obtained is complementary to the information provided by the short-time dynamics investigated in this article. An investigation in this direction is on the way. Furthermore, one may combine the interferometric detection schemes to separate the pathways as recently proposed [76, 77].
The theoretical description adopted for describing the dissipative polaritonic matter is formulated at the level of a quasi-particle approach. The quasi-particle Green’s function was chosen to describe the signal [78]. Over the last few years, a host of promising methods have been proposed that can potentially deal with the complexity of the quantum aggregates in near future. These methods have offered several different flavors of treating the matter and the cavity modes. We have, namely, quantum electrodynamics based hybrid (density) functional formulation [79–82], cluster-expansion [21], potential-energy surface-based dynamical calculations combined with the trajectory-based propagation for the cavity quadrature modes [83, 84], path integral based unified framework for nuclear modes and the idealized cavity modes [85–87], first-principles simulations [88–92]. Incorporating these methods to describe the phenomenology described in this article will require a qualified description of exciton formation, a description of exciton-exciton scattering in the presence of a dielectric environment, and nuclear propagation.
We also note that the numerical simulation adopted the analytical expressions which have been derived under the assumption of generalized time-translational invariance. For systems driven out-of-equilibrium via additional laser pulses, one may expect to see more correlation features in the signal. The scope to add additional pulses and using the biphotons as probes, although offers a more complicated scenario, is a promising avenue for future study of nonlinear response in correlated quantum materials [93–97].
Statements
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Author contributions
AD developed the project in consultation with AR. AD performed the research. AD and AR analyzed the results. AD wrote the manuscript with inputs from AR.
Acknowledgments
We acknowledge financial support from the European Research Council (ERC-2015-AdG-694097), by the Cluster of Excellence “Advanced Imaging of Matter” (AIM), Grupos Consolidados (IT1249-19) and SFB925 “Light induced dynamics and control of correlated quantum systems.” The Flatiron Institute is a division of the Simons Foundation.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1.
BoulaisÉSawayaNPDVenezianoRAndreoniABanalJLKondoTet alProgrammed Coherent Coupling in a Synthetic DNA-Based Excitonic Circuit. Nat Mater (2018) 17:159–66. 10.1038/nmat5033
2.
RomeroEAugulisRNovoderezhkinVIFerrettiMThiemeJZigmantasDet alQuantum Coherence in Photosynthesis for Efficient Solar-Energy Conversion. Nat Phys (2014) 10:676–82. 10.1038/nphys3017
3.
ScholesGDFlemingGRChenLXAspuru-GuzikABuchleitnerACokerDFet alUsing Coherence to Enhance Function in Chemical and Biophysical Systems. Nature (2017) 543:647–56. 10.1038/nature21425
4.
Schlau-CohenGSIshizakiACalhounTRGinsbergNSBallottariMBassiRet alElucidation of the Timescales and Origins of Quantum Electronic Coherence in LHCII. Nat Chem (2012) 4:389–95. 10.1038/nchem.1303
5.
CaoJCogdellRJCokerDFDuanH-GHauerJKleinekathöferUet alQuantum Biology Revisited. Sci Adv (2020) 6:eaaz4888. 10.1126/sciadv.aaz4888
6.
ScholesGDPolaritons and Excitons: Hamiltonian Design for Enhanced Coherence. Proc R Soc A (2020) 476:20200278. 10.1098/rspa.2020.0278
7.
BasovDAsenjo-GarciaASchuckPJZhuXRubioAPolariton Panorama. Nanophotonics (2021) 10:549.
8.
RuggenthalerMTancogne-DejeanNFlickJAppelHRubioAFrom a Quantum-Electrodynamical Light–Matter Description to Novel Spectroscopies. Nat Rev Chem (2018) 2:1. 10.1038/s41570-018-0118
9.
MewesLWangMIngleRABörjessonKCherguiMEnergy Relaxation Pathways between Light-Matter States Revealed by Coherent Two-Dimensional Spectroscopy. Commun Phys (2020) 3:1. 10.1038/s42005-020-00424-z
10.
Garcia-VidalFJCiutiCEbbesenTWManipulating Matter by strong Coupling to Vacuum fields. Science (2021) 373:eabd0336. 10.1126/science.abd0336
11.
CaoJGeneralized Resonance Energy Transfer Theory: Applications to Vibrational Energy Flow in Optical Cavities. arXiv preprint arXiv:2201.12117 (2022). 10.48550/ARXIV.2201.12117
12.
EngelhardtGCaoJUnusual Dynamical Properties of Disordered Polaritons in Micocavities. arXiv preprint arXiv:2112.04060 (2021) 105:064205.
13.
GroenhofGToppariJJCoherent Light Harvesting through Strong Coupling to Confined Light. J Phys Chem Lett (2018) 9:4848–51. 10.1021/acs.jpclett.8b02032
14.
LatiniSShinDSatoSASchäferCDe GiovanniniUHübenerHet alThe Ferroelectric Photo Ground State of SrTiO 3: Cavity Materials Engineering. Proc Natl Acad Sci (2021) 118. 10.1073/pnas.2105618118
15.
LengersFKuhnTReiterDEPhonon Signatures in Spectra of Exciton Polaritons in Transition Metal Dichalcogenides. Phys Rev B (2021) 104:L241301. 10.1103/physrevb.104.l241301
16.
AutryTMNardinGSmallwoodCLSilvermanKBajoniDLemaîtreAet alExcitation Ladder of Cavity Polaritons. Phys Rev Lett (2020) 125:067403. 10.1103/physrevlett.125.067403
17.
SidlerDSchäferCRuggenthalerMRubioAPolaritonic Chemistry: Collective Strong Coupling Implies Strong Local Modification of Chemical Properties. J Phys Chem Lett (2020) 12:508–16. 10.1021/acs.jpclett.0c03436
18.
YangP-YCaoJQuantum Effects in Chemical Reactions under Polaritonic Vibrational Strong Coupling. J Phys Chem Lett (2021) 12:9531–8. 10.1021/acs.jpclett.1c02210
19.
SalijATempelaarRMicroscopic Theory of Cavity-Confined Monolayer Semiconductors: Polariton-Induced valley Relaxation and the prospect of Enhancing and Controlling valley Pseudospin by Chiral strong Coupling. Phys Rev B (2021) 103:035431. 10.1103/physrevb.103.035431
20.
LatiniSDe GiovanniniUSieEJGedikNHübenerHRubioAPhonoritons as Hybridized Exciton-Photon-Phonon Excitations in a Monolayer H -BN Optical Cavity. Phys Rev Lett (2021) 126:227401. 10.1103/physrevlett.126.227401
21.
HauglandTSRoncaEKjønstadEFRubioAKochHCoupled Cluster Theory for Molecular Polaritons: Changing Ground and Excited States. Phys Rev X (2020) 10:041043. 10.1103/physrevx.10.041043
22.
SidlerDRuggenthalerMSchäferCRoncaERubioA,A Perspective on Ab Initio Modeling of Polaritonic Chemistry: The Role of Non-equilibrium Effects and Quantum Collectivity. arXiv preprint arXiv:2108.12244 (2021).
23.
SchäferCFlickJRoncaENarangPRubioAShining Light on the Microscopic Resonant Mechanism Responsible for Cavity-Mediated Chemical Reactivity. arXiv preprint arXiv:2104.12429 (2021).
24.
ZhangZScullyMOAgarwalGSQuantum Entanglement between Two Magnon Modes via Kerr Nonlinearity Driven Far from Equilibrium. Phys Rev Res (2019) 1:023021. 10.1103/physrevresearch.1.023021
25.
ZhangZAgarwalGSScullyMOQuantum Fluctuations in the Fröhlich Condensate of Molecular Vibrations Driven Far from Equilibrium. Phys Rev Lett (2019) 122:158101. 10.1103/physrevlett.122.158101
26.
WangDKelkarHMartin-CanoDUtikalTGötzingerSSandoghdarVCoherent Coupling of a Single Molecule to a Scanning Fabry-Perot Microcavity. Phys Rev X (2017) 7:021014. 10.1103/physrevx.7.021014
27.
MaserAGmeinerBUtikalTGötzingerSSandoghdarVFew-photon Coherent Nonlinear Optics with a Single Molecule. Nat Photon (2016) 10:450–3. 10.1038/nphoton.2016.63
28.
Campos-Gonzalez-AnguloJAYuen-ZhouJGeneralization of the Tavis-Cummings Model for Multi-Level Anharmonic Systems: Insights on the Second Excitation Manifold (2022) 156:194308. 10.1063/5.0087234
29.
DebnathARubioAEntangled Photon Assisted Multidimensional Nonlinear Optics of Exciton-Polaritons. J Appl Phys (2020) 128:113102. 10.1063/5.0012754
30.
DelPoCAKudischBParkKHKhanS-U -ZFassioliFFaustiDet alPolariton Transitions in Femtosecond Transient Absorption Studies of Ultrastrong Light-Molecule Coupling. J Phys Chem Lett (2020) 11:2667–74. 10.1021/acs.jpclett.0c00247
31.
RibeiroRFCampos-Gonzalez-AnguloJAGiebinkNCXiongWYuen-ZhouJEnhanced Optical Nonlinearities under Collective strong Light-Matter Coupling. Phys Rev A (2021) 103:063111. 10.1103/physreva.103.063111
32.
RenkenSPandyaRGeorgiouKJayaprakashRGaiLShenZet alUntargeted Effects in Organic Exciton-Polariton Transient Spectroscopy: A Cautionary Tale. J Chem Phys (2021) 155:154701. 10.1063/5.0063173
33.
ZhangZPengTNieXAgarwalGSScullyMODirectional Superradiant Emission from Statistically Independent Incoherent. arXiv preprint arXiv:2106.10988 (2021). 10.48550/ARXIV.2106.10988
34.
ZhangZWangKYiZZubairyMSScullyMOMukamelSPolariton-Assisted Cooperativity of Molecules in Microcavities Monitored by Two-Dimensional Infrared Spectroscopy. J Phys Chem Lett (2019) 10:4448–54. 10.1021/acs.jpclett.9b00979
35.
ZhangZSaurabhPDorfmanKEDebnathAMukamelSMonitoring Polariton Dynamics in the LHCII Photosynthetic Antenna in a Microcavity by Two-Photon Coincidence Counting. J Chem Phys (2018) 148:74302. 10.1063/1.5004432
36.
LeeDAlbrechtACThe Vibronic Theory of Resonance Hyper-Raman Scattering. Adv infrared Raman Spectrosc (1985) 12:179.
37.
TempelaarRJansenTLCKnoesterJExciton-Exciton Annihilation Is Coherently Suppressed in H-Aggregates, but Not in J-Aggregates. J Phys Chem Lett (2017) 8:6113–7. 10.1021/acs.jpclett.7b02745
38.
Gutiérrez-MezaEMalatestaRLiHBargigiaISrimath KandadaARValverde-ChávezDAet alFrenkel Biexcitons in Hybrid HJ Photophysical Aggregates. Sci Adv (2021) 7:eabi5197. 10.1126/sciadv.abi5197
39.
MukamelSPrinciples of Nonlinear Optical Spectroscopy, 6. New York: Oxford University Press on Demand (1999).
40.
LomsadzeBCundiffSTLine-shape Analysis of Double-Quantum Multidimensional Coherent Spectra. Phys Rev A (2020) 102:043514. 10.1103/physreva.102.043514
41.
GaoFCundiffSTLiHProbing Dipole-Dipole Interaction in a Rubidium Gas via Double-Quantum 2D Spectroscopy. Opt Lett (2016) 41:2954. 10.1364/ol.41.002954
42.
KimJMukamelSScholesGDTwo-Dimensional Electronic Double-Quantum Coherence Spectroscopy. Acc Chem Res (2009) 42:1375–84. 10.1021/ar9000795
43.
SvozilíkJPeřinaJLeón-MontielRd. JVirtual-state Spectroscopy with Frequency-Tailored Intense Entangled Beams. J Opt Soc Am B (2018) 35:460. 10.1364/josab.35.000460
44.
León-MontielRd. JSvozilíkJTorresJPU’RenABTemperature-Controlled Entangled-Photon Absorption Spectroscopy. Phys Rev Lett (2019) 123:023601. 10.1103/physrevlett.123.023601
45.
KangGNasiri AvanakiKMosqueraMABurdickRKVillabona-MonsalveJPGoodsonTIIIet alEfficient Modeling of Organic Chromophores for Entangled Two-Photon Absorption. J Am Chem Soc (2020) 142:10446–58. 10.1021/jacs.0c02808
46.
SchlawinFDorfmanKEMukamelSEntangled Two-Photon Absorption Spectroscopy. Acc Chem Res (2018) 51:2207–14. 10.1021/acs.accounts.8b00173
47.
SalehBEAJostBMFeiH-BTeichMCEntangled-Photon Virtual-State Spectroscopy. Phys Rev Lett (1998) 80:3483–6. 10.1103/physrevlett.80.3483
48.
RichterMMukamelSUltrafast Double-Quantum-Coherence Spectroscopy of Excitons with Entangled Photons. Phys Rev A (2010) 82:013820. 10.1103/physreva.82.013820
49.
RoslyakOMukamelSMultidimensional Pump-Probe Spectroscopy with Entangled Twin-Photon States. Phys Rev A (2009) 79:063409. 10.1103/physreva.79.063409
50.
MukamelSFreybergerMSchleichWBelliniMZavattaALeuchsGet alRoadmap on Quantum Light Spectroscopy. J Phys B: Mol Opt Phys (2020) 53:072002. 10.1088/1361-6455/ab69a8
51.
LandesTRaymerMGAllgaierMMerkoucheSSmithBJMarcusAHQuantifying the Enhancement of Two-Photon Absorption Due to Spectral-Temporal Entanglement. Opt Express (2021) 29:20022. 10.1364/oe.422544
52.
BittnerERLiHPiryatinskiASrimath KandadaARSilvaCProbing Exciton/exciton Interactions with Entangled Photons: Theory. J Chem Phys (2020) 152:071101. 10.1063/1.5139197
53.
SchlawinFBuchleitnerATheory of Coherent Control with Quantum Light. New J Phys (2017) 19:013009. 10.1088/1367-2630/aa55ec
54.
LiHPiryatinskiASrimath KandadaARSilvaCBittnerERPhoton Entanglement Entropy as a Probe of many-body Correlations and Fluctuations. J Chem Phys (2019) 150:184106. 10.1063/1.5083613
55.
RichterMSinghRSiemensMCundiffSTDeconvolution of Optical Multidimensional Coherent Spectra. Sci Adv (2018) 4:eaar7697. 10.1126/sciadv.aar7697
56.
AgranovichVToshichBCollective Properties of Frenkel Excitons. Sov Phys JETP (1968) 26:104.
57.
ChernyakVZhangWMMukamelSMultidimensional Femtosecond Spectroscopies of Molecular Aggregates and Semiconductor Nanostructures: The Nonlinear Exciton Equations. J Chem Phys (1998) 109:9587–601. 10.1063/1.477621
58.
AbramaviciusDPalmieriBVoronineDVŠandaFMukamelSCoherent Multidimensional Optical Spectroscopy of Excitons in Molecular Aggregates; Quasiparticle versus Supermolecule Perspectives. Chem Rev (2009) 109:2350–408. 10.1021/cr800268n
59.
NovoderezhkinVIPalaciosMAVan AmerongenHVan GrondelleRExcitation Dynamics in the LHCII Complex of Higher Plants: Modeling Based on the 2.72 Å Crystal Structure. J Phys Chem B (2005) 109:10493–504. 10.1021/jp044082f
60.
NovoderezhkinVIPalaciosMAVan AmerongenHVan GrondelleREnergy-Transfer Dynamics in the LHCII Complex of Higher Plants: Modified Redfield Approach. J Phys Chem B (2004) 108:10363–75. 10.1021/jp0496001
61.
NovoderezhkinVMarinAvan GrondelleRIntra- and Inter-monomeric Transfers in the Light Harvesting LHCII Complex: the Redfield–Förster Picture. Phys Chem Chem Phys (2011) 13:17093.
62.
van AmerongenHvan GrondelleRUnderstanding the Energy Transfer Function of Lhcii, the Major Light-Harvesting Complex of green Plants (2001). 10.1021/jp0028406
63.
De LiberatoSLight-Matter Decoupling in the Deep Strong Coupling Regime: The Breakdown of the Purcell Effect. Phys Rev Lett (2014) 112:016401. 10.1103/physrevlett.112.016401
64.
SchäferCRuggenthalerMRokajVRubioARelevance of the Quadratic Diamagnetic and Self-Polarization Terms in Cavity Quantum Electrodynamics. ACS Photon (2020) 7:975–90. 10.1021/acsphotonics.9b01649
65.
KeelingJCoulomb Interactions, Gauge Invariance, and Phase Transitions of the Dicke Model. J Phys Condens Matter (2007) 19:295213. 10.1088/0953-8984/19/29/295213
66.
VukicsAGrießerTDomokosPElimination of the A-Square Problem from Cavity QED. Phys Rev Lett (2014) 112:073601. 10.1103/physrevlett.112.073601
67.
DorfmanKEMukamelSMultidimensional Spectroscopy with Entangled Light: Loop vs Ladder Delay Scanning Protocols. New J Phys (2014) 16:033013. 10.1088/1367-2630/16/3/033013
68.
RaymerMGLandesTMarcusAHEntangled Two-Photon Absorption by Atoms and Molecules: A Quantum Optics Tutorial. J Chem Phys (2021) 155:081501. 10.1063/5.0049338
69.
CutipaPChekhovaMVBright Squeezed Vacuum for Two-Photon Spectroscopy: Simultaneously High Resolution in Time and Frequency, Space and Wavevector. Opt Lett (2022) 47:465. 10.1364/ol.448352
70.
AndersenULGehringTMarquardtCLeuchsG30 Years of Squeezed Light Generation. Phys Scr (2016) 91:053001. 10.1088/0031-8949/91/5/053001
71.
ArzaniFFabreCTrepsNVersatile Engineering of Multimode Squeezed States by Optimizing the Pump Spectral Profile in Spontaneous Parametric Down-Conversion. Phys Rev A (2018) 97:033808. 10.1103/physreva.97.033808
72.
SchlawinFDorfmanKEFingerhutBPMukamelSSuppression of Population Transport and Control of Exciton Distributions by Entangled Photons. Nat Commun (2013) 4:1. 10.1038/ncomms2802
73.
KellerTERubinMHTheory of Two-Photon Entanglement for Spontaneous Parametric Down-Conversion Driven by a Narrow Pump Pulse. Phys Rev A (1997) 56:1534–41. 10.1103/physreva.56.1534
74.
MalỳPMančalTSignatures of Exciton Delocalization and Exciton–Exciton Annihilation in Fluorescence-Detected Two-Dimensional Coherent Spectroscopy. J Phys Chem Lett (2018) 9:5654. 10.1021/acs.jpclett.8b02271
75.
MalỳPBrixnerTFluorescence-Detected Pump–Probe Spectroscopy. Angew Chem Int Edition (2021) 60:18867. 10.1002/anie.202102901
76.
DorfmanKEAsbanSGuBMukamelSHong-Ou-Mandel Interferometry and Spectroscopy Using Entangled Photons. Commun Phys (2021) 4:1. 10.1038/s42005-021-00542-2
77.
AsbanSMukamelSDistinguishability and “Which Pathway” Information in Multidimensional Interferometric Spectroscopy with a Single Entangled Photon-Pair. Sci Adv (2021) 7:eabj4566. 10.1126/sciadv.abj4566
78.
AbrikosovAAGorkovLPDzyaloshinskiIEMethods of Quantum Field Theory in Statistical Physics. New York: Courier Corporation (2012).
79.
RuggenthalerMFlickJPellegriniCAppelHTokatlyIVRubioAQuantum-electrodynamical Density-Functional Theory: Bridging Quantum Optics and Electronic-Structure Theory. Phys Rev A (2014) 90:012508. 10.1103/physreva.90.012508
80.
SchäferCBuchholzFPenzMRuggenthalerMRubioAMaking Ab Initio QED Functional(s): Nonperturbative and Photon-free Effective Frameworks for strong Light-Matter Coupling. Proc Natl Acad Sci (2021) 118. 10.1073/pnas.2110464118
81.
FlickJNarangPCavity-Correlated Electron-Nuclear Dynamics from First Principles. Phys Rev Lett (2018) 121:113002. 10.1103/physrevlett.121.113002
82.
YangJOuQPeiZWangHWengBShuaiZet alQuantum-electrodynamical Time-dependent Density Functional Theory within Gaussian Atomic Basis. J Chem Phys (2021) 155:064107. 10.1063/5.0057542
83.
HoffmannNMLacombeLRubioAMaitraNTEffect of many Modes on Self-Polarization and Photochemical Suppression in Cavities. J Chem Phys (2020) 153:104103. 10.1063/5.0012723
84.
ZhangYNelsonTTretiakSNon-adiabatic Molecular Dynamics of Molecules in the Presence of strong Light-Matter Interactions. J Chem Phys (2019) 151:154109. 10.1063/1.5116550
85.
LiXMandalAHuoPCavity Frequency-dependent Theory for Vibrational Polariton Chemistry. Nat Commun (2021) 12:1. 10.1038/s41467-021-21610-9
86.
ChowdhurySNMandalAHuoPRing Polymer Quantization of the Photon Field in Polariton Chemistry. J Chem Phys (2021) 154:044109. 10.1063/5.0038330
87.
MandalALiXHuoPTheory of Vibrational Polariton Chemistry in the Collective Coupling Regime. J Chem Phys (2022) 156:014101. 10.1063/5.0074106
88.
SchäferCRuggenthalerMAppelHRubioAModification of Excitation and Charge Transfer in Cavity Quantum-Electrodynamical Chemistry. Proc Natl Acad Sci U.S.A (2019) 116:4883–92. 10.1073/pnas.1814178116
89.
FlickJRuggenthalerMAppelHRubioAAtoms and Molecules in Cavities, from Weak to strong Coupling in Quantum-Electrodynamics (QED) Chemistry. Proc Natl Acad Sci U.S.A (2017) 114:3026–34. 10.1073/pnas.1615509114
90.
FlickJWelakuhDMRuggenthalerMAppelHRubioALight-Matter Response in Nonrelativistic Quantum Electrodynamics. ACS Photon (2019) 6:2757–78. 10.1021/acsphotonics.9b00768
91.
WelakuhDMFlickJRuggenthalerMAppelHRubioAFrequency-dependent Sternheimer Linear-Response Formalism for Strongly Coupled Light-Matter Systems. arXiv preprint arXiv:2201.08734 (2022). 10.48550/ARXIV.2201.08734
92.
SvendsenMKKurmanYSchmidtPKoppensFKaminerIThygesenKSCombining Density Functional Theory with Macroscopic QED for Quantum Light-Matter Interactions in 2D Materials. Nat Commun (2021) 12:1. 10.1038/s41467-021-23012-3
93.
RostamiHKatsnelsonMIVignaleGPoliniMGauge Invariance and Ward Identities in Nonlinear Response Theory. Ann Phys (2021) 431:168523. 10.1016/j.aop.2021.168523
94.
ParameswaranSAGopalakrishnanSAsymptotically Exact Theory for Nonlinear Spectroscopy of Random Quantum Magnets. Phys Rev Lett (2020) 125:237601. 10.1103/physrevlett.125.237601
95.
RostamiHCappellutiEDominant Role of Two-Photon Vertex in Nonlinear Response in Two-Dimensional Dirac Systems. npj 2D Mater Appl (2021) 5:1. 10.1038/s41699-021-00217-0
96.
ChoiWLeeKHKimYBTheory of Two-Dimensional Nonlinear Spectroscopy for the Kitaev Spin Liquid. Phys Rev Lett (2020) 124:117205. 10.1103/physrevlett.124.117205
97.
GrankinAHafeziMGalitskiVMEnhancement of Superconductivity with External Phonon Squeezing. Phys Rev B (2021) 104:L220503. 10.1103/physrevb.104.l220503
Appendix
In this appendix, we provide two sets of supplementary information that aid the main text. It is composed of a relevant figure of merit for entanglement in the biphoton source and a description of parameters regarding the phonon spectral density.
A1 Singular Value Decomposition of the Biphoton Spectral Function
The singular value decomposition of the joint spectral amplitude of the biphoton sources provides an estimate for the number of effective modes that participated in mode-squeezing. It can be analyzed by obtaining the singular value decomposition of the said function at different parameter regimes. It can be seen in Figure A1 that the cases involving the shorter temporal entanglement parameters increase the number of effectively squeezed modes. Also, a decrease in the pump bandwidth decreases the modal amplitudes for a case with comparable temporal entanglement parameters. The classical two-photon pulse, unsurprisingly, presents a plot reminiscent of the uncorrelated feature.
FIGURE A1
Singular values for the three sets of entangled biphoton sources corresponding to a similar parameter regime of the ones used in the simulations of Figure 3. For the first two plots [i.e., in (A,B)] the values are truncated at nsvd = 50, and for the third one [i.e., in (C)], at n = 25 (after normalization). The classical two-photon pulse yields, as expected, nsvd = 1.
A2 Phonon Parameters
The site-independent spectral function is composed of discrete frequencies which primarily induce multiple timescales Markovian dissipation. Here we enlist the parameter values for the 48 structured phonon modes used in the simulation [59–62]. Along with these values, the values of other parameters are given as, γj (cm−1) = 30.0 for all the multimode Brownian oscillators. Corresponding Figure A2 illustrates the spectral density distribution. The parameter values for the overdamped oscillator are given by, λ0 (cm−1) = 37.0 and γ0 (cm−1) = 30.0.
FIGURE A2
(A,B) The list of the υj (cm−1) parameters corresponding to the multi-mode Brownian oscillators.The list of the λj (cm−1) parameters corresponding to the multi-mode Brownian oscillators. (C) The spectral function of the structured phonon modes (corresponding to the multi-mode Brownian oscillators) used in the simulation.
Summary
Keywords
entanglement, exciton-polaritons, photosynthetic complex, nonlinear optics, multidimensional spectroscopy, quantum cavity electrodynamics, dissipation
Citation
Debnath A and Rubio A (2022) Entangled Biphoton Enhanced Double Quantum Coherence Signal as a Probe for Cavity Polariton Correlations in Presence of Phonon Induced Dephasing. Front. Phys. 10:879113. doi: 10.3389/fphy.2022.879113
Received
18 February 2022
Accepted
06 May 2022
Published
04 July 2022
Volume
10 - 2022
Edited by
Juan Torres, The Institute of Photonic Sciences (ICFO), Spain
Reviewed by
Che-Ming Li, National Cheng Kung University, Taiwan
Bin Fang, University of Texas at Austin, United States
Updates
Copyright
© 2022 Debnath and Rubio.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Arunangshu Debnath, arunangshu.debnath@desy.de; Angel Rubio, angel.rubio@mpsd.mpg.de
This article was submitted to Quantum Engineering and Technology, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.