Home | About | Journals | Submit | Contact Us | Français |

**|**HHS Author Manuscripts**|**PMC2964889

Formats

Article sections

- Abstract
- I. INTRODUCTION
- II. INCOHERENT OPTICAL SIGNALS
- III. DIAGRAMMATIC REPRESENTATION FOR SPONTANEOUS LIGHT EMISSION (SLE)
- IV. THE PUMP-PROBE SIGNAL
- V. COHERENT VS. INCOHERENT NONLINEAR OPTICAL PROCESSES
- VI. HETERODYNE-DETECTION AS A STIMULATED WAVE MIXING
- VII. CONCLUSIONS
- References

Authors

Related links

Phys Rev A. Author manuscript; available in PMC 2010 October 27.

Published in final edited form as:

PMCID: PMC2964889

NIHMSID: NIHMS212364

Department of Chemistry, University of California, Irvine, CA 92697

The publisher's final edited version of this article is available at Phys Rev A

See other articles in PMC that cite the published article.

Nonlinear optical signals from an assembly of *N* noninteracting particles consist of an incoherent and a coherent component, whose magnitudes scale ~ *N* and ~ *N*(*N* − 1), respectively. A unified microscopic description of both types of signals is developed using a quantum electrodynamical (QED) treatment of the optical fields. Closed nonequilibrium Green’s function expressions are derived that incorporate both stimulated and spontaneous processes. General (*n* + 1)-wave mixing experiments are discussed as an example of spontaneously generated signals. When performed on a single particle, such signals cannot be expressed in terms of the *n*th order polarization, as predicted by the semiclassical theory. Stimulated processes are shown to be purely incoherent in nature. Within the QED framework, heterodyne-detected wave mixing signals are simply viewed as *incoherent stimulated* emission, whereas homodyne signals are generated by *coherent spontaneous* emission.

Nonlinear optical processes in bulk materials are traditionally classified as either coherent or incoherent [1]–[4]. This classification refers to the relation of the macroscopic signal intensity of the sample to the nonlinear response of the individual molecules. Processes where the detected *signal* is obtained by simply adding up the contributions of the individual particles in the sample, are termed incoherent. Such signals hence scale as the number of particles *N* in the interaction volume. With coherent processes, on the other hand, the contributions of the constituent particles add up on the level of *amplitudes* rather than intensities, and the detected signal is proportional to *N*^{2}.

Coherent optical signals are usually calculated semiclassically by adopting a two step procedure [2] – [4]. First a (linear or nonlinear) polarization *P*^{(}^{n}^{)} induced in the system by interactions with classical external fields is calculated microscopically. In the second step, *P*^{(}^{n}^{)} is used as a source in the macroscopic Maxwell equations to generate the signal. Frequency-domain signals are given by the absolute square of an amplitude related to the susceptibilities *χ*^{(}^{n}^{)} which, in turn, are proportional to *N*. This gives rise to the ~ *N*^{2} dependence of coherent signals [5],[6]. Examples are Rayleigh scattering (*n* = 1), sum frequency generation (*n* = 2), four-wave mixing (*n* = 3). In contrast, signals such as Raman, fluorescence, two photon fluorescence etc. are incoherent and may not be calculated by the above semiclassical approach. Instead, they require a quantum description of the field to account for spontaneous emission.

Treating both types of processes by a unified approach is of fundamental interest [1], particularly when they coexist. For example, two photon fluorescence and second harmonic generation have been observed in the same sample [7]. There is a considerable current effort to perform nonlinear optical measurements on small samples and even single nano particles or molecules [8]–[11]. Recent experiments include 4-wave mixing spectroscopy of single semiconductor quantum dots [12] or nonlinear optical measurements on single electrons in F centers[13]. CARS microscopy [14] is being extended to small interaction volumes. The traditional macroscopic formulation does not apply for small samples and a quantum electrodynamical (QED) formulation is called for.

A QED treatment is essential for predicting observables related to the quantum nature of the field. A powerful quantum master equation approach has been successfully used in quantum optics to describe fluctuations of the laser field, photon statistics, entanglement and squeezing [15]–[25]. In this paper we consider simple observables that are usually calculated by the semiclassical approach. However a quantum description of the fields allows a more precise definition of the signal, provides new insights and shows the limitations of the semiclassical approach. Our final results are expressed in terms of multipoint correlation functions of the material system. These can be evaluated for more complicated models than the Bloch-type equations ordinarily used in quantum optics. Details of molecular complexity, coupling to a solvent with arbitrary time scale and excitonic effects may be readily described by these expressions [3].

The present approach is based on the many-body non-equilibrium Green’s functions (NEGF) technique [26] – [32]. It provides a single microscopic definition of the signal that includes both spontaneous and stimulated emission and generates both incoherent and coherent signals. We shall use a superoperator in Liouville space representation that is commonly applied to study electrical conduction in quantum junctions [33]–[35]. In Sec. II we give a general definition of the signal (Eq. (11)) and express it in terms of Green’s functions of the system and the field. The formalism will be demonstrated by first considering the signal generated by a single molecule. Applications are made to two types of incoherent experiments, spontaneous light emission (SLE) (Sec. III) and to pump-probe spectroscopy (Sec. IV). A purely diagrammatic derivation of the final expressions for the SLE signal using the technique of Keldysh Schwinger loops [30] is developed. We present rules that directly translate the diagrams to partially-time-ordered expressions for the signal, both in the time-and in the frequency-domain. This representation, has been shown to be most useful for the study of frequency-domain experiments within the semiclassical approximation [36]. We extend it to quantum fields, and show that it yields more compact expressions than the fully-time-ordered double-sided Feynman diagrams.

In Sec. V, we turn to the response of an assembly of noninteracting particles. The expressions of Sec. II are generalized by simply replacing the dipole operator of a single molecule by that of an assembly of particles. The signal is separated into an incoherent and a coherent part. The latter, which is described by semiclassical optical susceptibilities, is dominant in macroscopic samples. However for small *N* both have comparable magnitudes and coexist. The coherent signal does not exist in the limit of a single molecule. This indicates a breakdown of the semiclassical theory which does predict such a signal. Finally, in Sec. VI we show that within QED heterodyne detection emerges as an incoherent stimulated emission process in the detection mode. This is in contrast with semiclassical treatment where we first generate a coherent macroscopic signal, interfere it with a local oscillator and then obtain the signal. All of these steps are avoided by the microscopic QED definition given here. Sec. VII provides a summary of our results.

We first consider a single molecule driven by an optical field *E*(*t*). This is all we need for calculating incoherent signals. The extension to coherent signals from molecular ensembles will be done in Sec. V.

The total Hamiltonian is given by

$$\widehat{H}={\widehat{H}}_{0}+{\widehat{H}}_{F}+{\widehat{H}}_{\mathit{int}},$$

(1)

where *Ĥ*_{0} represents the free molecule, and

$${\widehat{H}}_{F}=\sum _{s}\hslash {\omega}_{s}{\widehat{a}}_{s}^{\u2020}{\widehat{a}}_{s},$$

(2)

is the radiation field Hamiltonian. In Eq. (2)
${\widehat{a}}_{s}({\widehat{a}}_{s}^{\u2020})$ denotes the destruction (creation) operator for the *s*th field mode and
${\widehat{a}}_{s}({\widehat{a}}_{s}^{\u2020})$ satisfy the boson commutation relation
$[{\widehat{a}}_{s},{\widehat{a}}_{{s}^{\prime}}^{\u2020}]={\delta}_{s{s}^{\prime}}$. For clarity, we shall omit the unit vector describing the field polarization of a particular mode, and treat the field as a scalar. The generalization to include the full tensorial expressions is straightforward.

We shall treat the emitted photon modes with frequency *ω _{s}* quantum mechanically. () are the positive (negative) frequency components of the quantized electric field

$$\widehat{E}(\mathbf{r},t)=\widehat{\mathcal{E}}(\mathbf{r},t)+{\widehat{\mathcal{E}}}^{\u2020}(\mathbf{r},t),$$

(3)

with

$$\widehat{\mathcal{E}}(\mathbf{r},t)=\sum _{s}{\left(\frac{2\pi \hslash {\omega}_{s}}{\mathrm{\Omega}}\right)}^{1/2}{\widehat{a}}_{s}\phantom{\rule{0.16667em}{0ex}}{\text{e}}^{i{\mathbf{k}}_{s}\xb7\mathbf{r}-i{\omega}_{s}t},$$

(4)

and Ω is the quantization volume.

In the rotating wave approximation (RWA) the molecule-field interaction is given by

$${\widehat{H}}_{\mathit{int}}(t)=\widehat{\mathcal{E}}(\mathbf{r},t){\widehat{V}}^{\u2020}+{\widehat{\mathcal{E}}}^{\u2020}(\mathbf{r},t)\widehat{V},$$

(5)

where we have partitioned the dipole operator as = + ^{†}. Here ^{†}() are the creation (annihilation) operators for excitations. The time-dependence in Eq. (5) is in the interaction picture with respect to *Ĥ _{F}*. Consequently, the time dependence of the total (field+molecule) density operator

$$\widehat{H}(t)={\widehat{H}}_{0}+{\widehat{H}}_{\mathit{int}}(t).$$

(6)

We assume that the field is initially in a coherent state,

$$\mid {\mathrm{\Psi}}_{F}\rangle ={A}_{0}exp\left\{\sum _{s}{\widehat{a}}_{s}^{\u2020}{\alpha}_{s}\right\}\mid 0\rangle ,$$

(7)

where

$${A}_{0}=exp(\sum _{s}{\mid {\alpha}_{s}\mid}^{2}),$$

(8)

is the normalization such that Ψ* _{F}*|Ψ

$$\langle {\mathrm{\Psi}}_{F}\mid \widehat{E}(\mathbf{r},t)\mid {\mathrm{\Psi}}_{F}\rangle =\mathcal{E}(\mathbf{r},t)+c.c.,$$

(9)

where

$$\mathcal{E}(\mathbf{r},t)=\sum _{s}{\left(\frac{2\pi \hslash {\omega}_{s}}{\mathrm{\Omega}}\right)}^{1/2}{\mathbf{e}}^{i{\mathbf{k}}_{s}\xb7\mathbf{r}-i{\omega}_{s}t}{\alpha}_{s},$$

(10)

is the field amplitude at space-point **r**.

We shall define the time and frequency resolved signal *S*(*t*) as the rate of change of the photon occupation, i.e.

$$S(t):=\frac{\text{d}}{\text{d}t}{\left(\stackrel{\sim}{\mathcal{N}}\right)}_{\text{D},\text{T}},$$

(11)

where

$$\widehat{\mathcal{N}}:=\sum _{s}{\widehat{a}}_{s}^{\u2020}{\widehat{a}}_{s},$$

(12)

denotes the photon number operator. This definition applies to both spontaneous and stimulated process, as will be demonstrated later. In Eq. (12) the sum extends over all modes of interest in the experiment under consideration.

We shall adopt of the following convention for the various ensemble averages, one of which has been used in Eq. (11). In Eq. (11), (…)_{D,T} stands for a trace with respect to the density operator of the molecule where *Ĥ _{int}* includes

The expectation value in Eq. (11) may be conveniently evaluated by switching to the Heisenberg picture, using the basic identity

$$\begin{array}{l}\frac{\text{d}}{\text{d}t}{\left(\widehat{\mathcal{N}}\right)}_{\text{D},\text{T}}\equiv \langle \frac{\text{d}}{\text{d}t}{\widehat{\mathcal{N}}}_{H}\rangle \\ =\langle \sum _{s}\frac{i}{\hslash}[{\widehat{H}}_{\mathit{int}}(t),{\widehat{a}}_{s,H}^{\u2020}{\widehat{a}}_{s,H}]\rangle ,\end{array}$$

(13)

and the canonical commutation relations for ${\widehat{a}}_{s}({\widehat{a}}_{s}^{\u2020})$. In Eq. (13) we denote operators in the Heisenberg picture by a subscript “H”.

Evaluating the commutator in Eq. (13) and transforming back to the original interaction picure, we finally obtain for the signal, Eq. (11),

$$S(t)=-\frac{2}{\hslash}Im\left\{{\left(\widehat{\mathcal{E}}(\mathbf{r},t){\widehat{V}}^{\u2020}\right)}_{\text{D},\text{T}}\right\}.$$

(14)

Taking into account the time dependence of the total (system + field) density operator, which is determined by the Hamiltonian of Eq. (6), the trace on the rhs of Eq. (14) can be computed perturbatively in the field *E*(**r**, *t*) by switching to the interaction picture with respect to the Hamiltonian *Ĥ*_{0}. We shall compute this trace using superoperators in Liouville space. To this end we associate with each operator *Â* in Hilbert space, a left (L) and a right (R) operation in Liouville space, by defining

$$\begin{array}{l}{\widehat{A}}_{L}\widehat{X}:=\widehat{A}\widehat{X},\\ {\widehat{A}}_{R}\widehat{X}:=\widehat{X}\widehat{A},\end{array}$$

(15)

where denotes an arbitrary Hilbert space operator. Furthermore, we introduce the following linear combinations of *L*/*R* operations, which will be referred to as +/− operations

$${\widehat{A}}_{\pm}:=\frac{1}{\sqrt{2}}[{\widehat{A}}_{L}\pm {\widehat{A}}_{R}].$$

(16)

A key tool in the following manipulations is the time-ordering operator in Liouville space, ; when acting on a product of superoperators, it reorders them such that their time arguments increase from right to left. In the interaction picture Eq. (14) can be expressed as

$$S(t)=-\frac{2}{\hslash}Im\left[\langle \mathcal{T}{\widehat{\mathcal{E}}}_{L}(\mathbf{r},t){\widehat{V}}_{L}^{\u2020}(t)exp\left\{-\frac{i}{\hslash}{\int}_{-\infty}^{t}\text{d}\tau \sqrt{2}{\mathcal{H}}_{\mathit{int}-}(\tau )\right\}\rangle \right],$$

(17)

where is the superoperator corresponding to *Ĥ _{int}*, Eq. (5).

Note that *Ĥ _{int}* contains two interaction terms, one with the incoming field modes
${\widehat{H}}_{\mathit{int}}^{\prime}$ and the other with the signal modes
${\widehat{H}}_{\mathit{int}}^{\u2033}$. As before, the explicit time dependence represents the interaction picture

$$\begin{array}{l}{\widehat{a}}_{sX}(t):={\text{e}}^{{\scriptstyle \frac{i}{\hslash}}{\mathcal{H}}_{FX}t}{\widehat{a}}_{sX}{\text{e}}^{-{\scriptstyle \frac{i}{\hslash}}{\mathcal{H}}_{FX}t},\\ {\widehat{V}}_{X}(t):={\text{e}}^{{\scriptstyle \frac{i}{\hslash}}{\mathcal{H}}_{0X}t}{\widehat{V}}_{X}{\text{e}}^{-{\scriptstyle \frac{i}{\hslash}}{\mathcal{H}}_{0X}t},\end{array}$$

(18)

where *X* = *L*, *R* and () is the superoperator corresponding to *Ĥ _{F}*(

To lowest order in c ${\widehat{H}}_{\mathit{int}}^{\u2033}$ Eq. (17) gives

$$S(t)=\frac{2\sqrt{2}}{{\hslash}^{2}}Re\sum _{s}{\left(\frac{2\pi \hslash {\omega}_{s}}{\mathrm{\Omega}}\right)}^{1/2}{\text{e}}^{i{\mathbf{k}}_{s}\xb7\mathbf{r}}\times {\int}_{-\infty}^{t}\text{d}\tau {\left(\mathcal{T}{\widehat{a}}_{sL}(t){\widehat{V}}_{L}^{\u2020}(t){\mathcal{H}}_{\mathit{int}-}^{\u2033}(\tau )\right)}_{\text{D}}.$$

(19)

The superoperator expression corresponding to Eq. (5) gives

$$\sqrt{2}{\mathcal{H}}_{\mathit{int}-}(\tau )={\widehat{\mathcal{E}}}_{L}(\mathbf{r},\tau ){\widehat{V}}_{L}^{\u2020}(\tau )+{\widehat{\mathcal{E}}}_{L}^{\u2020}(\mathbf{r},\tau ){\widehat{V}}_{L}(\tau )-{V}_{R}^{\u2020}(\tau ){\mathcal{E}}_{R}(\mathbf{r},\tau )-{V}_{R}(\tau ){\mathcal{E}}_{R}^{\u2020}(\mathbf{r},\tau ).$$

(20)

Substituting in Eq. (19) and factorizing the correlation functions into products of field and system parts, we get

$$S(t)=\frac{4\pi \hslash}{\mathrm{\Omega}}Re{\sum}_{s{s}^{\prime}}{({\omega}_{s}{\omega}_{{s}^{\prime}})}^{1/2}{\text{e}}^{i({\mathbf{k}}_{s}-{\mathbf{k}}_{{s}^{\prime}})\xb7\mathbf{r}}\times {\int}_{-\infty}^{t}d\tau \phantom{\rule{0.16667em}{0ex}}[{d}_{LR}^{s{s}^{\prime}}(t,\tau ){D}_{RL}(\tau ,t)-{d}_{LL}^{s{s}^{\prime}}(t,\tau ){D}_{LL}(\tau ,t)],$$

(21)

where we have introduced the field and system nonequilibrium Green’s functions, *d _{XY}* and

$$\begin{array}{l}{d}_{XY}^{s{s}^{\prime}}(t,\tau ):=-\frac{i}{\hslash}{\left(\mathcal{T}{\widehat{a}}_{sX}(t){\widehat{a}}_{{s}^{\prime}Y}^{\u2020}(\tau )\right)}_{\text{D}},\\ {D}_{XY}(t,\tau ):=-\frac{i}{\hslash}{\left(\mathcal{T}{\widehat{V}}_{X}(t){\widehat{V}}_{Y}^{\u2020}(\tau )\right)}_{\text{D}}.\end{array}$$

(22)

In Eq. (21) we assume that the system is in its ground state for *t* → −∞, hence the following system correlation functions vanish

$$\begin{array}{l}{\left(\mathcal{T}{\widehat{V}}_{L}^{\u2020}(t){\widehat{V}}_{L}^{\u2020}(\tau )\right)}_{\text{D}}=0,\\ {\left(\mathcal{T}{\widehat{V}}_{L}^{\u2020}(t){\widehat{V}}_{R}^{\u2020}(\tau )\right)}_{\text{D}}=0,\phantom{\rule{0.38889em}{0ex}}t>\tau .\end{array}$$

(23)

From the basic definitions of the Green’s functions we have

$$\begin{array}{l}{D}_{LL}(\tau ,t)=\theta (\tau -t){D}_{RL}(\tau ,t)+\theta (t-\tau ){D}_{LR}(\tau ,t),\\ {d}_{LL}(t,\tau )=\theta (t-\tau ){d}_{RL}(t,\tau )+\theta (\tau -t){d}_{LR}(t,\tau ).\end{array}$$

(24)

Substituting Eq. (24) in (21) gives

$$S(t)=\frac{4\pi \hslash}{\mathrm{\Omega}}Re{\sum}_{s{s}^{\prime}}{({\omega}_{s}{\omega}_{{s}^{\prime}})}^{1/2}{\text{e}}^{i({\mathbf{k}}_{s}-{\mathbf{k}}_{{s}^{\prime}})\xb7\mathbf{r}}\times {\int}_{-\infty}^{t}\text{d}\tau \phantom{\rule{0.16667em}{0ex}}[{d}_{LR}^{s{s}^{\prime}}(t,\tau ){D}_{RL}(\tau ,t)-{d}_{RL}^{s{s}^{\prime}}(t,\tau ){D}_{LR}(\tau ,t)].$$

(25)

Using Eqs. (7) and (22), the Green’s functions for the field are given by

$$\begin{array}{l}{d}_{LR}^{s{s}^{\prime}}(t,\tau )=-\frac{i}{\hslash}{\left({\widehat{a}}_{{s}^{\prime}}^{\u2020}(\tau ){\widehat{a}}_{s}(t)\right)}_{\text{D}}\\ =-\frac{i}{\hslash}{\text{e}}^{-i{\omega}_{s}t}{\text{e}}^{i{\omega}_{{s}^{\prime}}\tau}{\alpha}_{s}{\alpha}_{{s}^{\prime}}^{\ast},\end{array}$$

(26)

$$\begin{array}{l}{d}_{RL}^{s{s}^{\prime}}(t,\tau )=-\frac{i}{\hslash}{\left({\widehat{a}}_{s}(t){\widehat{a}}_{{s}^{\prime}}^{\u2020}(\tau )\right)}_{\text{D}}\\ =-\frac{i}{\hslash}{\text{e}}^{-i{\omega}_{s}t}{\text{e}}^{i{\omega}_{{s}^{\prime}}\tau}({\delta}_{s{s}^{\prime}}+{\alpha}_{s}{\alpha}_{{s}^{\prime}}^{\ast}).\end{array}$$

(27)

The signal, Eq. (25), then becomes

$$S(t)=\frac{4\pi \hslash}{\mathrm{\Omega}}Im\sum _{s{s}^{\prime}}{({\omega}_{s}{\omega}_{{s}^{\prime}})}^{1/2}{\text{e}}^{i({\mathbf{k}}_{s}-{\mathbf{k}}_{{s}^{\prime}})\xb7\mathbf{r}}{\text{e}}^{-i{\omega}_{s}t}\times {\int}_{-\infty}^{t}d\tau {\text{e}}^{i{\omega}_{{s}^{\prime}}\tau}\left[{\alpha}_{s}{\alpha}_{{s}^{\prime}}^{\ast}{D}_{RL}(\tau ,t)-({\delta}_{s{s}^{\prime}}+{\alpha}_{s}{\alpha}_{{s}^{\prime}}^{\ast}){D}_{LR}(\tau ,t)\right].$$

(28)

In Eq. (28), the term with *δ _{ss}*

Equation (28) can also be written as

$$S(t)=2Im{\int}_{-\infty}^{t}\text{d}\tau \phantom{\rule{0.16667em}{0ex}}[{D}_{RL}(\tau ,t)-{D}_{LR}(\tau ,t)]\mathcal{E}(\mathbf{r},t){\mathcal{E}}^{\ast}(\mathbf{r},\tau )-\frac{4\pi \hslash}{\mathrm{\Omega}}{\sum}_{s}{\omega}_{s}Im{\int}_{-\infty}^{t}d\tau {\text{e}}^{-i{\omega}_{s}(t-\tau )}{D}_{LR}(\tau ,t),$$

(29)

or equivalently

$$S(t)=2Im{\int}_{-\infty}^{t}d\tau \phantom{\rule{0.16667em}{0ex}}[{D}_{RL}(\tau ,t)\mathcal{E}(\mathbf{r},t){\mathcal{E}}^{\ast}(\mathbf{r},\tau )-{D}_{LR}(t,\tau ){\mathcal{E}}^{\ast}(\mathbf{r},t)\mathcal{E}(\mathbf{r},\tau )]-\frac{4\pi \hslash}{\mathrm{\Omega}}{\sum}_{s}{\omega}_{s}Im{\int}_{-\infty}^{t}d\tau {\text{e}}^{i{\omega}_{s}(t-\tau )}{D}_{LR}(t,\tau ).$$

(30)

In Eq. (30) we made use of the symmetry relation

$${D}_{LR}^{\ast}(\tau ,t)=-{D}_{LR}(t,\tau ),$$

(31)

to arrive at an expression, where interactions at time *t* always occur from the left (i.e. on the ket). This choice will be adopted in the diagrammatic representation, which will be introduced in Sec. III.

Equations (29) or equivalently (30) constitute our basic NEGF expression for incoherent optical signals. The first term in the brackets describes the creation of excitations in the system by absorbing photons from the field and the second term represents the reverse process. Both processes are stimulated. The signal is thus given by the net photon flux. The last term represents spontaneous emission. In the coming two sections we shall apply these expressions to compute different signals.

To describe spontaneous light emission (SLE), which is an incoherent process, we assume that the scattered field is initially in its vacuum state and is generated by spontaneous emission. Thus (**r**, *t*) = 0 and Eq. (30), reduces to

$${S}_{\text{SLE}}(t)=-\frac{4\pi \hslash}{\mathrm{\Omega}}{\omega}_{s}Im{\int}_{-\infty}^{t}d\tau {\text{e}}^{i{\omega}_{s}(t-\tau )}{D}_{LR}(t,\tau ),$$

(32)

where it is furthermore assumed that the scattered field has only one mode *ω _{s}* thus dropping the sum over

Expanding the Green’s function *D _{LR}* to second order in the incoming field yields

$${S}_{\text{SLE}}(t)=-\frac{4\pi}{\mathrm{\Omega}{\hslash}^{2}}{\omega}_{s}Re{\int}_{-\infty}^{t}\text{d}\tau \int \text{d}{\tau}_{1}\int \text{d}{\tau}_{2}{\text{e}}^{i{\omega}_{s}(t-\tau )}\times \langle \mathcal{T}{\widehat{V}}_{L}(t){\widehat{V}}_{R}^{\u2020}(\tau ){\mathcal{H}}_{\mathit{int}-}({\tau}_{1}){\mathcal{H}}_{\mathit{int}-}({\tau}_{2})\rangle .$$

(33)

In order to go any further, we need to expand the four-point correlation function of Eq. (33) in *L*/*R* operations. A systematic, diagrammatic technique for selecting the terms in Eq. (33) that contribute within the RWA, will be introduced next.

