Origin of emission from square-shaped organic microlasers
The emission from open cavities with non-integrable features remains a challenging problem of practical as well as fundamental relevance. Square-shaped dielectric microcavities provide a favorable case study with generic implications for other polygonal resonators. We report on a joint experimental and theoretical study of square-shaped organic microlasers exhibiting a far-field emission that is strongly concentrated in the directions parallel to the side walls of the cavity. A semiclassical model for the far-field distributions is developed that is in agreement with even fine features of the experimental findings. Comparison of the model calculations with the experimental data allows the precise identification of the lasing modes and their emission mechanisms, providing strong support for a physically intuitive ray-dynamical interpretation. Special attention is paid to the role of diffraction and the finite side length.
pacs:42.55.Sa, 03.65.Sq, 05.45.Mt
Semiclassical physics emerged during the development of quantum mechanics, almost one century ago, to account for the transition from wave physics in the quantum regime to classical mechanics Brack and Bhaduri (2003). Today, semiclassical physics plays an essential role in quantum chaos Stöckmann (2000), and its methods are applied in virtually any field that features wave dynamics, including acoustics Tanner and Søndergaard (2007), electromagnetism Stöckmann (2000), hydrodynamics Kudrolli et al. (2001), and loop quantum gravity Bianchi and Haggard (2011). We demonstrate the power of these methods when applied to optical microresonators by studying square-shaped polymer-based microlasers. Their emission properties are little understood for two reasons. Firstly, the square resonator with Dirichlet boundary conditions is a standard example of a separable system, whereas the dielectric square resonator is nonseparable and hence nonintegrable due to the diffraction at the dielectric corners. This remains an open problem in mathematical physics with tremendous impact on radar or telecommunication applications Gennarelli et al. (2015). Secondly, it was observed that these lasers emit very narrow lobes in only a few directions. Directional emission from microlasers has been intensely investigated due to possible applications Wiersig et al. (2011) and observed for, e.g., limaçon- and stadium-shaped microlasers. The underlying mechanism is understood in terms of their chaotic ray dynamics Schwefel et al. (2004); Wiersig and Hentschel (2008). However, these explanations cannot be applied to polygonal resonators since they do not exhibit chaotic ray dynamics.
While square and hexagonal microresonators and -lasers have been extensively studied, their far-field distributions have been scarcely investigated Wiersig (2003a); Yang and Huang (2007); Chen et al. (2009); Dai et al. (2009). This is partly because rounded corners significantly influence the resonance frequencies and field distributions Boriskina (2005); Wiersig (2003a); Dietrich et al. (2012); Kudo et al. (2013), and it is technologically challenging to fabricate sharp corners. Here we investigate the far-field emission of microlasers with sharp corners (i.e., with a radius of curvature much smaller than the wavelength). A semiclassical model for the dielectric square Bittner et al. (2013, 2014) is used to predict the far-field distributions, which are in very good agreement with numerical calculations. Careful comparison of the measured far-field distributions with the model allows to identify the lasing modes and understand the mechanism for their high directionality.
The microlasers consisted of a PMMA
The microlasers were pumped by a pulsed, frequency-doubled Nd:YAG laser ( nm, ns, Hz). The short pulses and low duty cycle help to avoid problems from heating and quenching due to dark states. The pump beam impinged vertically and it covered the whole cavity uniformly The lasing emission was collected in the sample plane by a lens and transferred to a spectrometer by an optical fiber. It was polarized parallel to the sample plane. The far-field intensity distributions were measured as a function of the azimuthal angle by rotating the cavity.
A typical spectrum is shown in Fig. 2. It features a sequence of equidistant peaks. Its Fourier transform (right inset of Fig. 2) shows that the free spectral range (FSR) of the spectrum corresponds to the optical length of the diamond orbit Lebental et al. (2007); Bogomolny et al. (2011). Therefore it is expected that the observed lasing modes are localized on trajectories with an angle of incidence close to .
A typical far-field intensity distribution of a single resonance is presented in Fig. 3(a). It exhibits four emission lobes in the directions parallel to the side walls. The emission lobes are very narrow with a full width at half maximum of about [see Fig. 3(b)]. Figure 1(b) shows an image of the lasing cavity taken by a camera with a high-magnification zoom lens. It evidences that the lasing emission (red) is emitted from the side walls parallel to the emission direction, whereas the side walls perpendicular to it do not emit but only scatter the green pump laser. Since the emission is at a grazing angle, it must stem from rays impinging on the cavity side walls with an angle of incidence very close to the critical one, , where is the refractive index Lebental et al. (2007). This is in agreement with the previous observation that the ray trajectories on which the lasing modes are based have an angle of incidence close to .
Iii Model calculations
The far-field distributions are calculated using the semiclassical model introduced in Ref. Bittner et al. (2013, 2014) to identify the observed lasing modes. The model is based on the observation that the wave functions are localized on classical tori and are thus composed of eight plane waves . Their directions are related to each other by the symmetry operations of the point group. The corresponding ray trajectories in Fig. 4 are given by the momentum vectors and their rotations by 90, .
To quantize the spectrum, the simplest approach is to separate the system along the and directions, and write two phase loop conditions analogous to that of a Fabry-Pérot cavity Bittner et al. (2013, 2014),
The main difference from a Fabry-Pérot cavity are the Fresnel coefficients corresponding to an incidence with angles instead of . The momentum vector components can thus be formally written as
The resonance wave number is where is the free-space wavelength. Note that are in general complex-valued. All wave functions are obtained by adding up the plane waves with their correct momentum vectors and signs. For transverse magnetic (TM) [transverse electric (TE)] polarization, the electric (magnetic) field is parallel to the plane of the cavity, and corresponds to the component of the electric (magnetic) field. The resonances [i.e., solutions of Eq. (1)] are labeled by their quantum numbers and their symmetries () with respect to the diagonal () as , where () means a symmetric (antisymmetric) wave function. There are six symmetry classes corresponding to the different combinations of and the quantum numbers that are labeled by the Mulliken symbols , and (see Table 1).
|Diagonal||Horizontal/vertical||Parity of||Parity of||Mulliken||Model wave function|
Green’s identity is used to infer the far-field distribution Yang and Huang (2007). The derivation is detailed in the Supplementary Materials along with the comparison to numerical simulations (see Figs. S1 and S2). The full expression for the far-field distribution of an mode is given in Eq. (S6). It consists of terms like
where () for TM (TE) polarization and . Each such term corresponds to emission in the two directions defined by the roots of the argument of the function. If the angle of incidence of the corresponding rays inside the resonator is smaller than , then is simply given by Snell’s law (see Fig. 4), and the emission lobe is called refractive. Then the term is in fact the diffraction pattern of a plane wave going through a slit with a width that is the projection of the side of the square on the emission direction (see Fig. 4). When is complex, corresponding to rays with angle of incidence larger than , the lobe is called nonrefractive and is emitted parallel to the side wall. The nonrefractive lobes are jagged and much broader. They have a much smaller amplitude than the refractive lobes due to the finite imaginary part of . While no energy is transmitted when a plane wave impinges on an infinite dielectric interface with an angle larger than the critical one (i.e., it is totally reflected), this is not the case for an interface of finite length as illustrated by the existence of the nonrefractive emission lobes Wiersig (2003a), which leads to a modification of the reflection and transmission coefficients.
Two different cases of far-field diagrams can appear for , depending on . If , the plane waves impinging on a side wall with angle escape refractively, whereas those impinging with are totally reflected. Hence refractive and nonrefractive lobes are observed, where the latter are typically negligible compared to the former. If , all plane waves are totally reflected, and nonrefractive lobes are observed, though their directions are four-fold degenerate.
Iv Comparison with experiment
We now compare the model with experiments. A spectrum generated by the model and corresponding to the experimental conditions is shown in Fig. 5. The resonances are arranged in branches. Each of them consists of modes with identical longitudinal quantum number and increasing transverse quantum number as decreases Bittner et al. (2014). The horizontal distance between these branches agrees well with the experimentally observed FSR, corresponding to the length of the diamond orbit. Therefore we conclude that the observed lasing resonances belong to different branches, i.e., each has a different .
We use the far-field intensity distributions to infer and thus their transverse quantum numbers. The model far-field patterns of three representative modes are presented in Fig. 6. They are labeled by (a), (b), and (c) in Fig. 5. Mode (c) is located just below the critical angle () and features two refractive emission lobes between and , while its nonrefractive lobes are not visible due to their negligible amplitude. Since the measured far-field patterns feature only lobes in the total range of , such modes with can be excluded. Mode (b) with is located just above the critical angle. The refractive lobes have become nonrefractive and merged into four narrow lobes parallel to the sides. This is precisely the kind of far-field pattern observed experimentally. Mode (a) finally is well above the critical angle (). The nonrefractive lobes have become very broad and jagged unlike what was observed experimentally. Also their amplitude has decreased significantly. This development continues with increasing . The qualitative comparison with the measured far-field intensity distribution shows that the experimental lasing modes are only consistent with mode (b), i.e., an angle of incidence just above the critical angle. This is confirmed by the fits described in the following.
The experimental far-field distribution measured with an angular resolution of 1 as in Fig. 3 is compared to the envelope of the model far-field distributions. This is detailed in the Supplementary Materials. The only two independent parameters are the amplitude and , which is considered a continuous variable due to the high resonance density. The fit of the lobe yields and is plotted in Fig. 3(b). The agreement is excellent. Similar results are obtained for the other lobes and the other resonances as well as squares with various sizes, yielding values of from to . As expected from the qualitative considerations in the previous paragraph and the photograph in Fig. 1(b), the angle of incidence is slightly above the critical angle.
In fact, the model predicts far-field distributions with a very quickly oscillating substructure that stems from the terms in the square brackets in Eq. (3), namely the and functions with the argument . The actual model far-field intensity distribution for the fitted parameters is shown as thin green line in Fig. 3(b). To confirm this prediction, the angular resolution of the setup was improved by putting a slit in front of the collection lens. The far-field intensity distribution measured this way is shown in the inset of Fig. 3(b). It features oscillations with a period of , while the expected period is according to Eq. (3). The agreement is excellent for different wavelengths and cavity sizes.
Systematic measurements with increased angular resolution will allow for a better determination of the symmetry of the lasing resonances. Furthermore, the symmetry of the far-field distributions can be changed by the polarization of the pump laser due to the fluorescence anisotropy of the laser dye Gozhyk et al. (2012): when the square microlaser is pumped with linear polarization parallel to the axis [instead of circular polarization as in Fig. 3(a)], the emission is almost bidirectional with strong emission lobes at and , whereas the emission lobes at and are strongly suppressed (not shown). This effect is a consequence of the lower degree of symmetry induced by the pump polarization. It can be explained by the formation of coherent superpositions of degenerate mode pairs like the and modes. It should be noted that the lasing modes and hence far-field emission of microlasers can also be controlled by localized or selective pumping Ben-Messaoud and Zyss (2005); Bachelard et al. (2012); Liew et al. (2014), whereas controlling it only with the pump polarization is an experimentally simpler approach. These topics are, however, beyond the scope of this article and will be explored in a future publication.
The far-field intensity distributions of square-shaped organic microlasers were investigated and revealed highly directional emission in the four directions parallel to the cavity sides. Sharp corners were essential for obtaining this high directionality because rounded corners significantly change the emission behavior Wiersig (2003a). Analytical formulas for the far-field distributions based on a semiclassical model Bittner et al. (2013, 2014) were developed and showed excellent agreement both with experimental data and with numerical simulations. They could be easily adapted to other polygonal cavities like hexagon or rectangles Yang and Huang (2007) and other refractive indices. The comparison between the model and experiments evidences that the lasing modes are based on trajectories with an angle of incidence almost equal to the critical angle that hence leave the resonator at a grazing angle. Interestingly, these modes do not feature the highest quality factors; these would be modes with larger than . Maybe these modes actually lase, but are not observed because of their lower amplitude in the far-field compared to other modes. Or they may actually have higher losses due to diffraction at the substrate interface Lozenko et al. (2012).
The semiclassical model does not take into account diffraction by the corners, whereas it is precisely this effect which renders the system nonintegrable. While the model calculations agree very well with the observed spectrum and far-field distributions, there are some minor observations that cannot be explained by them. For example, some light is emitted into other directions than the four main emission lobes due to diffraction at the corners, but it is much less intense
Acknowledgements.The numerical calculations are based on a code developed by C. Schmit. S. B. gratefully acknowledges funding from the European Union Seventh Framework Programme (FP7/2007-2013) under Grant No. 246.556.10. This work was supported by a public grant from the Laboratoire d’Excellence Physics Atom Light Matter (LabEx PALM) overseen by the French National Research Agency (ANR) as part of the Investissements d’Avenir program (Reference No. ANR-10-LABX-0039).
Appendix A Supplementary Material: Calculation of the far-field distributions
We use Green’s identity,
to infer the far-field distribution (cf. Ref. Yang and Huang (2007)), where is the derivative with respect to the surface normal of the domain and
is the Green’s function in two dimensions with being the wave number. We choose the boundary of the resonator as the integration path and assume continuity of and at the resonator boundaries where () for TM (TE) polarization Lebental et al. (2007). Then the wave function outside of the resonator, , is related to that inside, , by
We start with a single plane wave inside the resonator and calculate the corresponding outside wave function in the far field. The integral can be simplified in this case where and with by using for . Furthermore, in this limit . The integration yields
It is instructive to replace and and look at the maxima of the two terms in the large square brackets. The first term stems from the integration along the side walls parallel to the axis and is maximal when its denominator vanishes, i.e., for . This is simply Snell’s law for a ray refracted at a side wall parallel to the axis, and analogously the maximum of the second term corresponds to a ray refracted at a side wall parallel to the axis.
The results for the full wave functions of the dielectric square are obtained by adding up the plane waves with their correct momentum vectors and relative amplitudes. For the symmetry class , this yields
Appendix B Supplementary Material: Comparison with numerical calculations
To validate this semiclassical model, numerical simulations were performed using the boundary element method Wiersig (2003b); Lebental et al. (2007). The complete spectrum of TM modes as well as the near- and far-field distributions were calculated in the range of – for . Figure S1 shows the numerically calculated spectrum. For each resonance, we inferred the quantum numbers and , and then the corresponding , by comparing the near-field distributions to the model predictions. The very regular structure of the spectrum and the possibility to quantize it through and confirm that the dielectric square resonator is nearly an integrable system Peres (1984); Bittner et al. (2014).
Figure S2 presents typical far-field distributions, one with well below [panel (a)], and the other with just above [panel (b)]. In Fig. SS2, the mode exhibits at total of narrow refractive emission lobes. Their directions agree very well with those expected from ray optics indicated by the arrows. The analytically calculated far-field distribution (red line) agrees almost perfectly with the numerical calculation (black line). The nonrefractive emission lobes cannot be recognized since their amplitudes are negligible compared to those of the refractive lobes. It should be noted that the lobes of this example are much broader than those of the corresponding example in Fig. 6(c) of the main text for larger since the ratio between the effective slit size and the wavelength is much smaller here. In Fig. SS2, the directions of the various lobes are very well reproduced by the model although it sometimes fails to predict their amplitude correctly. In general, the ray-based model is less accurate for modes with since the Fresnel reflection coefficients for an infinite interface that are used in the quantization condition are no longer a good approximation close to and above the critical angle. This is evidenced by the prediction of zero losses for such modes by the model, whereas all resonances naturally have a finite in reality. The nonrefractive lobes are no longer narrow but broad and jagged. The model has only been tested for moderately sized cavities , and not for – like those investigated experimentally, because simulations in this size region are extremely cumbersome. However, thanks to the semiclassical nature of the model, it seems reasonable that it works well for very large cavities.
Appendix C Supplementary Material: Envelope of the far-field distributions
The fast oscillations of the far-field intensity distribution are not visible with an angular resolution of , and the measured far-field patterns are hence compared with the envelope of the model far-field distributions. The terms in the square brackets in Eq. (S6), for example the first term of the equation,
feature these very fast oscillations due to the and functions with the argument . The envelope of is calculated by considering as an independent variable and replacing it by the value
for which becomes extremal, yielding
Analogous formulas are obtained for the other terms. The function that is fitted to the experimental data can be further simplified since for the parameter range of interest only one of the terms contributes significantly to each emission lobe. So for example the lobe around is fitted by the formula
where the only two fit parameters are the amplitude and the angle of incidence which is related to the momentum vector components via and .
Each lobe of the measured far-field intensity distribution was fitted separately in the range where by Eq. (S10) or a corresponding expression. TE polarization was used. The envelopes of the far-field distributions are almost identical for the different symmetry classes and hence fitting the formula for another symmetry class to the experimental data yields nearly the same parameters and an equally good agreement. The effective refractive index of was assumed to be a constant since the shape of the envelope depends very little on it, and the fits are insofar independent of as the fit yields practically the same value for for a wide range of values of .
- Poly(methyl methacrylate) PMMA A6 resist by Microchem.
- 4-(Dicyanomethylene)-2-methyl-6-(4-dimethylaminostyryl)-4Hpyran by Exciton.
- For a video of the rotating square microlaser please see the Supplementary Material at http://dx.doi.org/10.1209/0295-5075/113/54002
- M. Brack and R. K. Bhaduri, Semiclassical Physics (Westview Press, Oxford, 2003).
- H.-J. Stöckmann, Quantum Chaos: An Introduction (Cambridge University Press, Cambridge, UK, 2000).
- G. Tanner and N. Søndergaard, J. Phys. A 40, R443 (2007).
- A. Kudrolli, M. C. Abraham, and J. P. Gollub, Phys. Rev. E 63, 026208 (2001).
- E. Bianchi and H. M. Haggard, Phys. Rev. Lett. 107, 011301 (2011).
- G. Gennarelli, M. Frongillo, and G. Riccio, IEEE Trans. Antennas Propag. 63, 374 (2015).
- J. Wiersig, J. Unterhinninghofen, Q. Song, H. Cao, M. Hentschel, and S. Shinohara, “Trends in nano- and micro-cavities,” (Bentham Science Publishers, Sharjah, United Arab Emirates, 2011) Chap. 4, pp. 109–152.
- H. G. L. Schwefel, N. B. Rex, H. E. Tureci, R. K. Chang, A. D. Stone, T. Ben-Messaoud, and J. Zyss, J. Opt. Soc. Am. B 21, 923 (2004).
- J. Wiersig and M. Hentschel, Phys. Rev. Lett. 100, 033901 (2008).
- J. Wiersig, Phys. Rev. A 67, 023807 (2003a).
- Y.-D. Yang and Y.-Z. Huang, IEEE J. Quantum Electron. 43, 497 (2007).
- R. C. C. Chen, Y. T. Yu, Y. J. Huang, C. C. Chen, Y. F. Chen, and K. F. Huang, Opt. Lett. 34, 1810 (2009).
- J. Dai, C. X. Xu, K. Zheng, C. G. Lv, and Y. P. Cui, Appl. Phys. Lett. 95, 241110 (2009).
- S. V. Boriskina, IEEE J. Quantum Electron. 41, 857 (2005).
- C. P. Dietrich, M. Lange, T. Böntgen, and M. Grundmann, Appl. Phys. Lett. 101, 141116 (2012).
- H. Kudo, R. Suzuki, and T. Tanabe, Phys. Rev. A 88, 023807 (2013).
- S. Bittner, E. Bogomolny, B. Dietz, M. Miski-Oglu, and A. Richter, Phys. Rev. E 88, 062906 (2013).
- S. Bittner, E. Bogomolny, B. Dietz, M. Miski-Oglu, and A. Richter, Phys. Rev. E 90, 052909 (2014).
- S. Lozenko, N. Djellali, I. Gozhyk, C. Delezoide, J. Lautru, C. Ulysse, J. Zyss, and M. Lebental, J. Appl. Phys. 111, 103116 (2012).
- M. Lebental, N. Djellali, C. Arnaud, J.-S. Lauret, J. Zyss, R. Dubertrand, C. Schmit, and E. Bogomolny, Phys. Rev. A 76, 023830 (2007).
- E. Bogomolny, N. Djellali, R. Dubertrand, I. Gozhyk, M. Lebental, C. Schmit, C. Ulysse, and J. Zyss, Phys. Rev. E 83, 036208 (2011).
- I. Gozhyk, G. Clavier, R. Méallet-Renault, M. Dvorko, R. Pansu, J.-F. Audibert, A. Brosseau, C. Lafargue, V. Tsvirkun, S. Lozenko, S. Forget, S. Chénais, C. Ulysse, J. Zyss, and M. Lebental, Phys. Rev. A 86, 043817 (2012).
- T. Ben-Messaoud and J. Zyss, Appl. Phys. Lett. 86, 241110 (2005).
- N. Bachelard, J. Andreasen, S. Gigan, and P. Sebbah, Phys. Rev. Lett. 109, 033903 (2012).
- S. F. Liew, B. Redding, L. Ge, G. S. Solomon, and H. Cao, Appl. Phys. Lett. 104, 231108 (2014).
- J. Wiersig, J. Opt. A 5, 53 (2003b).
- A. Peres, Phys. Rev. Lett. 53, 1711 (1984).