OSA's Digital Library

Journal of the Optical Society of America B

Journal of the Optical Society of America B

| OPTICAL PHYSICS

  • Editor: Henry van Driel
  • Vol. 28, Iss. 12 — Dec. 1, 2011
  • pp: 2947–2955
« Show journal navigation

Nonlinear effects in random lasers

Jonathan Andreasen, Patrick Sebbah, and Christian Vanneste  »View Author Affiliations


JOSA B, Vol. 28, Issue 12, pp. 2947-2955 (2011)
http://dx.doi.org/10.1364/JOSAB.28.002947


View Full Text Article

Acrobat PDF (847 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

Recent numerical and theoretical studies have demonstrated that the modes at threshold of a random laser are in direct correspondence with the resonances of the same system without gain, a feature which is well known in conventional lasers but not known until recently for random lasers. This paper presents numerical results of the multimode regime that takes place when the pumping rate is progressively increased above threshold. Behavior that is already known in standard lasers, such as mode competition and nonlinear wave mixing, are shown to also take place in random lasers thus reinforcing their recent modal description. However, due to the complexity of the laser modes and to the openness of such lasers, which require large external pumping to compensate for strong loss, one observes that these effects are systematic and can be more pronounced than in a conventional laser.

© 2011 Optical Society of America

1. INTRODUCTION

Since their prediction by Lethokov [1

V. S. Letokhov, “Generation of light by a scattering medium with negative resonance absorption,” Sov. Phys. JETP 26, 835–840 (1968).

], random lasers have been the subject of numerous studies. Unlike conventional lasers they have no cavity, such as a Fabry–Pérot resonator or ring cavity. Instead, they are made of a scattering random medium like a semiconductor powder or a suspension of scattering particles in a laser dye and are excited by an optical [2

H. Cao, “Lasing in random media,” Waves Random Media 13, R1–R39 (2003) and references therein. [CrossRef]

] or electrical [3

Y. Li, X. Ma, M. Xu, L. Xiang, and D. Yang, “Remarkable decrease in threshold for electrically pumped random ultraviolet lasing from ZnO fim by incorporation of Zn 2 TiO 4 nanoparticles,” Opt. Express 19, 8662–8669 (2011). [CrossRef] [PubMed]

] pump which introduces gain. Multiple scattering of light in the random medium provides optical feedback. An important feature of random lasers is that they are open systems, typically with strong leakage. Because of these unusual features, the nature of lasing modes of a random laser is at the center of a several years’ debate. Random lasing has been described as light diffusion with gain [4

S. John and G. Pang, “Theory of lasing in a multiple-scattering medium,” Phys. Rev. A 54, 3642–3652 (1996). [CrossRef] [PubMed]

, 5

D. S. Wiersma and A. Lagendijk, “Light diffusion with gain and random lasers,” Phys. Rev. E 54, 4256–4265 (1996). [CrossRef]

], in terms of well-localized modes inside the scattering medium [6

H. Cao, J. Y. Xu, S.-H. Chang, and S. T. Ho, “Transition from amplified spontaneous emission to laser action in strongly scattering media,” Phys. Rev. E 61, 1985–1989 (2000). [CrossRef]

, 7

X. Jiang and C. M. Soukoulis, “Time dependent theory for random lasers,” Phys. Rev. Lett. 85, 70–73 (2000). [CrossRef] [PubMed]

, 8

C. Vanneste and P. Sebbah, “Selective excitation of localized modes in active random media,” Phys. Rev. Lett. 87, 183903 (2001). [CrossRef]

, 9

R. C. Polson, M. E. Raikh, and Z. V. Vardeny, “Universality in unintentional laser resonators in π-conjugated polymer films,” C. R. Acad. Sci. Ser. IV A, 509–521 (2002).

, 10

V. M. Apalkov, M. E. Raikh, and B. Shapiro, “Random resonators and prelocalized modes in disordered dielectric films,” Phys. Rev. Lett. 89, 016802 (2002). [CrossRef] [PubMed]

, 11

A. Tulek, R. C. Polson, and Z. V. Vardeny, “Naturally occurring resonators in random lasing of π-conjugated polymer films,” Nat. Phys. 6, 303–310 (2010). [CrossRef]

] and by random walks of photons along exceptionally long paths [12

S. Mujumdar, M. Ricci, R. Torre, and D. S. Wiersma, “Amplified extended modes in random lasers,” Phys. Rev. Lett. 93, 053903 (2004). [CrossRef] [PubMed]

]. Such approaches suffer from different drawbacks. On the one hand, the diffusion equation approach, which is successful in describing random lasers with nonresonant feedback [2

H. Cao, “Lasing in random media,” Waves Random Media 13, R1–R39 (2003) and references therein. [CrossRef]

], and the random walk of photons do not take into account the wave aspect of light propagation, thus neglecting the part played by interference. On the other hand, experimental characterization of the scattering media which exhibited random lasing in the presence of gain showed that they were far from the localization regime, thus making the occurrence of well-localized modes unlikely.

Recently, important progress has been made in theoretical and numerical studies by recognizing that even the bad resonances of leaky systems play a role in random lasing, similar to the part played by cavity modes of a conventional laser [13

C. Vanneste, P. Sebbah, and H. Cao, “Lasing with resonant feedback in weakly scattering random systems,” Phys. Rev. Lett. 98, 143902 (2007). [CrossRef] [PubMed]

, 14

H. E. Türeci, L. Ge, S. Rotter, and A. D. Stone, “Strong interactions in multimode random lasers,” Science 320, 643–646 (2008). [CrossRef] [PubMed]

, 15

H. E. Türeci, A. D. Stone, L. Ge, S. Rotter, and R. J. Tandy, “ Ab initio self-consistent laser theory and random lasers,” Nonlinearity 22, C1–C18 (2009). [CrossRef]

, 16

O. Zaitsev, L. Deych, and V. Shuvayev, “Statistical properties of one-dimensional random lasers,” Phys. Rev. Lett. 102, 043906 (2009). [CrossRef] [PubMed]

, 17

O. Zaitsev and L. Deych, “Recent developments in the theory of multimode random lasers,” J. Opt. 12, 024001 (2010). [CrossRef]

, 18

O. Zaitsev and L. Deych, “Diagrammatic semiclassical laser theory,” Phys. Rev. A 81, 023822 (2010). [CrossRef]

]. In particular, in such a modal description of laser action, it was shown theoretically and numerically [19

J. Andreasen, A. Asatryan, L. Botten, M. Byrne, H. Cao, L. Ge, L. Labonté, P. Sebbah, A. D. Stone, H. E. Türeci, and C. Vanneste, “Modes of random lasers,” Adv. Opt. Photon. 3, 88–127 (2011). [CrossRef]

] that the first lasing mode at threshold is not strictly identical but very close to a system resonance (or equivalently, to a decaying quasimode, which generalizes the concept of modes to leaky systems [20

S. M. Dutra and G. Nienhuis, “Quantized modes of a leaky cavity,” Phys. Rev. B 62, 063805 (2000).

]). Despite the complex nature of these resonances, the understanding of random lasing is therefore greatly simplified in the single-mode lasing regime, just above threshold. As the pumping rate is increased, however, nonlinear effects come into play. These effects might be even more pronounced in weakly scattering active media where high external pumping is required to compensate for strong loss. This is particularly true in the multimode regime above the lasing threshold, with the onset of competition between different lasing modes. Detailed investigations of the multimode regime in random lasers are scarce [21

H. Cao, X. Jiang, Y. Ling, J. Y. Xu, and C. M. Soukoulis, “Mode repulsion and mode coupling in random lasers,” Phys. Rev. B 67, 161101(R) (2003). [CrossRef]

, 22

X. Jiang, S. Feng, C. M. Soukoulis, J. Zi, J. D. Joannopoulos, and H. Cao, “Coupling, competition, and stability of modes in random lasers,” Phys. Rev. B 69, 104202 (2004). [CrossRef]

, 23

C. Conti, M. Leonetti, A. Fratalocchi, L. Angelani, and G. Ruocco, “Condensation in disordered lasers: Theory, 3d+1 simulations, and experiments,” Phys. Rev. Lett. 101, 143901 (2008). [CrossRef] [PubMed]

]. The background of standard laser physics is not directly transposable to random lasers, where the complexity of the spatial and spectral properties of random lasers need to be taken into account. The recent steady-state ab initio laser theory (SALT) [24

H. E. Türeci, A. D. Stone, and B. Collier, “Self-consistent multimode lasing theory for complex or random lasing media,” Phys. Rev. A 74, 043822 (2006). [CrossRef]

, 25

L. Ge, R. J. Tandy, A. D. Stone, and H. E. Türeci, “Quantitative verification of ab initio self-consistent laser theory,” Opt. Express 16, 16895–16902 (2008). [CrossRef] [PubMed]

], which considers the openness of novel laser systems, such as random lasers, is the first theory to give predictions about nonlinear phenomena and mode interaction [14

H. E. Türeci, L. Ge, S. Rotter, and A. D. Stone, “Strong interactions in multimode random lasers,” Science 320, 643–646 (2008). [CrossRef] [PubMed]

, 15

H. E. Türeci, A. D. Stone, L. Ge, S. Rotter, and R. J. Tandy, “ Ab initio self-consistent laser theory and random lasers,” Nonlinearity 22, C1–C18 (2009). [CrossRef]

]. This theory relies, however, on several assumptions. In particular, it assumes the existence of a steady-state multiperiodic solution of the laser field, population inversion, and polarization of the atomic medium. However, the full nonlinear dynamics certainly play a role in determining behavior in such complex and highly nonlinear media. Including them would not only bring a more complete theoretical approach, but also bring theory closer to experiment.

In this paper, we numerically investigate the full dynamics of one-dimensional (1D) and two-dimensional (2D) random lasers using steady external pumping by progressively increasing the external pump intensity. We report several manifestations of optical nonlinearities, both in the single-mode and multimode regime, in the steady-state as well as in the transient regime. Strong relaxation oscillations, mode competition, and mode suppression are reported in the transient regime, which reveal the complexity of the laser dynamics in random lasers. Third-harmonic generation, four-wave mixing and sum-frequency generation are observed, which have never been reported before. When progressively increasing the pumping rate above threshold for lasing, the so-called “second threshold” [26

E. Roldán, G. J. de Valcárcel, F. Prati, F. Mitschke, and T. Voigt, “Multilongitudinal mode emission in ring cavity class B lasers,” in “Trends in Spatiotemporal Dynamics in Lasers. Instabilities, Polarization Dynamics, and Spatial Structures ,” O. Gomez-Calderon and J. M. Guerra, eds. (Research Signpost, 2005), pp. 1–80.

] is eventually reached, where the steady-state becomes unstable. Note that this “second threshold” is not the threshold of the second lasing mode which will be discussed later. A coherent instability manifests itself as temporal oscillations of the field intensity, atomic population inversion, and medium polarization. This issue has been addressed in a different paper [27

J. Andreasen, P. Sebbah, and C. Vanneste, “Coherent instabilities in random lasers,” Phys. Rev. A 84, 023826 (2011). [CrossRef]

] and will not be discussed here, although it is another manifestation of the nonlinear dynamics expected in random lasers. If these effects are not new, their observation in random lasers is interesting since they not only challenge the theoretical understanding of these systems but some of them are more pronounced than in conventional lasers. Our effort here is to stress how these observations are related to the particular nature of the modes of random lasers.

The paper is organized as follows. Section 2 describes the 1D and 2D random structures that have been studied and the numerical methods that have been used. Section 3 shows nonlinear effects like the buildup of laser oscillation and interactions between lasing modes, which include mode competition and nonlinear wave mixing. Conclusions are given in Section 4.

2. NUMERICAL APPROACHES

2A. Random Structures

The 1D random structures we consider are similar to those studied in [28

J. Andreasen and H. Cao, “Creation of new lasing modes with spatially nonuniform gain,” Opt. Lett. 34, 3586–3588 (2009). [CrossRef] [PubMed]

, 29

J. Andreasen and H. Cao, “Numerical study of amplified spontaneous emission and lasing in random media,” Phys. Rev. A 82, 063835 (2010). [CrossRef]

]. They are composed of 41 layers and shown in Fig. 1a. Dielectric material with optical index n1=1.25 separated by air gaps ( n2=1) results in a spatially modulated index n(x). Outside the random medium, the index is 1. The system is randomized by specifying thicknesses for each layer as d 1,2= d 1,2(1+ηζ), where d1=100nm and d2=200nm are the average thicknesses of the layers, η=0.9 represents the degree of randomness, and ζ is a random number in (1,1). The length of the random structure L is normalized to L=6.1μm. These parameters give a localization length ξ11μm within the wavelength range considered here. Here, the localization length ξ is calculated from the dependence of ensemble-averaged transmittance T on the system lengths L as ξ 1=dlnT/dL and averaged over the wavelength range of interest.

The 2D random structures considered [Fig. 1b] are made of circular dielectric particles with radius r=60nm, optical index n1=1.25, and surface filling fraction Φ=40%, which are randomly distributed in a background medium of size L2=5×5 μm 2 and index n2=1. The optical index of the domain that surrounds the random medium is n=1. The scattering mean free path s2μm and the localization length ξ12μm. Such systems, which are similar to those considered in [13

C. Vanneste, P. Sebbah, and H. Cao, “Lasing with resonant feedback in weakly scattering random systems,” Phys. Rev. Lett. 98, 143902 (2007). [CrossRef] [PubMed]

, 30

P. Sebbah and C. Vanneste, “Random laser in the localized regime,” Phys. Rev. B 66, 144202 (2002). [CrossRef]

], are in a weakly scattering regime close to the diffusive regime. We also consider a second set of smaller 2D systems with a size L2=1×1 μm 2 having a density of states comparable to the 1D systems. The optical index of the particles is chosen as n1=1.5. In this case, the scattering mean free path s0.3μm and the localization length ξ5μm. Such small 2D systems are in a quasi-ballistic regime.

2B. Methods to Study Lasing

In 1D and 2D systems the background medium of index 1 has been chosen as the active part of the system. Two different methods have been used to study lasing in 1D systems. First, it is possible to employ the transfer matrix (TM) method similar to that used in [31

J. Andreasen and H. Cao, “Spectral behavior of partially pumped weakly scattering random lasers,” Opt. Express 19, 3418–3433 (2011). [CrossRef] [PubMed]

]. However, this is a frequency-domain method, which cannot describe transient behavior in the time-domain and lasing mode competition in the multimode regime. For this reason, we also solve Maxwell’s equations using the finite-difference time-domain (FDTD) method [32

A. Taflove and S. Hagness, Computational Electrodynamics (Artech House, 2005), 3rd ed.

]. The two methods can be compared easily [28

J. Andreasen and H. Cao, “Creation of new lasing modes with spatially nonuniform gain,” Opt. Lett. 34, 3586–3588 (2009). [CrossRef] [PubMed]

], yet bring different advantages. In this case, the gain medium in FDTD is modeled by a four-level atomic system [7

X. Jiang and C. M. Soukoulis, “Time dependent theory for random lasers,” Phys. Rev. Lett. 85, 70–73 (2000). [CrossRef] [PubMed]

, 30

P. Sebbah and C. Vanneste, “Random laser in the localized regime,” Phys. Rev. B 66, 144202 (2002). [CrossRef]

, 33

A. S. Nagra and R. A. York, “FDTD analysis of wave propagation in nonlinear absorbing and gain media,” IEEE Trans. Antennas Propag. 46, 334–340 (1998). [CrossRef]

].

In 2D systems, we only use the FDTD method to solve Maxwell’s equations and consider the case of transverse magnetic polarization. We use perfectly matched layer (CPML) absorbing boundaries [32

A. Taflove and S. Hagness, Computational Electrodynamics (Artech House, 2005), 3rd ed.

] to approximate open boundary conditions. Gain in the system is again introduced by coupling Maxwell’s equations with the population equations of a four- level atomic system. The corresponding equations are given in Appendix A. Population inversion between the levels corresponding to the laser transition is created by an external pump which transfers atoms from the ground state (level 1) to the upper level (level 4) of the four-level atomic system. With the parameters used in our simulations, the gain curve is centered at λa=446.9nm and has a spectral width δ λa=11nm (see Appendix A). In the calculations presented in the paper, the control parameter is the pumping rate Pr of atoms from the ground state to the upper level. Note that spontaneous emission is not taken into account in our calculations.

3. RESULTS

Our results here are devoted to random laser modes above threshold and to their interactions in the multimode regime. We not only discuss the final stationary state but also the transient buildup of the laser field. In both situations, we observe mode competition. The large gain which must be introduced in order to counterbalance the large loss of the system leads to a significant evolution of the lasing modes as a function of the pumping rate. This effect is much larger than its counterpart in conventional lasers. Finally, mode interaction reveals itself by the systematic observation of nonlinear wave mixing, which appears as soon as multiple lasing modes coexist.

3A. Laser Field Buildup

Each FDTD calculation of a given random system with a fixed value of the pumping rate Pr starts in the same way. First, the initial atomic populations are set at the stationary values in absence of stimulated emission, i.e., for ( Ez/ ωa)dP/dt=0 (Appendix A). At the initial time t=0, a small-intensity broadband light pulse is launched in the system. Since there is no spontaneous emission in our model, this short pulse, which propagates and is scattered by the particles of the system, provides the nonzero initial field necessary for lasing action to begin. If Pr is larger than the threshold value, the laser field builds up in the system until it reaches a value where the population inversion starts to be depleted by increasing stimulated emission. The steady-state regime is reached when the value of the population inversion density corresponds to the gain, which exactly compensates the loss through the open boundaries of the system. This is the standard behavior of any laser oscillator.

However, random lasers are open systems usually with short “cavity” lifetimes, i.e., shorter than the atomic lifetime T1. This is the well-known condition [34

A. E. Siegman, Lasers (University Science Books, 1986).

] for the occurrence of relaxation oscillations or spiking, i.e., temporal oscillations of the field and of the population inversion before the laser oscillation converges to a steady-state. An example is given in Fig. 2 where relaxation oscillations are seen in the population inversion and laser intensity for a 1D random laser. As the pumping rate increases, the intensity and the frequency of the relaxation oscillations increase. Indeed, relaxation oscillations in random lasers have been reported by several groups [35

C. M. Soukoulis, X. Jiang, J. Y. Xu, and H. Cao, “Dynamic response and relaxation oscillations in random lasers,” Phys. Rev. B 65, 041103(R) (2002). [CrossRef]

, 36

M. A. Noginov, G. Zhu, A. A. Frantz, J. Novak, S. N. Williams, and I. Fowlkes, “Dependence of Nd Sc 3( BO 3 ) 4 random laser parameters on particle size,” J. Opt. Soc. Am. B 21, 191–200 (2004) and references cited therein. [CrossRef]

, 37

K. L. van der Molen, A. P. Mosk, and A. Lagendijk, “Relaxation oscillations in long-pulsed random lasers,” Phys. Rev. A 80, 055803 (2009). [CrossRef]

] and we have systematically observed them in our simulations thus stressing the importance of the openness of such systems.

3B. Mode Alterations Near Threshold

As mentioned earlier, even the bad resonances of open random systems play a role in lasing similar to the part played by the cavity modes of a conventional laser. In particular, the first lasing mode at threshold is not strictly identical but very close to a resonance of the passive system, a result which was established only recently for random lasers [19

J. Andreasen, A. Asatryan, L. Botten, M. Byrne, H. Cao, L. Ge, L. Labonté, P. Sebbah, A. D. Stone, H. E. Türeci, and C. Vanneste, “Modes of random lasers,” Adv. Opt. Photon. 3, 88–127 (2011). [CrossRef]

]. However, this conclusion relies on several conditions that are not always fulfilled in real experiments, for instance, uniform pumping of the total system. In real random lasers, pumping takes place in the active elements, which usually do not fill up the total volume, such as particles suspended in laser dye where there is no gain in the particles. This situation contrasts strongly with that of conventional lasers, where pumping is usually distributed uniformly in the gain domain.

Three modes of a random laser near threshold are considered in Fig. 3. One mode is the first lasing mode at threshold and the two other modes are second lasing modes above threshold for two different realizations of the disorder. In order to obtain a quantitative measure of the similarity of two intensity distributions Ei2(x) and Ej2(x), we introduce the difference
D ij= 0L Ei2(x) Ej2(x)dx,
(1)
where the distributions E i,j2(x) are normalized to one. The first column in Fig. 3a compares the intensity distributions of the lasing modes with uniform gain and those with partial gain, i.e., gain located in the gaps between the high index particles [see Fig. 1a]. Noticeable differences are clearly visible, up to the value D=0.16 for the first lasing mode. Previous investigation of partial pumping in random lasers [38

J. Andreasen, C. Vanneste, L. Ge, and H. Cao, “Effects of spatially nonuniform gain on lasing modes in weakly scattering random systems,” Phys. Rev. A 81, 043818 (2010). [CrossRef]

] has already demonstrated that lasing modes can become very different from resonances of the passive system. Here, we are in an intermediate situation where gain is not placed everywhere but homogeneously distributed across the system. The impact of partial gain varies from one lasing mode to another lasing mode, as shown by the reduced D values for the second lasing modes. This illustrates the role of randomness in these lasers and the effects of partial pumping that arises in such multiple scattering systems.

Another approximation used in several theoretical and numerical studies is a flat gain curve, i.e., gain that does not depend on frequency. The second column in Fig. 3a compares the intensity distributions of the first lasing mode with uniform and flat gain and that with partial and frequency- dependent gain. The frequency dependence of the gain does not contribute significantly to the modification of the spatial distribution of the first lasing mode, which is close to the atomic transition wavelength, but does impact the second lasing modes. The threshold for the second lasing mode is larger, which results in a larger modification of the medium susceptibility due to the gain and further modifies the lasing mode.

The large alteration of the medium susceptibility due to the intense pumping that is required to achieve lasing is again a particular feature of random lasers due to their openness and associated strong leakage. For comparison, we consider a Fabry–Pérot laser with mirrors of reflectivity 80%. The length between the mirrors is L=6.1μm, the same as the random lasers. The effective index of the random lasers, n eff=1.1, is used between the mirrors to obtain similar mode spacing. The corresponding differences for this Fabry–Pérot laser are D= 10 4 and 10 3 for the first and second lasing modes, respectively, compared to the smallest value D= 10 1 observed in the random laser. This illustrates that, although similar effects may be observed in random and conventional lasers, the inhomogeneity and the openness of random lasers make these effects more pronounced, causing more dramatic changes to steady-state laser behavior.

Finally, different from the first lasing mode at threshold, lasing modes in the multimode regime compete for gain with the occurrence of saturation effects. Here, we take advantage of the TM and FDTD methods described in Subsection 2B. The third column in Fig. 3a compares the intensity distributions of the first lasing mode with uniform and flat gain and that with gain saturation computed via the FDTD method. As expected, saturation effects at the threshold of the first lasing mode are negligible, but strongly influence the second lasing mode. Figure 3b compares the spatial profiles of the second lasing mode at threshold for these two situations. The openness of the structure is demonstrated by the fact that the lasing mode energy is peaked at the boundary of the system [38

J. Andreasen, C. Vanneste, L. Ge, and H. Cao, “Effects of spatially nonuniform gain on lasing modes in weakly scattering random systems,” Phys. Rev. A 81, 043818 (2010). [CrossRef]

, 39

X. Wu, W. Fang, A. Yamilov, A. A. Chabanov, A. A. Asatryan, L. C. Botten, and H. Cao, “Random lasing in weakly scattering systems,” Phys. Rev. A 74, 053812 (2006). [CrossRef]

, 40

C. Vanneste and P. Sebbah, “Complexity of two-dimensional quasimodes at the transition from weak scattering to Anderson localization,” Phys. Rev. A 79, 041802(R) (2009). [CrossRef]

]. For comparison, the corresponding differences for the Fabry–Pérot laser are D= 10 2 and 3· 10 2 for the first and second lasing modes, respectively. Once again such effects in random lasers are more pronounced compared to their conventional laser counterpart.

In addition to the example case above, we have performed a systematic investigation of lasing thresholds over ten different realizations of disorder in 1D random lasers. We compare results from the TM and FDTD methods, where the essential difference is that mode competition effects are included in FDTD. Figure 4a displays the lasing thresholds of the first and second lasing modes of the 10 random lasers. The thresholds for the first lasing mode are identical, whatever the method. This result is expected since, at the threshold of the first lasing mode, there is neither competition with other lasing modes nor saturation effects. To the contrary, there is a noticeable difference between the two methods for the thresholds of the second lasing modes. This indicates that competition effects are important. The threshold value calculated via the FDTD method is always larger than the value calculated via the TM method. This demonstrates that the first lasing mode delays the onset of the second lasing mode, an effect that the TM method does not incorporate. One consistently notices that the larger the threshold difference calculated via the TM method between lasing mode 1 and 2 [filled and open circles in Fig. 4a], the larger the difference for lasing mode 2 between the threshold calculated via the FDTD method and that calculated via the TM method [open triangles and circles in Fig. 4a]. This effect is clearly illustrated in Fig. 4b, which displays the ratio of the second lasing threshold over the first lasing threshold without (TM method) and with (FDTD method) mode interaction included. The larger the first ratio, even larger is the second ratio. Realizations 7 and 10 exhibit the largest threshold ratios via the TM method. The second lasing thresholds were never reached with the FDTD method, even when the pumping rate was extended almost two orders of magnitude past the second lasing thresholds of other realizations. In these systems, the single-mode regime is very robust. The resonances underlying the lasing mode in realizations 7 and 10 are examined in Appendix B.

3C. Mode Evolution Above Threshold

Further above the lasing threshold, mode competition results in an irregular increase of the intensities of the lasing modes as a function of the pumping rate and, possibly, in mode suppression. We confirm these predictions here. It is important to point out that in all cases, increasing the pumping rate eventually leads all the random lasers we have studied to reach the threshold of a coherent instability (called the “second threshold”) thus limiting the range of pumping rates for observing the stationary states we describe below. This is not surprising since such instabilities are known to occur in very leaky lasers [26

E. Roldán, G. J. de Valcárcel, F. Prati, F. Mitschke, and T. Voigt, “Multilongitudinal mode emission in ring cavity class B lasers,” in “Trends in Spatiotemporal Dynamics in Lasers. Instabilities, Polarization Dynamics, and Spatial Structures ,” O. Gomez-Calderon and J. M. Guerra, eds. (Research Signpost, 2005), pp. 1–80.

]. Details concerning such instabilities in random lasers have been discussed [27

J. Andreasen, P. Sebbah, and C. Vanneste, “Coherent instabilities in random lasers,” Phys. Rev. A 84, 023826 (2011). [CrossRef]

], but their effect on lasing mode interactions are beyond the scope of this paper. Thus, we focus now on laser behavior between the first threshold and this so-called second threshold.

Figure 5 shows an interesting case of a bimodal 1D random laser emission spectrum evolution as the pump rate is progressively increased. As the second lasing mode turns on at Pr=0.30 ns 1, the progression of the first mode is stopped and its intensity progressively reduces until it turns off at Pr=0.44 ns 1 [Fig. 5a]. Figure 5b displays the wavelengths of the two lasing modes as a function of the pumping rate. Frequency drift of lasing mode 2 toward the central frequency of the gain curve is a manifestation of frequency pulling and allows the lasing mode to experience more amplification, which certainly helps lasing mode 2 to overcome lasing mode 1 [41

B. Liu, A. Yamilov, Y. Ling, J. Y. Xu, and H. Cao, “Dynamic nonlinear effect on lasing in a random medium,” Phys. Rev. Lett. 91, 063903 (2003). The surprising drift of mode 1 across the maximum of the gain curve, toward mode 2 [Fig. 5b], is attributed to the nonlinear Kerr effect. [CrossRef] [PubMed]

]. However, the main reason for the extinction of the first lasing mode is self-saturation and cross-saturation effects between the two lasing modes in agreement with the recent predictions of SALT [14

H. E. Türeci, L. Ge, S. Rotter, and A. D. Stone, “Strong interactions in multimode random lasers,” Science 320, 643–646 (2008). [CrossRef] [PubMed]

].

Figure 6 displays color-coded wavelength power spectra (vertical axis) versus time (horizontal axis) for successive values of the pumping rate. Spectra are obtained by Fourier transform over a sliding Welch time window of 3.3ps and a time delay of 0.3ps. Relaxation oscillations are readily seen in all spectrograms at early times as vertical wavelets, which decay rapidly over time. Mode competition manifests itself already in Fig. 6b which is just below the threshold of lasing mode 2. Lasing mode 2 shows up for a very short duration before the first mode starts oscillating. One observes the same phenomenon in Figs. 6c, 6d but now the second lasing mode comes back after having been suppressed by the onset of the first lasing mode. The situation is finally reversed in Fig. 6f where only lasing mode 2 survives, after a brief appearance of lasing mode 1. It is interesting to note that the spectral drift of the lasing lines described in Fig. 5b occurs in the time-domain [Fig. 6f]. This is a signature of the dynamics of the cross-saturation between lasing modes.

We consider now the evolution of the spatial distribution of the lasing modes as a function of the pumping rate using the FDTD method. Figure 7 displays the difference D ii( Pr), as defined by Eq. (1), between one lasing mode at threshold ( Pr= Pr th) and the same lasing mode at pumping rate Pr. The difference increases monotonically with Pr, though not at a constant rate. In this example, the evolution of lasing mode 1 slows down abruptly at the onset of lasing mode 2. This indicates that self-saturation effects that first dominate the evolution of lasing mode 1 become less efficient at the onset of lasing mode 2. This is consistent with the evolution of the intensities displayed in Fig. 5a, which shows that lasing mode 1 slowly decays when lasing mode 2 appears. Figure 7 also shows that values of the differences become significant when Pr increases. Hence, the lasing mode distributions depart significantly from the distribution of the associated resonances. This shows that in open systems like random lasers, the conventional description of lasing modes in terms of resonances of the passive cavity is limited to values of the pumping rate close to lasing threshold and is not valid for larger pumping rates, a result in agreement with the recent predictions of SALT.

The previous discussion concerned 1D systems. We have also examined 10 different 2D random lasers of size 5×5 μm 2 and 10 different 2D random lasers of size 1×1 μm 2. For the smaller size, the density of states is comparable to the 1D systems so that the distinction of the main lasing peaks from other peaks like those due to four-wave mixing was easier than in large 2D systems. Figure 8 shows a representative sample of the evolution of lasing modes for one of those 2D lasers, size 1×1 μm 2. We observe irregular increase of the lasing mode intensities with increasing pump rate, which again is the signature of mode self-saturation and cross-saturation. The evolution of lasing mode intensities as a function of Pr is quite similar in the more realistic 2D random lasers of larger size. The difference is that more modes reach the lasing threshold in the larger 2D system, but similar irregular increase of the mode intensities is observed as in Fig. 8.

3D. Nonlinear Wave Mixing

Until now, we have described nonlinear effects involving gain saturation and gain competition, as well as nonlinear refraction. In our numerical simulations, in both the single- and multimode regimes, we also observe the effect of χ (3) nonlinearities, i.e., frequency-mixing processes including third- harmonic generation and four-wave mixing (because noise is not incorporated into our simulations, we did not observe such mixing below the threshold for lasing due to the lack of large fields in the system).

First, we observe that each lasing mode of wavelength λ has a peak at the wavelength λ/3 associated with it (Fig. 9). We found that third-harmonic generation is a consistent process in random lasers, which occurs whatever the value of the pumping rate at or above threshold.

Next, at the onset of a second lasing mode, we observe new peaks in the emission spectrum. Examples for 1D and 2D systems are displayed in Figs. 10a, 10c. On top of the two lasing modes with wavelengths λ1 and λ2, narrow lateral peaks appear at λ3 and λ4. We attribute the origin of these new lines to four-wave mixing. This can be checked as follows:
1/ λ31/ λ2=1/ λ21/ λ1,
(2a)
1/ λ41/ λ1=1/ λ11/ λ2.
(2b)
(2)
This signature of four-wave mixing has been verified for all couples of neighboring lasing modes in 1D and 2D lasers. The four-wave mixing peaks in Fig. 10 are very weak. We chose the spectra in Fig. 10 in order to show that spectra in 1D and 2D could be very similar. However, we found that the efficiency of four-wave mixing is strongly dependent on the realization of the disorder. Some of them have efficiencies up to four orders of magnitude below the intensities of the main peaks instead of 10 to 12 orders of magnitude as in Fig. 10. We cannot exclude that some realizations could give better efficiencies. This is a statistical investigation to be done in the future.

Normally, four-wave mixing requires not only energy conservation, such as Eqs. (2a, 2b), but is also improved by phase-matching conditions. In standard optical homogeneous systems, phase-matching is discussed in terms of wave vectors
k1 + k2 k3 + k4 .
(3)
Obviously, one cannot associate a unique wavevector to the lasing modes of a random system. Hence for a random laser, Eq. (3) is meaningless. Nevertheless, the systematic observation of four-wave mixing shows that randomness in random lasers naturally provides random quasi-phase-matching conditions.

Associated with four-wave mixing, we observe again third-harmonic generation and sum-frequency generation involving the two lasing modes. This is shown in Figs. 10b, 10d. One can check the following equalities:
1/ λ 3t=2/ λ2+1/ λ1,
(4a)
1/ λ 4t=2/ λ1+1/ λ2.
(4b)
(4)

To our knowledge, this is the first report of such nonlinear effects in random lasers. An important result of this section is that nonlinear wave mixing effects are systematic and can be significant depending on the realization of the disorder in such random systems. Similar conclusions for three-wave mixing in random nonlinear materials have been reported in the past [42

R. C. Miller, “Optical harmonic generation in single crystal BaTiO 3 ,” Phys. Rev. 134, A1313–A1319 (1964). [CrossRef]

, 43

C. F. Dewey Jr. and L. O. Hocker, “Enhanced nonlinear optical effects in rotationally twinned crystals,” Appl. Phys. Lett. 26, 442–444 (1975). [CrossRef]

] and more recently [44

E. Y. Morozov and A. S. Chirkin, “Stochastic quasi-phase matching in nonlinear-optical crystals with an irregular domain structure,” Quantum Electron. 34, 227–232 (2004). [CrossRef]

, 45

V. A. Mel’nikov, L. A. Golovan, S. O. Konorov, D. A. Muzychenko, A. B. Fedotov, A. M. Zheltikov, V. Y. Timoshenko, and P. K. Kashkarov, “Second-harmonic generation in strongly scattering porous gallium phosphide,” Appl. Phys. B 79, 225–228 (2004). [CrossRef]

, 46

M. Baudrier-Raybaut, R. Haïdar, P. Kupecek, P. Lemasson, and E. Rosencher, “Random quasi-phase-matching in bulk polycrystalline isotropic nonlinear materials,” Nature 432, 374–376 (2004). [CrossRef] [PubMed]

, 47

S. E. Skipetrov, “Disorder is the new order,” Nature 432, 285–286 (2004). [CrossRef] [PubMed]

]. Such studies reported that phase-matching conditions, which can be difficult to achieve in pure and homogeneous materials, can be met more easily in a random structure although efficiencies are not optimized. Our results concerning random lasers confirm this effect of so-called random quasi-phase-matching [46

M. Baudrier-Raybaut, R. Haïdar, P. Kupecek, P. Lemasson, and E. Rosencher, “Random quasi-phase-matching in bulk polycrystalline isotropic nonlinear materials,” Nature 432, 374–376 (2004). [CrossRef] [PubMed]

] in nonlinear random media.

4. CONCLUSION

In this paper, we have presented numerical results about random lasing when the external pumping rate is progressively increased above threshold. Strong mode competition and significant alteration of the lasing modes due to the large pumping powers that are required to compensate the loss have been observed. Moreover nonlinear wave mixing was shown to take place consistently, where the efficiency depends on the realization of the disorder. Such results confirm and extend the modal description of random lasers that has been recently proposed after the long debate about their nature. Our results show that the complexity and the openness of random lasers lead to the systematic observation of nonlinear effects, some of them being larger than in conventional lasers.

Appendices

APPENDIX A: LASER EQUATIONS

We write below Maxwell’s equations for a 2D system in the case of transverse magnetic polarization and the population equations of the four-level atomic system,
μ0 Hx/t= Ez/y μ0 Hy/t= Ez/x ϵi ϵ0 Ez/t+P/t= Hy/x Hx/y,
(A1)
where ϵi= ni2 i=1,2, ϵ0 is the electric permittivity, and μ0 the magnetic permeability of vacuum. P is the polarization density, which acts as a source in Maxwell’s equations. The time evolution of the four-level atomic system is described by population equations [34

A. E. Siegman, Lasers (University Science Books, 1986).

].
d N1/dt= N2/ τ21 Pr N1 d N2/dt= N3/ τ32 N2/ τ21( Ez/ ωa)dP/dt d N3/dt= N4/ τ43 N3/ τ32+( Ez/ ωa)dP/dt d N4/dt= N4/ τ43+ Pr N1,
(A2)
Ni is the population density in level i, i=1 to 4. The electrons in the ground level 1 are transferred to the upper level 4 by an external pump at a fixed rate Pr. Electrons in level 4 flow downward to level 3 by means of nonradiative decay processes with a characteristic time τ43. This time is very short so that the electrons excited in level 4 quickly populate level 3. The intermediate levels 3 and 2 are the upper and lower levels, respectively, of the laser transition. The decay rate downward from level 3 is τ32. Stimulated transitions due to the electromagnetic field take place between these two levels. Electrons then decay nonradiatively from level 2 to level 1 with a characteristic time τ21. The stimulated transition rate is given as ( Ez/ ωa)dP/dt, where ωa=( E3 E2)/ is the transition frequency between levels 2 and 3. The quantities Ez, P, and Ni depend on the time t but also on the position r in the system. Eventually, the polarization obeys the equation
d2P/d t2+Δ ωadP/dt+ ωa2P=κΔN Ez,
(A3)
ΔN= N2 N3 is the population difference density between the populations in the lower and upper levels of the atomic transition.

Amplification takes place when the external pumping mechanism produces an inverted population difference ΔN<0. The linewidth of the atomic transition is Δ ωa=1/ τ32+2/ T2=1/ T1+2/ T2, where we have used the usual notation T1 for τ32. The collision time T2 is usually much smaller than the lifetime T1. The constant κ is given by
κ=6π ϵ0 c3/ ωa2 T1.
(A4)
We have chosen ωa=4.215× 1015Hz corresponding to a gain curve centered at λa=446.9nm and the following values T1=100ps and T2=20fs. Hence, the gain curve has a spectral width Δ λa=11nm.

APPENDIX B: QUASIMODE CONTRIBUTIONS TO LASING MODE SUPPRESSION

Figure 11 displays the quasimode quality factors (Q-factors) and the spectral distance from the gain center wavelength for 10 realizations of passive 1D random structures. The Q- factor of a quasimode is defined by Q=ν/δν, where ν is the quasimode frequency and δν its width equal to the reciprocal of its lifetime τ. One observes that the quasimodes associated with the first lasing modes are all located near the wavelength corresponding to the maximum of the gain curve while the quasimodes associated to the second lasing modes are all farther from λa. In this case, for a quasimode likely to be lasing, the proximity to λa is more important than the value of the Q-factor.

In Fig. 4, it was seen that two realizations of disorder for 1D random lasers exhibited an extremely robust single-mode lasing regime. The quasimode properties of those two realizations (7 and 10) are shown explicitly in Fig. 11. For both realizations, the quasimode associated with the first lasing mode is nearly coincident with λa. For realization 7, the quasimode associated with the second lasing mode is away from λa and its Q-factor is half of the first. For realization 10, the Q-factors of both quasimodes are the same but the second quasimode is farther from the gain center. Thus, both wavelength and Q-factor play a role in the suppression of the second lasing mode.

Our observation of robust single-mode lasing seems to be in conflict with experimental results for which single-mode random lasers are scarce. We believe that our observation of robust single-mode lasing is a consequence of the small size of our systems. Note that experimentally, preliminary evidence of robust single-mode lasing is described in [48

R. Bardoux, A. Kaneta, M. Funato, K. Okamoto, Y. Kawakami, A. Kikuchi, and K. Kishino, “Single mode emission and non-stochastic laser system based on disordered point-sized structures: toward a tuneable random laser,” Opt. Express 19, 9262–9268 (2011). [CrossRef] [PubMed]

].

ACKNOWLEDGMENTS

We thank O. Alibart and H. Cao for stimulating discussions. This work was supported by the Agence Nationale de la Recherche (ANR) under Grant No. ANR-08-BLAN-0302-01, the Provence-Alpes-Côte-d'Azur (PACA) region, the CG06, and the Groupement de Recherche 3219 MesoImage. JA acknowledges support from the the Embassy of France in the United States. This work was performed using high performance computing resources from GENCI-CINES (Grant 2010-99660).

References and links

1.

V. S. Letokhov, “Generation of light by a scattering medium with negative resonance absorption,” Sov. Phys. JETP 26, 835–840 (1968).

2.

H. Cao, “Lasing in random media,” Waves Random Media 13, R1–R39 (2003) and references therein. [CrossRef]

3.

Y. Li, X. Ma, M. Xu, L. Xiang, and D. Yang, “Remarkable decrease in threshold for electrically pumped random ultraviolet lasing from ZnO fim by incorporation of Zn 2 TiO 4 nanoparticles,” Opt. Express 19, 8662–8669 (2011). [CrossRef] [PubMed]

4.

S. John and G. Pang, “Theory of lasing in a multiple-scattering medium,” Phys. Rev. A 54, 3642–3652 (1996). [CrossRef] [PubMed]

5.

D. S. Wiersma and A. Lagendijk, “Light diffusion with gain and random lasers,” Phys. Rev. E 54, 4256–4265 (1996). [CrossRef]

6.

H. Cao, J. Y. Xu, S.-H. Chang, and S. T. Ho, “Transition from amplified spontaneous emission to laser action in strongly scattering media,” Phys. Rev. E 61, 1985–1989 (2000). [CrossRef]

7.

X. Jiang and C. M. Soukoulis, “Time dependent theory for random lasers,” Phys. Rev. Lett. 85, 70–73 (2000). [CrossRef] [PubMed]

8.

C. Vanneste and P. Sebbah, “Selective excitation of localized modes in active random media,” Phys. Rev. Lett. 87, 183903 (2001). [CrossRef]

9.

R. C. Polson, M. E. Raikh, and Z. V. Vardeny, “Universality in unintentional laser resonators in π-conjugated polymer films,” C. R. Acad. Sci. Ser. IV A, 509–521 (2002).

10.

V. M. Apalkov, M. E. Raikh, and B. Shapiro, “Random resonators and prelocalized modes in disordered dielectric films,” Phys. Rev. Lett. 89, 016802 (2002). [CrossRef] [PubMed]

11.

A. Tulek, R. C. Polson, and Z. V. Vardeny, “Naturally occurring resonators in random lasing of π-conjugated polymer films,” Nat. Phys. 6, 303–310 (2010). [CrossRef]

12.

S. Mujumdar, M. Ricci, R. Torre, and D. S. Wiersma, “Amplified extended modes in random lasers,” Phys. Rev. Lett. 93, 053903 (2004). [CrossRef] [PubMed]

13.

C. Vanneste, P. Sebbah, and H. Cao, “Lasing with resonant feedback in weakly scattering random systems,” Phys. Rev. Lett. 98, 143902 (2007). [CrossRef] [PubMed]

14.

H. E. Türeci, L. Ge, S. Rotter, and A. D. Stone, “Strong interactions in multimode random lasers,” Science 320, 643–646 (2008). [CrossRef] [PubMed]

15.

H. E. Türeci, A. D. Stone, L. Ge, S. Rotter, and R. J. Tandy, “ Ab initio self-consistent laser theory and random lasers,” Nonlinearity 22, C1–C18 (2009). [CrossRef]

16.

O. Zaitsev, L. Deych, and V. Shuvayev, “Statistical properties of one-dimensional random lasers,” Phys. Rev. Lett. 102, 043906 (2009). [CrossRef] [PubMed]

17.

O. Zaitsev and L. Deych, “Recent developments in the theory of multimode random lasers,” J. Opt. 12, 024001 (2010). [CrossRef]

18.

O. Zaitsev and L. Deych, “Diagrammatic semiclassical laser theory,” Phys. Rev. A 81, 023822 (2010). [CrossRef]

19.

J. Andreasen, A. Asatryan, L. Botten, M. Byrne, H. Cao, L. Ge, L. Labonté, P. Sebbah, A. D. Stone, H. E. Türeci, and C. Vanneste, “Modes of random lasers,” Adv. Opt. Photon. 3, 88–127 (2011). [CrossRef]

20.

S. M. Dutra and G. Nienhuis, “Quantized modes of a leaky cavity,” Phys. Rev. B 62, 063805 (2000).

21.

H. Cao, X. Jiang, Y. Ling, J. Y. Xu, and C. M. Soukoulis, “Mode repulsion and mode coupling in random lasers,” Phys. Rev. B 67, 161101(R) (2003). [CrossRef]

22.

X. Jiang, S. Feng, C. M. Soukoulis, J. Zi, J. D. Joannopoulos, and H. Cao, “Coupling, competition, and stability of modes in random lasers,” Phys. Rev. B 69, 104202 (2004). [CrossRef]

23.

C. Conti, M. Leonetti, A. Fratalocchi, L. Angelani, and G. Ruocco, “Condensation in disordered lasers: Theory, 3d+1 simulations, and experiments,” Phys. Rev. Lett. 101, 143901 (2008). [CrossRef] [PubMed]

24.

H. E. Türeci, A. D. Stone, and B. Collier, “Self-consistent multimode lasing theory for complex or random lasing media,” Phys. Rev. A 74, 043822 (2006). [CrossRef]

25.

L. Ge, R. J. Tandy, A. D. Stone, and H. E. Türeci, “Quantitative verification of ab initio self-consistent laser theory,” Opt. Express 16, 16895–16902 (2008). [CrossRef] [PubMed]

26.

E. Roldán, G. J. de Valcárcel, F. Prati, F. Mitschke, and T. Voigt, “Multilongitudinal mode emission in ring cavity class B lasers,” in “Trends in Spatiotemporal Dynamics in Lasers. Instabilities, Polarization Dynamics, and Spatial Structures ,” O. Gomez-Calderon and J. M. Guerra, eds. (Research Signpost, 2005), pp. 1–80.

27.

J. Andreasen, P. Sebbah, and C. Vanneste, “Coherent instabilities in random lasers,” Phys. Rev. A 84, 023826 (2011). [CrossRef]

28.

J. Andreasen and H. Cao, “Creation of new lasing modes with spatially nonuniform gain,” Opt. Lett. 34, 3586–3588 (2009). [CrossRef] [PubMed]

29.

J. Andreasen and H. Cao, “Numerical study of amplified spontaneous emission and lasing in random media,” Phys. Rev. A 82, 063835 (2010). [CrossRef]

30.

P. Sebbah and C. Vanneste, “Random laser in the localized regime,” Phys. Rev. B 66, 144202 (2002). [CrossRef]

31.

J. Andreasen and H. Cao, “Spectral behavior of partially pumped weakly scattering random lasers,” Opt. Express 19, 3418–3433 (2011). [CrossRef] [PubMed]

32.

A. Taflove and S. Hagness, Computational Electrodynamics (Artech House, 2005), 3rd ed.

33.

A. S. Nagra and R. A. York, “FDTD analysis of wave propagation in nonlinear absorbing and gain media,” IEEE Trans. Antennas Propag. 46, 334–340 (1998). [CrossRef]

34.

A. E. Siegman, Lasers (University Science Books, 1986).

35.

C. M. Soukoulis, X. Jiang, J. Y. Xu, and H. Cao, “Dynamic response and relaxation oscillations in random lasers,” Phys. Rev. B 65, 041103(R) (2002). [CrossRef]

36.

M. A. Noginov, G. Zhu, A. A. Frantz, J. Novak, S. N. Williams, and I. Fowlkes, “Dependence of Nd Sc 3( BO 3 ) 4 random laser parameters on particle size,” J. Opt. Soc. Am. B 21, 191–200 (2004) and references cited therein. [CrossRef]

37.

K. L. van der Molen, A. P. Mosk, and A. Lagendijk, “Relaxation oscillations in long-pulsed random lasers,” Phys. Rev. A 80, 055803 (2009). [CrossRef]

38.

J. Andreasen, C. Vanneste, L. Ge, and H. Cao, “Effects of spatially nonuniform gain on lasing modes in weakly scattering random systems,” Phys. Rev. A 81, 043818 (2010). [CrossRef]

39.

X. Wu, W. Fang, A. Yamilov, A. A. Chabanov, A. A. Asatryan, L. C. Botten, and H. Cao, “Random lasing in weakly scattering systems,” Phys. Rev. A 74, 053812 (2006). [CrossRef]

40.

C. Vanneste and P. Sebbah, “Complexity of two-dimensional quasimodes at the transition from weak scattering to Anderson localization,” Phys. Rev. A 79, 041802(R) (2009). [CrossRef]

41.

B. Liu, A. Yamilov, Y. Ling, J. Y. Xu, and H. Cao, “Dynamic nonlinear effect on lasing in a random medium,” Phys. Rev. Lett. 91, 063903 (2003). The surprising drift of mode 1 across the maximum of the gain curve, toward mode 2 [Fig. 5b], is attributed to the nonlinear Kerr effect. [CrossRef] [PubMed]

42.

R. C. Miller, “Optical harmonic generation in single crystal BaTiO 3 ,” Phys. Rev. 134, A1313–A1319 (1964). [CrossRef]

43.

C. F. Dewey Jr. and L. O. Hocker, “Enhanced nonlinear optical effects in rotationally twinned crystals,” Appl. Phys. Lett. 26, 442–444 (1975). [CrossRef]

44.

E. Y. Morozov and A. S. Chirkin, “Stochastic quasi-phase matching in nonlinear-optical crystals with an irregular domain structure,” Quantum Electron. 34, 227–232 (2004). [CrossRef]

45.

V. A. Mel’nikov, L. A. Golovan, S. O. Konorov, D. A. Muzychenko, A. B. Fedotov, A. M. Zheltikov, V. Y. Timoshenko, and P. K. Kashkarov, “Second-harmonic generation in strongly scattering porous gallium phosphide,” Appl. Phys. B 79, 225–228 (2004). [CrossRef]

46.

M. Baudrier-Raybaut, R. Haïdar, P. Kupecek, P. Lemasson, and E. Rosencher, “Random quasi-phase-matching in bulk polycrystalline isotropic nonlinear materials,” Nature 432, 374–376 (2004). [CrossRef] [PubMed]

47.

S. E. Skipetrov, “Disorder is the new order,” Nature 432, 285–286 (2004). [CrossRef] [PubMed]

48.

R. Bardoux, A. Kaneta, M. Funato, K. Okamoto, Y. Kawakami, A. Kikuchi, and K. Kishino, “Single mode emission and non-stochastic laser system based on disordered point-sized structures: toward a tuneable random laser,” Opt. Express 19, 9262–9268 (2011). [CrossRef] [PubMed]

Fig. 1 Spatially dependent index of refraction. (a)  n(x) of a 1D random structure. (b)  n(x,y) of a 2D random structure, 5×5 μm 2. The outer black border represents the absorbing boundary.
Fig. 2 Relaxation oscillations in (a) population inversion and (b) intensity for Pr=0.25 ns 1 (solid black line), Pr=0.26 ns 1 (solid gray line), Pr=0.27 ns 1 (dashed dark red line), and Pr=0.28 ns 1 (dotted red line) where single-mode lasing occurs. The intensity and oscillation frequency increases with Pr.
Fig. 3 (a) Differences between the intensity distributions of the threshold lasing mode with uniform and flat gain and that with successively including (PP) partial gain, (FD) frequency- dependent gain, and (GS) gain saturation. (Triangles) first lasing mode, (squares) second lasing mode. (Open and closed squares) two different realizations for second lasing modes to verify the effect of gain saturation. (b) Intensity distributions of second threshold lasing mode with (solid red line) uniform gain and (dashed black line) gain saturation with gain only in the air gaps.
Fig. 4 Lasing thresholds for 10 realizations of 1D random lasers. (a) First (filled symbols) and second (open symbols) lasing thresholds without mode interaction (circles) and with mode interaction (triangles) included. (b) Ratio of the second lasing threshold over the first lasing threshold without and with mode interaction included. The single-mode regime persisted for realizations 7 and 10, even for the largest ratio checked (9000).
Fig. 5 (a) Intensity and (b) wavelength versus pumping rate Pr of the (solid black line) first and (dotted gray line) second lasing mode. Their respective thresholds are at Pr=0.24 ns 1 and Pr=0.30 ns 1. Mode 1 is suppressed for Pr0.44 ns 1.
Fig. 6 Spectrograms of output intensity (color on a log scale ranging from 10 108) of a 1D random laser. (a)  Pr=0.24 ns 1, single-mode lasing. (b)  Pr=0.29 ns 1, just below the lasing threshold of the second mode. The second mode appears in the transient regime. (c–e)  Pr=0.30, 0.31, 0.43 ns 1, multimode lasing. (f)  Pr=0.47 ns 1, the first lasing mode is suppressed, though it appears in the transient regime.
Fig. 7 Differences D between (circles) lasing mode 1 and lasing mode 1 at threshold, (triangles) lasing mode 2 and lasing mode 2 at threshold, and (squares) lasing mode 1 and 2 at the same pumping rate Pr.
Fig. 8 Intensity versus pumping rate Pr of the (solid black line) first, (dotted gray line) second, and (dashed dark-gray line) third lasing mode of a 2D random laser, size 1×1 μm 2.
Fig. 9 Emission spectrum at the threshold of the first lasing mode and (inset) peak from third-harmonic generation for a (a) 1D random laser at Pr=0.24 ns 1 and (b) 2D random laser at Pr=2.50 ns 1.
Fig. 10 Emission spectrum for a (a–b) 1D random laser at Pr=0.30 ns 1 and (c–d) 2D random laser at Pr=3.00 ns 1. (a,c) Two lasing modes (labeled 1 and 2) and the peaks resulting from four-wave mixing (labeled 3 and 4). (b,d) Peaks resulting from third-harmonic generation (labeled 1t, 2t) and sum-frequency generation (labeled 3t, 4t).
Fig. 11 Quasimode Q-factors and spectral distance from the gain center wavelength λa for 10 realizations of passive random structures in 1D. Quasimodes correspond to the first lasing mode (red circles) and second lasing mode (blue triangles) found without mode interaction (TM); the first and second modes for the same realization of disorder are connected with a line. With mode interaction (FDTD), the two cases in which the first lasing mode strongly suppresses the second lasing mode are marked as 7 and 10.

OCIS Codes
(140.3460) Lasers and laser optics : Lasers
(290.4210) Scattering : Multiple scattering
(190.4223) Nonlinear optics : Nonlinear wave mixing
(260.2710) Physical optics : Inhomogeneous optical media

ToC Category:
Lasers and Laser Optics

History
Original Manuscript: July 29, 2011
Revised Manuscript: September 14, 2011
Manuscript Accepted: September 30, 2011
Published: November 18, 2011

Citation
Jonathan Andreasen, Patrick Sebbah, and Christian Vanneste, "Nonlinear effects in random lasers," J. Opt. Soc. Am. B 28, 2947-2955 (2011)
http://www.opticsinfobase.org/josab/abstract.cfm?URI=josab-28-12-2947


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. V. S. Letokhov, “Generation of light by a scattering medium with negative resonance absorption,” Sov. Phys. JETP 26, 835–840(1968).
  2. H. Cao, “Lasing in random media,” Waves Random Media 13, R1–R39 (2003) and references therein. [CrossRef]
  3. Y. Li, X. Ma, M. Xu, L. Xiang, and D. Yang, “Remarkable decrease in threshold for electrically pumped random ultraviolet lasing from ZnO fim by incorporation of Zn2TiO4 nanoparticles,” Opt. Express 19, 8662–8669 (2011). [CrossRef] [PubMed]
  4. S. John and G. Pang, “Theory of lasing in a multiple-scattering medium,” Phys. Rev. A 54, 3642–3652 (1996). [CrossRef] [PubMed]
  5. D. S. Wiersma and A. Lagendijk, “Light diffusion with gain and random lasers,” Phys. Rev. E 54, 4256–4265 (1996). [CrossRef]
  6. H. Cao, J. Y. Xu, S.-H. Chang, and S. T. Ho, “Transition from amplified spontaneous emission to laser action in strongly scattering media,” Phys. Rev. E 61, 1985–1989 (2000). [CrossRef]
  7. X. Jiang and C. M. Soukoulis, “Time dependent theory for random lasers,” Phys. Rev. Lett. 85, 70–73 (2000). [CrossRef] [PubMed]
  8. C. Vanneste and P. Sebbah, “Selective excitation of localized modes in active random media,” Phys. Rev. Lett. 87, 183903(2001). [CrossRef]
  9. R. C. Polson, M. E. Raikh, and Z. V. Vardeny, “Universality in unintentional laser resonators in π-conjugated polymer films,” C. R. Acad. Sci. Ser. IV A, 509–521 (2002).
  10. V. M. Apalkov, M. E. Raikh, and B. Shapiro, “Random resonators and prelocalized modes in disordered dielectric films,” Phys. Rev. Lett. 89, 016802 (2002). [CrossRef] [PubMed]
  11. A. Tulek, R. C. Polson, and Z. V. Vardeny, “Naturally occurring resonators in random lasing of π-conjugated polymer films,” Nat. Phys. 6, 303–310 (2010). [CrossRef]
  12. S. Mujumdar, M. Ricci, R. Torre, and D. S. Wiersma, “Amplified extended modes in random lasers,” Phys. Rev. Lett. 93, 053903(2004). [CrossRef] [PubMed]
  13. C. Vanneste, P. Sebbah, and H. Cao, “Lasing with resonant feedback in weakly scattering random systems,” Phys. Rev. Lett. 98, 143902 (2007). [CrossRef] [PubMed]
  14. H. E. Türeci, L. Ge, S. Rotter, and A. D. Stone, “Strong interactions in multimode random lasers,” Science 320, 643–646 (2008). [CrossRef] [PubMed]
  15. H. E. Türeci, A. D. Stone, L. Ge, S. Rotter, and R. J. Tandy, “Ab initio self-consistent laser theory and random lasers,” Nonlinearity 22, C1–C18 (2009). [CrossRef]
  16. O. Zaitsev, L. Deych, and V. Shuvayev, “Statistical properties of one-dimensional random lasers,” Phys. Rev. Lett. 102, 043906(2009). [CrossRef] [PubMed]
  17. O. Zaitsev and L. Deych, “Recent developments in the theory of multimode random lasers,” J. Opt. 12, 024001 (2010). [CrossRef]
  18. O. Zaitsev and L. Deych, “Diagrammatic semiclassical laser theory,” Phys. Rev. A 81, 023822 (2010). [CrossRef]
  19. J. Andreasen, A. Asatryan, L. Botten, M. Byrne, H. Cao, L. Ge, L. Labonté, P. Sebbah, A. D. Stone, H. E. Türeci, and C. Vanneste, “Modes of random lasers,” Adv. Opt. Photon. 3, 88–127 (2011). [CrossRef]
  20. S. M. Dutra and G. Nienhuis, “Quantized modes of a leaky cavity,” Phys. Rev. B 62, 063805 (2000).
  21. H. Cao, X. Jiang, Y. Ling, J. Y. Xu, and C. M. Soukoulis, “Mode repulsion and mode coupling in random lasers,” Phys. Rev. B 67, 161101(R) (2003). [CrossRef]
  22. X. Jiang, S. Feng, C. M. Soukoulis, J. Zi, J. D. Joannopoulos, and H. Cao, “Coupling, competition, and stability of modes in random lasers,” Phys. Rev. B 69, 104202 (2004). [CrossRef]
  23. C. Conti, M. Leonetti, A. Fratalocchi, L. Angelani, and G. Ruocco, “Condensation in disordered lasers: Theory, 3d+1 simulations, and experiments,” Phys. Rev. Lett. 101, 143901 (2008). [CrossRef] [PubMed]
  24. H. E. Türeci, A. D. Stone, and B. Collier, “Self-consistent multimode lasing theory for complex or random lasing media,” Phys. Rev. A 74, 043822 (2006). [CrossRef]
  25. L. Ge, R. J. Tandy, A. D. Stone, and H. E. Türeci, “Quantitative verification of ab initio self-consistent laser theory,” Opt. Express 16, 16895–16902 (2008). [CrossRef] [PubMed]
  26. E. Roldán, G. J. de Valcárcel, F. Prati, F. Mitschke, and T. Voigt, “Multilongitudinal mode emission in ring cavity class B lasers,” in “Trends in Spatiotemporal Dynamics in Lasers. Instabilities, Polarization Dynamics, and Spatial Structures,” O.Gomez-Calderon and J.M.Guerra, eds. (Research Signpost, 2005), pp. 1–80.
  27. J. Andreasen, P. Sebbah, and C. Vanneste, “Coherent instabilities in random lasers,” Phys. Rev. A 84, 023826 (2011). [CrossRef]
  28. J. Andreasen and H. Cao, “Creation of new lasing modes with spatially nonuniform gain,” Opt. Lett. 34, 3586–3588 (2009). [CrossRef] [PubMed]
  29. J. Andreasen and H. Cao, “Numerical study of amplified spontaneous emission and lasing in random media,” Phys. Rev. A 82, 063835 (2010). [CrossRef]
  30. P. Sebbah and C. Vanneste, “Random laser in the localized regime,” Phys. Rev. B 66, 144202 (2002). [CrossRef]
  31. J. Andreasen and H. Cao, “Spectral behavior of partially pumped weakly scattering random lasers,” Opt. Express 19, 3418–3433(2011). [CrossRef] [PubMed]
  32. A. Taflove and S. Hagness, Computational Electrodynamics (Artech House, 2005), 3rd ed.
  33. A. S. Nagra and R. A. York, “FDTD analysis of wave propagation in nonlinear absorbing and gain media,” IEEE Trans. Antennas Propag. 46, 334–340 (1998). [CrossRef]
  34. A. E. Siegman, Lasers (University Science Books, 1986).
  35. C. M. Soukoulis, X. Jiang, J. Y. Xu, and H. Cao, “Dynamic response and relaxation oscillations in random lasers,” Phys. Rev. B 65, 041103(R) (2002). [CrossRef]
  36. M. A. Noginov, G. Zhu, A. A. Frantz, J. Novak, S. N. Williams, and I. Fowlkes, “Dependence of NdSc3(BO3)4 random laser parameters on particle size,” J. Opt. Soc. Am. B 21, 191–200 (2004) and references cited therein. [CrossRef]
  37. K. L. van der Molen, A. P. Mosk, and A. Lagendijk, “Relaxation oscillations in long-pulsed random lasers,” Phys. Rev. A 80, 055803 (2009). [CrossRef]
  38. J. Andreasen, C. Vanneste, L. Ge, and H. Cao, “Effects of spatially nonuniform gain on lasing modes in weakly scattering random systems,” Phys. Rev. A 81, 043818 (2010). [CrossRef]
  39. X. Wu, W. Fang, A. Yamilov, A. A. Chabanov, A. A. Asatryan, L. C. Botten, and H. Cao, “Random lasing in weakly scattering systems,” Phys. Rev. A 74, 053812 (2006). [CrossRef]
  40. C. Vanneste and P. Sebbah, “Complexity of two-dimensional quasimodes at the transition from weak scattering to Anderson localization,” Phys. Rev. A 79, 041802(R) (2009). [CrossRef]
  41. B. Liu, A. Yamilov, Y. Ling, J. Y. Xu, and H. Cao, “Dynamic nonlinear effect on lasing in a random medium,” Phys. Rev. Lett. 91, 063903 (2003). The surprising drift of mode 1 across the maximum of the gain curve, toward mode 2 [Fig. ], is attributed to the nonlinear Kerr effect. [CrossRef] [PubMed]
  42. R. C. Miller, “Optical harmonic generation in single crystal BaTiO3,” Phys. Rev. 134, A1313–A1319 (1964). [CrossRef]
  43. C. F. Dewey Jr. and L. O. Hocker, “Enhanced nonlinear optical effects in rotationally twinned crystals,” Appl. Phys. Lett. 26, 442–444 (1975). [CrossRef]
  44. E. Y. Morozov and A. S. Chirkin, “Stochastic quasi-phase matching in nonlinear-optical crystals with an irregular domain structure,” Quantum Electron. 34, 227–232 (2004). [CrossRef]
  45. V. A. Mel’nikov, L. A. Golovan, S. O. Konorov, D. A. Muzychenko, A. B. Fedotov, A. M. Zheltikov, V. Y. Timoshenko, and P. K. Kashkarov, “Second-harmonic generation in strongly scattering porous gallium phosphide,” Appl. Phys. B 79, 225–228 (2004). [CrossRef]
  46. M. Baudrier-Raybaut, R. Haïdar, P. Kupecek, P. Lemasson, and E. Rosencher, “Random quasi-phase-matching in bulk polycrystalline isotropic nonlinear materials,” Nature 432, 374–376 (2004). [CrossRef] [PubMed]
  47. S. E. Skipetrov, “Disorder is the new order,” Nature 432, 285–286 (2004). [CrossRef] [PubMed]
  48. R. Bardoux, A. Kaneta, M. Funato, K. Okamoto, Y. Kawakami, A. Kikuchi, and K. Kishino, “Single mode emission and non-stochastic laser system based on disordered point-sized structures: toward a tuneable random laser,” Opt. Express 19, 9262–9268 (2011). [CrossRef] [PubMed]

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.


« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited