PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Phys Rev A. Author manuscript; available in PMC 2010 July 2.
Published in final edited form as:
PMCID: PMC2895952
NIHMSID: NIHMS212370

Multidimensional pump-probe spectroscopy with entangled twin-photon states

Abstract

We show that entangled photons may be used in coherent multidimensional nonlinear spectroscopy to provide information on matter by scanning photon wave function parameters (entanglement time and delay of twin photons), rather than frequencies and time delays, as is commonly done with classical pulses. Signals are expressed and interpreted intuitively in terms of products of matter and field correlation functions using a diagrammatic close time path loop formalism which reveals the entangled quantum pathways of the fields and matter. The pump-probe signal measured when the pump and the probe are in a twin entangled state shows two-photon resonant contributions which scale linearly rather than quadratically with the incident beam intensity and reveal frequencies of off-resonant transitions. Two-dimensional spectrograms obtained by double Fourier transform of the signal with respect to the entanglement time and delay of the twins could provide detailed information on correlations among states and dynamical processes with high temporal resolution. The analogy with multidimensional time-domain optical techniques which use sequences of short classical pulses and pulse shaping algorithms is pointed out.

I. INTRODUCTION

Multidimensional techniques have long been used in NMR to study the structure and dynamics of complex molecules by measuring the response of spins to sequences of magnetic pulses [1]. Correlation plots obtained by varying several time delays spread the information in multiple dimensions thus improving the resolution and obtaining information unavailable from ordinary one-dimensional (1D) spectroscopy (e.g., absorption). Over the past 15 years these techniques have been extended to the femtosecond regime using infrared and visible pulses [24] and were applied to many types of systems and relaxation processes including vibrations in proteins, protein folding, hydrogen bonded liquids, and exciton dynamics in electronically excited photo-synthetic complexes [5]. Pulse shaping techniques make it possible to manipulate and control these signals [611]. All of these measurements employ classical fields and use multiple pulse delays or frequencies as control parameters.

In this paper we propose multidimensional signals which exploit the quantum nature of the field and are based on the manipulation of photon wave functions. These can now be engineered to create entangled photons which depend on multiple parameters that can be varied to generate two-dimensional (2D) correlation spectra. Entanglement is a purely quantum type of correlation among particles which exists even when they are spatially well separated [12,13]. Many applications to quantum computing, information and communication protocols, cryptography, and teleportation are based on the manipulation of entangled particles [1416]. Nonlocal quantum correlations brought about by entanglement have been used to test the foundations of quantum mechanics by demonstrating violation of Bell's inequalities and resolving the Einstein-Podolsky-Rosen (EPR) paradox [1719]. It is easier to maintain entanglement of photons than of massive particles. Entangled photons have been generated by a variety of schemes. These include atomic cascades [2022], parametric down conversion (PDC) [23,24], or hyperparametric scattering (HPS) [25,26] in nonlinear crystals and, most recently, biexciton decay in semiconductors [27,28]. The photons may be entangled with respect to their polarization, momentum-position, energy-time, their orbital angular momentum, or number state. Photons may be hyperentangled when the entanglement is encoded in several degrees of freedom all mentioned above [29,30].

Entangled photons have been used recently in nonlinear spectroscopic applications, including two-photon absorption (TPA) [31,32,53], sum frequency generation (SFG) [3335], two-photon-induced fluorescence (TPIF) [3640], and two-photon-induced transparency (TPIT) [41]. The simplest signature of entanglement is that TPA signals scale linearly rather than quadratically with the pump intensity [42]. This reflects the fact that in these experiments, the entangled photons act as a single object. This effect has been predicted and demonstrated in molecular systems. As a consequence, power levels required for TPA excitation can be dramatically reduced by many orders of magnitude, becoming closer to those for single-photon excitation [36]. This could be used to minimize damage in two-photon microscopy. The spatial resolution can also be improved since TPA only takes place in regions where correlated photon pairs overlap in space and time [43,44].

In a recent study we calculated pump-probe (PP) signals using polarization entangled photons generated by type II PDC and postselected by the Mach-Zehnder interferometer. We showed that only the TPA pathways scale linearly with pump intensity and at low pump intensity dominate the other, stimulated Raman (SR) contributions, which scale quadratically. Thus the entangled photons enable quantum pathway selectivity and can help disentangle complex spectra with overlapping features [45]. Photon-entanglement-induced strong correlations of quantum pathways of matter have also been predicted for the difference frequency generation (DFG) technique [46].

Here we focus on a different type of entangled two-photon states known as “twin photons” [40,41,47,54]. The twins are generated by PDC of a single pump photon and correlated both in time and space. Their TPA which occurs within the entanglement time T12 and area A12 is strongly enhanced since it does not rely on the accidental simultaneous arrival of two uncorrelated photons. TPA with twin photons had also been used to detect off-resonant states. The sample becomes transparent for certain values of the entanglement time T12, thereby revealing off-resonant state frequencies. Alternatively, one can fix the entanglement time and introduce a time delay |τ|≤T12 between the twin photons. The Fourier transform of the TPA signal with respect to τ can also reveal the structure of the virtual states. Since the delay can be varied with 0.25 fs steps, the technique can monitor intermediate-state levels separated by ~eV and may therefore be used for atomic spectroscopy. Current theory is limited to the regime where the entangled time is shorter than the dephasing time [32]. We develop a microscopic theory of nonlinear optical signals with twin photons using close time path loop (CTPL) diagrams [48,49] and apply it to predict pump-probe signals in a three-level ladder |gright angle bracket, |eright angle bracket, |fright angle bracket (Fig. 1). We focus on resonances in PP spectroscopy where the sum of the pump and probe frequencies is resonant with one of the molecular states |fright angle bracket i.e., ω12fg. Difference frequency ω1−ω2 (stimulated Raman) type resonances are interesting as well, but will not be considered here. In Sec. II we show that the signals for general quantum optical fields are given by sums over products of matter and light quantum pathways. The pathways can be divided into two groups. The first, representing TPA, is expressed as a modulus square of a two-photon scattering amplitude (Tmatrix)~Tfg(T12,τ)2. It may be calculated as two-photon counting rate using Glauber's theory [23,39,50]. The second component [pump modulated probe (PMP)] represents a process where the pump acts as a “catalyst.” It modulates the probe absorption but returns to its initial state at the end. This group is given by the imaginary part of an amplitude product ~TgfTfg (not modulus square) and cannot be calculated as a photon counting rate. The TPA and PMP can be measured separately by fluorescence from the doubly excited |fright angle bracket and the singly excited states |eright angle bracket, respectively. Each group has a characteristic resonances pattern. In Sec. III we show that the TPA and PMP contributions to the signal are identical when the fields are classical and the intermediate state |eright angle bracket is off resonant. In Sec. IV we study the variation of the signal with the entanglement time T12. We show that the TPA signal with twin photons scales linearly with the field intensity and is inversely proportional to the entanglement time. In Sec. V we introduce the time delay τ between the twins and present 2D correlation spectra obtained by a double Fourier transform with respect to the entanglement time T12 and τ. Various groups of resonances can be clearly separated by these spectra which contain a rich pattern of cross peaks. In Sec. VI we discuss the connection with 2D signals obtained with shaped classical pulses. We show that two mutually delayed rectangularly shaped optical fields may be used to reproduce the TPA spectra obtained with twins. A general discussion is finally given in Sec. VII.

FIG. 1
(Color online)(a) Pump-probe experiment with the twin photons. (b) One of the possible schemes for controlling the external time delay τ between the twins [40]. (c) The twin photon transition amplitude in Eq. (22) is nonzero in the diagonal strip ...

II. CTPL FORMALISM OF PUMP-PROBE SIGNALS

We consider an assembly of N three-level |gright angle bracket, |eright angle bracket, fright angle bracket molecules [Fig. 1(d)] which are initially in the ground state |gright angle bracket. We focus on one component of a general PP signal obtained with two optical modes k1 and k2 where the sum of the two-photon frequencies coincides with one of the transitions ω12gf. The stimulated Raman part with ω1−ω2 resonances is neglected. The frequency domain (cw) PP signal with classical modes is defined as the change in the transmitted intensity of the probe with and without the pump [48]. We treat both beams on the same footing as shown in Fig. 1(a). One of the detectors (d2) measures the change in intensity of the “probe” mode E2 with and without the sample. The second detector (d1) does the same for the “pump” mode E1. We thus define two signals:

SPPd2(ω2;ω1)=ImNTTTdt4E2(t4)P(3)(k2,t4),
(1)

SPPd1(ω1;ω2)=ImNTTTdt4E1(t4)P(3)(k1,t4).
(2)

Here P(3)(k, t) is the third-order polarization and 2T is the measurement time. The nonlinear polarization is induced by sequence of interactions with the optical fields. When the last interaction occurs with k2 (k1) mode it gives P(3)(k2, t4) [P(3)(k1, t4)].

We further define the following symmetric and asymmetric combinations of the two signals:

SPPSym(ω1;ω2)=12[SPPd1(ω1;ω2)+SPPd2(ω2;ω1)],

SPPAsym(ω1;ω2)=12[SPPd1(ω2;ω1)SPPd2(ω1;ω2)].
(3)

In the rotating wave approximation the matter field coupling is given by

Hint=V(t)E(t)+V(t)E(t).

Here the dipole operator is partitioned as V(t)+V(t), where

V(t)=jk>jμjkeiωjgtjkeiωkgt
(4)

is the positive frequency (de-excitation) part; {|kright angle bracket} represents the eigenstates of the bare material system labeled according to increasing energy. V represents the excitation part of the dipole operator.

The field of mode j=1,2 is similarly partitioned as Ej(t)+Ej(t), where the positive frequency component is

Ej(t)=(2πωjΩ)12ajexp(iωjt).
(5)

Here, aj (aj) is the annihilation (creation) boson operator for the j mode and Ω is the quantization volume.

Using the superoperator nonequilibrium Green's function (SNGF) formalism [48,51] and the loop diagrams of Fig. 2 we can express Eq. (1) in terms of four-point correlation function of the matter and the field and obtain

SPPd2(ω2;ω1)=2NTTTdt4Imi3dt3dt2dt1[θ(t4t2)θ(t2t3)θ(t3t1)V(t1)V(t3)V(t2)V(t4)×E1(t1)E2(t3)E1(t2)E2(t4)+θ(t4t2)θ(t2t1)θ(t1t3)V(t3)V(t1)V(t2)V(t4)E2(t3)E1(t1)E1(t2)E2(t4)+θ(t4t2)θ(t4t3)θ(t3t1)V(t1)V(t3)V(t4)V(t2)E1(t1)E2(t3)E2(t4)E1(t2)+θ(t4t2)θ(t4t1)θ(t1t3)V(t3)V(t1)V(t4)V(t2)E2(t3)E1(t1)E2(t4)E1(t2)].
(6)

Here the step functions θ(t−τ) maintain the time orderings of interactions along the loop (but not in real time). The interactions are time ordered within each branch of the diagram, but the two branches are not time ordered with respect to each other. By convention, the chronologically last interaction is an absorption and occurs at the left branch of the loop. The four terms in Eq. (6) correspond to diagrams i–iv of Fig. 2. SPPd1(ω1;ω2) is given by diagrams v–vii and can be obtained from SPPd2(ω1;ω2) [Eq. (6)] by interchanging {ω1ω2}. To proceed further we must calculate the optical field correlation functions by specifying its initial state.

FIG. 2
CTPL diagrams for the contributions to the pump-probe signal which contain ω12fg resonances. By convention, the last interaction is an absorption and occurs at the left branch of the loop. Diagrams i–iv represent k ...

III. PUMP-PROBE SPECTROSCOPY WITH CLASSICAL OPTICAL FIELDS

To set the stage and introduce the relevant building blocks of the signal we shall first consider conventional frequency domain PP performed with classical modes. In this case all optical field correlation functions in Eq. (6) may be factorized into the same product of the field intensities: left angle bracketEEEEright angle bracket=|ε1|22|2. The time arguments in the optical field correlation functions transformed to the loop times (the arguments of the step functions) provide the Fourier components of the material correlation functions and yield the signal

SPPd1(ω1;ω2)=4Nε12ε22Imχ(3)(ω1;ω1,ω2,ω2),

SPPd2(ω2;ω1)=4Nε12ε22Imχ(3)(ω2;ω2,ω1,ω1),

where the explicit CTPL expression for nonlinear susceptibility χ(3) is given in Refs. [48,49]. For comparison with the entangled photon generated PP signal we recast the above expressions in terms of elements of the symmetric two-photon scattering T matrix,

T(ω1,ω2)=(V+V)[G(ω1+ωg)+G(ω2+ωg)](V+V),
(7)

where G(ω) is the Fourier transform of the retarded Green's function operator G(τ)=θ(τ)exp(−iHτ),

G(ω)=idτG(τ)eiωτ,

where H is the matter Hamiltonian. We shall need the following matrix elements of the two-photon operator (two-photon transition amplitudes):

fT(ω1,ω2)g=Tfg(ω1,ω2)=μfeIeg(ω1)μeg+μfeIeg(ω2)μeg,

fT(ω1,ω2)g=Tfg(ω1,ω2)=μfeIeg(ω1)μeg+μfeIeg(ω2)μeg,

gT(ω1,ω2)f=Tgf(ω1,ω2)=μgeIge(ω1)μef+μgeIge(ω2)μef,

gT(ω1,ω2)f=Tgf(ω1,ω2)=μgeIge(ω1)μef+μgeIge(ω2)μef.
(8)

Here we have introduced the auxiliary quantity

Ieg(ω)=Gee(ω+ωg)=1ωωeg+iγ,
(9)

where Gee(τ)=θ(τ)exp(−iωeτ− γτ), γ is a dephasing rate, and exp(iωgτ) describes the free evolution of the molecule in the ground state |gright angle bracket.

Tfg(ω1,ω2) describes retarded TPA on the left branches of the loop diagrams in Fig. 2 where the loop time runs forward. Tfg(ω1,ω2) represents advanced TPA on the right branches where the loop time runs backward. Similarly, Tgf(ω1,ω2)[Tgf(ω1,ω2)] describes retarded (advanced) TP emission. The system evolution after the absorption of the two photons is described by the doubly excited state Green's function Gff=θ(τ)exp(−iωfτ) and we define

Ifg(ω1+ω2)=Gff(ω1+ω2+ωg)=1ω1+ω2ωfg+iγ.
(10)

Using the loop diagram rules [48] and Eq. (9) we obtain

SPPd2(ω2;ω1)~Imi6eeμgeμeμefμfeIfg(ω1+ω2)×[Ieg(ω2)Ieg(ω1)+Ieg(ω2)Ieg(ω2)+Ieg(ω1)Ieg(ω1)+Ieg(ω1)Ieg(ω2)].

The four terms correspond to pathways i–iv in Fig. 2, respectively. The other signal SPPd1(ω1;ω2) can be obtained from the above equation by interchanging ω1 and ω2.

The symmetric signal (3) assumes a compact form in terms of the two-photon transition amplitudes (8),

SPPSym(ω1,ω2)=SPMP(ω1;ω2)+STPA(ω1;ω2)2Nε12ε22ImIfg(ω1+ω2)×[Tgf(ω1,ω2)Tfg(ω1,ω2)+Tfg(ω1,ω2)2].
(11)

The two terms in the brackets represent the PMP and TPA groups of pathways, respectively, as shown in Fig. 2. To interpret these terms we recall that PP is accompanied by spontaneous light-emission process. Each partially time-ordered CTPL diagram may be dissected into a set of fully time-ordered Feynman diagrams, where chronologically the last interaction occurs at time t4 [52]. In the TPA pathways the molecule ends up in the doubly excited state |fright angle bracketleft angle bracketf| and two photons are absorbed. The fluorescence from this state can be detected and described by the TP counting [50] formalism ~Tfg2. The fundamental difference between the single-photon absorption and the PMP is the involvement of the second photon. The state of the second mode at time t4 is the same as its initial state, but it provides intermediate coherence between the ground and double excited states, i.e., catalyzes the absorption of the other mode. The PMP pathways are related to fluorescence from the intermediate states |eright angle bracketleft angle brackete|. Since only one of the photons is absorbed, their contribution is not described by the single-photon transition amplitude as in the photon counting formalism but by ~Tgf(ω1,ω2)Tfg(ω1,ω2).

The asymmetric signal vanishes in the degenerate case (ω12) where the two detectors give the same PP signal. When the intermediate states |eright angle bracket are off resonant one has Tfg(ω1,ω2)=Tfg(ω1,ω2). Equation (11) then assumes the two-photon counting form

SPPSym(ω1,ω2)=4Nε12ε22Tgf(ω1,ω2)2δ(ω1+ω2ωfg).

However, as will be shown in Sec. IV, these results do not hold for quantum optical fields, where the two photons are correlated. Additional terms coming from the PMP group of pathways contribute to all PP signals. As a result, spontaneous emission from the final |fright angle bracket and the intermediate states |eright angle bracket become comparable.

IV. PUMP-PROBE WITH TWIN PHOTONS

We now assume that modes k1 and k2 are initially prepared in a twin state [41] representing an entangled photon pair. This can be produced by spontaneous parametric down conversion in a birefringent second-order nonlinear crystal χ(2). The photon pair is described by a nonseparable wave function even when the photons are physically separated [47],

ψ(2)(0)=Ck1k2tIexp(iΔωtI2)sinc(ΔωtI2)×exp(i(Δk)zLz2)sinc((Δk)zLz2)k1,k2.
(12)

The Fock state |k1,k2| has one photon in each mode k1 and k2; Δω=ωp−ω1−ω2k=kpk1k2; tI and Lz are the interaction time and length within the PDC crystal. The normalization factor C=χ(2)εpA12 is proportional to the nonlinearity of the PDC crystal, the pump electric field amplitude εp, and the entanglement area A12 [24]. The detectors d1 and d2 can distinguish between the two modes by their polarization (vertical for k1 and horizontal for k2).

For the twin state each of the optical field correlation functions left angle bracketEEEEright angle bracket in Eq. (6) is factorized into a product of two field transition amplitudes [47],

0,0E2(t4)E1(t2)ψ(2)(0)=θ(t4t2)F(t4,t2,T12)+θ(t2t4)F(t2,t4,T12)×ψ(2)(0)E1(t3)E2,(t1)0,0=θ(t3t1)F(t3,t1,T12)+θ(t1t3)F(t1,t3,T12),
(13)

where

F(t4,t2,T12)=2πχ(2)εpΩω1ω2A12T12×exp[i(ω1t2+ω2t4)]rect(t4t2T12),
(14)

rect(t) is the rectangular function equal to 1 for 0≤t≤1 and 0 otherwise, and T12=(1/v1−1/v2)Lz is the entanglement time determined by the group velocities v1(v2) of the mode k1(k2) within the PDC crystal.

By substituting Eq. (13) into Eq. (6), we get

SPPd2(ω2;ω1,T12)=4N2TTTdt4Imi3dt3dt2dt1×[θ(t4t2)θ(t2t3)θ(t3t1)V(t1)V(t3)V(t2)V(t4)F(t3,t1,T12)F(t4,t2,T12)+θ(t4t2)θ(t2t1)θ(t1t3)V(t3)V(t1)V(t2)V(t4)F(t1,t3,T12)F(t4,t2,T12)+θ(t4t2)θ(t4t3)θ(t3t1)V(t1)V(t3)V(t4)V(t2)F(t3,t1,T12)F(t2,t4,T12)+θ(t4t2)θ(t4t1)θ(t1t3)V(t3)V(t1)V(t4)V(t2)F(t1,t3,T12)F(t2,t4,T12)].
(15)

We now introduce the Fourier transform of the joined matter and entangled photon Green's functions,

Ieg(ω,T12)=iGee(τ)rect(τT12)exp[i(ωωg)τ]dτ=ei(ωωeg+iγ)T121ωωeg+iγ.
(16)

By substituting Eqs. (4) and (13) into Eq. (15) and changing the time variables to the loop times, we can express the TPA signal in terms of the entangled photons and matter Green's functions as follows:

SPPd2(ω2;ω1,T12)=16π2Nχ(2)2εp2ω1ω2A12T12Ω2Imf,e,eμegμegμefμefIfg(ω1+ω2)[Ieg(ω2,T12)Ieg(ω1,T12)+Ieg(ω2,T12)Ieg(ω2,T12)+Ieg(ω1,T12)Ieg(ω1,T12)+Ieg(ω1,T12)Ieg(ω2,T12)],
(17)

where Ifg12) is the material doubly excited state Green's function (10). The signal SPPd1(ω1;ω2,T12) can be obtained from Eq. (17) by interchanging ω1 and ω2.

Introducing the entanglement modified two-photon transition amplitude

Tgf(ω1,ω2,T12)=μgeμefIeg(ω1,T12)+μgeμefIeg(ω2,T12),

Tfg(ω1,ω2,T12)=μgeμefIeg(ω1,T12)+μgeμefIeg(ω2,T12),
(18)

we obtain the symmetric TPA signal induced by the twin photons,

SPPSym(ω1,ω2,T12)=8π2Nχ(2)2εp2ω1ω2A12T12Ω2×ImIfg(ω1+ω2)[Tgf(ω1,ω2,T12)×Tfg(ω1,ω2,T12)+Tfg(ω1,ω2,T12)2].
(19)

As in Eq. (11), the two terms correspond to PMP and TPA pathways. It is interesting to note that the symmetric PP signal generated by twin photons (19) may not be described by the absolute square of the transition amplitude.

We have carried out simulations for the following model material system. The single excited state frequencies ωeg={1.5, 1.6, 1.8} are off resonant (ωeg−ωp/2)T12[dbl greater than]γT12. The doubly excited state at ωfg=1 is resonant with the pump photon ωgfp used to generate the degenerate twins ω12p/2. All dipole moments are real and identical.

We have varied the entanglement time T12 from Tmin=212ΔT12 to Tmax=2Tmin with the time stem ΔT12=2π/[max(ωegeg−ωp)]=0.025. That is if the unit of energy is ωp=1 eV then Δτ=1.3 fs. This can be done by tilting the PDC crystal [36] which increases its size Lz.

In Fig. 3 we display the entanglement time spectra defined as follows:

SPPν(ωp2,ωT)=NTminTmaxdT12exp(iωT12T12)SPPν(ωp2,T12),
(20)

where superscript ν can be {d1,d2,Sym,TPA,PMP} and N=1STPASym(ωp2,T12) is a normalization factor. The spectra show several groups of resonances. Group a correspond to the intraband resonances of the intermediate-state manifold 0≤ωT12eg−ωeg≤0.3, group b corresponds to the 1≤ωT12eg−ωp/2≤1.3 resonances, and group c describes 2≤ωT12egp≤2.6 resonances. Note that the resonances of group b can be detected by a conventional pump probe with short well-separated classical pulses.

FIG. 3
(Color online) The pump-probe spectra log10STPAν(ωp2,ωT) [Eq. (25)]. The signals recorded by the two detectors are identical Sd1=Sd2=SSym. Left (right) column correspond to large (small) dephasing rate, ...

All pathways (Fig. 2) which contribute to the signals interfere constructively. The left column in Fig. 3 uses a larger dephasing rate 2.68<T12γ<5.36 than the right column 0.03<γT12<0.05. The larger dephasing rate in the left column quenches the resonances in regions a and c. PMP pathways contribute to spectral regions b and c, while TPA pathways show resonances in regions i and ii. We note that for this technique Sd1=Sd2=SSym and the asymmetric signal SAsym vanishes. In both columns the STPA signal has no resonances in region c whereas the SPMP signal does not show resonances in region a.

V. PUMP-PROBE WITH MUTUALLY DELAYED TWIN PHOTONS

We next turn to a more general twin photon PP setup as proposed by Saleh et al. [32]. The twins remain degenerate ω12p/2 but we introduce a relative time delay 0<τ<T12 between them by shifting the phase of the k1 mode [12,40],

E1(t)=(πωpΩ)12a1exp[iωp(tτ2)2],

E2(t)=(πωpΩ)12a2exp[iωp(t+τ2)2].
(21)

Using Eq. (21) we obtain the optical field transition amplitudes for diagrams i–iv,

0,0E2(t4)E1(t2)ψ(2)(0)=θ(t4t2)F(t4,t2,T12)+θ(t2t4)F(t2,t4,T12),

ψ(2)(0)E1(t3)E2(t1)0,0=θ(t3t1)F(t3,t1,T12)+θ(t1t3)F(t1,t3,T12),
(22)

where

F(t4,t2,T12)=πχ(2)εpωpΩA12T12exp[iωp2(t4+t2)]rect(t4t2T12τ).

For diagrams v–viii one has to interchange t4t2,t3t1 in the above expression.

The signal in Eq. (15) now becomes

SPPd2(ωp2,T12,τ)=16π2Nχ(2)2εp2ωp2A12T12Ω2Imf,e,eμegμegμefμefIfg(ωp)[Ieg(ωp2,T12τ)Ieg(ωp2,T12τ)+Ieg(ωp2,T12τ)Ieg(ωp2,T12+τ)+Ieg(ωp2,T12τ)Ieg(ωp2,T12τ)+Ieg(ωp2,T12τ)Ieg(ωp2,T12+τ)].
(23)

The signal from the other detector is SPPd1(ωp2,T12,τ)=SPPd2(ωp2,T12,τ). Here Ieg(ωp2,T12τ) describes the Fourier transform of the system (twin+matter) Green's function which depends on the delay τ between the entangled photons. This delay breaks the symmetry of the process making it dependent on the order in which the twins are absorbed or emitted even in the degenerate case [see Fig. 1(c)]. If mode k1 precedes k2, the system's evolution is described by the following Green's function: Ieg(ω,T12−τ). When k2 comes first we have Ieg(ω,T12+τ). The broken symmetry implies that SPPd2(ωp2,T12,τ)SPPd1(ωp,2T12,τ).

Using the field transition amplitudes (22) the two-photon transition amplitude assumes the form

Tgf(ωp2,T12,τ)=eμgeμefIeg(ωp2,T12τ)+μgeμefIeg(ωp2,T12+τ),

Tfg(ωp2,T12,τ)=eμgeμefIeg(ωp2,T12τ)+μgeμefIeg(ωp2,T12+τ).
(24)

The symmetrized signal (3) can be written in terms of the matter field transition amplitudes (24) as

SPPSym(ωp2,T12,τ)=12[SPPPMP(ωp2,T12,τ)+STAPSym(ωp2,T12,τ)],
(25)

where

SPPPMP(ωp2,T12,τ)=16π2Nχ(2)2εp2ωp2A12T12Ω2×ImfIfg(ωp)Tfg(ωp2,T12,τ)×Tgf(ωp2,T12,τ),
(26)

SPPTPA(ωp2,T12,τ)=16π2Nχ(2)2εp2ωp2A12T12Ω2×fTfg(ωp2,T12,τ)2δ(ωpωfg).
(27)

Equation (27) is identical to the result of Ref. [32]. For the degenerate case ω12, the symmetry of the above signal implies the following symmetry with respect to the time delay SPPSym(ωp2,T12,τ)=SPPSym(ωp,2T12,τ).

We next turn to the antisymmetric signal [the difference between the detector d1 and d2 in Fig. 1(a)],

SPPAsym(ωp2,T12,τ)=8π2Nχ(2)2εp2ωp2A12T12Ω2ImfIfg(ωp)×[Tgf(ωp2,T12,τ)Tgf(ωp2,T12,τ)+Tgf(ωp2,T12,τ)Tfg(ωp2,T12,τ)].
(28)

Here we have introduced the antisymmetric two-photon transition amplitude

Tgf(ωp2,T12,τ)=eμgeμefIeg(ωp2,T12τ)μgeμefIeg(ωp2,T12+τ).
(29)

This signal vanishes in frequency-domain experiments with classical fields where the order in which the modes interact is not important,

Tgf(ω1,ω2)=μgeμefIeg(ω1,ω2)μgeμefIeg(ω1,ω2)=0.

The asymmetric signal also vanishes at zero twin photon delay since Tgf(ωp2,T12,0)=0. Absorption of k1 is always accompanied by absorption of k2 and the same holds for emission. For the twin photons with delay τ the Green's functions do depend on the order of interaction and the asymmetric signal does not vanish. This signal is thus induced by the delay between twin photons.

In Fig. 4 we compare the individual detector readings (23), the symmetrized signal (25), and the TPA contribution (27) to the PP signal. First, we hold the entanglement time T12 fixed and plot the Fourier transforms with respect to the time delay τ:

SPPν(ωp2,T12,ωτ)=N0T12dτexp(iωττ)SPPν(ωp2,T12,τ).
(30)

The entanglement time increases from zero with a step ΔτΔT12 up to the value of the entanglement time T12=Tmin. As in Fig. 3 the resonances of the entanglement time PP spectra may be classified by the same three spectral groups a,b,c. Recall that the spectra with respect to the entanglement time associate the region a with PMP and region c with TPA. However their meanings for the spectra given by Eq. (30) is quite different. For PMP pathways region a resonances are produced by the pathways i, v, iii, and vii where emitted and absorbed photons follow the same chronological order; the region b resonances are given by the pathways ii, vi, iv, and viii where the chronological photon order for the emission is opposite to those of the absorption. For TPA the meaning of the regions of a and c is exactly opposite to those of PMP. This peculiar feature is a result of broken symmetry of the field transition amplitude. For small dephasing rate γT12=0.03 (right column) PMP and TPA contribute to the same regions in all spectra and spectrally nonseparable by the pump-probe technique. At large dephasing rate γT12=2.68 (left column) pathways v–viii do not contribute to regions a and c. PMP pathways i and ii contribute to the spectral regions b and c while TPA pathways iii and iv contribute to regions a and b. We note that the spectra obtained from the two detectors are different Sd1Sd2. Only the Sd2 signal shows resonances in the nonclassical regions a and c. However, both SSym and SAsym signals contain all resonances. The TPA signal misses the region c resonances and PMP misses resonances in region a; thus, one can separate PMP from TPA if focuses on the regions a and c only.

FIG. 4
(Color online) The pump-probe spectra log10STPAν(ωp2,T12,ωτ) vs ωτ and a fixed entanglement time T12. The two columns differ by the dephasing rate, as indicated.

In Fig. 5 we display correlation 2D spectra between the entanglement and the delay times,

STPAν(ωp2,ωT12,ωτ)=NTminTmaxdT12exp(iωT12T12)×0T12dτexp(iωττ)STPAν(ωp2,T12,τ).
(31)

The resonances are classified into five regions in the (ωT12τ) plane. The PMP pathways contribute to regions (b,b), (c,a), and (c,c); the TPA pathways create resonances at (a,a), (b,b), and (a,c). Therefore regions a,c and pathways c,a directly show the correlation between the PMP and TPA. We find that for both values of the dephasing rate the TPA and PMP signals miss the (c,c) and (a,a) resonances correspondingly. The correlation spectra from the two detectors are different. Unlike the delay spectra, detector d1 reveals all spectral regions, while resonances (a,a) and (c,c) from detector d2 are suppressed. Finally the off-diagonal, as well as (b,b), resonances are sensitive to the dephasing rate and vanish for large dephasing. The diagonal (a,a) and (c,c) resonances overlap as the dephasing rate is increased. In summary, the conventional frequency-domain PP signal when ω1 and ω2 are off resonant can be described as a two-photon counting process. The twin photons characterized by their frequencies, polarizations, and the entanglement time T12 generate the symmetric PP signal [Eq. (19)] obtained from CTPL formalism. The signal consists of TPA~Tfg(T12)2 as well as PMP~Tgf(T12)Tfg(T12) terms. The two terms may be observed separately by measuring two-photon-induced fluorescence from the doubly and singly excited states, respectively. Alternatively, they can be spectrally separated by focusing on certain regions in SPPSym(ωp2,ωT12). When the time delay τ is introduced the two detectors yield two different signals SPPd1,2(ωp2,ωτ) giving rise to the asymmetric signal SPPAsym(ωp2,ωτ) or spectra. For SPPSym(ωp2,ωτ) the separation between PMP and TPA is possible only at large dephasing. The correlation between PMP and TPA pathways can be seen in the off-diagonal regions on the 2D correlation spectra SPPν(ωp2,ωT12,ωτ). Spectral separation and correlation between PMP and TPA is important if the direct observation of the fluorescence from |eright angle bracket or |fright angle bracket is not possible due to selection rules or fast nonradiative relaxation of the |fright angle bracket state.

FIG. 5
(Color online) 2D correlation pump-probe spectra log10STPAν(ωp2,ωT12,ωτ) with entangled twin photons. The left and right columns represent different dephasing rates.

VI. PUMP-PROBE WITH RECTANGULARLY SHAPED CLASSICAL OPTICAL FIELDS

We now demonstrate how similar spectra can be obtained with shaped classical pulses. The key ingredients of the TPA signals are the optical field transition amplitudes, which can be shaped and manipulated. Let us assume that modes k1 and k2 are created from a rectangular (width T12) pump beam with frequency ωp by means of 50:50 beam splitter. Both modes are classical and initially in a coherent state |βright angle bracket. The beam-splitter-induced π/2 phase shift between the reflected k1 and transmitted k2 modes can be compensated by difference in their optical path from the splitter to the sample. We further introduce relative time delay τ in k1 beam path. The optical field operators are

E1(t)=(2πωpΩ)12a1exp[iωp(tτ2)2]×rect(tτ2T122T12),

E2(t)=(2πωpΩ)12a2exp[iωp(t+τ2)2]×rect(t+τ2T122T12).
(32)

Making use of the following identity:

rect(t4τ2T122T12)rect(t2+τ2T122T12)=rect(t4t2+τT12),

we obtain the following field transition amplitudes for diagrams i–iv:

β1,β2E2(t4)E1(t2)β1,β2=θ(t4t2)E(t4,t2,T12)+θ(t2t4)E(t2,t4,T12)=2πωpΩexp[iωp2(t4+t2)]×[θ(t4t2)rect(t4t2T12τ)+θ(t2t4)rect(t2t4T12+τ)],

β1,β2E1(t3)E2,(t1)β1,β2=θ(t3t1)E(t3,t1,T12)+θ(t1t3)E(t1,t3,T12)=2πωpΩexp[iωp2(t3+t1)]×[θ(t3t1)rect(t3t1T12τ)+θ(t1t3)rect(t1t3T12+τ)].
(33)

Substituting Eq. (33) into Eq. (15) gives

STPASym(ωp2,T12,τ)=16π2Nωp2Ω32ImIfg(ωp)[Tfg,2(ωp2,T12,τ)+Tfg(ωp,T12,τ)2]
(34)

Equation (34) derived for shaped classical modes coincides with the entanglement-induced TPA signal [Eq. (25)] up to an overall scaling factor. Equation (25) scales linearly with the pump photon intensity, while the classical signal scales quadratically. The scaling factor is independent on the τ parameter. Therefore their Fourier spectra with respect to τ are identical.

Pump-probe spectroscopy performed with shaped classical pulses where the pump is split into two copropagating pulses separated by a delay τ gives information equivalent to four-wave mixing carried out with four pulses. 2D spectra can be obtained by, e.g., varying τ and the T12. The twin delay τ plays a similar role but enters in a different way than with classical fields. The precise connection between the two approaches to multidimensional spectroscopy based controlling either the delays between classical pulses or the parameters of entangled photon wave function that determine their relative timing statistics will be an interesting topic for future investigations.

VII. CONCLUSIONS

In this paper we demonstrated how quantum degrees of freedom of the radiation field may be used as multidimensional spectroscopic probes of matter. Mutually delayed twin photons tuned to two-photon resonances ω12fg can provide direct information on off-resonant transitions. The twin photon state is a wave packet of two groups of modes that can be distinguished and detected separately by two detectors d1 and d2 by exploiting their different polarizations or propagation directions. The two detectors measure the intensity change in one of the beams with and without the sample. An intuitive microscopic formalism based on close time path loop (CTPL) diagrams is developed and employed to calculate the pump-probe signals generated with the entangled photons SPPd1 and SPPd2. Symmetric (asymmetric) signals are further defined as the sum (difference) of the detector readings (SPPd1±SPPd2). For general quantum optical fields the signal is recast in terms of quantum pathways of matter, each scaled with a corresponding four-point field correlation functions. Such signals may not be expressed in terms of ordinary causal response functions and susceptibilities as is the case for classical fields since they depend on correlated spontaneous fluctuations of field and matter. For classical fields the signal is given by a product of a susceptibility and field amplitude. In the entangled case the signal is given by a sum of such products. The symmetrized pump-probe signal can be divided into components which represent two groups of pathways SPPPMP+SPPTPA. TPA may be measured by fluorescence from the doubly excited states |fright angle bracket and may be calculated as a modulus square of a two-photon counting amplitude. PMP on the other hand represents a process where one of the pump photon creates an intermediate coherence between the ground and double excited states and “catalyzes” the probe absorption, without changing its state at the end of the process. This component can be measured by fluorescence from the singly excited |eright angle bracket state and is given by the imaginary part of a product of two-photon complex transition amplitudes (not modulus square). For classical off-resonant fields SPPPMP=SPPTPA, and the antisymmetric signal SPPAsym vanishes. However it becomes finite with twin photons. The distinction between the pump-probe and TPA signals has not been clearly made in the literature, where the two terms are used interchangeably. 2D signals obtained by varying the delay time τ and the entanglement time T12 reveal the off-resonant singly excited state frequencies and allow us to distinguish between various groups of resonances. Focusing on different regions of the spectra one can separate PMP and TPA contributions or measure the correlation between the two. This spectral resolution is especially important if direct measurement of the fluorescence from the |fright angle bracket is impossible. We further showed that one can mimic the entangled-photon pump-probe signals by using two mutually delayed rectangularly shaped classical optical fields. The present formalism allows us to describe pulse shaping and photon entanglement within the same framework. This may allow the design of other types of time-domain spectroscopies which make use of the broad band width of entangled photons to probe relaxation processes in the resonantly excited states.

ACKNOWLEDGMENTS

This work was supported by the National Science Foundation Grant No. CHE-0745892 and the National Institutes of Health Grant No. GM59230.

References

[1] Ernst R, Bodenhausen G, Wokaun A. Principles of Nuclear Magnetic Resonance in One and Two Dimensions. Clanderon Press; Oxford: 1998.
[2] Mukamel S. Annu. Rev. Phys. Chem. 2000;51:691. [PubMed]
[3] Asplund M, Zanni M, Hochstrasser R. Proc. Natl. Acad. Sci. U.S.A. 2000;97:8219. [PubMed]
[4] Mukamel S, Abramavicius D, Yang Lijun, Zhuang Wei, Schweigert IV, Voronine DV. Acc. Chem. Res. (to be published)
[5] Engel GS, Calhoun TR, Read EL, Ahn TK, Mancal T, Cheng YC, Blankenship RE, Fleming GR. Nature (London) 2007;446:782. [PubMed]
[6] Oron D, Dudovich N, Silberberg Y. Phys. Rev. Lett. 2003;90:213902. [PubMed]
[7] Li H, Ahmasi D, Xu HB, Wrzesinski PJ, Lozovoy VV, Dantus M. Opt. Photonics News. 2008;19:46.
[8] Postma S, Van Rhijn AC, Korterik JP, Gross P, Herek JL, Offerhaus HL. Opt. Express. 2008;16:7985. [PubMed]
[9] DeFlores LP, Nicodemus RA, Tokmakoff A. Opt. Lett. 2007;32:2966. [PubMed]
[10] Grumstrup EM, Shim SH, Montgomery MA, Damrauer NH, Zanni MT. Opt. Express. 2007;15:16681. [PubMed]
[11] Myers JA, Lewis KL, Tekavec PF, Ogilvie JP. Opt. Express. 2008;16:17420. [PubMed]
[12] Edamatsu K. Jpn. J. Appl. Phys. 2007;46:7175.
[13] Hong CK, Mandel L. Phys. Rev. A. 1985;31:2409. [PubMed]
[14] Furusawa A, Sørensen JL, Braunstein SL, Fuchs CA, Kimble HJ, Polzik ES. Science. 1998;282:706. [PubMed]
[15] Gisin N, Ribordy G, Tittel W, Zbinden H. Rev. Mod. Phys. 2002;74:145.
[16] Scully M, Zubairy M. Quantum Optics. Cambridge University Press; Cambridge, UK: 1997.
[17] Aspect A, Grangier P, Roger G. Phys. Rev. Lett. 1982;49:91.
[18] Haroche S, Raimond J. Exploring the Quantum: Atoms, Cavities, and Photons. Oxford University Press; New York: 2006.
[19] Kwiat PG, Mattle K, Weinfurter H, Zeilinger A, Sergienko AV, Shih Y. Phys. Rev. Lett. 1995;75:4337. [PubMed]
[20] Fry ES, Thompson RC. Phys. Rev. Lett. 1976;37:465.
[21] Aspect A, Grangier P, Roger G. Phys. Rev. Lett. 1981;47:460.
[22] Yu S, Eberly JH. Science. 2009;323:598. [PubMed]
[23] Mandel L, Wolf E. Optical Coherence and Quantum Optics. Cambridge University Press; Cambridge, England: 1995.
[24] Gerry C, Knight P. Introductory Quantum Optics. Cambridge University Press; Cambridge, England: 2005.
[25] Fiorentino M, Voss PL, Sharping JE, Kumar P. IEEE Photonics Technol. Lett. 2002;14:983.
[26] Inoue K, Shimizu K. Jpn. J. Appl. Phys. 2004;43:8048.
[27] Edamatsu K, Oohata G, Shimizu R, Itoh T. Nature (London) 2004;431:167. [PubMed]
[28] Stevenson RM, Hudson AJ, Young RJ, Atkinson P, Cooper K, Ritchie DA, Shields AJ. Opt. Express. 2007;15:6507. [PubMed]
[29] Barreiro JT, Langford NK, Peters NA, Kwiat PG. Phys. Rev. Lett. 2005;95:260501. [PubMed]
[30] Barbieri M, Cinelli C, Mataloni P, De Martini F. Phys. Rev. A. 2005;72:052110.
[31] Dayan B, Peer A, Friesem AA, Silberberg Y. Phys. Rev. Lett. 2004;93:023005. [PubMed]
[32] Saleh BEA, Jost BM, Fei HB, Teich MC. Phys. Rev. Lett. 1998;80:3483.
[33] Dayan B, Pe'er A, Friesem AA, Silberberg Y. Phys. Rev. Lett. 2005;94:043602. [PubMed]
[34] Pe'er A, Bromberg Y, Dayan B, Silberberg Y, Friesem AA. Opt. Express. 2007;15:8760. [PubMed]
[35] Pe'er A, Dayan B, Friesem AA, Silberberg Y. Phys. Rev. Lett. 2005;94:073601. [PubMed]
[36] Harpham MR, Suzer O, Ma C-Q, Bauerle P, Good-son T. J. Am. Chem. Soc. 2009;131:973. [PubMed]
[37] Lee D, Goodson T., III Proc. SPIE. 2007;6653:66530V.
[38] Lee D, Goodson T. J. Phys. Chem. B. 2006;110:25582. [PubMed]
[39] Süzer Ö, Goodson T. Opt. Express. 2008;16:20166. [PubMed]
[40] Teich M, Saleh B. U.S. Patent. 5,796,477 Aug 18, 1998.
[41] Fei HB, Jost BM, Popescu S, Saleh BEA, Teich MC. Phys. Rev. Lett. 1997;78:1679.
[42] Javanainen J, Gould PL. Phys. Rev. A. 1990;41:5088. [PubMed]
[43] Abouraddy AF, Saleh BEA, Sergienko AV, Teich MC. Phys. Rev. Lett. 2001;87:123602. [PubMed]
[44] Nasr MB, Saleh BEA, Sergienko AV, Teich MC. Phys. Rev. Lett. 2003;91:083601. [PubMed]
[45] Roslyak O, Marx C, Mukamel S. Phys. Rev. A. 2009;79:033832. [PMC free article] [PubMed]
[46] Roslyak O, Mukamel S. Opt. Express. 2009;17:1093. [PMC free article] [PubMed]
[47] Peřina J, Jr., Saleh BEA, Teich MC. Phys. Rev. A. 1998;57:3972.
[48] Marx CA, Harbola U, Mukamel S. Phys. Rev. A. 2008;77:022110. [PMC free article] [PubMed]
[49] Roslyak O, Marx C, Mukamel S. Mol. Phys. (to be published)
[50] Glauber R. Quantum Theory of Optical Coherence: Selected Papers and Lectures. Wiley-VCH; New York: 2007.
[51] Harbola U, Mukamel S. Phys. Rep. 2008;465:191.
[52] Mukamel S. Principles of Nonlinear Optical Spectroscopy. Oxford University Press; New York: 1995.
[53] Georgiades N.Ph., Polzik ES, Kimble HJ. Phys. Rev. A. 1999;59:676.
[54] Smithey DT, Beck M, Belsley M, Raymer HG. Phys. Rev. Lett. 1992;69:2650. [PubMed]