OSA's Digital Library

Optics Express

Optics Express

  • Editor: Michael Duncan
  • Vol. 14, Iss. 23 — Nov. 13, 2006
  • pp: 11418–11432
« Show journal navigation

Noise figure and photon statistics in coherent anti-Stokes Raman scattering

D. Dimitropoulos, D. R. Solli, R. Claps, and B. Jalali  »View Author Affiliations


Optics Express, Vol. 14, Issue 23, pp. 11418-11432 (2006)
http://dx.doi.org/10.1364/OE.14.011418


View Full Text Article

Acrobat PDF (220 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

Coherent anti-Stokes Raman scattering (CARS) is a well-known Raman scattering process that occurs when Stokes, anti-Stokes and pump waves are properly phase-matched. Using a quantum-field approach with Langevin noise sources, we calculate the noise figure for wavelength conversion between the Stokes and anti-Stokes waves in CARS and show its dependence on phase mismatch. Under phase matched conditions, the minimum noise figure is approximately 3 dB, with a correction that depends on the pump frequency, Stokes shift, refractive indices, and nonlinear susceptibilities. We calculate the photon statistics of CARS and show that the photon number distribution is non-Gaussian. Our findings may be significant for currently pursued applications of CARS including wavelength conversion in data transmission and spectroscopic detection of dilute biochemical species.

© 2006 Optical Society of America

1. Introduction

Coherent anti-Stokes Raman scattering (CARS) is a nonlinear optical process in which energy is parametrically transferred between phase-matched Stokes and anti-Stokes fields via interaction with vibrational modes of the material in the presence of a strong pump wave. Since its discovery by Maker and Terhune in the early days of nonlinear optical research [1

1. P.D. Maker and R.W. Terhune, “Study of optical effects due to an induced polarization third order in the electric field strength,” Phys. Rev. 137, A801 (1965). [CrossRef]

], it has found numerous applications in spectroscopy and substance identification [2

2. H. Vogt, “Coherent and hyper-Raman techniques” in Topics in Applied Physics (eds. M. Cardona and G. Guntherodt) vol. 50, p. 207, Springer-Verlag (1982).

,3

3. A.M. Zheltikov and P. Radi, “Non-linear Raman spectroscopy 75 years after the Nobel Prize for the discovery of Raman scattering and 40 years after the first CARS experiments,” J. of Raman Spectrosc. 36, 92–94 (2005). [CrossRef]

] and continues as an active field of applied research [4–9

4. G. Beadie, Z.E. Sariyanni, Y. V. Rostovtsev, T. Opatrny, J. Reintjes, and M.O. Scully, “Towards a FAST CARS anthrax detector: coherence preparation using simultaneous femtosecond laser pulses,” Opt. Commun. 244, 423–430 (2005). [CrossRef]

]. CARS was also recently demonstrated in silicon and has been proposed as means of wavelength conversion between the technologically important 1300 nm and the 1550 nm wavelength bands [10

10. R. Claps, V. Raghunathan, D. Dimitropoulos, and B. Jalali, “Anti-Stokes Raman conversion in Silicon waveguides,” Opt. Express 11, 2862–2872 (2003). [CrossRef] [PubMed]

].

An important topic that has been discussed is the fundamental noise in the CARS process. To our knowledge, the first investigation to address the issue of noise in CARS was conducted by Perina et al [11–13

11. J. Perina, “Photon statistics in Raman scattering with frequency mismatch,” Optica Acta 28, 1529 (1981). [CrossRef]

] with more recent investigations due to Voss et al [14–15

14. P.L. Voss and P. Kumar, “Raman-noise-induced noise-figure limit for χ(3) parametric amplifiers,” Opt. Lett. 29, 445 (2004). [CrossRef] [PubMed]

]. Scully et al [16

16. K.K. Das, G. S. Agarwal, Y.M. Golubev, and M.O. Scully, “Langevin analysis of fundamental noise limits in coherent anti-Stokes Raman spectroscopy,” Phys. Rev. A 71, 013802 (2005). [CrossRef]

] have also addressed this issue in the context of an ultra-fast spectroscopic technique they proposed for detection of rarified bacterial spores [4

4. G. Beadie, Z.E. Sariyanni, Y. V. Rostovtsev, T. Opatrny, J. Reintjes, and M.O. Scully, “Towards a FAST CARS anthrax detector: coherence preparation using simultaneous femtosecond laser pulses,” Opt. Commun. 244, 423–430 (2005). [CrossRef]

]. In References [11–13

11. J. Perina, “Photon statistics in Raman scattering with frequency mismatch,” Optica Acta 28, 1529 (1981). [CrossRef]

], the photon statistics were calculated for the CARS process, but the treatment is in the time-domain. The treatment in [16

16. K.K. Das, G. S. Agarwal, Y.M. Golubev, and M.O. Scully, “Langevin analysis of fundamental noise limits in coherent anti-Stokes Raman spectroscopy,” Phys. Rev. A 71, 013802 (2005). [CrossRef]

] is also in the time-domain. Voss et al used spatial propagation equations for the optical fields and calculated the noise figure for Raman amplification and CARS [14–15

14. P.L. Voss and P. Kumar, “Raman-noise-induced noise-figure limit for χ(3) parametric amplifiers,” Opt. Lett. 29, 445 (2004). [CrossRef] [PubMed]

], but did not calculate the photon statistics.

In this paper, we formulate the problem using spatial propagation equations and obtain both the noise figure and, for the first time, the photon statistics using such an approach. Although the time-domain approach is appropriate for cavity dynamics, this spatial treatment is more suitable for treating waveguide propagation where the effects of waveguide and material dispersion, as well as phase mismatch are explicitly included. We obtain the propagation equations for Stokes and anti-Stokes fields in a different manner than that employed in [14–15

14. P.L. Voss and P. Kumar, “Raman-noise-induced noise-figure limit for χ(3) parametric amplifiers,” Opt. Lett. 29, 445 (2004). [CrossRef] [PubMed]

]. We demonstrate that the minimum noise figure obtained at perfect phase matching is close to 3 dB, but deviates from this value by a factor that depends on the contribution of both linear and nonlinear dispersion, as well as the Stokes shift and pump frequency. There are two distinct contributions to the noise in the output signal: (1) noise originating from the zero point fluctuations of the inputs and (2) noise originating from the coupling of the Stokes and anti-Stokes fields with the damped material vibrations. We also calculate the full photon probability distribution at the output wavelength, and show that its tails deviate significantly from those of the Gaussian distribution. These findings may be important for determining the minimum bit error rate (BER) in CARS data transfer, as well as the fundamental detection limit of CARS spectroscopy.

We start with the classical equations for CARS and convert them into propagation equations for the quantum field operators for the Stokes and anti-Stokes modes. Using these equations, we derive the noise sources for the Stokes and the anti-Stokes fields by imposing the requirement that the commutators of their field operators must be conserved. Throughout the paper we deal with an ideal medium that does not have optical losses or permit other nonlinear optical processes. Optical losses (linear or nonlinear) can be easily included numerically, but are not discussed here because they are not fundamental to the CARS process.

2. Equations of motion

The CARS process increases the noise in both the Stokes and anti-Stokes fields. These modes are coupled together through interaction with the damped vibrational phonons of the material and since CARS involves dissipative interactions, it must introduce noise into the optical signals. The propagation of the Stokes and anti-Stokes field operators (denoted âS and âAS respectively), including the effect of this excess noise (using the Langevin noise operators (x), + (x) discussed below) are described by the following two equations:

dâsdx=gsâs+rgsexp(j(Δβx+2φP))âAS+j2gSN̂+
(1a)
dâAS+dx=r2gSâAS+rgSexp(j(Δβ·x+2φP))âS
+jexp(j(Δβ·x+2φP))2r2gSN̂+.
(1b)

These equations are derived in Appendix A starting from the classical field-amplitude equations of CARS, replacing the field amplitudes with quantum field operators, and invoking the principle of “commutator conservation.” In the above equations the operators that appear obey the following commutation relationships: [âS , a^S+ ]=1, [âAS , a^AS+ ]=1, [âS , a^AS+ ]=0 and, since the noise is uncorrelated at each spatial point, [(x), +(x′)]=δ(x-x′). The quantity Δβ=2βP -βS -βAS is the wavevector mismatch between the waves (βj is the propagation constant of the j-th wave, with j=P, S or AS), gS is equal to one half of the Raman (power) gain coefficient, φP is the phase of the pump wave and r=ωASns(χ*ASχS)12ωSnAS , where ωS ,ωAS are Stokes and anti-Stokes frequencies, nS , nAS the indices of refraction and χS ,χ*AS are the Raman susceptibilities. Since we will be interested in a narrow frequency band around the center of the Raman resonance, gS and r are real quantities.

The above equations (1a) and (1b) have the property that they conserve the commutators of the Stokes and anti-Stokes fields along propagation in the x-direction, so that for example, [âS (x), a^S+(x)]=[âS (x′), a^S+(x′)]. This is achieved in the equations through the presence of the terms involving the “Langevin noise source operators” (x), +(x), which introduce fluctuations in the Stokes and anti-Stokes fields, as described in more detail below.

The general solution to the above equations is:

âS(x)=A(x)âS(0)+B(x)âAS+(0)+0xdxNS(xx)N̂+(x)
(2a)
âAS+(x)=C(x)âAS+(0)+D(x)âS(0)+0xdxNAS(xx)N̂+(x).
(2b)

The details to the above solution along with general expressions for the A, B, C and D coefficients are presented in full in Appendix B.

3. Wavelength conversion noise figure

We can now calculate the fluctuations of the Stokes (anti-Stokes) field at the output x=L of the wavelength converter, when there is anti-Stokes (Stokes) field present at the input x=0. Let us consider, for example, the case where there are no photons at the Stokes frequency and the anti-Stokes wave at the input is a coherent state. The assumption of zero photons in the Stokes input is a good one, since at room temperature the mean photon number of the thermal optical radiation is practically zero. In our notation, |〉S,|〉AS denotes the states of the Stokes and anti-Stokes inputs, respectively.

We will calculate the noise figure first for the case of anti-Stokes to Stokes conversion when the anti-Stokes input is a coherent state of amplitude a. This situation corresponds to the state |0〉S |aAS |0〉R and the mean input photon number in the anti-Stokes field is |a|2.

Referring to Eqs. (2a) and (2b), we can express the mean photon rate at the Stokes output as:

âS+(x)âS(x)=B(x)2a2Convertedsignal+B(x)2+NS,TOTd.c.term,
(3)

where NS,TOT =0x dx′|NS (x)|2. For the fluctuations in photon number we find:

Δ(âS+(x)âS(x))2=B(x)2(NS,TOT+A(x)2+B(x)2)a2
+A(x)2(B(x)2+NS,TOT).
(4)

So far we have discussed single frequency signals; if the photons are sent in a pulse of spectral width Δf, then |a|2=Rf, where R is the input photon rate. The optical (signal-tonoise ratio) SNR for an ideal photodetector is:

SNR(x)=(aS+(x)âS(x)SIGNAL)2Δ(âS+(x)âS(x))2
=B(x)2R2B(x)2(NS,TOT+A(x)2+B(x)2)R·Δf+A(x)2(B(x)2+NS,TOT)(Δf)2.
(5)

The noise figure F of the wavelength conversion process is:

FASS=SNR(0)SNR(L)=RΔfSNR(L)
=1+NS,TOT+A(L)2B(L)2+A(L)2(B(L)2+NS,TOT)B(L)2(RΔf),
(6)

where the anti-Stokes input wave shot noise has been used as the reference. In the limit of a large input signal the minimum noise figure reduces to:

FASS,MIN=1+NS,TOT+A(L)2B(L)2.
(7)

On the other hand, the noise figure of Stokes to anti-Stokes conversion is given by:

FSAS=1+C(L)2+NAS,TOTD(L)2(1+1RΔf),
(8)

where NAS,TOT =0x dx′|NAS (x)|2. In the limit of a large input signal, this expression reduces to:

FSAS,MIN=1+C(L)2+NAS,TOTD(L)2.
(9)

Equations (7) and (9) provide a physical picture of the origin of the noise in the CARS conversion process. The excess noise figure has two parts. In the anti-Stokes to Stokes conversion, for example, the term proportional to |A(L)|2 originates from the zero-point fluctuations of the Stokes input, which propagates from x=0 to x=L while mixing with the converted signal. On the other hand, the term proportional to NAS,TOT quantifies the noise from the material vibrations that are coupled to the Stokes wave and the converted signal.

FASS,MIN=1+(1r2)2,
(10)

which is close to the minimum noise figure of an optical amplifier.

The second case is when gS (r 2-1)x≫1 and the conversion efficiency saturates. In this case, A(x)2r4(r21)2 , B(x)2r2(r21)2 and NS,TOT1(r21) , so the minimum noise figure reduces to:

FASS,MIN=2+r2(1r2)2.
(11)

These expressions deviate from the ideal amplifier by a factor which depends on the ratio of the photon energies of the Stokes and anti-Stokes waves, as well as the difference in their indices of refraction and their nonlinear susceptibilities. In the case of Δβ=0 for Stokes to anti-Stokes conversion, we have the following limiting values for the noise figure:

FSAS,MIN={2+1(rgSx)2,gS(r21)x12+(r21)2r2,gS(r21)x1.
(12)

In Fig. 1, we display the noise figure vs. pump intensity for the Stokes to anti-Stokes and anti-Stokes to Stokes conversion processes for perfect phase matching with an interaction length L of 2 cm, 2gS=(30 cm/GW)×IP and r=1.2 (typical parameters for semiconductor Raman media). We note that the curves for the two processes are qualitatively similar, converging in the weak pump limit, and asymptotically approaching the limiting value near 3 dB in the strong pump regime. When the intensity is low the noise figure is heavily impaired due to the low conversion gain, but at high intensity the quantum noise approaches a limit similar to the 3 dB minimum in amplifiers.

Fig. 1. The noise figure for signal conversion with perfect phase-matching is shown for interaction length L=2 cm, a gain coefficient 2gS=(30 cm/GW)×IP and r=1.2.

We have also plotted the noise figure for both Stokes to anti-Stokes and anti-Stokes to Stokes conversion as a function of the phase-mismatch, with r=1.2, gS = 0.25 cm -1 (which translates to a gain coefficient of 0.5 cm-1) and L=2 cm (see Fig. 2). Interestingly, the noise figure rapidly becomes large as the phase mismatch is increased, yet also shows periodic undulations corresponding with the well-known relationship between CARS efficiency and phase mismatch.

Fig. 2. The noise figure variation with the phase-mismatch for both Stokes and anti-Stokes scattering with L=2 cm, r=1.2, a gain coefficient 30 cm/GW and a pump intensity of 17 MW/cm2 (gS =0.25 cm -1).

4. Photon number distribution

We have also calculated the photon-number distribution for anti-Stokes (Stokes) photons at the output when the Stokes (anti-Stokes) field at the input is a coherent state. To perform the calculation we assume the “lumped” input-output relationships:

âS(x)=AâS(0)+BâAS+(0)+NSN̂+
(13a)
âAS+(x)=DâS(0)+CâAS+(0)+NASN̂+,
(13b)

with the commutator [, +]=1. The quantity |NS |2 equals the quantity NS,TOT in Section 4 of the paper and |NAS |2 equals NAS,TOT . The quantities A,B,C,D are the same quantities that arppear in equations (2a) and (2b). Formulated in this way, the lumped description of (13a) and (13b) gives the same results as the description of (1a) and (1b).

To derive the photon statistics, we assume that the reservoir is in a thermal state with the following density matrix:

ρ̂=(1q)n=0qnnRRn,
(14)

where q=〈nR/(1+〈nR) and 〈nR is the thermal occupation number of the reservoir. The probability distribution for both Stokes to anti-Stokes conversion and anti-Stokes to Stokes conversion has the following form:

p(no)=eμa2(1λ)λnoLno(μa2(λ1)λ),
(15)

where Ln0(x)=n=0n0n0n(x)nn! is the Laguerre polynomial of order n 0 and |a|2 is the mean photon number of the input coherent state. For Stokes to anti-Stokes conversion µ/λ=|D|2/(u+|D|2) and λ=(u+|D|2)/(1+u+|D|2) with u=|NAS |2 q/(1-q), whereas, for anti-Stokes to Stokes conversion µ/λ=|B|2/(u′+|B|2) and λ=(u′+|B|2)/(1+u′+|B|2) with u′=|NS |2/(1-q). The distribution has the mean value: 〈n 0〉=(λ+µ|a|2)/(1-λ). A derivation of this result for the simplified case where the reservoir is in the ground state is presented in Appendix C. Notice that when there is no input signal (a=0) the photon distribution in Eq. (15) becomes a thermal distribution; the spontaneously emitted Stokes and anti-Stokes fields in CARS have thermal distributions with a temperature that is determined by the parameter λ.

In Fig. 3, we show a representative plot of the probability distribution for anti-Stokes to Stokes conversion in the case of a coherent state input with a mean photon number of 400. This input yields a mean output photon number 45; for comparison, we have also plotted a Gaussian distribution with this mean value, illustrating that the tails of the CARS distribution deviate significantly from the perfect Gaussian. In particular, we see that the Gaussian approximation is adequate near the peak of the distribution, but becomes completely invalid for improbable events. This deviation may be very significant for applications in which the tails of the distribution are critical, such as data transmission and spectroscopic detection of dilute species. For example, knowledge of the tails may be important for determining the minimum bit error rate (BER) in a CARS wavelength converter or the concentration threshold for detection and identification in CARS spectroscopy. In Fig. 3, we also show the probability distribution for the Stokes output in the case of no input signal (spontaneous emission). In all the plots, the reservoir is taken to be in the ground state.

Fig. 3. Photon number probability distribution for the Stokes output with an anti-Stokes input of 400 photons (mean number) and a gain parameter gS ×L=0.33, assuming perfect phasematching. The mean photon number at the Stokes output is 45; the best fit Gaussian for this mean value is plotted for comparison. The Stokes at the output with no anti-Stokes input is also shown (only spontaneous emission present). The noise reservoir is taken to be in the ground state.

5. Conclusions

We have determined the noise in the CARS process by employing Langevin sources to keep the commutators of the field operators of the Stokes and anti-Stokes waves constant. We have calculated the noise figure and photon statistics for the wavelength conversion processes of Stokes to anti-Stokes and vice versa. The noise in the output Stokes and anti-Stokes fields has two contributions. One part is due to the input zero-point fluctuations coupling to the output and the other part is due to the coupling of the optical fields to the damped material vibrations. The best noise figure achievable is close to the 3 dB noise figure of the optical amplifier with a correction that depends on the ratio of Stokes to anti-Stokes frequencies, the ratio of the indices of refraction at the two frequencies, as well as the ratio of their nonlinear susceptibilities. The photon probability distribution has non-Gaussian tails, which may be significant for applications in data transmission and spectroscopic detection of highly rarified species.

Appendix A

We start with the classical propagation equations for the Stokes and anti-Stokes fields, which can be found in standard textbooks of nonlinear optics (see for example [17

17. R.W. Boyd, Nonlinear Optics, 3rd ed. (Academic Press, 1992).

]):

dASdx=gS(ωSnS)AP2AS+κ(ωSnS)AP2ejΔβxAAS*
(A1)
dAAS*dx=αAS(ωASnAS)AP2AAS*κ(ωASnAS)(AP*)2ejΔβxAS,
(A2)

To obtain the corresponding quantum-field equations, we replace the classical amplitudes with quantum mechanical operators by the substitution Aj →(ωj /nj )1/2âj . This replacement is similar in purpose to substitutions found in standard textbooks (see for example [18

18. A. Yariv, Quantum Electronics, 3rd ed. (Wiley: New York1989), Chapter 5.

]), but in this case the modes are enumerated in frequency rather than wavevector. This substitution yields the following differential equations for the field operators:

dâSdx=gSâS+κe2jφPejΔβxâAS+
(A3)
dâAS+dx=αASâAS+κe2jφPejΔβxâS,
(A4)

where φP is the phase of the pump wave, gS=gSωSnSAP2 , κ=κ(ωASωSnASnS)12AP2 and αAS=αASωASnASAP2 . These operators must obey the usual commutation relationships [âS , a^S+ ]=1, [âAS , a^AS+]=1 and [âS , a^AS+]=0.

It is not difficult to see however, that the propagation equations are not consistent with the uncertainty principle. For example, in the absence of coupling between the Stokes and anti-Stokes waves, it is well known that these modes experience only gain and loss, respectively. In this situation, however, the propagation equations displayed above predict that the commutators of the field operators are not conserved: d[âj , a^j+]/dx≠0. But this is inconsistent with basic quantum theory. The commutator of the field operators determines the noise limit imposed by the uncertainty principle. More specifically, 〈Δa^j,I2〉〈Δa^j,Q2〉≥(1/4)〈[âj,I , âj,Q ]〉2, where âj,I =(1/2)(âj +a^j+) and âj,Q =(1/2i)(âj -a^j+) represent the in-phase and quadrature components of the electromagnetic field. We must modify the propagation equations to include a so-called Langevin noise source operator, which restores the commutator to a constant value, i.e., d[âj,I , âj,Q ]/dx=0. The analysis presented here is inspired by the analysis of an optical amplifier in [19

19. H.A. Haus, Electromagnetic Noise and Quantum Optical Measurement (Springer2000), p. 217–223.

]. Since all the terms in the propagation equations are due to coupling of the optical fields with the damped vibrational modes in the material, we need to introduce only one ‘noise reservoir’ into the propagation equations. As we will show we can obtain consistent results with this approach.

Consider first the equation describing the evolution of the Stokes mode (decoupled from the anti-Stokes) and add a noise source operator G :

dâSdx=gSâS+F̂G.
(A5)

The commutator conservation equation gives the condition: [G ,a^S+]+[âS , F^G+]=-2gS . The noise operator and the field operator do not commute. The inhomogeneous solution of the propagation equation gives:

âSINH=xdxexp(gS(xx))F̂G(x).
(A6)

Since the noise is uncorrelated along the propagation, the commutator of the noise source with its conjugate should be proportional to a delta function; thus,

[âSINH(x),F̂G+(x)]=(12)C,
(A7)

where C is a proportionality constant. In order to satisfy Eq. (A5), we must have C=-2gS and the noise sources must obey the following commutation relation:

[F̂G(x),F̂G+(x)]=2gSδ(xx).
(A8)

Notice that the commutator has a negative value, implying that G (x) acts as a creation operator and F^G+(x) as an annihilation operator for “noise reservoir” excitations. To make the correspondence more explicit we rename the noise source operator: F̂G(x)=exp(jθS)2gSN̂+(x) , where (x) and +(x) are now standard annihilation and creation operators with [(x), +(x′)]=δ(x-x′) and exp(S ) is an undetermined phase factor whose purpose will become clear below. In a similar manner, we can also show that for the decoupled anti-Stokes mode, the noise source F̂L(x)=exp(jθAS)2αASN̂(x) is required.

κej(Δβx+2φP)[âAS+,âAS]κej(Δβx+2φP)[âS,âS+]
+exp(jθS)2gS[N̂+,âAS]+exp(jθAS)2αAS[âS,N̂]=0.
(A9)

We can evaluate the last two commutators in this expression using the inhomogeneous solution of the propagation equations:

âS(x)=exp(jθS)2gS0xdxPS(xx)N̂+(x)
(A10)
âAS+(x)=exp(jθAS)2αAS0xdxPAS(xx)N̂+(x),
(A11)

where in the above expressions PS (0)=PAS (0)=1. The functions PS (x) and PAS (x) are calculated explicitly in Appendix B (see equations (B9) and (B10)). Inserting (A10) and (A11) in (A9) we obtain:

θSθAS=Δβx+2φP+π,
(A12)

where we have incorporated the relations, αAS =r 2 gS and κ=rgS in the last step, with r=ωASnS(χ*ASχS)12ωSnAS . In what follows, we choose θS =-π/2 and θAS =-Δβ·x-2ψP -(3π/2) without loss of generality. Our final equations are:

dâSdx=gSâS+rgSexp(j(Δβx+2φP))âAS+j2gSN̂+
(A13)
dâAS+dx=r2gSâAS+rgSexp(j(Δβ·x+2φP))âS
+jexp(j(Δβ·x+2φP))2r2gSN̂+.
(A14)

Appendix B

Here we derive the solution to the coupled differential equations (1a) and (1b) for the Stokes and anti-Stokes field operators. In Eqs. (1a) and (1b) we make the substitutions:

âS=b̂Sexp(jΔβx2),âAS+=b̂AS+exp(jΔβx2).
(B1)

This yields:

ddxb̂s=(gSj(Δβ2))b̂s+rgsθb̂AS+j2gsejΔβx2N̂+(x)
(B2)
ddxb̂AS+=(j(Δβ2)r2gs)b̂AS+rg2θ*b̂s
(B3)

where θ=exp(j2ψP ). It is easy to see that the linear combinations:

Û(+)=b̂S+C(+)b̂AS+,Û()=b̂S+C()b̂AS+,
(B4)

with

C(±)=((r2+1)gsiΔβ)±((r2+1)gsiΔβ)24r2gs22rgsθ,
(B5)

obey the following equations:

ddxÛ(±)=Λ(±)Û(±)+j2gs(rθ*C(±)1)ejΔβ·x2N̂+(x),
(B6)

where C(±)=(gS -iβ/2)-Λ(±) θ/(rgS ) or

Λ(±)=gs(1r2)(gs(1+r2)iΔβ)24r2gs22.
(B7)

The solution of Eq. (B6) is:

Û(±)(x)=exp(Λ(±)x)Û(±)(0)
+j2gs(rθ*C(±)1)0xdxexp(Λ(±)(xx))ejΔβ·x2N̂+(x).
(B8)

Using Eq. (B4) we find the final expressions for the Stokes and anti-Stokes operators:

b̂s(x)=C()exp(Λ(+)x)C(+)exp(Λ()x)C()C(+)b̂s(0)
+C()C(+)exp(Λ(+)x)exp(Λ()x)C()C(+)b̂AS+(0)
+j2gs0xdx[(rθ*C(+)1)C()exp(Λ(+)(xx))C()C(+)
(rθ*C()1)C(+)exp(Λ()(xx))C()C(+)]ejΔβ·x2N̂+(x),
(B9)
b̂AS+(x)=exp(Λ(+)x)exp(Λ()x)C(+)C()b̂s(0)
+C(+)exp(Λ(+)x)C()exp(Λ()x)C(+)C()b̂AS+(0)
+j2gs0xdx[(rθ*C(+)1)exp(Λ(+)(xx))C(+)C()
(rθ*C()1)exp(Λ()(xx))C(+)C()]ejΔβ·x2N̂+(x).
(B10)

Substituting Eqs. (B9) and (B10) into Eqs. (B1) we find the solution for the field operators âS and a^AS+. The parameters A, B, C and D that appear throughout the main text are then given by the respective coefficients of the equations.

Appendix C

Here we derive the expression given above for the photon number probability distribution in CARS wavelength conversion. To calculate the probability distribution, we must determine its characteristic function defined as the Fourier series of the probability distribution:

P˜(k)=n=0p(n)ejkn=ejkn,
(C1)

which is also the expectation value of e -jkn where n is the photon number. The probability distribution is expressed in terms of the characteristic function as:

p(no)=12π02πdkejknoejkn.
(C2)

We present the calculation of the characteristic function for anti-Stokes to Stokes conversion; the derivation for Stokes to anti-Stokes conversion proceeds in a similar fashion. Also we take the reservoir to be in the ground state. The case where the reservoir is in a thermal state is slightly more complicated to compute but similar to the calculation below.

The characteristic function to be evaluated is:

exp(jkâS+(x)âS(x))=R0ASaS0exp(jkâS+(x)âS(x))0SaAS0R.
(C3)

We will need the expansion of the exponential with the operators put in normal and in anti-normal order:

exp(ξâ+â)=n=0(eξ1)nn!(â+)nân,
(C4)
exp(ξâ+â)=eξn=0(1eξ)nn!ân(â+)n.
(C5)

These results can be found for example in [20

20. W.H. Louisell, Quantum Statistical Properties of Radiation (John Wiley & Sons1973).

]. Next, we compute 〈(âS (x))n(a^S+ (x))n〉.The expectation value is given by:

R0ASaS0(âS(x))n(âS+(x))n0SaAS0R=l=0n(nl)n!l!Ba2lA2(nl), ,

l=0nnln(n1)...(l+1)(ξ)lηn=ηnl=0nnlnlξn1(ξ)n=ηn(1+ξ)n(ξ)n. .

The characteristic function can now be expressed as:

exp(ξâS+(x)âS(x))=eξn=0(1eξ)nn!âS(x)n(âS+(x))n
=eξn=0(1eξ)nn!ηn(1+ξ)n(ξ)n.

. Next, considering 1+/∂ξ′ and ξ′ as operators we note that [1+/∂ξ′,ξ′]=1 so that they have the same commutation relationship as â and â +. We define a new variable ξ″ as η-ηeξ =1-eξ ″, which we use to rewrite the characteristic function:

exp(ξâS+(x)âS(x))=eξ+η1ηn=0(1eξ)nn!(1+ξ)n(ξ)n
=1+(η1)eξηexp(ξξ(1+ξ)).

exp(ξâS+(x)âS(x))=1+(η1)eξηn=0(eξ1)nn!(ξ)n(1+ξ)n
=1+(η1)eξηn=0(eξ1)nn!(ξ)n
=1+(η1)eξηexp(ξ(1eξ)).

Finally, substituting ξ″ and setting ξ=jk we obtain the characteristic function in closed form:

exp(jkâS+(x)âS(x))=eξejk1+η(ejk1)exp(ξ1+η(ejk1)).
(C6)

As previously stated, the probability that the photon number a^S+(x)âS (x) has a value no given a coherent state with amplitude a as the input is given by the Fourier transform of the characteristic function. To compute the transform, we rewrite the integral obtained from Eqs. (C2) and (C6) by changing to the variable z=ejk :

p(noa)=eξdz2πizno(zzo)exp(ξηzzo),
(C7)

where Zo =(η-1)/η and the integrand has a pole at z=zo . To evaluate this integral, we expand the exponential and carry out the integration in the complex plane of every term in the expansion. The result is:

p(noa)=eξn=0(ξη)nηn!znodz2πi(zzo)n+1=eξn=0no(non)(ξη)nηn!zonon.
(C8)

This is the distribution we are seeking.

Acknowledgements

The authors would to like to thank Dr. J.D. Shah of DARPA/MTO for his support.

References and links

1.

P.D. Maker and R.W. Terhune, “Study of optical effects due to an induced polarization third order in the electric field strength,” Phys. Rev. 137, A801 (1965). [CrossRef]

2.

H. Vogt, “Coherent and hyper-Raman techniques” in Topics in Applied Physics (eds. M. Cardona and G. Guntherodt) vol. 50, p. 207, Springer-Verlag (1982).

3.

A.M. Zheltikov and P. Radi, “Non-linear Raman spectroscopy 75 years after the Nobel Prize for the discovery of Raman scattering and 40 years after the first CARS experiments,” J. of Raman Spectrosc. 36, 92–94 (2005). [CrossRef]

4.

G. Beadie, Z.E. Sariyanni, Y. V. Rostovtsev, T. Opatrny, J. Reintjes, and M.O. Scully, “Towards a FAST CARS anthrax detector: coherence preparation using simultaneous femtosecond laser pulses,” Opt. Commun. 244, 423–430 (2005). [CrossRef]

5.

K.P. Knutsen, J.C. Johnson, A.E. Miller, P.B. Petersen, and R.J. Saykally, “High-spectral-resolution multiplex CARS spectroscopy using chirped pulses,” SPIE Proc. 5323, 230–239 (2004). [CrossRef]

6.

K. Ishii and H. Hamaguchi, “Picosecond time-resolved multiplex CARS spectroscopy using optical Kerr gating,” Chem. Phys. Lett. 367, 672–677 (2003). [CrossRef]

7.

E.O. Potma, C.L. Evans, and X.S. Xie, “Heterodyne coherent anti-Stokes Raman scattering (CARS) imaging,” Opt. Lett. 31, 241–243 (2006). [CrossRef] [PubMed]

8.

F. Vestin, M. Afzelius, C. Brackmann, and P-E. Bengtsson, “Dual-broadband rotational CARS thermometry in the product gas of hydrocarbon flames,” Proc. of the Combustion Inst. 30, 1673–1680 (2004). [CrossRef]

9.

N. Djaker, P-F. Lenne, and H. Rigneault, “Vibrational imaging by coherent anti-Stokes Raman scattering (CARS) microscopy,” SPIE Proc. 5463, 133–139 (2004). [CrossRef]

10.

R. Claps, V. Raghunathan, D. Dimitropoulos, and B. Jalali, “Anti-Stokes Raman conversion in Silicon waveguides,” Opt. Express 11, 2862–2872 (2003). [CrossRef] [PubMed]

11.

J. Perina, “Photon statistics in Raman scattering with frequency mismatch,” Optica Acta 28, 1529 (1981). [CrossRef]

12.

J. Perina, “Photon statistics in Raman scattering of intense coherent light,” Optica Acta 28, 325 (1981). [CrossRef]

13.

M. Karska and J. Perina, “Photon statistics in stimulated Raman scattering of squeeze light,” J. Mod. Optics 37, 195 (1990). [CrossRef]

14.

P.L. Voss and P. Kumar, “Raman-noise-induced noise-figure limit for χ(3) parametric amplifiers,” Opt. Lett. 29, 445 (2004). [CrossRef] [PubMed]

15.

R. Tang, P.L. Voss, J. Lasri, P. Devgan, and P. Kumar, “Noise-figure of fiber-optical parametric amplifiers and wavelength converters: experimental investigation,” Opt. Lett. 29, 2372 (2004). [CrossRef] [PubMed]

16.

K.K. Das, G. S. Agarwal, Y.M. Golubev, and M.O. Scully, “Langevin analysis of fundamental noise limits in coherent anti-Stokes Raman spectroscopy,” Phys. Rev. A 71, 013802 (2005). [CrossRef]

17.

R.W. Boyd, Nonlinear Optics, 3rd ed. (Academic Press, 1992).

18.

A. Yariv, Quantum Electronics, 3rd ed. (Wiley: New York1989), Chapter 5.

19.

H.A. Haus, Electromagnetic Noise and Quantum Optical Measurement (Springer2000), p. 217–223.

20.

W.H. Louisell, Quantum Statistical Properties of Radiation (John Wiley & Sons1973).

OCIS Codes
(190.2620) Nonlinear optics : Harmonic generation and mixing
(190.4410) Nonlinear optics : Nonlinear optics, parametric processes
(190.5650) Nonlinear optics : Raman effect
(270.5290) Quantum optics : Photon statistics

ToC Category:
Quantum Optics

History
Original Manuscript: August 7, 2006
Revised Manuscript: October 13, 2006
Manuscript Accepted: October 13, 2006
Published: November 13, 2006

Virtual Issues
Vol. 1, Iss. 12 Virtual Journal for Biomedical Optics

Citation
D. Dimitropoulos, D. R. Solli, R. Claps, and B. Jalali, "Noise figure and photon statistics in coherent anti-Stokes Raman scattering," Opt. Express 14, 11418-11432 (2006)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-14-23-11418


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. P. D. Maker and R. W. Terhune, "Study of optical effects due to an induced polarization third order in the electric field strength," Phys. Rev. 137, A801 (1965). [CrossRef]
  2. H. Vogt, "Coherent and hyper-Raman techniques" in Topics in Applied Physics, M. Cardona and G. Guntherodt, eds., (Springer-Verlag, 1982) Vol. 50, p. 207.
  3. A. M. Zheltikov and P. Radi, "Non-linear Raman spectroscopy 75 years after the Nobel Prize for the discovery of Raman scattering and 40 years after the first CARS experiments," J. Raman Spectrosc. 36, 92-4 (2005). [CrossRef]
  4. G. Beadie, Z. E. Sariyanni, Y. V. Rostovtsev, T. Opatrny, J. Reintjes, and M. O. Scully, "Towards a FAST CARS anthrax detector: coherence preparation using simultaneous femtosecond laser pulses," Opt. Commun. 244, 423-430 (2005). [CrossRef]
  5. K. P. Knutsen, J. C. Johnson, A. E. Miller, P. B. Petersen, and R. J. Saykally, "High-spectral-resolution multiplex CARS spectroscopy using chirped pulses," in Multiphoton Microscopy in the Biomedical Sciences IV, A. Periasamy and P. T. So, eds, Proc. SPIE 5323, 230-239 (2004). [CrossRef]
  6. .K. Ishii and H. Hamaguchi, "Picosecond time-resolved multiplex CARS spectroscopy using optical Kerr gating," Chem. Phys. Lett. 367, 672-7 (2003). [CrossRef]
  7. .E.O. Potma, C.L. Evans, X.S. Xie, "Heterodyne coherent anti-Stokes Raman scattering (CARS) imaging," Opt. Lett. 31, 241-3 (2006). [CrossRef] [PubMed]
  8. F. Vestin, M. Afzelius, C. Brackmann, P-E. Bengtsson, "Dual-broadband rotational CARS thermometry in the product gas of hydrocarbon flames," Proc. of the Combustion Inst. 30, 1673-80 (2004). [CrossRef]
  9. N. Djaker, P-F. Lenne, H. Rigneault, "Vibrational imaging by coherent anti-Stokes Raman scattering (CARS) microscopy," in Femtosecond Laser Applications in Biology, S. Avrillier, and J.-M. Tualle, eds., SPIE Proc. 5463, 133-139 (2004). [CrossRef]
  10. R. Claps, V. Raghunathan, D. Dimitropoulos, and B. Jalali, "Anti-Stokes Raman conversion in Silicon waveguides," Opt. Express 11, 2862 - 2872 (2003). [CrossRef] [PubMed]
  11. J. Perina, "Photon statistics in Raman scattering with frequency mismatch," Optica Acta 28, 1529 (1981). [CrossRef]
  12. J. Perina, "Photon statistics in Raman scattering of intense coherent light," Optica Acta 28, 325 (1981). [CrossRef]
  13. M. Karska and J. Perina, "Photon statistics in stimulated Raman scattering of squeeze light," J. Mod. Optics 37, 195 (1990). [CrossRef]
  14. P. L. Voss and P. Kumar, "Raman-noise-induced noise-figure limit for χ(3) parametric amplifiers," Opt. Lett. 29, 445 (2004). [CrossRef] [PubMed]
  15. R. Tang, P. L. Voss, J. Lasri, P. Devgan, and P. Kumar, "Noise-figure of fiber-optical parametric amplifiers and wavelength converters: experimental investigation," Opt. Lett. 29, 2372 (2004). [CrossRef] [PubMed]
  16. K. K. Das, G. S. Agarwal, Y. M. Golubev, and M. O. Scully, "Langevin analysis of fundamental noise limits in coherent anti-Stokes Raman spectroscopy," Phys. Rev. A 71, 013802 (2005). [CrossRef]
  17. R. W. Boyd, Nonlinear Optics, 3rd ed. (Academic Press, 1992).
  18. A. Yariv, Quantum Electronics, 3rd ed. (Wiley: New York 1989), Chap. 5.
  19. H. A. Haus, Electromagnetic Noise and Quantum Optical Measurement (Springer 2000), p. 217 - 223.
  20. W. H. Louisell, Quantum Statistical Properties of Radiation (John Wiley & Sons 1973).

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Figures

Fig. 1. Fig. 2. Fig. 3.
 

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited