# Fourier Analysis of Slow and Fast Image Propagation through Single and Coupled Image Resonators

Authors: Parvin Sultana; Takahiro Matsumoto; Makoto Tomita
Journal: International Journal of Optics
 (2012)

---

## Abstract

We applied Fourier space analysis to a comprehensive study of the propagation of pulsed two-dimensional images through single and coupled image resonators.  The Fourier method shows that the image can propagate through the resonator successfully as long as the spatial and temporal Fourier components of the image are within the bandwidth of the amplitude and phase transfer functions. The relevant steep dispersion of the cavity can yield delayed or advanced images. The Fourier method reproduces characteristic aspects of the experimental observations of the image propagation, and also predicts new aspects, such as the spatial image profile dependence on the observation time and the coupling strength. To demonstrate the time evolution of the experiment, space- and time-resolved image propagations were performed using a streak camera.


---

## Body

## 1. Introduction

When propagating through media, the velocity of an optical wave packet depends on the variation of the refractive index. In recent years, there has been increasing interest in media in which the refractive index changes drastically with respect to changes in frequency, resulting in significant modification in the corresponding group velocities of the optical pulses. We refer to this phenomenon as “slow light” or “fast light.” Experimental demonstrations of slow and fast light have been realized based on stimulated Brillouin scattering in optical fibers [1], electromagnetically induced absorption (EIA) and transparency (EIT) in atomic vapors [2–5], structural dispersion in microsphere-optical fiber systems [6], coupled-resonator-induced transparency (CRIT) with classical light [7], and squeezed vacuum studies [8]. Recently, various groups have reported stopping light altogether, which leads to new methods of storing information [9, 10].It is challenging to apply slow and fast light techniques to two-dimensional (2D) image propagation. All-optical methods for delaying images have great potential in many applications, including holography and optical pattern correlation measurements or futuristic all-optical image routers for large capacity image information.  The first experimental demonstration of an all-optical image delay was reported in Cs atoms [11]. The storage and retrieval of a transverse image in an atomic vapor was demonstrated using a technique based on EIT [12, 13]. Recently, researchers succeeded in eliminating arbitrary paraxial image diffraction by manipulating the susceptibility in momentum space [14]. This work opens up new possibilities for utilizing EIT to realize high-resolution imaging and position-independent beam expansion.We reported an optical image propagation experiment using single and coupled image resonators. Two-dimensional images encoded on optical pulses were either advanced or delayed through a single image resonator using relevant dispersion [15]. Further, 2D tunable delayed images were demonstrated experimentally using CRIT. The coupled resonator is recognized as the all-optical analogue of EIT [16]. A steep positive slope induced by the transparency effect of the second cavity, in otherwise opaque absorption from the first cavity, allows the observation of the slow image. Because the image resonator was based on structural dispersion, it had many advantages compared with atomic systems.  In addition to the relative design simplicity, the advantages of the photonic structures are that the available frequencies are not restricted to intrinsic atomic transitions; an additional laser beam is not required to prepare the slow light medium; the system is not affected by Doppler motion. Although there have been studies demonstrating pulsed 2D image propagation, until now there has been no comprehensive theoretical analysis for image propagation through the single and coupled resonators, specifically from the dispersion point of view.  In contrast to real-space analysis [15, 16], in this study we applied Fourier space analysis to study the propagation of 2D images encoded on temporal pulses through single and coupled image resonators, comprised of Fabry-Perot resonators. Real-space analysis provided an intuitive interpretation; however, this approach is not easily applied to complicated images or coupled resonators. The Fourier method reproduces the characteristic aspects of the experimental observations on the image propagation reported to date and also predicts the time evolution of image propagation through the cavity.Fabry-Perot resonators are widely used experimentally and often discussed as a model cavity. The longitudinal domain effects on frequency filtering, the laser cavity, the phase-amplitude converter, and nonlinear characteristics have been studied extensively. However, little is known on the general transverse effect of a cavity on an image, particularly from a dispersion point of view. In this study, a planar cavity is treated as a parallel processing system with image transmission. Rays making a small angle with respect to the planar cavity axis escape from the cavity after a great number of reflections. However, as long as the spatial and temporal Fourier components of the image are within the bandwidth of the amplitude and phase transfer functions, the images can successfully propagate through the cavity. It has been discussed that particular types of cavities, such as confocal or hemiconfocal cavities with transverse degenerate cavity modes, could transmit a portion of the image [17, 18]; this provides interesting possibilities regarding the future of image resonators.


## 2. Fourier Analysis for Image Propagation

We applied Fourier space analysis to the wave vector (k)-frequency (ν) domain.  Figure 1(a) shows a schematic of a coupled Fabry-Perot cavity consisting of three plane parallel mirrors with reflectivity Ri=riri*=ri′ri′* and transitivity Ti=titi*=ti′ti′*, where i=1,2, and 3 correspond to the first, second, and third mirrors, respectively.  This configuration yielded the coupled-resonator-induced transparency.  The single cavity configuration can also be treated by setting R3=0. Figure 1(b) is an example of delayed pulsed 2D image propagation through the coupled-resonator-induced transparency. The insets (top) show the reflected image at the center of the transparency window (CRIT), indicated by the red arrow in the frequency spectrum (bottom). The pulse peak was delayed by 10.6 ns, due to the steep dispersion of the resonator. The incident electric field of the image with wave vectors kx,ky and frequency ν is denoted as A0(kx,ky,ν); the output electric field is given as A(kx,ky,ν). The notation B(kx,ky,ν), B0(kx,ky,ν), B1(kx,ky,ν), and B2(kx,ky,ν) refers to the electric fields in the first cavity; C(kx,ky,ν), C0(kx,ky,ν), C1(kx,ky,ν), and C2(kx,ky,ν) correspond to the electric fields in the second cavity, as shown in Figure 1. The electric fields of the cavities are given by
(1)A=A0r1+B0t1′,B=A0t1+B0r1′,B2=B1r2+C0t2′,C=B1t2+C0r2′,C2=C1r3,
where
(2)B1=Bexp[-ρ1+iϕ1],B0=B2exp[-ρ1+iϕ1],C1=Cexp[-ρ2+iϕ2],C0=Cr3exp[-2ρ2+2iϕ2].
Here, ρ1,ρ2 are the losses inside the cavity. The phase relation for the two cavities can be written as ϕj(kx,ky,ν)=iLj(k02-kx2-ky2)1/2, where i is the imaginary number, j=1,2, and k0=2πν/c, where c is the velocity of light, and ν is the frequency of the laser. One part of the incident light does not enter the medium, but instead is reflected from the boundary of the first mirror surface and appears directly at the output of the system; we refer to this part as the direct component. The other part of the incident light enters the cavities, where it is reflected multiple times before exiting through the first mirror as a multiple-reflected component. The total electric field at the output of this system is represented as the sum of the direct- and multiple-reflected components. Therefore, the resultant reflected image is given by
(3)T(kx,ky,ν)={r1-t1t1′T2e-2ρ1+2iϕ11-r1′T2e-2ρ1+2iϕ1}=|T(kx,ky,ν)|eiθ(kx,ky,ν),
where
(4)T2(kx,ky,ν)=r2-t2t2'r3e-2ρ2+2iϕ21-r2'r3e-2ρ2+2iϕ2.
Here, T(kx,ky,ν) is the amplitude transfer function of the coupled image resonator. The first term in (3) is the direct-reflected component, and the second term is the sum of the multiple-reflected components. For the single image resonator, the third mirror was removed (i.e., R3=0). Hence, multiple components reflected back and forth between the first and second mirrors that made up the single cavity.(a) Schematic illustration of the reflected and transmitted image propagation through a coupled Fabry-Perot image resonator. In this figure,R1, R2, and R3 are the reflectivities of the three mirrors that comprise the coupled cavity. (b) An example of delayed pulsed 2D image propagation through the coupled-resonator-induced transparency. The insets show the reflected image (top) at the center of the transparency window (CRIT), indicated by the red arrow in the frequency spectrum (bottom).
(a)
(b)Provided that the system is linear, the transfer function of (3) operates on the input image to produce an output image. We can write the Fourier transformation of the real space image encoded on optical pulses A0(x,y,τ) as
(5)A0(kx,ky,ν)=∭-∞∞A0(x,y,τ)ei(kxx+kyy)e-i2πτνdxdydτ.
By transforming A0(x,y,τ), we are basically decomposing the input image into the complex exponential functions of the various spatial and temporal frequencies (kx,ky,ν). Multiplication of the input spectrum by the transfer function then takes into account the effects of the system on each Fourier component. As a result, when an image propagates through the cavity, the amplitude and phase of each Fourier component will be modified as follows:
(6)A(kx,ky,ν)=A0(kx,ky,ν)T(kx,ky,ν).
Finally, taking the inverse Fourier transform of the output spectrum produces the real output image [19]:
(7)A(x,y,τ)=∭-∞∞A(kx,ky,ν)e-i(kxx+kyy)ei2πτνdkxdkydν.We calculated the image propagation through the single and coupled image resonators using (5)−(7). A variety of configurations can be considered for both single and coupled image resonators, depending on the coupling conditions of the cavities. Both slow and fast images come from the steep slope of the dispersion. In the reflection geometry, we obtained both anomalous and normal dispersion that depended on R1 and R2. When R1>R2, the coupling is weak compared with the loss; this is known as the under coupling condition. Under this condition, a negative gradient ∂θ/∂ν<0  (1) appears in the dispersion curve within the resonance dip. Such a slope results in anomalous dispersion, and a fast image is expected. When R1<R2, the coupling is strong compared with the loss; this is referred to as the over coupling condition. In this case, the phase increases monotonically as the frequency of the laser light is increased, that is, ∂θ/∂ν>0, and slow light is expected. When R1=R2, the coupling and the loss become equal, and the reflection minimum becomes zero. This is known as the critical coupling condition.Figures2(a)–2(c) show the intensity transfer function for single and coupled resonators in k-space and ν-space |T(k,ν)|2. The k-space cross-section of the intensity transfer function, that is, |T(k,0)|2, is shown along the vertical axis. Figures 2(d)–2(f) are the ν-space cross-sectional profiles, |T(0,ν)|2, and Figures 2(g)–2(i) are the reflected phase shift as a function of the detuning frequency. To discuss the image’s optical resolution limitations, we defined the cavity bandwidths of the transfer function |T(kx,ky,ν)|2 in k and ν-space as Δk and Δν, respectively. The Fourier components of the image outside the bandwidth were lost or phase-shifted, and the final images showed deterioration. Figure 2(b) shows the intensity transfer functions for the coupled resonator at Δν12=ν2-ν1=0, that is, when two cavities are adjusted to be on-resonant by tuning the cavity lengths. Here, ν1 and ν2 are the resonance frequencies of the first and second cavities, respectively. When Δν12=0, a narrow transmission peak appeared at δν=ν-ν2=0 due to the destructive interference between the two resonators, exhibiting CRIT [16]. Consider the situation in which the two cavities have different resonance frequencies, that is, Δν12≠0, as shown in Figure 2(c). For Δν12≠0, a sharp asymmetric structure appeared, due to the Fano effect. For Δν12>0, this asymmetry increased the relevant k-space bandwidth; this may be useful for fine image propagation. Using the parameters from Figures 2(b) and 2(c), the k-space bandwidth (Δk) can be increased from Δk=9.8 mm−1 to Δk=13.3 mm−1. It is depicted in Figures 2(b) and 2(c) by the black double-arrow line. The k-space bandwidth for the single resonator in Figure 2(a) is Δk=10.1 mm−1.Figure 2
(a)−(c) Calculated intensity transfer functions|T(k,ν)|2 for (a) a single image resonator, (b) coupled image resonator with Δν12=0, and (c) coupled image resonator with Δν12≠0. The black solid curves along the k-axis in Figures (a)–(c) are the k-space cross-sections |T(k,0)|2. (d)−(f) ν-space cross-sections |T(0,ν)|2. (g)−(i) Calculated phase shift as a function of the detuning frequency. Note that, in (a), the reflectances of the mirrors were R1=0.9 and R2=0.7, and the cavity length was L1=10 cm; in (b) and (c), the reflectances were R1=0.5, R2=0.9, and R3=0.99, and the cavity lengths were L1=5 cm and L2=10 cm, respectively. The k-space cavity bandwidths in (a), (b), and (c) were Δk=10.1 mm−1, Δk=9.8 mm−1, and Δk=13.3 mm−1, respectively, as indicated by the solid black double-arrow lines.


## 3. Stationary Image Propagation through  Single and Coupled Resonators

Using the transfer function shown in Figure2, we first examined the propagation of 2D images E(x,y) through single and coupled image resonators. We used a continuous wave input beam for the transmission of a monochromatic image in the x-y domain. Figure 3(a) is the original image injected into the cavity. Figure 3(b) is the Fourier-transformed input image. The full width at half maximum of the Fourier image was δk=4.9 mm−1. Figures 3(c)–3(e) are the reflected images obtained from the single resonator using a continuous wave input beam for the under, critical, and over coupling conditions, respectively. Several characteristic aspects of the 2D image propagation that was obtained experimentally [15, 16] were also reproduced in the Fourier analysis. We observed that the edge region of the image was enhanced for the under coupling condition (Figure 3(c)). This edge-enhanced effect was most prominent for the critical coupling condition (Figure 3(d)), in which we observed edge filtering. For the over coupling condition, the hollow effect was observed, that is, a superimposed outline structure appeared around the original image (Figure 3(e)). Figures 3(f) and 3(g) are the output reflected images observed from the coupled image resonator (R1=0.5, R2=0.9, R3=0.99, and L1=5 cm, L2=10 cm, Δk=9.0 mm−1, Δν=80.6 MHz) using a continuous wave input beam and temporal pulses of 25 ns (δν=17.6 MHz), respectively. It is clear from Figures 3(f) and 3(g) that the pulsed image can propagate through the image cavity at almost the same quality as the steady-state image, provided that the spatial and temporal Fourier components of the images are within the bandwidth of the cavity.Figure 3
(a) Original image injected into the cavity. (b) Fourier-transformed input image. The full width at half maximum of the Fourier image isδk=4.9 mm−1. (c)–(e) Reflected images obtained using continuous wave input light for under (R1=0.7, R2=0.5, Δk=13.0 mm−1), critical (R1=0.5, R2=0.5, Δk=9.0 mm−1), and over coupling (R1=0.5, R2=0.7, Δk=13.0 mm−1) conditions, respectively. (f) and (g) Reflected output images from the coupled image resonator (R1=0.5, R2=0.9, R3=0.99, L1=5 cm, L2=10 cm, Δk=9.0 mm−1, Δν=80.6MHz) using a continuous wave input beam and temporal pulses of 25 ns (δν=17.6 MHz), respectively.It is convenient to separate the direct- and multiple-reflected components to interpret the spatial and temporal characteristics of the propagated image. The total output image is the sum of the direct component and theπ(rad) phase-shifted and time-delayed multiple-reflected components [15]. Because the first mirror reflects all of the k-space Fourier components, the k-space bandwidth is effectively infinite for the direct component. On the other hand, for multiple components, the bandwidth is restricted by the cavity resonance mode; hence, the higher order components are lost.  As a consequence of the limited bandwidth, the edge region of the resultant multiple-reflected components is broadened. This was apparent in the single resonator system, as shown in Figures 3(c)–3(e). Multiple-reflected components were the major components for the over coupling condition (R1<R2); hence, the spatially broadened multiple component formed the hollow structure around the original image (Figure 3(e)). In contrast, for the under coupling condition (R1>R2), the direct component played the major role. As a result, the broadening effect on the edge region of the multiple components was weak. The destructive interference between the direct and multiple components in the edge region was not perfect, resulting in enhancement of the edge structure (Figure 3(c)). This effect was most evident under the critical coupling condition (Figure 3(d)).


## 4. Time Evolution of the Pulsed Image Propagation

We next considered the time evolution of the 2D image encoded on the pulses,E(x,y,t), using the transfer function shown in Figure 2. Figure 4(f) illustrates the cross-section that was considered for time evolution analysis. Although we were mostly interested in the time evolution of 2D image propagation, it was convenient to show certain cross-sections of the image over the course of the measurement, to observe the other sections of the image, which behaved in a similar manner. Figures 4(a)–4(d) show the Fourier analysis for the time evolution of the cross-section in reflected and transmitted images in the space-time (x-t) domain, using a single resonator. Figure 4(e) shows the input image of width 1.4 mm (the horizontal axis) encoded on the temporal Gaussian pulse of 25 ns (the vertical axis). Figures 4(a), 4(b), and 4(c) show the reflected images from the single resonator for under, critical, and over coupling conditions, respectively. The vertical profiles on the left side of Figures 4(a)–4(d) show the temporal profiles of the input and output images at the center position, indicated by the white and yellow curves, respectively. For the under coupling condition (Figure 4(a)), the temporal profile of the output image was advanced compared with the input image, thus, showing a manifestation of the fast image. Figure 4(c) shows the time evolution calculated for the over coupling condition. The hollow structures around the original image consisted of diffracted multiple components only, simply delayed by the effective optical path length. A similar effect was also shown in Figure 3(e). The center of the cross-section was more delayed as a result of the destructive interference between the direct and multiple components, which is a manifestation of the slow image. Under the critical coupling condition (Figure 4(b)), the temporal profile of the output image was split as a result of the destructive interference effect of the direct- and multiple-reflected components.Figure 4
(a)−(d) Numerical simulation of the reflected and transmitted images in the space-time (x-t) domain for a single resonator. (a)–(c) Reflected images from the single resonator for under, critical, and over coupling conditions, respectively. (d) Transmitted image from the single resonator. (e) Input image of width 1.4 mm (the horizontal axis) encoded on the temporal 25 ns Gaussian pulse (the vertical axis). (f) Illustration of the cross-section of the 2D image that was considered for time evolution analysis. The earlier, intermediate, and later time profiles are indicated by the notations 1, 2, and 3, respectively, in Figures 4(a)–4(d).In Figures4(a)–4(d), the spatial profiles at earlier, intermediate, and later times are plotted on the right hand side, indicated by the notations 1, 2, and 3, respectively. It is interesting to note that the time-resolved spatial profiles transformed continuously from under to critical to over coupling conditions.  The spatial profile at the earlier time for the under coupling condition (Figure 4(a), profile 1) showed almost no diffraction; the diffraction effect increased from intermediate to later time profiles (Figure 4(a), profiles 2 and 3). This is reasonable because the profiles at the later time traveled effectively longer propagation distances. The later time profile for the under coupling condition (Figure 4(a), profile 3) was followed by the intermediate time profile for critical coupling condition (Figure 4(b), profile 2). The image further developed from intermediate to later time profiles with increasing diffraction effect (Figure 4(b), profiles 2 and 3). Finally, the later time profile for the critical coupling condition (Figure 4(b), profile 3) was followed by the intermediate time profile for the over coupling condition (Figure 4(c), profile 2). This time development is reasonable because the diffraction effect from the multiple-reflected components increased as the observation time or coupling strength increased.A similar Fourier analysis approach was applied to the coupled image resonator. Figures5(a) and 5(b) are simulations of the reflected cross-sectional images of different sizes from the coupled resonator. In Figures 5(a) and 5(b), images with a cross-section of 1.4 mm and 2.8 mm were used as inputs, respectively. The cavity parameters were the same as those shown in Figures 3(f) and 3(g). The full widths of the Fourier images were δk=5.5 mm−1 and δk=2.7 mm−1, respectively, and the k-space bandwidth of the cavity was Δk=9.0 mm−1. The frequency bandwidths of the image and cavity were δν=17.6 MHz and Δν=80.6 MHz in both cases. The image that propagated through the coupled resonator had aspects of both under and over coupling conditions: the edge-enhancement effect at the earlier time due to the under coupling condition in the first cavity and the delayed time profile due to the over coupling condition in the second cavity.  In Figure 5(a), we observed a striped structure. When the image size was reduced, the higher Fourier components of the image overlapped the higher cavity resonance modes, which led to the striped effect.Figure 5
(a) and (b) Reflected images from the coupled resonator for different image sizes. The sizes of the images used in (a) and (b) are 1.4 mm and 2.8 mm, respectively, the full widths of the Fourier images areδk=5.5 mm−1 and δk=2.7 mm−1, respectively, and the spatial bandwidth of the cavity is Δk=9.0 mm−1. The notations 1, 2, and 3 correspond to the earlier, intermediate, and later time profiles, respectively, in both figures. In both cases, the cavity parameters (R1=0.5, R2=0.9,  R3=0.99,  L1=5 cm, L2=10 cm) are the same as those used for the simulations shown in Figures 3(f) and 3(g).


## 5. Experimental Demonstration

The present Fourier analysis not only reproduces the reported experimental observations [15, 16] but also predicts new aspects of the time evolution of image propagation through the cavity; the propagated image profiles are time dependent. To demonstrate the new aspect of the time evolution,we performed time- and space-resolved image propagation experiments. The experimental setup is illustrated in Figure 6(a). While various configurations are possible in the image resonators, we studied a particular case using an over coupled single resonator system corresponding to the simulation of Figure 4(c), because the time evolution effect appeared clearly in this system. The light source for the experiment was a continuous wave Nd+3:yttrium aluminum garnet (Nd:YAG) laser. The Nd:YAG laser, operated with a linewidth of 1 kHz, produced second-harmonic 532 nm light. The laser beam was expanded using two lenses (not shown in Figure 6(a)) to a diameter of 3 cm. An image mask with the shape of the symbol “ψ” was inserted into the collimated beam. The parameters used in this experiment were R1=0.7,R2=0.9, and L=30cm, which corresponded to Δk=5.9 mm−1. An image mask of width 1.6 mm (δk=4.7 mm−1) was used for this experiment, and the injected temporal pulse width was 40 ns (δν=11.0 MHz). While in the previous experiments [15, 16], the delay time was examined at a specific spatial point, here, we used a streak camera with a time resolution of 10 ps. The streak camera was a suitable device for the present experiment because it could analyze both the space and time domains simultaneously. The variable slit at the streak camera entrance was used to select the vertical height of the image.(a) Experimental setup for time- and space-resolved image propagation through the image resonator. Light pulses of 40 ns duration passed through an image mask and an image resonator system. The reflected images were sent to the streak camera using a half mirror (HF) to analyze time- and space-resolved reflected images. (b) (A) and (b) (B) show the reflected image for the off-resonance and on-resonance conditions, respectively. (b) (C) is the experimental observation of the temporal profiles of the reflected image. The dashed black curve is at the center of the off-resonance image; the solid blue and red curves are at the edge region and center region of the on-resonance image, respectively. (b) (D) is the cross-sectional profile of image (b) (B). The notations 1, 2, and 3 indicate the earlier, intermediate, and later times of the image, respectively.
(a)
(b)Figures6(b) (A) and 6(b) (B) show the experimental observation of the reflected off-resonance and on-resonance images, respectively. The experimental result shown in Figure 6(b) (B) corresponds to the simulation in Figure 4(c). In Figure 6(b) (C), the solid blue and red profiles are the results observed from the edge region and central region of the on-resonance image, respectively. The black-dashed line is the temporal profile obtained from the off-resonance image of Figure 6(b) (A). Profiles (1), (2), and (3) are the cross-sectional profiles of the earlier, intermediate, and later times of the image in Figure 6(b) (B), respectively. The observed propagated spatial patterns were time dependent. It was also noticeable from Figure 6(b) (B) that the central region of the image was more delayed (upshifted), compared with the edge and off-resonance image. This result showed good agreement with the discussion for the over coupling condition (Figure 4(c)) and experimentally verified the analysis based on the Fourier approach. The off-resonance image (Figure 6(b) (A)) showed some fringe structures. This structure may be attributed to free space propagation from the image mask to the resonator and from the resonator to the streak camera.


## 6. Conclusion

In summary, we analyzed the propagation of 2D images encoded on temporal pulses through single and coupled resonators in Fourier space, using amplitude and phase transfer functions. The resonators were capable of propagating images, either delayed or advanced, as long as the spatial and temporal Fourier components of the image were within the bandwidth of the resonator. We also studied the observation time and coupling strength-dependent spatial image profiles. It has proven difficult to this point to increase thek-space bandwidth for realistic image applications. We can consider several possible solutions for this dilemma. One such possibility would be the use of a transverse degenerate paraxial stable Fabry-Perot cavity, consisting of confocal mirrors [18]. Previous studies have shown that such a cavity can transmit the self-transformed portion of the image [17]. Another possibility would be to introduce image guide photonic structures, such as those used in endoscopes or fiber image plates. If we can introduce a cavity structure into these image guide photonic structures, then this may improve the k-space bandwidth significantly, allowing the transmission of fine images. It should be noted that, in this case, the steep dispersion of the cavities could also yield fast and slow images, as discussed in this paper.



---

*Source: 960985-2012-11-13.xml*