We start with Eq. (33), where the correlation function originates from expanding the Green’s function *D _{LR}*(

In order to derive an explicit expression for this correlation function in terms of superoperators, we first note that the interactions at times *t* and *τ* both represent photon emission, one on the ket at time *t*, the other on the bra at time *τ*. Since the system is assumed to be initially in its ground state, both emission events have to be preceded by two absorptions on either side. Although this implies a certain order of interactions on both the ket and the bra individually, their relative ordering in physical time remains a priori undetermined. This *partial time ordering* can be expressed diagrammatically by arranging the interactions on a loop. We adopt the following general rules for constructing and reading the diagrams:

- Time runs along the loop clockwise from bottom left to bottom right.
- The left strand of the loop represents the ket, the right corresponds to the bra.
- Each interaction with a field mode is represented by a wavy line on either the right (R-operators) or the left (L-operators).
- The field is indicated by dressing the wavy lines with arrows, where an arrow pointing to the right represents the field annihilation operator (
**r**,*t*), which involves the term e^{i(ks·r−ωst)}(see Eq. (4)). Conversely, an arrow pointing to the left corresponds to the field creation operator (**r**,*t*), being associated with e^{−i(ks·r−ωst)}. This is made explicit by adding the wavevectors ±**k**to the arrows._{s} - Within the RWA (Eq. (5)), each interaction with (
**r**,*t*) is accompanied by applying the operator*V*^{†}, which leads to excitation of the state represented by ket and deexcitation of the state represented by the bra, respectively. Arrows pointing “inwards” (i.e. pointing to the right on the ket and to the left on the bra) consequently cause absorption of a photon by exciting the system, whereas arrows pointing “outwards” (i.e. pointing to the left on the bra and to the right on the ket) represent deexciting the system by photon emission. - The interaction at the observation time
*t*, is fixed and is always the last. As a convention, it is chosen to occur from the left. This can always be achieved by a reflection of all interactions through the center line between the ket and the bra, which corresponds to taking the complex conjugate of the original correlation function. - Interactions within each strand are time-ordered, but interactions on different strands are not. Each loop can be further decomposed into several fully-time-ordered diagrams (double sided Feynman diagrams). These can be generated from the loop by simply shifting the arrows along each strand, thus changing their position relative to the interactions on the other strand. Each of these relative positions then gives rise to a particular fully time-ordered diagram.
- The overall sign of the correlation function is given by (−1)
^{NR}, where*N*stands for the number of interactions from the right, at times_{R}*τ*_{1}and*τ*_{2}.

We note that loop diagrams drawn according to the rules presented above, lead to double sided Feynman diagrams that follow the standard conventions employed within the semiclassical theory of nonlinear optics [3]. Using these rules, the SLE is represented by the single loop diagram displayed in Fig. 1. We denote the incoming field as ±**k**_{1} and the signal field as ±**k**_{2}. The loop translates Eq. (33) into the following expression for the signal,

Loop diagram for SLE. **k**_{1} is the incoming field, and **k**_{2} is the signal field. Note that the interactions are time ordered within each strand, but not between strands. *s*_{1}, *s*_{2}, *s*_{3} are the delay times between the interactions along the loop.

$$\begin{array}{r}{S}_{\text{SLE}}(t)=\frac{4\pi}{\mathrm{\Omega}{\hslash}^{2}}{\omega}_{2}Re{\int}_{-\infty}^{t}\text{d}\tau \int \text{d}{\tau}_{1}\int \text{d}{\tau}_{2}{\text{e}}^{i{\omega}_{2}(t-\tau )}\\ \langle \mathcal{T}{\widehat{V}}_{L}(t){\widehat{V}}_{R}^{\u2020}(\tau ){\widehat{V}}_{L}^{\u2020}({\tau}_{2}){\widehat{V}}_{R}({\tau}_{1})\rangle \langle \mathcal{T}{\widehat{\mathcal{E}}}_{L}(\mathbf{r},{\tau}_{2}){\widehat{\mathcal{E}}}_{R}^{\u2020}(\mathbf{r},{\tau}_{1})\rangle ,\end{array}$$

(34)

where an additional prefactor of 2 was added to take into account the symmetry with respect to the dummy variables *τ*_{1} and *τ*_{2}.

If desired, this can further be decomposed into fully time-ordered terms using the double sided diagrams of Fig. 2. As indicated above we can now e.g. shift the arrow in Fig. 1 corresponding to time *τ*_{2} along the ket, thus obtaining 3 possible relative positions with respect to the interactions at times *τ* and *τ*_{1}. The loop diagram consequently splits into 3 double sided Feynman diagrams. These are depicted in Fig. 2.

The resulting fully-time-ordered expression for the signal reads

$$\begin{array}{l}{S}_{\text{SLE}}(t)=\frac{4\pi}{\mathrm{\Omega}{\hslash}^{2}}{\omega}_{2}Re{\int}_{-\infty}^{t}\text{d}\tau \int \text{d}{\tau}_{1}\int \text{d}{\tau}_{2}{\mathcal{E}}_{1}^{\ast}(\mathbf{r},{\tau}_{1}){\mathcal{E}}_{1}(\mathbf{r},{\tau}_{2})\\ {\text{e}}^{i{\omega}_{s}(t-\tau )}\left[\theta (t\tau )\theta (\tau {\tau}_{1})\theta ({\tau}_{1}{\tau}_{2})\langle {V}_{L}(t){V}_{R}^{\u2020}(\tau ){V}_{R}({\tau}_{1}){V}_{L}^{\u2020}({\tau}_{2})\rangle +\theta (t{\tau}_{2})\theta ({\tau}_{2}\tau )\theta (\tau {\tau}_{1})\langle {\widehat{V}}_{L}(t){\widehat{V}}_{L}^{\u2020}({\tau}_{2}){\widehat{V}}_{R}^{\u2020}(\tau ){\widehat{V}}_{R}({\tau}_{1})\rangle +\theta (t\tau )\theta (\tau {\tau}_{2})\theta ({\tau}_{2}{\tau}_{1})\langle {V}_{L}(t){V}_{R}^{\u2020}(\tau ){V}_{L}^{\u2020}({\tau}_{2}){V}_{R}({\tau}_{1})\rangle \right].\end{array}$$

(35)

Here, we have made use of Eq. (7) to evaluate the correlation functions for the incoming field, since

$$\begin{array}{c}\langle \mathcal{T}{\widehat{\mathcal{E}}}_{L}(\mathbf{r},{\tau}_{2}){\widehat{\mathcal{E}}}_{R}^{\u2020}(\mathbf{r},{\tau}_{1})\rangle =\frac{2\pi \hslash}{\mathrm{\Omega}}{\omega}_{1}{\mid {\alpha}_{1}\mid}^{2}{\text{e}}^{-i{\omega}_{1}({\tau}_{1}-{\tau}_{2})}\\ ={\mathcal{E}}_{1}^{\ast}(\mathbf{r},{\tau}_{1}){\mathcal{E}}_{1}(\mathbf{r},{\tau}_{2}).\end{array}$$

(36)

Note that, in Eq. (35) full time-ordering is expressed by adding a product of step functions to each of the contributing terms. We apply the short hand notation *θ*(*τ*_{1}*τ*_{2}) = *θ*(*τ*_{1} − *τ*_{2}).

Next, we calculate the signal for a frequency-domain experiment with stationary beams. This will also illustrate how to extend the diagrammatic rules to frequency-domain experiments. We adopt the following convention for the Fourier transform of a function *f*

$$\begin{array}{l}\stackrel{\sim}{f}(\omega )=\int \text{d}t{\text{e}}^{i\omega t}f(t),\\ f(t)=\frac{1}{2\pi}\int \text{d}\omega {\text{e}}^{-i\omega t}\stackrel{\sim}{f}(\omega ).\end{array}$$

(37)

When the incoming field is stationary, the integrals in Eq. (34) can be evaluated directly, which by application of Eq. (37) leads to a multiple Fourier transform of the partially ordered 4-point system correlation function in Eq. (34). To simplify this calculation, it is advantageous to switch to a new set of time variables (“*s*-variables”), *s*_{1}:= *t* − *τ*_{2}, *s*_{2}:= *t* − *τ*, *s*_{3}:= *τ* − *τ*_{1}, which represent the time intervals of interactions *along* the loop. This choice becomes intuitive when looking at the diagram, Fig. 1.

We can now show that the correlation function in Eq. (35) depends on the *s*-variables in a simple way. To this end, we transform the system correlation function in Eq. (34) to a form that only involves *L*-superoperators, hence

$$\langle \mathcal{T}{\widehat{V}}_{L}(t){\widehat{V}}_{R}^{\u2020}(\tau ){\widehat{V}}_{R}({\tau}_{1}){\widehat{V}}_{L}^{\u2020}({\tau}_{2})\rangle =\langle {\widehat{V}}_{L}({\tau}_{1}){\widehat{V}}_{L}^{\u2020}(\tau ){\widehat{V}}_{L}(t){\widehat{V}}_{L}^{\u2020}({\tau}_{2})\rangle $$

(38)

In terms of the *s*-variables, Eq. (38) becomes

$$\langle \mathcal{T}{\widehat{V}}_{L}(t){\widehat{V}}_{R}^{\u2020}(\tau ){\widehat{V}}_{R}({\tau}_{1}){\widehat{V}}_{L}^{\u2020}({\tau}_{2})\rangle =\langle {\widehat{V}}_{L}{\mathcal{G}}^{\u2020}({s}_{3}){\widehat{V}}_{L}^{\u2020}{\mathcal{G}}^{\u2020}({s}_{2}){\widehat{V}}_{L}\mathcal{G}({s}_{1}){\widehat{V}}_{L}^{\u2020}\rangle ,$$

(39)

where we define the Liouville space propagator

$$\mathcal{G}(s):=\theta (s)exp\left(-\frac{i}{\hslash}\sqrt{2}{\mathcal{H}}_{0-}s\right).$$

(40)

Note that, when associating all interactions with the ket, one needs to propergate twice backwards in time. This amounts to applying the retarded propagator twice.

Using Eq. (39) the integration in Eq. (34), leads to a multiple Fourier transform of the correlation function in terms of the loop (s)- variables. This is readily calculated, giving the following expression for the frequency-domain SLE signal

$$\begin{array}{c}{S}_{\text{SLE}}({\omega}_{1},{\omega}_{2})=\frac{4\pi}{{\hslash}^{2}\mathrm{\Omega}}{\mid {\mathcal{E}}_{1,0}\mid}^{2}{\omega}_{2}\\ Im\langle {\widehat{V}}_{L}{\mathcal{G}}^{\u2020}({\omega}_{1}){\widehat{V}}_{L}^{\u2020}{\mathcal{G}}^{\u2020}({\omega}_{1}-{\omega}_{2}){\widehat{V}}_{L}\mathcal{G}({\omega}_{1}){\widehat{V}}_{L}^{\u2020}\rangle .\end{array}$$

(41)

Since all interactions in Eq. (41) are “left”-operations that act on the ket, it is possible to express the final result for the signal in Hilbert space, giving

$$\begin{array}{c}{S}_{\text{SLE}}({\omega}_{1},{\omega}_{2})=\frac{4\pi}{{\hslash}^{2}\mathrm{\Omega}}{\mid {\mathcal{E}}_{1,0}\mid}^{2}{\omega}_{2}\\ Im\langle \widehat{V}{\widehat{G}}^{\u2020}({\omega}_{g}+{\omega}_{1}){\widehat{V}}^{\u2020}{\widehat{G}}^{\u2020}({\omega}_{g}+{\omega}_{1}-{\omega}_{2})\widehat{V}\widehat{G}({\omega}_{g}+{\omega}_{1}){\widehat{V}}^{\u2020}\rangle .\end{array}$$

(42)

Upon conversion of Eq. (41), the frequency-domain Liouville space propagators,

$$\mathcal{G}(\omega )=\frac{\hslash}{\hslash \omega -\sqrt{2}{\mathcal{H}}_{0-}+i\eta},$$

(43)

are replaced by Hilbert-space propagators

$$\widehat{G}(\omega +{\omega}_{g})=\frac{\hslash}{\hslash \omega +\hslash {\omega}_{g}-{\widehat{H}}_{0}+i\eta}.$$

(44)

Here, *ω _{g}* is the material frequency of the ground state, which accounts for the free evolution of the bra. In both Eqs. (43) and (44) the infinitesimal

Equation (42) can be recast in the form

$$\begin{array}{c}{S}_{\text{SLE}}({\omega}_{1},{\omega}_{2})=-\frac{2\pi i}{{\hslash}^{2}\mathrm{\Omega}}{\mid {\mathcal{E}}_{1,0}\mid}^{2}{\omega}_{2}\\ \langle \widehat{V}{\widehat{G}}^{\u2020}({\omega}_{g}+{\omega}_{1}){\widehat{V}}^{\u2020}\left[{\widehat{G}}^{\u2020}({\omega}_{g}+{\omega}_{1}-{\omega}_{2})-\widehat{G}({\omega}_{g}+{\omega}_{1}-{\omega}_{2})\right]\widehat{V}\widehat{G}({\omega}_{g}+{\omega}_{1}){\widehat{V}}^{\u2020}\rangle ,\end{array}$$

(45)

which by using the level scheme of Fig. 3(a) yields the Kramers-Heisenberg formula

$${S}_{\text{SLE}}({\omega}_{1},{\omega}_{2})=-\frac{4{\pi}^{2}}{{\hslash}^{2}\mathrm{\Omega}}{\omega}_{2}{\mid {\mathcal{E}}_{1,0}\mid}^{2}{\mid {\chi}_{ca}\mid}^{2}\frac{{\mathrm{\Gamma}}_{ca}/\pi}{{({\omega}_{1}-{\omega}_{2}-{\omega}_{ca})}^{2}+{\mathrm{\Gamma}}_{ca}^{2}},$$

(46)

where

$${\chi}_{ca}=\frac{{\mu}_{ab}{\mu}_{bc}}{{\omega}_{1}-{\omega}_{ba}+i{\mathrm{\Gamma}}_{ba}}.$$

(47)

In Eqs. (46) and (47) *ω _{ca}* =

By proceeding along the same line for an arbitrary loop diagram, one can establish the following rules, that allow to translate a given diagram into its frequency-domain expression. These complement the rules given earlier for time-domain expressions.

- In the frequency-domain the loop translates into an alternating product of interactions (arrows) and periods of free evolutions (vertical solid lines) along the loop.
- Since the loop time goes clockwise along the loop, periods of free evolution on the left amount to propagating forward in real time (
*Ĝ*), whereas evolution on the right corresponds to backward propergation (*Ĝ*^{†}). - Each
*Ĝ*adds the multiplicative factor*i*, whereas each*Ĝ*^{†}results in a multiplication by (−*i*). - Frequency arguments of each propagator are cumulative, i.e. they are given by the sum of all “earlier” interactions along the loop. The ground state frequency
*ω*must be added to all Green’s functions._{g}

Equation (42) can be immediately generated from Fig. 4 by applying these rules.

The simplest nonlinear optical technique involves two fields, **k**_{1} (the pump) and **k**_{2} (the probe). The signal is defined as the difference of the probe transmitted intensity with and without the pump. We assume that the probe intensity is high so that spontaneous emission is negligible compared to stimulated emission. The second term in Eqs. (29) and (30), which describes spontaneous emission, is thus neglected. Using Eq. (29) and expanding to second order in the incoming field, the signal then reduces to

$$\begin{array}{l}{S}_{\text{PP}}(t)=2\frac{1}{{\hslash}^{3}}Re{\int}_{-\infty}^{t}d\tau \int d{\tau}_{1}\int d{\tau}_{2}{\mathcal{E}}_{2}(\mathbf{r},t){\mathcal{E}}_{2}^{\ast}(\mathbf{r},\tau )\\ \left[\langle \mathcal{T}{\widehat{V}}_{R}(\tau ){\widehat{V}}_{L}^{\u2020}(t){\mathcal{H}}_{\mathit{int}-}({\tau}_{1}){\mathcal{H}}_{\mathit{int}-}({\tau}_{2})\rangle -\langle \mathcal{T}{\widehat{V}}_{L}(\tau ){\widehat{V}}_{R}^{\u2020}(t){\mathcal{H}}_{\mathit{int}-}({\tau}_{1}){\mathcal{H}}_{\mathit{int}-}({\tau}_{2})\rangle \right].\end{array}$$

(48)

Using the identity

$${D}_{RL}(\tau ,t)-{D}_{LR}(\tau ,t)=\frac{i}{\hslash}{\left(\mathcal{T}{\widehat{V}}_{+}^{\u2020}(t){\widehat{V}}_{-}(\tau )\right)}_{\text{D}},$$

(49)

we can express the pump-probe signal in terms of +/− operators as

$$\begin{array}{c}{S}_{\text{PP}}(t)=-\frac{2}{{\hslash}^{3}}Re{\int}_{-\infty}^{t}d\tau \int \text{d}{\tau}_{1}\int \text{d}{\tau}_{2}{\mathcal{E}}_{2}(\mathbf{r},t){\mathcal{E}}_{2}^{\ast}(\mathbf{r},\tau )\\ \langle \mathcal{T}{\widehat{V}}_{+}^{\u2020}(t){V}_{-}(\tau ){\mathcal{H}}_{\mathit{int}-}({\tau}_{1}){\mathcal{H}}_{\mathit{int}-}({\tau}_{2})\rangle .\end{array}$$

(50)

This is equivalent to the classical expression for the signal in terms of *χ*^{(3)}. This becomes even more apparent if Eq. (50) is rewritten in the form

$$\begin{array}{c}{S}_{\text{PP}}(t)=-\frac{2}{{\hslash}^{3}}Re{\int}_{-\infty}^{t}\text{d}\tau \int \text{d}{\tau}_{1}\int \text{d}{\tau}_{2}{\mathcal{E}}_{2}(\mathbf{r},t){\mathcal{E}}_{2}^{\ast}(\mathbf{r},\tau )\\ {\mathcal{E}}_{1}(\mathbf{r},{\tau}_{1}){\mathcal{E}}_{1}^{\ast}(\mathbf{r},{\tau}_{2})\langle \mathcal{T}{\widehat{V}}_{+}^{\u2020}(t){\widehat{V}}_{-}(\tau ){\widehat{V}}_{-}^{\u2020}({\tau}_{1}){\widehat{V}}_{-}({\tau}_{2})\rangle ,\end{array}$$

(51)

which is equivalent to Eq. (50), if only terms proportional to the field intensities are taken into account. The details are given in appendix A.

The forgoing example of SLE illustrated the use of Keldysh-Schwinger loops to compute the signal. In particular, by using this *partially time ordered* diagrammatic representation, we were able to reduce the number of diagrams from three to one. The merits of this compact notation will become even more apparent for the pump-probe signal, as will be shown next.

We start with an expression for the signal using Eq. (30)

$$\begin{array}{c}{S}_{\text{PP}}(t)=2\frac{1}{{\hslash}^{3}}Re{\int}_{-\infty}^{t}d\tau \int d{\tau}_{1}\int d{\tau}_{2}\\ \left[{\mathcal{E}}_{2}(\mathbf{r},t){\mathcal{E}}_{2}^{\ast}(\mathbf{r},\tau )\langle \mathcal{T}{\widehat{V}}_{R}(\tau ){\widehat{V}}_{L}^{\u2020}(t){\mathcal{H}}_{\mathit{int}-}({\tau}_{1}){\mathcal{H}}_{\mathit{int}-}({\tau}_{2})\rangle -{\mathcal{E}}_{2}^{\ast}(\mathbf{r},t){\mathcal{E}}_{2}(\mathbf{r},\tau )\langle \mathcal{T}{\widehat{V}}_{L}(t){\widehat{V}}_{R}^{\u2020}(\tau ){\mathcal{H}}_{\mathit{int}-}({\tau}_{1}){\mathcal{H}}_{\mathit{int}-}({\tau}_{2})\rangle \right],\end{array}$$

(52)

where in analogy to Eq. (48), we neglect the spontaneous emission term.

We note that the second term in Eq. (52) coincides with the SLE correlation function, Eq. (33). What remains therefore is to resolve the first term in terms of *L*/*R* operations. For simplicity we assume a three level ladder with sequential transition dipole moments as shown in Fig. 3(b).

We again start our analysis by looking at the two interactions at times *t* and *τ*. Contrary to the second term in Eq. (52), both of these correspond to absorptions rather than emissions. For the dipole selection rules shown in Fig. 3(b), this leaves two possibilities for placing the remaining interactions along the loop: We either have three interactions on one side, or two interactions on either side (see Fig. 4).

For the former case and the model of Fig. 3(b), these can only correspond to two absorbtions and one emission. By applying the same argument used for SLE, the earliest of these three interactions must necessarily be an absorption. To proceed further, one then needs to distinguish whether these interactions occur on the ket or the bra. Since the last interaction is fixed to time *t*, the former case correspondingly allows only one possible loop diagram, which is displayed in panel (b) of Fig. 4.

When placing the three integrations on the right strand, one has to allow for all possible orderings of their associated times *τ*, *τ*_{1} and *τ*_{2}. At this point we reiterate that when drawing a loop diagram, even though the relative time ordering of the ket and the bra interactions is unspecified, the temporal order within each strand is fixed. Hence having three interactions with the bra, generates 4 loop diagrams, which arise upon performing permutions of the arrows on the right side. Panels (e) – (h) of Fig. 4 show these diagrams. Note that in this case the loop diagrams are actually fully time ordered.

We now turn to the diagrams with two left and two right interactions. Since the last interaction occurs on the ket at time *t*, for a system starting off in the ground state, the other interaction on the ket also must be an absorption. Finally, we note that when calculating the correlation functions Eq. (52), one needs to take a trace in the end, which restricts both interactions on the bra to be absorptions as well. Arguing along this line leads to the two loop diagrams shown in panels (c) and (d) of Fig. 4. The two loops again take into account the two possible permutations of placing interactions on the bra.

In summary, the pump-probe signal can be represented by eight loop diagrams, seven stem from the first correlation function in Eq. (52), whereas the remaining one is SLE-type. The corresponding expression for the signal derived using the time-domain rules is given in appendix B.

For the sake of completness, we shall compare the QED signal with the semiclassical result, where the pump-probe signal is calculated as a third order response in the direction +**k**_{1} − **k**_{1} + **k**_{2}. in Note that the diagrams in panels (b) – (h) in Fig. 4, correspond to +**k**_{1} − **k _{1}** −

$${S}_{\text{PP}}({\omega}_{1},{\omega}_{2})=-\frac{4}{{\hslash}^{3}}{\mid {\mathcal{E}}_{1,0}\mid}^{2}{\mid {\mathcal{E}}_{2,0}\mid}^{2}Im\{-\langle \widehat{V}{\widehat{G}}^{\u2020}({\omega}_{g}+{\omega}_{1}){\widehat{V}}^{\u2020}{\widehat{G}}^{\u2020}({\omega}_{g}+{\omega}_{1}-{\omega}_{2})\widehat{V}\widehat{G}({\omega}_{g}+{\omega}_{1}){\widehat{V}}^{\u2020}\rangle +\langle \widehat{V}{G}^{\u2020}({\omega}_{g}+{\omega}_{2}){\widehat{V}}^{\u2020}\widehat{G}({\omega}_{g}+{\omega}_{1}-{\omega}_{1})\widehat{V}\widehat{G}({\omega}_{g}+{\omega}_{1}){\widehat{V}}^{\u2020}\rangle +\langle \widehat{V}{\widehat{G}}^{\u2020}({\omega}_{g}+{\omega}_{1})\widehat{V}{\widehat{G}}^{\u2020}({\omega}_{g}+{\omega}_{1}+{\omega}_{2}){\widehat{V}}^{\u2020}\widehat{G}({\omega}_{g}+{\omega}_{1}){\widehat{V}}^{\u2020}\rangle +\langle \widehat{V}{\widehat{G}}^{\u2020}({\omega}_{g}+{\omega}_{2})\widehat{V}{\widehat{G}}^{\u2020}({\omega}_{g}+{\omega}_{1}+{\omega}_{2}){\widehat{V}}^{\u2020}\widehat{G}({\omega}_{g}+{\omega}_{1}){\widehat{V}}^{\u2020}\rangle +\langle \widehat{V}{\widehat{G}}^{\u2020}({\omega}_{g}+{\omega}_{1}){\widehat{V}}^{\u2020}{\widehat{G}}^{\u2020}({\omega}_{g}+{\omega}_{2}-{\omega}_{2})\widehat{V}{\widehat{G}}^{\u2020}({\omega}_{g}+{\omega}_{2}){\widehat{V}}^{\u2020}\rangle +\langle \widehat{V}{\widehat{G}}^{\u2020}({\omega}_{g}+{\omega}_{2}){\widehat{V}}^{\u2020}{\widehat{G}}^{\u2020}({\omega}_{g}+{\omega}_{2}-{\omega}_{1})\widehat{V}{\widehat{G}}^{\u2020}({\omega}_{g}+{\omega}_{2}){\widehat{V}}^{\u2020}\rangle +\langle \widehat{V}{\widehat{G}}^{\u2020}({\omega}_{g}+{\omega}_{1})\widehat{V}{\widehat{G}}^{\u2020}({\omega}_{g}+{\omega}_{1}+{\omega}_{2}){\widehat{V}}^{\u2020}\widehat{G}({\omega}_{g}+{\omega}_{2}){\widehat{V}}^{\u2020}\rangle +\langle \widehat{V}{\widehat{G}}^{\u2020}({\omega}_{g}+{\omega}_{2})\widehat{V}{\widehat{G}}^{\u2020}({\omega}_{g}+{\omega}_{1}+{\omega}_{2}){\widehat{V}}^{\u2020}{\widehat{G}}^{\u2020}({\omega}_{g}+{\omega}_{2}){\widehat{V}}^{\u2020}\rangle \}.$$

(53)

The eight terms correspond respectively to the eight diagrams in Fig. 4.

So far we have focused on incoherent processes. We could thus consider a single molecule and simply multiply the signal by *N* in the end. To describe coherent signals we consider in Eq. (14) as the dipole moment of a collection of noninteracting molecules at positions **r*** _{α}*, i.e.

$$\widehat{V}=\sum _{\alpha}{\widehat{V}}_{\alpha}\delta (\mathbf{r}-{\mathbf{r}}_{\alpha}),$$

(54)

where * _{α}* now denotes the dipole operator of a single molecule. Subsequent interactions of the detecting field can take place with the same (

$${S}_{\text{T}}=\frac{4\pi}{\mathrm{\Omega}}{\omega}_{s}Re\left\{{\int}_{-\infty}^{+\infty}\text{d}t{\int}_{-\infty}^{t}\text{d}\tau {\sum}_{\alpha ,\beta}{\text{e}}^{-i{\omega}_{s}(t-\tau )}{\text{e}}^{i{\mathbf{k}}_{s}\xb7({\mathbf{r}}_{\beta}-{\mathbf{r}}_{\alpha})}{\left(\mathcal{T}{\widehat{V}}_{\alpha ,L}(\tau ){\widehat{V}}_{\beta ,R}^{\u2020}(t)\right)}_{\text{D}}\right\},$$

(55)

where we have further added a time-integration to describe a frequency-domain experiment [37]. Equation (55) is a general expression for spontaneously generated signals. The *α* = *β* and *α* ≠ *β* terms in Eq. (55) give rise to incoherent (*S _{I}*) and coherent (

$${S}_{\text{T}}({\omega}_{s})=N\xb7{S}_{I}({\omega}_{s})+N(N-1)\xb7{S}_{C}({\omega}_{s}),$$

(56)

where the incoherent signal is given by

$${S}_{I}({\omega}_{s})=\frac{4\pi}{\mathrm{\Omega}}{\omega}_{s}Re\left\{{\int}_{-\infty}^{+\infty}\text{d}t{\int}_{-\infty}^{t}\text{d}\tau {\text{e}}^{-i{\omega}_{s}(t-\tau )}{\left(\mathcal{T}{\widehat{V}}_{\alpha ,L}(\tau ){\widehat{V}}_{\alpha ,R}^{\u2020}(t)\right)}_{\text{D}}\right\}.$$

(57)

For uncorrelated particles, we have [* _{α}, _{β}*] = 0. The

$$\begin{array}{l}{\left({\widehat{V}}_{\alpha ,L}(\tau ){\widehat{V}}_{\beta ,R}^{\u2020}(t)\right)}_{\text{D}}={\left({\widehat{V}}_{\alpha ,L}(\tau )\right)}_{\text{D}}{\left({\widehat{V}}_{\beta ,R}^{\u2020}(t)\right)}_{\text{D}}\\ ={\left({\widehat{V}}_{L}(\tau )\right)}_{\text{D}}{\left({\widehat{V}}_{L}^{\u2020}(t)\right)}_{\text{D}}.\end{array}$$

(58)

This gives rise to the coherent term in Eq. (56),

$$\begin{array}{l}{S}_{C}({\omega}_{s})=\frac{2\pi}{\mathrm{\Omega}}{\omega}_{s}{\mid P({\omega}_{s})\mid}^{2},\\ P({\omega}_{s})={\int}_{-\infty}^{+\infty}\text{d}\tau {\text{e}}^{i{\omega}_{s}\tau}{\left({\widehat{V}}_{L}^{\u2020}(\tau )\right)}_{\text{D}}.\end{array}$$

(59)

For this we have also made use of

$${\int}_{-\infty}^{+\infty}\text{d}t{\int}_{-\infty}^{t}\text{d}\tau =\frac{1}{2}{\int}_{-\infty}^{+\infty}\text{d}t{\int}_{-\infty}^{+\infty}\text{d}\tau ,$$

(60)

in Eq. (57).

In Eq. (59) we assumed that the sample is much smaller than the optical wavelength or that we have exact phase matching Δ**k** = 0, where Δ**k** is the difference between the wavevector of the signal mode and the sum of wavevectors of the *n* incoming fields. More generally, we should replace the *N*(*N* − 1) factor in Eq. (56) by

$$F(\mathrm{\Delta}\mathbf{k})=\sum _{\alpha ,\beta}{\text{e}}^{-i\mathrm{\Delta}\mathbf{k}\xb7({\mathbf{r}}_{\alpha}-{\mathbf{r}}_{\beta})}.$$

(61)

For macroscopic samples, the number of terms with *α* ≠ *β* in the double sum of Eq. (61) will by far exceed the terms with equal indices. Evaluating this double sum in the continuous limit and letting the sample volume go to infinity, gives the standard phase matching condition *F*(Δ**k**) → *δ*(Δ**k**), which yields a directed signal.

For small samples both coherent and incoherent terms need to be considered. It is interesting to note that for a single molecule we never get *S _{C}*, we only have

Below we compare the expressions for
${S}_{I}^{(n)}$ and *P*^{(}^{n}^{)} obtained by expanding Eqs. (57) and (59) for *n* = 1, 2, 3, where *n* denotes the number of incoming modes. We demonstrate that the incoherent signal for *n* incoming modes, relates to a 2(*n* + 1)-point correlation function in the dipole moment, compared to an (*n* + 1)-point correlation function contained in *P*^{(}^{n}^{)}.

For *n* = 1

$${P}^{(1)}({\omega}_{s})=\int \text{d}t{\text{e}}^{i{\omega}_{s}t}{\int}_{-\infty}^{t}\text{d}\tau \langle \mathcal{T}{\widehat{V}}_{L}(t){\widehat{V}}_{-}^{\u2020}(\tau )\rangle \mathcal{E}(\text{r},\tau ),$$

gives Rayleigh scattering, whereas

$${S}_{I}^{(1)}({\omega}_{s})=\int \text{d}t{\int}_{-\infty}^{t}\text{d}\tau {\text{e}}^{i{\omega}_{s}(t-\tau )}\int \text{d}{\tau}_{1}\int \text{d}{\tau}_{2}\langle \mathcal{T}{\widehat{V}}_{L}(t){\widehat{V}}_{R}^{\u2020}(\tau ){\widehat{V}}_{-}({\tau}_{1}){\widehat{V}}_{-}^{\u2020}({\tau}_{2})\rangle {\mathcal{E}}^{\ast}(\mathbf{r},{\tau}_{1})\mathcal{E}(\mathbf{r},{\tau}_{2}),$$

is responsible for SLE (Raman and Fluorescence) as discussed in Sec. III.

For *n* = 2

$${P}^{(2)}({\omega}_{s})=\int \text{d}t{\text{e}}^{i{\omega}_{s}t}{\int}_{-\infty}^{t}\text{d}\tau \int \text{d}{\tau}_{1}\langle \mathcal{T}{\widehat{V}}_{L}(t){\widehat{V}}_{-}^{\u2020}(\tau ){\widehat{V}}_{-}({\tau}_{1})\rangle \mathcal{E}(\text{r},\tau ){\mathcal{E}}^{\ast}(\mathbf{r},{\tau}_{1})$$

describes 3-wave mixing e.g. second harmonic generation [41]–[43], and

$$\begin{array}{c}{S}_{I}^{(2)}({\omega}_{s})=\int \text{d}t{\int}_{-\infty}^{t}\text{d}\tau {\text{e}}^{i{\omega}_{s}(t-\tau )}\int \text{d}{\tau}_{1}\int \text{d}{\tau}_{2}\int \text{d}{\tau}_{3}\int \text{d}{\tau}_{4}{\mathcal{E}}^{\ast}(\mathbf{r},{\tau}_{1})\mathcal{E}(\mathbf{r},{\tau}_{2}){\mathcal{E}}^{\ast}(\mathbf{r},{\tau}_{3})\mathcal{E}(\mathbf{r},{\tau}_{4})\\ \langle \mathcal{T}{\widehat{V}}_{L}(t){\widehat{V}}_{R}^{\u2020}(\tau ){\widehat{V}}_{-}({\tau}_{1}){\widehat{V}}_{-}^{\u2020}({\tau}_{2}){\widehat{V}}_{-}({\tau}_{3}){\widehat{V}}_{-}^{\u2020}({\tau}_{4})\rangle \end{array}$$

may represent many possible processes, depending on the level scheme and off-resonant detunings, e.g. two photon fluorescence [38]–[40].

For *n* = 3

$$\begin{array}{c}{P}^{(3)}({\omega}_{s})=\int \text{d}t{\text{e}}^{i{\omega}_{s}t}{\int}_{-\infty}^{t}\text{d}\tau \int \text{d}{\tau}_{1}\int \text{d}{\tau}_{2}\mathcal{E}(\text{r},\tau )\mathcal{E}(\mathbf{r},{\tau}_{1}){\mathcal{E}}^{\ast}(\mathbf{r},{\tau}_{1})\\ \langle \mathcal{T}{\widehat{V}}_{L}(t){\widehat{V}}_{-}^{\u2020}(\tau ){\widehat{V}}_{-}({\tau}_{1}){\widehat{V}}_{-}^{\u2020}({\tau}_{2})\rangle \end{array}$$

represents 4-wave mixing e.g. third harmonic generation and

$$\begin{array}{c}{S}_{I}^{(3)}({\omega}_{s})=\int \text{d}t{\int}_{-\infty}^{t}\text{d}\tau {\text{e}}^{i{\omega}_{s}(t-\tau )}\int \text{d}{\tau}_{1}\int \text{d}{\tau}_{2}\int \text{d}{\tau}_{3}\int \text{d}{\tau}_{4}\int \text{d}{\tau}_{5}\int \text{d}{\tau}_{6}{\mathcal{E}}^{\ast}(\mathbf{r},{\tau}_{1})\mathcal{E}(\mathbf{r},{\tau}_{2}){\mathcal{E}}^{\ast}(\mathbf{r},{\tau}_{3})\\ \mathcal{E}(\mathbf{r},{\tau}_{4}){\mathcal{E}}^{\ast}(\mathbf{r},{\tau}_{5})\mathcal{E}(\mathbf{r},{\tau}_{6})\langle \mathcal{T}{\widehat{V}}_{L}(t){\widehat{V}}_{R}^{\u2020}(\tau ){\widehat{V}}_{-}({\tau}_{1}){\widehat{V}}_{-}^{\u2020}({\tau}_{2}){\widehat{V}}_{-}({\tau}_{3}){\widehat{V}}_{-}^{\u2020}({\tau}_{4}){\widehat{V}}_{-}({\tau}_{4}){\widehat{V}}_{-}^{\u2020}({\tau}_{6})\rangle \end{array}$$

can represent many possible processes, e.g. three photon fluorescence.

A nanocrystal is made out of many unit cells. Its dipole operator is given by = Σ* _{α}^{α}*, where

When detecting nonlinear optical signals in macroscopic samples, one can either measure the intensity (homodyne detection) or the magnitude and phase of the signal field (heterodyne detection). The later technique uses an intense external field , usually referred to as “local oscillator”, which interferes with the created nonlinear signal in the same direction, but is assumed not to interact with the molecule. Using the standard semiclassical approach to nonlinear optics, the quantity detected in a heterodyne measurement is [3]

$${S}_{\mathit{HET}}({\mathbf{k}}_{s})\sim {\int}_{-\infty}^{+\infty}\text{d}tIm\{{\mathcal{E}}_{LO}^{\ast}(t){P}_{s}({\mathbf{k}}_{s},t)\},$$

(62)

where *P _{s}* is the polarization in the signal direction

Equation (55) describes the signal generated from an assembly of particles by spontaneous emission. Applying Eq. (54) to Eq. (29), we can give a similar expression for the stimulated signal

$${S}_{\text{stim}}(t)=-\frac{2}{\hslash}Re{\int}_{-\infty}^{t}\text{d}\tau \left[{\left(\mathcal{T}{\widehat{V}}_{\alpha ,R}(\tau ){\widehat{V}}_{\beta ,L}^{\u2020}(t)\right)}_{\text{D}}-{\left(\mathcal{T}{\widehat{V}}_{\alpha ,L}(\tau ){\widehat{V}}_{\beta ,R}^{\u2020}(t)\right)}_{\text{D}}\right]\mathcal{E}(\mathbf{r},t){\mathcal{E}}^{\ast}(\mathbf{r},\tau ).$$

(63)

In analogy to Sec. (V), we can now split Eq. (63) into a coherent (*α* ≠ *β*) and an incoherent (*α* = *β*) contribution. Since the dipole operators of different, uncorrelated particles *α* ≠ *β* commute, the coherent contribution vanishes identically. The present formalism hence yields the general result that stimulated processes are always incoherent in nature. In particular this applies to heterodyne experiments where the local oscillator serves as the stimulating field. It is hence sufficient to look only at a single particle and use Eq. (14) derived in Sec. II.

Equation (14) is reminiscent of the expression for heterodyne detection (Eq. (62)), if one only looks at one particular mode “s” in the definition of the signal (see Eqs. (11) and (12)), thus

$${S}_{\text{SWM}}(t)=-\frac{2}{\hslash}Im\left[{\left({\widehat{\mathcal{E}}}_{s}({\mathbf{r}}_{\alpha},t){\widehat{V}}^{\u2020}\right)}_{\text{D},\text{T}}\right].$$

(64)

Here, by adding a subscript *α*, we make explicit that Eq. (64) is the signal generated by *one* molecule at position **r*** _{α}*. For an assembly of noninteracting particles, this has then to be summed over

Heterodyne detection can be viewed as an incoherent stimulated emission process in the detected mode. Since Eq. (64) is already first order in the mode *s*, in a subsequent perturbative expansion of the density operator the detecting mode will only contribute to higher order. Hence, to first order in the coupling with the detecting mode, the time-dependence of the density operator is only given by all other modes *s*′ ≠ *s*. The field operator therfore acts directly on the state of the system (molecule + field) at *t* → −∞, which when assuming a coherent state of the field in the same limit (see Eq. (7)) yields

$$S(t)=\frac{2}{\hslash}Im\left[{\mathcal{E}}_{s}^{\ast}({\mathbf{r}}_{\alpha},t)\phantom{\rule{0.16667em}{0ex}}{\left(\widehat{V}\right)}_{\text{D}}\right].$$

(65)

In Eq. (65) we additionally took the complex conjugate within the imaginary part.

The expectation value in Eq. (65) can be further expanded in the incoming modes. For *n* modes this gives the classical *n*th order polarization

$${P}^{(n)}(\mathbf{r},t)={P}^{(n)}(t){\text{e}}^{i\left({\sum}_{{s}^{\prime}\ne s}{\mathbf{k}}_{{s}^{\prime}}\right)\xb7\mathbf{r}}.$$

(66)

Consequently, integrating over all times, we obtain the following generic result for an (*n*+1)-wave mixing process (i.e. 1 detected mode “s” + *n* incoming modes “*s*′ ≠ *s*”)

$$\begin{array}{l}S(\mathrm{\Delta}\mathbf{k})=\frac{2}{\hslash}{\int}_{-\infty}^{+\infty}\text{d}t\text{Im}\left[{\mathcal{E}}_{s,0}^{\ast}{P}^{(n)}(t)f(\mathrm{\Delta}\mathbf{k})\right],\\ f(\mathrm{\Delta}\mathbf{k})=\sum _{\alpha}{\text{e}}^{-i\mathrm{\Delta}\mathbf{k}\xb7{\mathbf{r}}_{\alpha}},\end{array}$$

(67)

where we explicitly display the spatial dependence of the fields, e.g.

$${\mathcal{E}}_{s}(\mathbf{r},t)={\mathcal{E}}_{s,0}(t){\text{e}}^{i{\mathbf{k}}_{s}\xb7\mathbf{r}}.$$

(68)

In Eq. (67) we have added the sum over **r*** _{α}* to make the transition from the signal generated by single molecule to one of several molecules (incoherent process).

Under exact phase-matching conditions (Δ**k** = 0), or if the sample is much smaller than the optical wavelength, we have *f*(Δ**k**) = *N*. Rather than creating a coherent signal, propagating it and interfering as done in the semiclassical approach, we can simply describe it microscopically as an incoherent signal.

Equation (67), hence allows us to interpret heterodyne detection as follows: *All* fields (the *n* fields (*ω _{s}*

A microscopic QED treatment of nonlinear optical processes induced by a weak quantized field is developed in this paper. In the traditional, widely-used approaches to nonlinear optical response, the matter-field interaction is treated semiclassically, hence the signal field is obtained by simultaneously solving macroscopic (Maxwell-) and microscopic (quantum Liouville-) equations. Here we treat the entire process as a single event.

Based on a microscopic definition of the spontaneous and stimulated signal, general expressions are derived for the nonlinear optical signal involving non-equilibrium Green’s functions in the incoming field. Depending on the experiment under consideration, these can then be expanded order by order in the incoming field modes. Application is made to spontaneous light emission and pump-probe spectroscopy. Coherent and incoherent processes involving both spontaneous and stimulated emission are treated using a unified framework.

A diagrammatic derivation of the contributing terms within the RWA using the technique of Keldysh Schwinger loop is developed. This respresentation reflects the partially time-ordered nature of the Green’s functions and hence yields more compact expressions than the well established fully time ordered double sided Feynman diagrams. The loop diagrams are particularly useful for frequency-domain measurements where the bookeeping of time ordering is not necessary anymore. If needed for time-domain experiments, the fully time ordered expressions can be obtained from the loop diagrams using a simple prescription. Rules for constructing and reading these diagrams are laid out, making it possible to intuitively derive the expressions for the signals.

The practical merits of this partially time-ordered approach are illustrated for the pump-probe technique where the number of diagrams is reduced from 16 to 8. This example also demonstrates an interesting fundamental merit of the QED treatment of nonlinear optical processes, even though in pump-probe spectroscopy quantum effects of the field (i.e. terms stemming from spontaneous emission) are usually negligible: Semiclassically, the pump-probe signal is obtained by calculating a third order response induced by two fields in the direction **k**_{1} − **k**_{1} +**k**_{2}. The resulting signal field is then obtained using heterodyne detection, where mode 2 acts as its own local oscillator field. Even though the final expression for the signal is identical to a QED result, the semiclassical approach creats an artificial asymmetry between the local oscillator (mode 2), which by definition does not interact with the molecule, and the remaining fields (mode 1 and 2), which give rise to the nonlinear polarization. This asymmetry, which is eliminated here, is a direct consequence of the limitations imposed by separately treating matter and field within the semiclassical approach.

Finally, using the QED formalism, we are able to treat both coherent and incoherent processes in a unified way. This distinction comes about by calculating the signal for an assembly of molecules rather than for a single particle. The signal then splits into two parts, a coherent term scaling ~ *N*(*N* − 1) and an incoherent term ~*N*, where *N* denotes the number of particles. As the number of particles in the sample is increased, the coherent term becomes more pronounced. For experiments carried out with only a few molecules however, both coherent and incoherent terms will generally make comparable contributions to the signal. In single molecule experiments only the incoherent term survives.

Since incoherent signals may not be calculated semiclassically, this leads to a striking consequence. For an (*n* +1)-wave mixing processes carried out on a single molecule we show that the signal is not given by an (*n* + 1)-point response function of the dipole operator related to *P*^{(}^{n}^{)}, as predicted by the semiclassical approach. It is rather given by a different 2(*n* + 1)-point combination of correlation functions. This doubling arises since the signal may not be recast as an amplitude square and we must calculate the signal itself, not an amplitude. The present formalism allows a microscopic calculation of nonlinear optical experiments on single molecules which have only become feasible recently.

Finally, we address this apparent limitation of the semiclassical description from a more general viewpoint and show that heterodyne detection for an (*n*+1)-wave mixing experiment can be simply viewed as an incoherent stimulated emission process in the detected mode. In contrast, homodyne-detected *n*-wave mixing is a coherent spontaneous emission process. For a macroscopic sample, by summing over the molecules in the interaction volume, we recover the phase-matching condition, i.e. the heterodyne-signal is finite only if the stimulating mode (i.e. the local oscillator) matches this condition.

We wish to thank Prof. Sunney Xie for useful discussions. The support of the National Science Foundation (Grant No. CHE-0446555, CBC-0533162), NIRT (Grant No. EEC 0303389)) and the National Institutes of Health (GM59230) is gratefully acknowledged.

In this appendix we illustrate the equivalence of the QED expression for the pump-probe signal with the classical signal expressed in terms of *χ*^{(3)}. Starting from Eq. (50) and applying the definition of Eq. (16), we can write

$${\mathcal{H}}_{\mathit{int}-}(\tau )=\frac{1}{\sqrt{2}}\left[{\widehat{\mathcal{E}}}_{+}(\mathbf{r},\tau ){\widehat{V}}_{-}^{\u2020}+{\widehat{\mathcal{E}}}_{-}(\mathbf{r},\tau ){\widehat{V}}_{+}^{\u2020}+{\widehat{\mathcal{E}}}_{+}^{\u2020}(\mathbf{r},\tau ){\widehat{V}}_{-}+{\widehat{\mathcal{E}}}_{-}^{\u2020}(\mathbf{r},\tau ){\widehat{V}}_{+}\right].$$

(A1)

Making use of Eq. (A1), the correlation function in Eq. (50) can be recast in terms of correlations of the system and the field. Calculating the field correlations, only 6 of the possible 16 terms are nonzero, giving

$$2\langle \mathcal{T}{\widehat{V}}_{+}^{\u2020}(t){\widehat{V}}_{-}(\tau ){{\mathcal{H}}^{\prime}}_{\mathit{int}-}({\tau}_{1}){{\mathcal{H}}^{\prime}}_{\mathit{int}-}({\tau}_{2})\rangle =\langle \mathcal{T}{\widehat{V}}_{+}^{\u2020}(t){\widehat{V}}_{-}(\tau ){\widehat{V}}_{-}^{\u2020}({\tau}_{1}){\widehat{V}}_{-}({\tau}_{2})\rangle \langle \mathcal{T}{\widehat{\mathcal{E}}}_{+}(\mathbf{r},{\tau}_{1}){\widehat{\mathcal{E}}}_{+}^{\u2020}(\mathbf{r},{\tau}_{2})\rangle +\langle \mathcal{T}{\widehat{V}}_{+}^{\u2020}(t){\widehat{V}}_{-}(\tau ){\widehat{V}}_{-}({\tau}_{1}){\widehat{V}}_{-}^{\u2020}({\tau}_{2})\rangle \langle \mathcal{T}{\widehat{\mathcal{E}}}_{+}^{\u2020}(\mathbf{r},{\tau}_{1}){\widehat{\mathcal{E}}}_{+}(\mathbf{r},{\tau}_{2})\rangle +\langle \mathcal{T}{\widehat{V}}_{+}^{\u2020}(t){\widehat{V}}_{-}(\tau ){\widehat{V}}_{-}^{\u2020}({\tau}_{1}){\widehat{V}}_{+}({\tau}_{2})\rangle \langle \mathcal{T}{\widehat{\mathcal{E}}}_{+}(\mathbf{r},{\tau}_{1}){\widehat{\mathcal{E}}}_{-}^{\u2020}(\mathbf{r},{\tau}_{2})\rangle +\langle \mathcal{T}{\widehat{V}}_{+}^{\u2020}(t){\widehat{V}}_{-}(\tau ){\widehat{V}}_{+}({\tau}_{1}){\widehat{V}}_{-}^{\u2020}({\tau}_{2})\rangle \langle \mathcal{T}{\widehat{\mathcal{E}}}_{-}^{\u2020}(\mathbf{r},{\tau}_{1}){\widehat{\mathcal{E}}}_{+}(\mathbf{r},{\tau}_{2})\rangle +\langle \mathcal{T}{\widehat{V}}_{+}^{\u2020}(t){\widehat{V}}_{-}(\tau ){\widehat{V}}_{-}({\tau}_{1}){\widehat{V}}_{+}^{\u2020}({\tau}_{2})\rangle \langle \mathcal{T}{\widehat{\mathcal{E}}}_{+}^{\u2020}(\mathbf{r},{\tau}_{1}){\widehat{\mathcal{E}}}_{-}(\mathbf{r},{\tau}_{2})\rangle +\langle \mathcal{T}{\widehat{V}}_{+}^{\u2020}(t){\widehat{V}}_{-}(\tau ){\widehat{V}}_{+}^{\u2020}({\tau}_{1}){\widehat{V}}_{-}({\tau}_{2})\rangle \langle \mathcal{T}{\widehat{\mathcal{E}}}_{-}(\mathbf{r},{\tau}_{1}){\widehat{\mathcal{E}}}_{+}^{\u2020}(\mathbf{r},{\tau}_{2})\rangle ,$$

(A2)

where the system is assumed to be in its ground state for *t* → −∞, hence

$$\langle \mathcal{T}{V}_{+}^{\u2020}(t){\widehat{V}}_{-}(\tau ){\widehat{V}}_{-}({\tau}_{1}){\widehat{V}}_{-}({\tau}_{2})\rangle =0.$$

(A3)

Note that since Eq. (50) is symmetric with respect to the variables *τ*_{1} and *τ*_{2}, Eq. (A2) simplifies further and only includes three distinct terms upon integration, i.e.

$$\int \text{d}{\tau}_{1}\int \text{d}{\tau}_{2}2\langle \mathcal{T}{\widehat{V}}_{R}(\tau ){\widehat{V}}_{L}^{\u2020}(t){{\mathcal{H}}^{\prime}}_{\mathit{int}-}(\tau ){{\mathcal{H}}^{\prime}}_{\mathit{int}-}({\tau}_{2})\rangle =2\int \text{d}{\tau}_{1}\int \text{d}{\tau}_{2}\times \langle \mathcal{T}{\widehat{V}}_{+}^{\u2020}(t){\widehat{V}}_{-}(\tau ){\widehat{V}}_{-}^{\u2020}({\tau}_{1}){\widehat{V}}_{-}({\tau}_{2})\rangle \langle \mathcal{T}{\widehat{\mathcal{E}}}_{+}(\mathbf{r},{\tau}_{1}){\widehat{\mathcal{E}}}_{+}^{\u2020}(\mathbf{r},{\tau}_{2})\rangle +\langle \mathcal{T}{\widehat{V}}_{+}^{\u2020}(t){\widehat{V}}_{-}(\tau ){\widehat{V}}_{-}^{\u2020}({\tau}_{1}){\widehat{V}}_{+}({\tau}_{2})\rangle \langle \mathcal{T}{\widehat{\mathcal{E}}}_{+}(\mathbf{r},{\tau}_{1}){\widehat{\mathcal{E}}}_{-}^{\u2020}(\mathbf{r},{\tau}_{2})\rangle +\langle \mathcal{T}{\widehat{V}}_{+}^{\u2020}(t){\widehat{V}}_{-}(\tau ){\widehat{V}}_{-}({\tau}_{1}){\widehat{V}}_{+}^{\u2020}({\tau}_{2})\rangle \langle \mathcal{T}{\widehat{\mathcal{E}}}_{+}^{\u2020}(\mathbf{r},{\tau}_{1}){\widehat{\mathcal{E}}}_{-}(\mathbf{r},{\tau}_{2})\rangle .$$

(A4)

For a classical incoming field, the correlation function in Eq. (A2) is equal to only the first term on the rhs, where

$$\langle \mathcal{T}{\widehat{\mathcal{E}}}_{+}(\mathbf{r},{\tau}_{1}){\widehat{\mathcal{E}}}_{+}^{\u2020}(\mathbf{r},{\tau}_{2})\rangle =\mathcal{E}(\mathbf{r},{\tau}_{1}){\mathcal{E}}^{\ast}(\mathbf{r},{\tau}_{2})+\frac{2\pi \hslash \omega}{\mathrm{\Omega}}{\text{e}}^{i\omega ({\tau}_{1}-{\tau}_{2})}.$$

(A5)

For a classical field the second term in Eq. (A5) will be negligible.

The last two terms in Eq. (A4) are due to the quantum character of the incoming field, since

$$\begin{array}{l}\langle \mathcal{T}{\widehat{\mathcal{E}}}_{+}(\mathbf{r},{\tau}_{1}){\widehat{\mathcal{E}}}_{-}^{\u2020}(\mathbf{r},{\tau}_{2})\rangle =\frac{2\pi \hslash \omega}{\mathrm{\Omega}}{\text{e}}^{-i\omega ({\tau}_{1}-{\tau}_{2})}\theta ({\tau}_{1}-{\tau}_{2}),\\ \langle \mathcal{T}{\widehat{\mathcal{E}}}_{+}^{\u2020}(\mathbf{r},{\tau}_{1}){\widehat{\mathcal{E}}}_{-}(\mathbf{r},{\tau}_{2})\rangle =\frac{2\pi \hslash \omega}{\mathrm{\Omega}}{\text{e}}^{-i\omega ({\tau}_{2}-{\tau}_{1})}\theta ({\tau}_{1}-{\tau}_{2}).\end{array}$$

(A6)

Hence, neglecting the last two terms in Eq. (A4) and making use of Eq. (A5), we obtain Eq. (51).

In this appendix we give the expressions for the pump-probe signal in the time-domain obtained by applying the rules summarized in Sec. III to the loop diagrams in Fig. 4. This complements the frequency-domain expression given in Eq. (53). Applying the time-domain rules for the loop diagrams in Fig. 4, we obtain

$${S}_{\text{PP}}(t)=\frac{4}{{\hslash}^{3}}Re\{{\int}_{-\infty}^{t}\text{d}\tau \int \text{d}{\tau}_{1}\int \text{d}{\tau}_{2}{\mathcal{E}}_{2}^{\ast}(\mathbf{r},t){\mathcal{E}}_{2}(\mathbf{r},\tau )\phantom{\rule{0ex}{0ex}}\langle \mathcal{T}{V}_{L}(t){V}_{R}^{\u2020}(\tau ){V}_{L}^{\u2020}({\tau}_{2}){V}_{R}({\tau}_{1})\rangle \langle \mathcal{T}{\widehat{\mathcal{E}}}_{1,L}(\mathbf{r},{\tau}_{2}){\widehat{\mathcal{E}}}_{1,R}^{\u2020}(\mathbf{r},{\tau}_{1})\rangle \phantom{\rule{0ex}{0ex}}+{\int}_{-\infty}^{t}\text{d}\tau \int \text{d}{\tau}_{1}\int \text{d}{\tau}_{2}{\mathcal{E}}_{2}(\mathbf{r},t){\mathcal{E}}_{2}^{\ast}(\mathbf{r},\tau )\phantom{\rule{0ex}{0ex}}+[\langle \mathcal{T}{V}_{L}^{\u2020}(t){V}_{R}(\tau ){V}_{1}({\tau}_{1}){V}_{L}^{\u2020}({\tau}_{2})\rangle \langle \mathcal{T}{\widehat{\mathcal{E}}}_{1,L}^{\u2020}(\mathbf{r},{\tau}_{1}){\widehat{\mathcal{E}}}_{1,L}(\mathbf{r},{\tau}_{2})\rangle \phantom{\rule{0ex}{0ex}}-\langle \mathcal{T}{V}_{L}^{\u2020}(t){V}_{R}(\tau ){V}_{L}^{\u2020}({\tau}_{1}){V}_{R}({\tau}_{2})\rangle \langle \mathcal{T}{\widehat{\mathcal{E}}}_{1,L}(\mathbf{r},{\tau}_{1}){\widehat{\mathcal{E}}}_{1,R}^{\u2020}(\mathbf{r},{\tau}_{2})\rangle \phantom{\rule{0ex}{0ex}}+\langle \mathcal{T}{V}_{L}^{\u2020}(t){V}_{R}(\tau ){V}_{R}^{\u2020}({\tau}_{1}){V}_{R}({\tau}_{2})\rangle \langle \mathcal{T}{\widehat{\mathcal{E}}}_{1,R}(\mathbf{r},{\tau}_{1}){\widehat{\mathcal{E}}}_{1,R}^{\u2020}(\mathbf{r},{\tau}_{2})\rangle ]\},$$

(B1)

where we kept the order of the diagrams in Fig. 4. Similar to Eq. (34), an additional factor 2 accounts for the symmetry with respect to *τ*_{1} and *τ*_{2}. This can again be broken into fully time-ordered terms using the double sided Feynman diagrams displayed in Fig. 5. We then get

$${S}_{\text{PP}}(t)=\frac{4}{{\hslash}^{3}}Re{\int}_{-\infty}^{t}\text{d}\tau \int \text{d}{\tau}_{1}\int \text{d}{\tau}_{2}{\mathcal{E}}_{2}^{\ast}(\mathbf{r},t){\mathcal{E}}_{2}(\mathbf{r},\tau ){\mathcal{E}}_{1}^{\ast}(\mathbf{r},{\tau}_{1}){\mathcal{E}}_{1}(\mathbf{r},{\tau}_{2})[\theta (t\tau )\theta (\tau {\tau}_{1})\theta ({\tau}_{1}{\tau}_{2})\langle {\widehat{V}}_{L}(t){\widehat{V}}_{R}^{\u2020}(\tau ){\widehat{V}}_{R}({\tau}_{1}){\widehat{V}}_{L}^{\u2020}({\tau}_{2})\rangle \phantom{\rule{0ex}{0ex}}+\theta (t{\tau}_{2})\theta ({\tau}_{2}\tau )\theta (\tau {\tau}_{1})\langle {\widehat{V}}_{L}(t){\widehat{V}}_{L}^{\u2020}({\tau}_{2}){\widehat{V}}_{R}^{\u2020}(\tau ){\widehat{V}}_{R}({\tau}_{1})\rangle \phantom{\rule{0ex}{0ex}}+\theta (t{\tau}_{2})\theta ({\tau}_{2}\tau )\theta (\tau {\tau}_{1})\langle {\widehat{V}}_{L}(t){\widehat{V}}_{L}^{\u2020}({\tau}_{2}){\widehat{V}}_{R}^{\u2020}(\tau ){\widehat{V}}_{R}({\tau}_{1})\rangle \phantom{\rule{0ex}{0ex}}+\theta (t\tau )\theta (\tau {\tau}_{2})\theta ({\tau}_{2}{\tau}_{1})\langle {\widehat{V}}_{L}(t){\widehat{V}}_{R}^{\u2020}(\tau ){\widehat{V}}_{L}^{\u2020}({\tau}_{2}){\widehat{V}}_{R}({\tau}_{1})\rangle ]\phantom{\rule{0ex}{0ex}}+{\mathcal{E}}_{2}(\mathbf{r},t){\mathcal{E}}_{2}^{\ast}(\mathbf{r},\tau )\{{\mathcal{E}}_{1}^{\ast}(\mathbf{r},{\tau}_{1}){\mathcal{E}}_{1}(\mathbf{r},{\tau}_{2})[\theta (t\tau )\theta (\tau {\tau}_{1})\theta ({\tau}_{1}{\tau}_{2})\langle {\widehat{V}}_{L}^{\u2020}(t){\widehat{V}}_{R}(\tau ){\widehat{V}}_{L}({\tau}_{1}){\widehat{V}}_{L}^{\u2020}({\tau}_{2})\rangle \phantom{\rule{0ex}{0ex}}+\theta (t{\tau}_{1})\theta ({\tau}_{1}\tau )\theta (\tau {\tau}_{2})\langle {\widehat{V}}_{L}^{\u2020}(t){\widehat{V}}_{L}({\tau}_{1}){\widehat{V}}_{R}(\tau ){\widehat{V}}_{L}^{\u2020}({\tau}_{2})\rangle \phantom{\rule{0ex}{0ex}}+\theta (t{\tau}_{1})\theta ({\tau}_{1}{\tau}_{2})\theta ({\tau}_{2}\tau )\langle {\widehat{V}}_{L}^{\u2020}(t){\widehat{V}}_{L}({\tau}_{1}){\widehat{V}}_{L}^{\u2020}({\tau}_{2}){\widehat{V}}_{R}(\tau )\rangle ]-{\mathcal{E}}_{1}(\mathbf{r},{\tau}_{1}){\mathcal{E}}_{1}^{\ast}(\mathbf{r},{\tau}_{2}))[\theta (t\tau )\theta (\tau {\tau}_{1})\theta ({\tau}_{1}{\tau}_{2})\langle {\widehat{V}}_{L}^{\u2020}(t){\widehat{V}}_{R}(\tau ){\widehat{V}}_{L}^{\u2020}({\tau}_{1}){\widehat{V}}_{R}({\tau}_{2})\rangle \phantom{\rule{0ex}{0ex}}+\theta (t{\tau}_{1})\theta ({\tau}_{1}\tau )\theta (\tau {\tau}_{2})\langle {\widehat{V}}_{L}^{\u2020}(t){\widehat{V}}_{L}^{\u2020}({\tau}_{1}){\widehat{V}}_{R}(\tau ){\widehat{V}}_{R}({\tau}_{2})\rangle \phantom{\rule{0ex}{0ex}}+\theta (t\tau )\theta (\tau {\tau}_{2})\theta ({\tau}_{2}{\tau}_{1})\langle {\widehat{V}}_{L}^{\u2020}(t){\widehat{V}}_{R}(\tau ){\widehat{V}}_{R}({\tau}_{2}){\widehat{V}}_{L}^{\u2020}({\tau}_{1})\rangle \phantom{\rule{0ex}{0ex}}+\theta (t{\tau}_{2})\theta ({\tau}_{2}{\tau}_{1})\theta ({\tau}_{1}\tau )\langle {\widehat{V}}_{L}^{\u2020}(t){\widehat{V}}_{R}^{\u2020}({\tau}_{2}){\widehat{V}}_{L}^{\u2020}({\tau}_{1}){\widehat{V}}_{R}(\tau )\rangle \phantom{\rule{0ex}{0ex}}+\theta (t{\tau}_{1})\theta ({\tau}_{1}{\tau}_{2})\theta ({\tau}_{2}\tau )\langle {\widehat{V}}_{L}^{\u2020}(t){\widehat{V}}_{L}^{\u2020}({\tau}_{1}){\widehat{V}}_{R}({\tau}_{2}){\widehat{V}}_{R}(\tau )\rangle \phantom{\rule{0ex}{0ex}}+\theta (t{\tau}_{2})\theta ({\tau}_{2}\tau )\theta (\tau {\tau}_{1})\langle {\widehat{V}}_{L}^{\u2020}(t){\widehat{V}}_{R}({\tau}_{2}){\widehat{V}}_{R}(\tau ){\widehat{V}}_{L}^{\u2020}({\tau}_{1})\rangle ]\phantom{\rule{0ex}{0ex}}+{\mathcal{E}}_{1}(\mathbf{r},{\tau}_{1}){\mathcal{E}}_{1}^{\ast}(\mathbf{r},{\tau}_{2})[\theta (t\tau )\theta (\tau {\tau}_{1})\theta ({\tau}_{1}{\tau}_{2})\langle {\widehat{V}}_{L}^{\u2020}(t){\widehat{V}}_{R}(\tau ){\widehat{V}}_{R}^{\u2020}({\tau}_{1}){\widehat{V}}_{R}({\tau}_{2})\rangle \phantom{\rule{0ex}{0ex}}+\theta (t{\tau}_{2})\theta ({\tau}_{2}{\tau}_{1})\theta ({\tau}_{1}\tau )\langle {\widehat{V}}_{L}^{\u2020}(t){\widehat{V}}_{R}({\tau}_{2}){\widehat{V}}_{R}^{\u2020}({\tau}_{1}){\widehat{V}}_{R}(\tau )\rangle \phantom{\rule{0ex}{0ex}}+\theta (t{\tau}_{1})\theta ({\tau}_{1}\tau )\theta (\tau {\tau}_{2})\langle {\widehat{V}}_{L}^{\u2020}(t){\widehat{V}}_{R}^{\u2020}({\tau}_{1}){\widehat{V}}_{R}(\tau ){\widehat{V}}_{R}({\tau}_{2})\rangle \phantom{\rule{0ex}{0ex}}+\theta (t{\tau}_{1})\theta ({\tau}_{1}{\tau}_{2})\theta ({\tau}_{2}\tau )\langle {\widehat{V}}_{L}^{\u2020}(t){\widehat{V}}_{R}^{\u2020}({\tau}_{1}){\widehat{V}}_{R}({\tau}_{2}){\widehat{V}}_{R}(\tau )\rangle \left]\right\}.$$

(B2)

Comparison of Eqs. (B1) and (B2) illustrates the tremendous simplification achieved by the more compact partially-time-ordered approach. Employing the diagrammatic Keldysh-Schwinger loop technique, the number of diagrams is reduced by half (16 vs. 8 diagrams).

PACS numbers: 42.65-k, 12.20.-m, 42.62. Fi, 42.50.-p, 42.50.Ct

1. Andrews D, Allcock P. Optical Harmonics in Molecular Systems; Quantum Electrodynamical Theory. Wiley-VCH; Weinheim: 2002.

2. Bloembergen N. Nonlinear optics. Benjamin; New York: 1965.

3. Mukamel S. Principles of Nonlinear Optical Spectroscopy. Oxford University Press; New York: 1995.

4. Scully M, Zubairy M. Quantum Optics. Cambridge University Press; Cambridge: 1997.

5. Hanamura E, Mukamel S. J Opt Soc Am B. 1986;3:1124–1126.

6. Mukamel S. Adv Chem Phys. 1988;70(Part I):165–230.

7. Bidault S, Brasselet S, Zyss J. Optics Letters. 2004;29:1242. [PubMed]

8. Van Dijk MA, Lippitz M, Orrit M. Accounts of chemical research. 2005;38:594–601. [PubMed]

9. Van Dijk MA, Lippitz M, Orrit M. Phys Rev Lett. 2005;95:267406. [PubMed]

10. Sanchez EJ, Novotny L, Holtom GR, Xie XS. Phys Chem A. 1997;101:7019.

11. Muskens OL, Del Fatti N, Vallee F. Nano Letters. 2006;6:552. [PubMed]

12. Langbein W, et al. J Phys¿ Condens Matter. 2007;19:295203. [PubMed]

13. Lukin MD, et al. Science. 2007;316:1312. [PubMed]

14. Potma EO, Xie XS. Chem Phys Chem. 2005;6:77. [PubMed]

15. Haroche S, Raimond JM. Exploring the Quantum: Atoms, Cavities, and Photons. Oxford University Press; Oxford: 2006.

16. Schleich WP. Quantum Optics in Phase Space. Wiley-VCH; Berlin: 2001.

17. Scully M, Lamb WE., Jr Physical Review. 1967;159:159.

18. Sarget M, III, Holm DA, Zubairy MS. Physical Review A. 1985;31:3112. [PubMed]

19. Stenholm S, Holm DA, Sarget M., III Physical Review A. 1985;31:3124. [PubMed]

20. Reid MD, Walls DF. Physical Review A. 1986;34:4929. [PubMed]

21. Agarwal GS. Physical Review A. 1986;34:4055. [PubMed]

22. Raymer MG, Walmsley IA. Progress in Optics. 1990;28:181.

23. Gorbachev VN, Zhiliba AI. J Opt B. 2000;2:764.

24. Panteleev AA, Rerikh VlK, Starostin AN. Journal of Experimental and Theoretical Physics. 2000;90:50.

25. Macovei M, Li GX. Physical Review A. 2007;76:023818.

26. Negele J, Orland H. Quantum Many Particle Physics. Westview Press; 1998.

27. Rammer J. Quantum field theory of non-equilibrium states. Cambridge University Press; 2007.

28. Keldysh LV. Sov Phys JETP. 1965;20:1018.

29. Gorkov LP, Dzyaloshinski IE, Abrikosov AA. Methods of quantum field theory. Prentice Hall; 1963.

30. Mills R. Propagators for many-particle systems; an elementary treatment. New York: Gordon and Breach; 1969.

31. Datta S. Electronic Transport in Mesoscopic Systems. Cambridge University Press; Cambridge: 1995.

32. Haug H, Jauho A-P. Quantum Kinetics in Transport and Optics of Semiconductors. Springer-Verlag; Berlin, Heidelberg: 1996.

33. Harbola U, Maddox JB, Mukamel S. Phys Rev B. 2006;73:205404.

34. Harbola U, Mukamel S. The Journal of Chemical Physics. 2006;124:044106. [PubMed]

35. Harbola U, Maddox JB, Mukamel S. Phys Rev B. 2006;73:075211.

36. Mukamel S. submitted to Phys. Rev. A.

37. Chernyak V, Wang N, Mukamel S. Physics Reports. 1995;263:213–309.

38. Callis PR. Annu Rev Phys Chem. 1997;48:271. [PubMed]

39. Pons T, Mertz J. Optics Communications. 2006;258:203.

40. Rehms AA, Callis PR. Chemical Physics Letters. 1993;208:276.

41. Campagnola PJ, Loew LM. Nature Biotechnology. 2003;21:1356. [PubMed]

42. Moreaux L, Pons T, Dambrin V, Blanchard-Desce M, Mertz J. Optics Letters. 2003;28:625. [PubMed]

43. Naujok RR, Paul HJ, Corn RM. The Journal of Physical Chemistry. 1996;100:10497.

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |