Abstract
We develop a method for finding the number and shapes of the independently squeezed or amplified modes of a spatially-broadband, travelling-wave, frequency- and polarization-degenerate optical parametric amplifier in the general case of an elliptical Gaussian pump. The obtained results show that for tightly focused pump only one mode is squeezed, and this mode has a Gaussian TEM00 shape. For larger pump spot sizes that support multiple modes, the shapes of the most-amplified modes are close to Hermite- or Laguerre-Gaussian profiles. These results can be used to generate matched local oscillators for detecting high amounts of squeezing and to design parametric image amplifiers that introduce minimal distortion.
©2011 Optical Society of America
1. Introduction
Phase-sensitive optical parametric amplifiers (PSAs), unlike any other optical amplifiers, can increase the magnitude of the signal without adding any noise [1]. This property, coupled with the wide temporal bandwidth of fiber-based parametric amplifiers, has lead to their use as nearly noiseless inline amplifiers for optical communication systems [2–6]. Similarly, the broad spatial bandwidth of the parametric amplifiers enables their use as noiseless amplifiers of faint optical images, as was theoretically proposed in [7–9] and experimentally demonstrated in [10–13]. The signal-to-noise ratio improvement provided by such image amplifiers can in turn lead to resolution improvement in the detection of faint images [14–16]. The same devices can also be used for the generation of spatially-broadband squeezed vacuum for quantum information processing applications. The design of practical parametric image amplifiers, however, represents a big challenge. Indeed, the traveling-wave nature of gain in these devices requires the use of a tightly focused pump beam (typically, the fundamental Gaussian TEM00 mode) of very high power (~1 kW per pixel of resolution) [17]. The resulting spatially-varying PSA gain, together with the limited spatial bandwidth of the PSA [determined from phase-matching conditions to be ~(kp/L)1/2, where kp is the pump propagation constant and L is the nonlinear crystal’s length] [7, 12, 18], couples and mixes up in both space and spatial-frequency domains the modes representing the information content of the image [17]. These spatial mode-mixing effects, known as gain-induced diffraction [19, 20], also make it difficult to detect the lowest-noise mode of a traveling-wave squeezer [21, 22], because the properly mode-matched homodyne detector requires exact knowledge of this mode’s spatial profile. The coupling between the various spatial modes of the signal was previously studied either in a travelling-wave PSA in which only the radially symmetric signal modes were considered [23, 24], or in a cavity-based sub-threshold optical-parametric oscillator [25], all pumped by a circular Gaussian beam. In this paper, we extend the mode-coupling analysis to the travelling-wave PSA with a more general elliptical Gaussian TEM00 pump and then develop a procedure to find, for the first time to our knowledge, the orthogonal set of independently squeezed (or amplified) eigenmodes of such a PSA. We presented the preliminary results of this approach in a recent conference paper [26].
The rest of this paper is organized as follows: Section 2 describes the theoretical model of mode coupling in the PSA and the procedure for finding the PSA’s eigenmodes, Section 3 describes the details of our computational approach, Section 4 reports and discusses the results, and Section 5 summarizes our work.
2. Theory of 2-D spatially-multimode PSA
2.1 Hermite-Gaussian mode representation
Our work builds upon the Laguerre-Gaussian (LG) expansion of the parametric amplifier equations, originally developed for cavities [25] and later extended to 1-D (radially symmetrical) traveling-wave PSAs [23]. We generalize this approach by replacing the LG expansion with a Hermite-Gaussian (HG) expansion to describe arbitrary 2-D images and a general elliptical Gaussian pump beam with potentially unequal 1/e intensity radii a0px and a0py in x- and y-dimensions, respectively. By expanding the PSA input over signal TEMmn HG modes having the same Rayleigh ranges zRx = kpa0px2 and zRy = kpa0py2 as the pump, we reduce the PSA propagation to a system of coupled ordinary differential equations for the HG-mode amplitudes Amn = Xmn + iYmn. After integrating this system of equations, we can obtain the PSA’s Green’s function and all the quantum correlators needed for finding the independently squeezed or amplified modes [27].
We start by considering the nonlinear paraxial wave equation of a degenerate optical parametric amplifier in the undepleted pump approximation [17, 27, 28]
We are looking for the solutions in the form , where is a slowly-varying field envelope, is a transverse vector with coordinates (x,y), the intensity is given by with index i taking value of either s or p, denoting the signal or pump field, respectively, and ωp = 2ωs. We expand the signal and pump beams over the HG basis with potentially unequal 1/e intensity radii in x- and y-dimensions:
where the one-dimensional HG modes gm are defined as (β = x or y)with the orthogonality conditionand the Rayleigh range zR, 1/e intensity radius a(z), Gouy phase shift θ(z), and beam’s radius of curvature R(z), respectively, given byThe x- and y-radii for the signal expansion basis are chosen to be 21/2 times greater than those of the pump, which ensures that the signal and the pump have the same wavefront curvature and the same Rayleigh ranges zRx = kpa0px2 and zRy = kpa0py2. The beam waists for both the pump and the signal HG basis occur at z = 0.
After substituting the HG expansions of the signal and the pump [Eq. (2)] into Eq. (1) and projecting the results onto the signal’s HG basis, we arrive at the following coupled-mode equations for the signal’s HG mode amplitudes Amn(z):
whereΔk is the wavevector mismatch, P0 is the pump power and θp is the initial pump phase. The overlap integral Bmm′ of the pump and two signal modes with indices m and m′ has a closed-form expressionwhere the double factorial (m + m′ – 1)!! = 1 for m + m′ = 0. Let us list several important properties of the coupling matrix Bmm′ comprising coefficients Bmm′. First of all,i.e., propagation away from the center of the crystal reduces the magnitude of Bmm′ and introduces a phase shift due to the Gouy phase mismatch between the pump and the two signal modes. Second, the magnitude of Bmm′(0) for large indices m and m′ can be approximated asi.e., for a fixed (m + m′), it exhibits a fast Gaussian decay as a function of (m – m′), whereas the decay along the main diagonal m = m′ is very slow and proportional to (m + m′)–1/2, as illustrated in Fig. 1 . Hence, the fast decay of Bmm′ versus (m – m′) serves as the selection rule favoring coupling between the signal modes with close indices. On the other hand, the slow decay versus (m + m′) means that the maximum range of the amplified signal modes is determined not by the magnitude of the coupling coefficient, but by its Gouy phase mismatch [numerator in Eq. (9)] that leads to fast oscillations at large (m + m′), which limits the PSA gain. The signal evolution in the PSA in Eq. (6) is governed by the outer product of matrices Bmm′ and Bnn′, yielding a 4th-rank tensor in which only 25% of the elements are not zero.2.2 Laguerre-Gaussian mode representation
In the case of a circular Gaussian pump beam with a0px = a0py = a0p, it is more convenient to use the LG expansion [25]:
where the LG modes fpl of radial index p and azimuthal index l have the formr, φ are the polar coordinates,and the orthogonality condition isThe fundamental LG mode LG00 coincides with the fundamental circular HG mode TEM00:
The coupled-mode equations in the LG representation take the form of
wherethe overlap integral of the pump mode LG00, with the signal modes LGpl and LGp′l′ is given byand the coupling coefficients between the modes LGpl and LGp′l′ are zero for l ≠ –l′. The latter selection rule means that the circular symmetry of the pump reduces the rank of the coupling tensor in the PSA propagation equation from 4 to 3. Apart from this simplification, other properties of the coupling coefficients in Eq. (16) are similar to those in the HG case. Namely, these coefficients can be expressed asand for large value of (p + p′ + |l|) they have an asymptotic dependencewhich imposes fast Gaussian decay when indices deviate from p = p′ and l = 0 condition, as well as a slow p –1/2 decay when they satisfy this condition. We also note that, as a consequence of Eq. (15), the first diagonal element of the coupling tensor is the same in both LG and HG representations:2.3 PSA Green’s functions and eigenmodes
In what follows, we will use index notation corresponding to the HG representation; results for the LG representation can be obtained in a similar manner. We will assume that the pump beam waist is located exactly at the center of a nonlinear crystal of length L, i.e., the crystal input is at z = –L/2 and crystal output is at z = L/2. Following our approach to solving the PSA equation and finding its eigenmodes, presented in [27], we re-write the complex signal mode amplitudes Amn = Xmn + iYmn in a matrix-vector form
where and are real-valued matrices. Then, we can express the solution aswhere the tensor is the PSA Green’s function in HG representation, relating the mode amplitudes at the input of the crystal to those at the output. Unlike the numerically challenging evaluation of the PSA Green’s function in xy-representation [29], can be easily computed by exciting a finite number of input HG modes , one at a time, and recording the output mode pattern in each case. Once the Green’s function is known, the amount of squeezing can be found using the Rayleigh quotientFrom the symplectic property of the tensor [27, 30–32] it follows that the output can be decomposed into a set of independently squeezed modes (a Karhunen-Loève expansion) determined from the eigenvalue equation
where λ is the gain (same as squeezing factor) of the eigenmode . Symplecticity of the tensor ensures that for each eigenmode with eigenvalue λ, there is an eigenmode with eigenvalue 1/λ that is obtained by –90° phase shift of (i.e., by substituting for , and for ):Equation (24) allows one to find the independently squeezed / amplified modes at the PSA output (at z = L/2). In order to obtain their shapes at the PSA input (at z = –L/2), one can invoke the reciprocity argument [23], which stems from the equivalence between z-reversal and conjugation of the coupling coefficients in Eq. (6) (for this, it is helpful to assume θp = –π/2 as the initial phase reference). The reciprocity principle states that the squeezing factor detected by a particular spatial profile of the local oscillator is the same as the classical PSA gain seen by the input signal whose profile is conjugate to that of the local oscillator. When applied to the local-oscillator shape that matches one of the PSA eigenmodes (and for the PSA eigenmode the gain and squeezing factors are the same), this principle means that the eigenmode profile at the PSA input is the conjugate of the eigenmode profile at the PSA output.
3. Computational considerations
We numerically integrate Eq. (6) with 4th-rank coupling tensor Bmm′Bnn′ by the 4th-order Runge-Kutta method (RK-4). The two critical parameters of the computation are the maximum HG order mmax, at which the tensor is truncated, and the step size Δz over the propagation distance z.
To estimate the truncation order mmax for each dimension, we recall that it is determined by the Gouy phase mismatch [the numerator in Eq. (9)]. When the mismatch reaches π, the coupling coefficient changes sign, i.e., switches from amplification to de-amplification. If this reversal occurs within the crystal, it means that the mode is outside of the crystal’s spatial bandwidth and can be ignored. Thus, at m = m′ = mmax the mismatch should reach π at the end of the crystal:
which results in the truncation-order conditionThe same condition can be obtained by noting that the spatial-frequency bandwidth of the mmax-order HG function, which is approximately , should be equal to the spatial bandwidth of the crystal (πkp/L)1/2 [17]. For a PPKTP crystal with L = 2 cm, ns ≈np ≈1.78, free-space λp = 780 nm, and a0p = 100 μm, Eq. (27) yields mmax = 22.5, which we have confirmed by observing the sufficient decay of the Green’s function versus (m + m′) at this index value. To provide additional accuracy margin, we use mmax = 32 for these PSA parameters. For a circular pump spot size, the x- and y-dimensions require 32 modes each.
As the pump waist increases, Eq. (27) demands significant increase in the number of considered modes in both the x- and y-dimensions, proportional to the square of the corresponding waist radii:
Thus, for a0p = 200 μm we need to use mmax = 128, i.e., mmax that is 4 times greater than that for a0p = 100 μm. Thus, for a 200-μm circular pump waist, Eq. (24) requires diagonalization of a tensor consisting of 2×2×(128×128)×(128×128) = 230 elements (only quarter of which are non-zero), leading to multi-gigabyte memory requirements for the computer. This unfortunately fast scaling with pump waist size takes place because our choice of the HG expansion basis for the signal is not optimal. More specifically, in order to eliminate the beam curvature from the overlap integral, the signal waist of the expansion basis has been chosen to be 21/2 times greater than the pump’s, whereas all the interesting dynamics takes place on a smaller transverse scale that varies between the inverse spatial bandwidth of the crystal on the low end and the pump waist size on the high end. The number of independently amplified (or squeezed) PSA modes can be estimated by the product of the pump beam area and the square of the crystal’s spatial bandwidth [17, 33]:
which is square root of the number of modes required for computation by Eq. (28). Significant potential computational memory savings [reduction of mxmaxnymax down to the value given by Eq. (29)] can be achieved by using a different HG expansion basis for the signal, that with a beam waist in the vicinity of the geometric average of the pump waist and the inverse spatial bandwidth of the crystal [34], as was previously discussed in the context of temporal modes [35]. These savings, however, are achieved at the expense of a considerably more demanding computation of the coupling coefficients at every z-step, and have not been attempted in the present work.The estimation of the step size Δz can be done by demanding both the coefficients Bmm′ and the solution to not change much over one step. The fastest variation in the coupling coefficients comes from the Gouy phase mismatch for the maximum order mmax, and we require this phase to change over one step by no more than π/150 ≈0.02 rad:
which leads to the step size that is independent of the pump waist radius:The second requirement on variation of the solution over one step can be met by requiring the magnitude of the most amplified eigenmode to change by no more than 2% (0.086 dB) over one step, which means that the step size given by Eq. (31) meets this requirement up to the gain of ≈300 × 0.086 dB = 25.8 dB for the most amplified mode.
4. Results and discussion
Using the approach discussed in Sections 2.1 and 2.3, we have numerically generated the Green’s functions of Eq. (6) in the HG representation and obtained the PSA eigenmodes for various pump powers P0 and spot sizes a0px × a0py for a PPKTP crystal of length L = 2 cm, deff = 8.7 pm/V, ns ≈np ≈1.78, and signal wavelength of 1560 nm. For pump spot sizes of 100×100 μm2 and smaller (100×50, 100×25, and 25×25 μm2) we have used 32×32 HG functions (i.e., mxmax = nymax = 32) as the signal expansion basis. For pump waists larger than 100 μm, we have increased the truncation order according to Eq. (28); e.g., for 400×100 μm2 pump spot size we have used 512×32 expansion basis. For 100×100 μm2 pump spot size, the Rayleigh range equals zR = kp(a0p)2 ≈143 mm, i.e., it is ~7 times longer than the crystal.
In the first set of results, whose eigenvalue spectra are shown in Fig. 2 , we have chosen the pump powers such that the gains (which are also the same as the squeezing factors) of the most amplified (fundamental) eigenmode for each spot size are approximately the same and equal to 15 (or 11.8 dB). These powers are shown in Fig. 3(b) and also listed in the legends of Figs. 4 , 5 . From Fig. 2 one can see that the number of modes amplified by the PSA grows with the pump spot size. For circular pump waists, some eigenvalues are degenerate, whereas for elliptical waists the degeneracy is lifted. The smallest pump spot size (25×25 μm2) seems to support the amplification of only the fundamental mode.
In the second set of results, whose eigenvalue spectra are shown in Fig. 3(a), we have adjusted the pump powers so that the peak intensity of the pump for all spot sizes is the same as that in the 100×100 μm2 case of the first set of results, i.e., the pump power is scaled proportionally to the beam area. The spectra shown in Fig. 3(a) indicate that, under the constant intensity condition, the gain of the fundamental eigenmode decreases with decreasing spot size. This point is further illustrated in Fig. 3(b) (red), which shows the gain, at constant pump intensity, as a function of the spot size a0px × a0py. In the same Fig. 3(b), we also plot the pump powers (blue) needed to achieve the fundamental-mode gain of ~15, as in the first set of data. The fact that this dependence is not linear indicates that higher intensity is required to achieve the same gain for the fundamental eigenmode when the tightly focused pump size starts to approach the inverse spatial bandwidth of the crystal (L/πkp)1/2 ≈21 μm.
All the remaining Figs. 4–9 present results from the first set of data (i.e., where the pump powers are adjusted to produce the same gain of ~15 for the fundamental eigenmode). Figures 4 and 5 show the (x,y) representation of the first few computed eigenmodes (i.e., the spatial profiles of their intensity), whereas Figs. 6 –9 show the HG representation of these eigenmodes (i.e., |Amn|2). The inset in Fig. 4 also shows the number of modes with gains within 3 dB of the gain of the fundamental eigenmode [i.e., the number of modes above the dashed black line in Fig. 2(a)] as a function of the pump power, exhibiting roughly linear dependence with the inverse slope of ~690 W / mode, which can be considered as the power efficiency of the PSA under consideration. The large required pump power justifies the use of the undepleted pump approximation in Eq. (1) for input signal powers up to ~1 W (i.e., as long as the amplified signal is much weaker than the pump).
From Figs. 4–9, one can make two main observations: 1) the fundamental eigenmode always exhibits a TEM00-like single-lobe (x,y) shape, even though its HG representation consists of a superposition of many modes (except in the 25×25 μm2 case, for which almost no HG modes other than the TEM00 are present); 2) although at small pump spot sizes the higher-order eigenmodes have complicated (x,y) shapes, at spot sizes of 100×100 μm2 and larger the eigenmodes’ amplitude and phase profiles, at least qualitatively, resemble those of the HG modes TEMmn (for elliptical pump waists) or the LG modes (for circular pump waists). We discuss these two observations separately below.
The first observation is further illustrated in Fig. 6, which shows the PSA’s fundamental eigenmode #0 in the HG representation for various pump spot sizes. The overlap of this mode with the TEM00 mode of our HG basis is 97% for the 25×25 μm2, 61% for the 100×25 μm2, 50% for the 100×50 μm2, 35% for the 100×100 μm2, 16.2% for the 800×50 μm2, 16.1% for the 400×100 μm2, and 15.7% for the 200×200 μm2 pump spot sizes. Thus, even though the eigenmode #0 is a single-lobe in all 7 cases, it requires a significant number of HG modes for its representation in all cases except the 25×25 μm2 pump spot size. Higher-order eigenmodes involve superpositions of even greater numbers of HG modes, as shown in Figs. 7 –9.
The second observation leads us to the conclusion that, for PSAs with pump spot sizes large enough to support amplification of many eigenmodes, it should be possible to represent each eigenmode by a superposition of just a few HG modes. The fact that our HG basis requires a great number of modes for representing the PSA eigenmodes simply indicates that it is not the optimal basis, as we have discussed in Section 3 after Eq. (28). We have recently obtained some preliminary results on finding the optimum (compact) basis for representation of the PSA modes [34], and the complete study of that will be the subject of a separate publication.
Let us also point out that for the 25×25 μm2 pump spot size, only the fundamental PSA eigenmode #0 sees any significant gain, and this mode has 97% overlap with the TEM00 HG mode. This indicates that, in spite of gain-induced diffraction [20], for tightly focused pumps or long crystals (zR / L ≈0.45 in this case) the crystal’s limited spatial bandwidth forces the PSA to produce squeezed vacuum in a single, well-defined fundamental Gaussian mode. This fact greatly simplifies mode matching of the local oscillator for subsequent homodyne detection and opens the possibilities for generation and observation of significant squeezing factors in such a regime of the PSA operation.
5. Conclusion
To summarize, we have developed and demonstrated an approach to rigorous calculation of the independently squeezed (or amplified) modes of the spatially-broadband PSA. While an order of magnitude estimate of the number of PSA modes is (pump waist × spatial bandwidth)2, our rigorous method determines the exact number and the shapes of the PSA modes supported by the pump beam of a given spot size. This method is the spatial version of the quantum Karhunen-Loève expansion [36] previously used in the temporal domain to study squeezing during quantum soliton propagation [32] and in four-wave-mixing processes [37] in optical fiber. The generated eigenmodes are the spatial analogs of the temporal “supermodes” found in optical parametric oscillators [38]. The PSA eigenmodes are also closely related to the Schmidt modes of spontaneous parametric down conversion [39], previously discussed in both temporal [40] and spatial [41] contexts: at very low PSA gains, the PSA modes correspond to frequency-degenerate transverse Schmidt modes.
The obtained results indicate that the PSA has two important regimes: a) when the pump is tightly focused (zR / L < 1), the limited spatial bandwidth forces the PSA to produce squeezing in a single, well-defined TEM00 Gaussian eigenmode; b) for larger pump spot sizes that support many PSA eigenmodes, the shapes of the most-amplified PSA modes are close to the first few Hermite-Gaussian (for elliptical pump waists) or Laguerre-Gaussian (for circular pump waists) modes. The exact amount of this resemblance (mode overlap) will be the subject of a separate study.
The mode-calculation procedure and the obtained eigenmode shapes are important for both squeezing (finding the maximally squeezed modes as well as maximizing the number of highly-squeezed modes) and noiseless-amplifier (boosting a faint image before detection) applications of the PSA. For the squeezing application, the mode shape gives the profile of the optimum matched local oscillator needed to measure the high-degree of quadrature noise suppression in a travelling-wave PSA. For the image-amplifier application, the optimum performance will be reached when the image modes are mapped onto the PSA eigenmodes, so that they can be amplified without cross-coupling to each other.
This material is based upon work funded by DARPA’s Quantum Sensor Program under AFRL Contract No. FA8750-09-C-0194. Any opinions, findings, conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of DARPA or the U.S. Air Force.
References and links
1. C. M. Caves, “Quantum limits on noise in linear amplifiers,” Phys. Rev. D Part. Fields 26(8), 1817–1839 (1982). [CrossRef]
2. D. Levandovsky, M. Vasilyev, and P. Kumar, “Amplitude squeezing of light by means of a phase-sensitive fiber parametric amplifier,” Opt. Lett. 24(14), 984–986 (1999). [CrossRef] [PubMed]
3. D. Levandovsky, M. Vasilyev, and P. Kumar, “Near-noiseless amplification of light by a phase-sensitive fibre amplifier,” PRAMANA–J. Phys. 56(2-3), 281–285 (2001). [CrossRef]
4. W. Imajuku, A. Takada, and Y. Yamabayashi, “Low-noise amplification under the 3 dB noise figure in high-gain phase-sensitive fibre amplifier,” Electron. Lett. 35(22), 1954–1955 (1999). [CrossRef]
5. P. L. Voss, K. G. Köprülü, and P. Kumar, “Raman-noise-induced quantum limits for χ(3) nondegenerate phase-sensitive amplification and quadrature squeezing,” J. Opt. Soc. Am. B 23(4), 598–610 (2006). [CrossRef]
6. Z. Tong, A. Bogris, C. Lundström, C. J. McKinstrie, M. Vasilyev, M. Karlsson, and P. A. Andrekson, “Modeling and measurement of the noise figure of a cascaded non-degenerate phase-sensitive parametric amplifier,” Opt. Express 18(14), 14820–14835 (2010). [CrossRef] [PubMed]
7. M. I. Kolobov and L. A. Lugiato, “Noiseless amplification of optical images,” Phys. Rev. A 52(6), 4930–4940 (1995). [CrossRef] [PubMed]
8. M. Kolobov, “The spatial behavior of nonclassical light,” Rev. Mod. Phys. 71(5), 1539–1589 (1999). [CrossRef]
9. K. Wang, G. Yang, A. Gatti, and L. Lugiato, “Controlling the signal-to-noise ratio in optical parametric image amplification,” J. Opt. B Quantum Semiclassical Opt. 5(4), S535–S544 (2003). [CrossRef]
10. S.-K. Choi, M. Vasilyev, and P. Kumar, “Noiseless optical amplification of images,” Phys. Rev. Lett. 83, 1938–1941 (1999) [erratum: Phys. Rev. Lett. 84, 1361–1361 (2000)].
11. A. Mosset, F. Devaux, and E. Lantz, “Spatially noiseless optical amplification of images,” Phys. Rev. Lett. 94(22), 223603 (2005). [CrossRef] [PubMed]
12. E. Lantz and F. Devaux, “Parametric amplification of images: from time gating to noiseless amplification,” IEEE J. Sel. Top. Quantum Electron. 14(3), 635–647 (2008). [CrossRef]
13. L. Lopez, N. Treps, B. Chalopin, C. Fabre, and A. Maître, “Quantum processing of images by continuous wave optical parametric amplification,” Phys. Rev. Lett. 100(1), 013604 (2008). [CrossRef] [PubMed]
14. P. Kumar, V. Grigoryan, and M. Vasilyev, “Noise-free amplification: towards quantum laser radar,” the 14th Coherent Laser Radar Conference, Snowmass, CO, July 2007. http://space.hsv.usra.edu/CLRC/presentations/Kumar.ppt
15. Z. Dutton, J. H. Shapiro, and S. Guha, “LADAR resolution improvement using receivers enhanced with squeezed-vacuum injection and phase-sensitive amplification,” J. Opt. Soc. Am. B 27(6), A63–A72 (2010). [CrossRef]
16. O.-K. Lim, G. Alon, Z. Dutton, S. Guha, M. Vasilyev, and P. Kumar, “Optical resolution enhancement with phase-sensitive preamplification,” in Conference on Lasers and Electro-Optics, OSA Technical Digest (CD) (Optical Society of America, 2010), paper CTuPP7.
17. M. Vasilyev, N. Stelmakh, and P. Kumar, “Phase-sensitive image amplification with elliptical Gaussian pump,” Opt. Express 17(14), 11415–11425 (2009). [CrossRef] [PubMed]
18. M. Vasilyev, N. Stelmakh, and P. Kumar, “Estimation of the spatial bandwidth of an optical parametric amplifier with plane-wave pump,” J. Mod. Opt. 56(18-19), 2029–2033 (2009). [CrossRef]
19. A. La Porta and R. E. Slusher, “Squeezing limits at high parametric gains,” Phys. Rev. A 44(3), 2013–2022 (1991). [CrossRef] [PubMed]
20. S.-K. Choi, R.-D. Li, C. Kim, and P. Kumar, “Traveling-wave optical parametric amplifier: investigation of its phase-sensitive and phase-insensitive gain response,” J. Opt. Soc. Am. B 14(7), 1564–1575 (1997). [CrossRef]
21. C. Kim and P. Kumar, “Quadrature-squeezed light detection using a self-generated matched local oscillator,” Phys. Rev. Lett. 73(12), 1605–1608 (1994). [CrossRef] [PubMed]
22. R.-D. Li, S.-K. Choi, C. Kim, and P. Kumar, “Generation of sub-Poissonian pulses of light,” Phys. Rev. A 51(5), R3429–R3432 (1995). [CrossRef] [PubMed]
23. K. G. Köprülü and O. Aytür, “Analysis of Gaussian-beam degenerate optical parametric amplifiers for the generation of quadrature-squeezed states,” Phys. Rev. A 60(5), 4122–4134 (1999). [CrossRef]
24. K. G. Köprülü and O. Aytür, “Analysis of the generation of amplitude-squeezed light with Gaussian-beam degenerate optical parametric amplifiers,” J. Opt. Soc. Am. B 18(6), 846–854 (2001). [CrossRef]
25. C. Schwob, P. F. Cohadon, C. Fabre, M. A. M. Marte, H. Ritsch, A. Gatti, and L. Lugiato, “Transverse effects and mode couplings in OPOs,” Appl. Phys. B 66(6), 685–699 (1998). [CrossRef]
26. M. Annamalai, N. Stelmakh, M. Vasilyev, and P. Kumar, “Spatial modes of phase-sensitive image amplifier with elliptical Gaussian pump,” in Laser Science, OSA Technical Digest (CD) (Optical Society of America, 2010), paper LTuB5.
27. M. Vasilyev, M. Annamalai, N. Stelmakh, and P. Kumar, “Quantum properties of a spatially-broadband traveling-wave phase-sensitive optical parametric amplifier,” J. Mod. Opt. 57(19), 1908–1915 (2010). [CrossRef]
28. Please note that the definition of deff in our prior work (Refs. 17, 18, and 27) is different from that in the present paper. The prior-work deff denotes the quantity that is more commonly known as the effective χ(2) and equals 2deff in the present paper’s notations. As a result, the nonlinear paraxial wave equation in Refs. 17, 18, and 27 does not have the factor of 2 in front of deff. One fallout of this unfortunate choice of notation in our prior work is that Ref. 17 assumes effective χ(2) = 8.7 pm/V for PPKTP crystal, which is about half of the actual value of that crystal’s nonlinearity, and the resulting pump powers listed in Refs. 17 and 26 are four times larger than those required for the same gain in a real PPKTP crystal. The present paper’s definitions rectify the previous inconsistencies.
29. E. Lantz and F. Devaux, “Numerical simulation of spatial fluctuations in parametric image amplification,” Eur. Phys. J. D 17(1), 93–98 (2001). [CrossRef]
30. H. P. Yuen, “Two-photon coherent states of the radiation field,” Phys. Rev. A 13(6), 2226–2243 (1976). [CrossRef]
31. H. P. Yuen, “Multimode two-photon coherent states and unitary representation of the symplectic group,” Nucl. Phys. B 6, 309–313 (1989). [CrossRef]
32. D. Levandovsky, “Quantum noise suppression using optical fibers,” Ph.D. thesis, Northwestern University, 1999.
33. L. Lopez, S. Gigan, N. Treps, A. Maître, C. Fabre, and A. Gatti, “Multimode squeezing properties of a confocal optical parametric oscillator: Beyond the thin-crystal approximation,” Phys. Rev. A 72(1), 013806 (2005). [CrossRef]
34. M. Annamalai, M. Vasilyev, N. Stelmakh, and P. Kumar, “Compact Representation of Spatial Modes of Phase-Sensitive Image Amplifier,” in Conference on Lasers and Electro-Optics, OSA Technical Digest (CD) (Optical Society of America, 2011), paper JThB77.
35. W. Wasilewski, A. I. Lvovsky, K. Banaszek, and C. Radzewicz, “Pulsed squeezed light: Simultaneous squeezing of multiple modes,” Phys. Rev. A 73(6), 063819 (2006). [CrossRef]
36. J. H. Shapiro and A. Shakeel, “Optimizing homodyne detection of quadrature-noise squeezing by local-oscillator selection,” J. Opt. Soc. Am. B 14(2), 232–249 (1997). [CrossRef]
37. C. J. McKinstrie, “Unitary and singular value decompositions of parametric processes in fibers,” Opt. Commun. 282(4), 583–593 (2009). [CrossRef]
38. G. Patera, N. Treps, C. Fabre, and G. J. de Valcárcel, “Quantum theory of synchronously pumped type I optical parametric oscillators: characterization of the squeezed supermodes,” Eur. Phys. J. D 56(1), 123–140 (2010). [CrossRef]
39. A. Ekert and P. L. Knight, “Entangled quantum systems and the Schmidt decomposition,” Am. J. Phys. 63(5), 415–423 (1995). [CrossRef]
40. C. K. Law, I. A. Walmsley, and J. H. Eberly, “Continuous frequency entanglement: effective finite hilbert space and entropy control,” Phys. Rev. Lett. 84(23), 5304–5307 (2000). [CrossRef] [PubMed]
41. C. K. Law and J. H. Eberly, “Analysis and interpretation of high transverse entanglement in optical parametric down conversion,” Phys. Rev. Lett. 92(12), 127903 (2004). [CrossRef] [PubMed]