|
1.IntroductionTwo-wave mixing, as an interesting phenomenon of great research and application value in nonlinear optics, has received growing interest in the past four decades. In a photorefractive crystal, two-wave mixing can be explained as the nonlinear interaction of intensity and phase between two incident beams through the photorefractive effect.1–3 The diffusion effects in the crystal are not negligible, since there is a strong possibility that beam shape distortion will occur in a diffusion-dominated photorefractive material, and this has also been observed in the previous experimental results of volume holographic recording and reconstruction using a Fe-doped crystal.4 Plenty of effective applications of two-wave mixing have been applied to different fields, such as real-time holography, self-pumped phase conjugation, coherent image amplification, optical storage, etc.5–7 For real-time beam shaping, laser beam profile may be monitored during pulse duration, researchers suggested adaptive method, based on dynamic holographic gratings, recorded in photosensitive materials.8 As for two-dimensional (2-D) image amplification, it is a practical way to make use of photorefractive crystals in two-wave mixing arrangements because of the relative convenient operation and the potential for high gain.9–11 To achieve high-quality image amplification, not only the value of gain needs to be considered but also the uniformity of the amplification for each pixel in the image, which takes account of the difference in spatial frequency in a beam. The gain nonuniformity will introduce serious wavefront distortion and apparent polarization state degradation to the output beam. Especially, the induced perturbation between the spatial intensity and phase of the beam, which increases rapidly with the propagation in nonlinear media, will cause small-scale self-focusing in the beam and even devastating damage to the elements in the light path. Therefore, it is of great significance to explore the uniformity of image amplification by two-wave mixing. However, most research on coherent image amplification by two-wave mixing in the previous literature only emphasized on the condition of a single spatial frequency.1,12–15 In this article, the uniformity of image amplification by two-wave mixing in photorefractive crystal, which takes the difference in spatial frequency in a beam into account, has been studied in detail. 2.Theoretical AnalysisThe model of two-wave mixing in the previous literatures, which are concentrated on the situation for a single spatial frequency, has been presented in Fig. 1. As seen in Fig. 1, the two beams and enter into a photorefractive crystal symmetrically. The interference intensity and refractive index distributing equations 10 can be expressed as where is the intensity modulation depth; is the saturation value of photoinduced refractive index change; is the value of index gratings wave vector; is the phase offset between the interference intensity pattern and refractive index gratings.In order to study the uniformity of amplification by two-wave mixing in photorefractive crystal, a model that takes account of the difference in spatial frequency in a beam has been built. And the geometrical configuration has been presented in Fig. 2. As seen in Fig. 2, is the optical axis of the crystal, which is parallel to -axis. Assuming the pump and signal beam, both of which are superposition of a series of plane waves, enter into the medium symmetrically with respect to the normal (-axis) from the left side (). And both of them are -polarized (electrical field vector is perpendicular to the incident plane). Sometimes, two lenses are being used to focus energy of the two beams to the crystal (Fourier plane of the lenses), the plane waves are converted to spherical waves. But, here we just consider the situation of plane wave incidence. The two beams couple with each other in the crystal, which leads to the amplification of the signal beam at the expense of the pump beam intensity decrease. For both the pump and signal beam, the optical field analytical expression on each pixel can be represented by and ; separately, where () indicates the spatial locations of light on individual pixel in the pump beam, and () similarly in the signal beam; ; ; ; , ; and are separately the wave vectors on each pixel in the pump and signal beam, which can be expressed with the coordinate system in Fig. 2. To explore the gain distribution in the signal beam, the difference among incident angles for each pixel has been taken into account for both the pump and signal beam. As a small divergence angle will be induced with the spatial propagation in Gaussian beam shot from a laser,16 the pump and signal waves span a rather small angular range, which results in different spatial frequencies in a beam. The offsets of the propagation direction on each pixel to the average propagation direction are represented by and separately for the pump and signal waves. Therefore, when entering into the crystal, the incident angle on each pixel in the pump and signal beam can be represented separately by and , where and indicate the average incident angle of the pump and signal beam, respectively. Ordinarily, is assumed. As is known, the far-field small divergence angle of laser beams can be figured out by the formula16 where is the resonant cavity length of the laser, is the wavelength of the beams. Consequently, for the pump and signal waves, and vary over a range of separately, i.e., the incident angle on each pixel in the signal beam range from to .The interference intensity of the pump and signal waves in the crystal, which has been deduced from Eq. (1), can be expressed as where is the total intensity of the incident pump and signal beam; c.c. indicates the complex conjugate of the previous term.The charge carriers resulted from interference intensity pattern are liberated and redistribute themselves with the effect of drift or diffusion. A charge imbalance, which results in a modulated electric field, occurs and thereby, alters the spatial distribution of the refraction index in the crystal through the electro-optic effect.17 Since the angular range of the incident beam is much smaller than the average incident angle for both the pump and signal beam [i.e., and ], the interference between the pump waves, as well as that between the signal waves, will induce index gratings which are much weaker than those written by the interference between the pump and signal waves.18 This is the scenario that can be approximately simplified as that, only the index gratings that written by the interference between the pump and signal waves make contribution to the energy transfer between the two beams. A phase offset between the interference intensity pattern and the index gratings may exists in the crystal. When no external electric field is applied (i.e., diffusion is the dominated effect for the transport of the charge carriers), the phase offset becomes , in which case the energy transfer reaches its peak whereas the phase coupling disappears.19,20 Therefore, based on Eq. (2), the refractive index distribution can then be approximated as where is the refractive index of the crystal when no light is present; is the modulation factor of refractive index, which depends on the spacing and direction of the grating, as well as on the material properties of the crystal, e.g., the electro-optic coefficient. Equation (5) is solved for the steady state so that the amplitudes and are taken to be time independent.With this and the slowly varying amplitude approximation,21 Maxwell’s scalar equations yield the following coupled amplitude equations where is the absorption coefficient of the crystal, which is regulated by the wavelength of the incident beam.Since the relation between the intensities and amplitudes of the incident beams are and , the coupled intensity equations are given by In the assumed diffusion driven scenario, the phases of the beams are decoupled so that the intensity Eqs. (8) and (9) describe the two-wave mixing process completely. With Eqs. (8) and (9), one can integrate to yield the intensity versus effectively acting distance on each pixel in the signal beam, which can be normalized by the initial intensity as Therefore, the gain, which is defined as the intensity ratio of the output signal wave in the presence of a pump beam to that in the absence of a pump beam, can be given by where is the thickness of the crystal, is the angle between the propagation direction and the normal (-axis) on each pixel in the signal beam inside the crystal, and is the effective interaction length on individual pixel in the signal beam.To calculate the gain, the expression of intensity coupling coefficient ,22 presented below, can be inserted into Eq. (11): In Eq. (12), the related parameters are included in coefficients and , which can be given by where is the effective electric–optic coefficient; is the refractive index; is the recombination constant, which is fixed when the value of is set; is the effective charge density, which is regulated by the wavelength of the incident beams; is the Boltzmann constant; is the absolute temperature; is the charge on the electron; is the wavelength of the incident beams.According to the analysis above, the gain nonuniformity in the signal beam, which indicates the unevenness for image amplification, is induced by the difference in spatial frequency in the beam. As seen in Eqs. (11)–(14), the factors, which impact on the uniformity of amplification in the signal beam, are the crystal thickness , the incident angle , and the wavelength . 3.Simulation and DiscussionTo explore the effects of these factors on the uniformity of amplification in the signal beam, which loads a 2-D image, one can vary the average incident angle , the crystal thickness and the wavelength of incident beam , respectively, and observe the behavior of the gain. The nonuniformity of gain can be measured with the relative standard deviation . Here, AVG and SD are the mean value and standard deviation of the gain in the signal image beam, respectively. The gain nonuniformity is negligible compared with the value of gain when reaches a low enough value. The usual requirement is . In this article, the numerical simulation for two-wave mixing is based on the model analyzed above. Assuming the signal beam and the pump beam are split from a laser beam at a ratio of 1:1120. Then, the two of them recombines at a photorefractive crystal from the same side with the average incident angle . Three scenarios with different wavelengths have been analyzed to obtain the gain distribution in the signal image beam. 3.1.Argon Laser at nmIn the first scenario, the light source is an argon laser with a cavity length of operating at wavelength of 488 nm, in which case the crystal has the parameters of and according to Ref. 22. Substituting and into Eq. (3) yields the far-field small divergence angle of the laser beam: . Arbitrary values within range can be assigned to , the offset of the propagation direction on individual pixel to the average propagation direction in the signal image beam, to work out the gain with Eqs. (11)–(14). Computer-generated plots of the main value of gain AVG and nonuniformity parameter versus average incident angle and crystal thickness have been presented in Fig. 3. The conditions that the value of gain is (i.e., the intensity of the signal image wave is reduced in the crystal23) have been ignored. As seen in Fig. 3(a), for different values of , AVG represents difference in fluctuating amplitude while similar in trend of variation with the change of . That is, the effects of angle on AVG are similar for different values of . The value of AVG increases with the increasing of , first reaches a maximum, then decreases. For each value of , AVG reaches its maximum with a same fixed value of , i.e., 17 deg. The similar trend of variation can be seen in Fig. 3(b) for value of . For arbitrary value of , reaches its minimum with the same value of , and this is due to the minimal gradient of AVG. When the average incident angle is 17 deg, the incident angle of signal and pump beam on each pixel are around the value of 17 deg, which leads to a minimal fluctuation of AVG. As for this, the image can achieve the maximum gain and the highest level of amplification uniformity when the incident angle of beams is set as 17 deg. The values of AVG, SD, and versus for the optimal average incident angle are presented in Fig. 4. The situation for should be ignored due to its reduction of signal wave intensity in the crystal. Within the 0- to 0.33-cm range of , increases with the increasing , then decreases to a minimum and afterward increases again. Particularly, the optimal thickness of the crystal , in which case reaches its minimum 0.00016%, is turned out to be 0.17 cm. Moreover, the maximum gain has been reached with simultaneously. 3.2.Argon Laser at nmThe second scenario is that the operating wavelength of the argon laser used in the first scenario is changed to 514.5 nm, in which case the crystal has the parameters of and according to Ref. 22. What is more, the far-field small divergence angle of the laser beam is turned into . Substituting , , and into the simulation carried on in the first scenario yields the gain distribution for each pixel in the signal image beam. The values of AVG and versus angle and are shown in Fig. 5. According to the results, the phenomenon that AVG and represent difference in fluctuating amplitude while similar in trend of variation with the change of for arbitrary value of , is the same as that seen in the first scenario. The optimal average incident angle , which will guarantee the maximum gain and the highest level of amplification uniformity regardless of the value of , is appeared to be 15 deg. Simulation results of AVG, SD, and versus for are presented in Fig. 6. Within the 0 to 0.43 cm range of , increases with the increasing , then basically remain unchanged. Since the signal beam has not been amplified when the gain is , the situation with should be ignored. 3.3.He–Ne Laser at nmThe last scenario is that, the argon laser used in the first scenario is changed for a He–Ne laser with a cavity length of at wavelength 632.8 nm, in which case the crystal has the parameters of and . Besides, the far-field small divergence angle of the laser beam is turned out to be . Similarly, bring , , and into the simulation carried on in the first scenario yields the gain distribution on each pixel in the output signal beam. The results of AVG and with varying angle and are presented in Fig. 7. Note that the fluctuation amplitude of is larger than that in the previous two scenarios. The optimal average incident angle is 4 deg. The values of AVG, SD, and versus for have been presented in Fig. 8. increases with the increasing , then decreases to a minimum and afterward increases again. The minimum value of (0.000184%) can be reached with . By comparing the three scenarios above, the influence of the wavelength on the uniformity of image amplification can be observed. For different wavelengths of the incident beams, the values of fluctuate within different ranges and the optimal values of nonuniformity that can be reached are different. 4.ConclusionsIn short, to explore the uniformity of image amplification by two-wave mixing in photorefractive crystals, a model in view of different spatial frequency in a beam has been built. According to the gain of signal beam deduced from the model, it is the difference in spatial frequency attributed to different incident angle on each pixel in the beam that induces the unevenness of image amplification. It has been proved that the factors which impact on the uniformity of image amplification are: the thickness of photorefractive crystal, the incident angle on each pixel in signal image beam, and the wavelength of incident beams. As an example, the effects of these factors in photorefractive crystal have been researched through simulation. It turns out that a considerable nonuniformity of can be reached when the wavelength of the incident beam is 488 nm (at , ). AcknowledgmentsThis work was supported by the National Natural Science Foundation of China (No. 61071190, No. 61171158) and the Fundamental Research Funds for the Central Universities of China (No. CDJZR12160010). ReferencesN. V. Kukhrarevet al.,
“Holographic storage in electro-optic crystals. I. Steady state,”
Ferroelectrics, 22
(1), 949
–953
(1978). http://dx.doi.org/10.1080/00150197908239450 FEROA8 0015-0193 Google Scholar
M. C. Sørenet al.,
“Two-wave mixing in a broad-area semiconductor amplifier,”
Opt. Express, 14
(25), 12373
–12379
(2006). OPEXFF 1094-4087 Google Scholar
M. C. Sørenet al.,
“Nonlinear gain amplification due to two-wave mixing in a broad-area semiconductor amplifier with moving gratings,”
Opt. Express, 16
(8), 5565
–5571
(2008). http://dx.doi.org/10.1364/OE.16.005565 OPEXFF 1094-4087 Google Scholar
P. P. Banerjeeet al.,
“Volume holographic recording and readout for 90-deg geometry,”
Opt. Eng., 43
(9), 2053
–2060
(2004). http://dx.doi.org/10.1117/1.1774195 OPEGAR 0091-3286 Google Scholar
Photorefractive Materials and their Applications. I, Springer-Verlag, Berlin
(1987). Google Scholar
Photorefractive Materials and their Applications. II, Springer-Verlag, Berlin
(1989). Google Scholar
N. Katyalet al.,
“The influence of pump beam polarization on the signal beam, in TWM in photorefractive crystals: transmission geometry,”
Optik, 122
(3), 207
–210
(2011). http://dx.doi.org/10.1016/j.ijleo.2009.11.028 OTIKAJ 0030-4026 Google Scholar
N. Kukhtarevet al.,
“Laser beam shaping by holographic optical elements,”
Proc. SPIE, 5257 152
–162
(2003). http://dx.doi.org/10.1117/12.545847 PSISDG 0277-786X Google Scholar
N. V. Kukhrarevet al.,
“Holographic storage in electro-optic crystals. II. beam coupling-light amplification,”
Ferroelectrics, 22
(1), 961
–964
(1978). http://dx.doi.org/10.1080/00150197908239451 FEROA8 0015-0193 Google Scholar
J. P. HuignardA. Marrakchi,
“Coherent signal beam amplification in two-wave mixing experiments with photorefractive crystals,”
Opt. Commun., 38
(4), 249
–254
(1981). http://dx.doi.org/10.1016/0030-4018(81)90392-8 OPCOB8 0030-4018 Google Scholar
Y. FainmanE. KlancnikS. H. Lee,
“Optimal coherent image amplification by two-wave coupling in photorefractive ,”
Opt. Eng., 25
(2), 228
–234
(1986). http://dx.doi.org/10.1117/12.7973810 OPEGAR 0091-3286 Google Scholar
P. TayebatiD. Mahgerefteh,
“Theory of the photorefractive effect for and with shallow traps,”
J. Opt. Soc. Am. B, 8
(5), 1053
–1064
(1991). http://dx.doi.org/10.1364/JOSAB.8.001053 JOBPDE 0740-3224 Google Scholar
M. H. Garettet al.,
“High beam-coupling gain and deep- and shallow-trap effects in cobalt-doped barium titanate,”
J. Opt. Soc. Am. B, 9
(8), 1407
–1415
(1992). http://dx.doi.org/10.1364/JOSAB.9.001407 JOBPDE 0740-3224 Google Scholar
K. Buseet al.,
“Photorefractive properties of tetragonal crystals and explanation by the three-valence charge-transport model,”
J. Opt. Soc. Am. B, 13
(11), 2644
–2651
(1996). http://dx.doi.org/10.1364/JOSAB.13.002644 JOBPDE 0740-3224 Google Scholar
B. Penget al.,
“Two-wave mixing of ion-implanted photorefractive waveguides in near-stoichiometric crystals,”
Opt. Mater., 33
(6), 773
–776
(2011). http://dx.doi.org/10.1016/j.optmat.2010.12.016 OMATET 0925-3467 Google Scholar
J. Chen, Principle and Application of Laser, Publishing House of Electronics Industry, Beijing
(2004). Google Scholar
P. G. Kazanskyet al.,
“Pockels effect in thermally poled silica optical fibers,”
Electron. Lett., 31
(1), 62
–65
(1995). http://dx.doi.org/10.1049/el:19950036 ELLEAK 0013-5194 Google Scholar
D. Z. AndersonR. Saxena,
“Theory of multimode operation of a unidirectional ring oscillator having photorefractive gain: weak field limit,”
J. Opt. Soc. Am. B., 4
(2), 164
–176
(1987). http://dx.doi.org/10.1364/JOSAB.4.000164 JOBPDE 0740-3224 Google Scholar
M. ChiJ. P. HuignardP. M. Petersen,
“A general theory of two-wave mixing in nonlinear media,”
J. Opt. Soc. Am. B, 26
(8), 1578
–1584
(2009). http://dx.doi.org/10.1364/JOSAB.26.001578 JOBPDE 0740-3224 Google Scholar
N. Koukourakiset al.,
“Photorefractive two-wave mixing for image amplification in digital holography,”
Opt. Express, 19
(22), 22004
–22023
(2011). http://dx.doi.org/10.1364/OE.19.022004 OPEXFF 1094-4087 Google Scholar
Z. Y. LiB. Y. GuG. Z. Yang,
“Slowly varying amplitude approximation appraised by transfer-matrix approach,”
Phys. Rev. B, 60
(15), 10644
–10647
(1999). http://dx.doi.org/10.1103/PhysRevB.60.10644 PRBMDO 1098-0121 Google Scholar
C. Yanget al., Photorefractive Nonliear Optical Materials:LiNbO3 Crystals, Science Press, Beijing
(2009). Google Scholar
T. K. YadavM. K. MauryaR. A. Yadav,
“Effect of photoconductivity and oscillation frequency shift on the signal beam intensity in two beam coupling in photorefractive materials,”
Optik, 122
(18), 1607
–1614
(2011). http://dx.doi.org/10.1016/j.ijleo.2010.10.011 OTIKAJ 0030-4026 Google Scholar
BiographyDanyu Chen received her BS degree in communication engineering from PLA Information Engineering University, Zhengzhou, China, in 2011. She is now a MS candidate in the College of Communication Engineering, Chongqing University, Chongqing, China. Her research focuses on optical information processing and image processing systems. Fengchun Tian received his BE, ME, and DE degrees in radio engineering, biomedical instruments and engineering, theoretical electric engineering from Chongqing University, Chongqing, P.R. China, in 1984, 1986, and 1996, respectively. Since 1984, he has been working in Chongqing University as a teacher. Since 2007, he is also an adjunct professor in the University of Guelph, Canada. His current research interests are image processing (including optical image and video processing), biomedical and modern signal processing technology. Youwen Hu received his BS degree in telecommunication engineering from Chongqing University, Chongqing, China, in 2004. He received his MS degree in signal and information processing in Chongqing University, Chongqing, China. He is now a PhD candidate at Chongqing University. His research interests are adaptive optics, optical negative feedback system, and optical information processing. Ying Liu received her BS degree in communication engineering from Qingdao University of Science and Technology, Qingdao, China, in 2012. Now, she is a graduate student in Chongqing University, Chongqing, China. Her research focuses on the optical image processing. Liang Han received her MS degree in communication and information systems in 2004 and her PhD degree in circuits and systems in 2008 from the College of Communication Engineering, Chongqing University, Chongqing, China. In 1997, he joined Chongqing University, where he is currently an associate professor in the College of Communication Engineering. His research focuses on imaging processing, especially in multiresolution analysis and its optical implementation method, image fusion, image denoising, and weak signal extraction methods. |