two - dimensional ( 2d ) billiards long have been studied as model systems with hamiltonian dynamics .
this is in great part due to their seeming simplicity that contrasts the wealth of different dynamical behaviors that they can exhibit , including integrable , chaotic , pseudointegrable , and mixed dynamics .
one interesting class of 2d billiards are polygons , of which triangles are the simplest case .
while some classes of triangular billiards are well understood , many unsolved problems remain for triangles of asymmetric shape .
these open questions concern , for example , the existence , number , and stability with respect to geometric perturbations of their periodic orbits ( pos ) @xcite . while the existence of at least one po , the so - called fagnano s orbit , is assured for acute triangles , and the existence of pos for obtuse triangles with no angle greater than @xmath0
has also been proven @xcite , it is not known whether any po exists at all in an arbitrary obtuse triangular billiard . but even for triangles that are known to have one or more pos , their actual construction is often nontrivial , and even the shortest pos can be quite complicated .
the search for pos in triangular billiards hence stays a field of active research .
two - dimensional billiards are also studied in the context of quantum and wave - dynamical chaos @xcite to understand the manifestation of ray dynamics in the properties of the corresponding wave - dynamical systems .
early experiments concentrated on microwave and acoustic resonators @xcite , and new interest has arisen with the advent of applications like optical microcavities and -lasers @xcite .
in particular , the influence of pos on the spectral and emission properties of microlasers is important in view both of a fundamental understanding of these devices and their applications @xcite .
while many microcavities have circular or deformed circular shape , different types of polygonal microresonators also have been investigated .
examples of such structures include semiconductor @xcite , organic @xcite and crystal microlasers @xcite , silicon and silica microresonators @xcite , vertical - cavity surface - emitting lasers ( vcsels ) @xcite , and hexagonal zinc oxide nanocavities and -rods @xcite . these studies , however , examined only equilateral polygons , while very few experimental investigations of nonequilateral polygonal resonators have been reported @xcite .
furthermore the scattering properties of triangular and other polygonal structures have been studied @xcite .
thus the properties of nonequilateral triangular and polygonal microcavities are poorly understood . in this article
we study the lasing characteristics of triangular organic microlasers of different shapes with an emphasis on the role of symmetries and their absence .
a laser resonator is usually designed to confine light on a specific periodic ray trajectory and the active medium is positioned to provide optimal overlap with the corresponding resonant states . in the experiments reported here ,
the point of view was reversed : the geometry of the cavity was given and defined the various possible pos and the distribution of the active medium , but it was not known _ a priori _ on which po , if any , the lasing modes would be based .
the laser was essentially left free to decide which po was being favored . while that po could often be guessed in advance for simple , symmetric triangles ,
the question of which specific po is chosen was particularly interesting in cases where even the shortest pos are long and complicated or the pos are even not known at all .
furthermore , the lasing modes need not necessarily be based on a particular po . the aim was therefore to understand which parameters determine the dominant po , to what extent the properties of the lasing modes can be explained by that orbit , and which features of triangular microlasers are beyond simple ray - optical explanations .    the article is organized as follows
. section [ sec : triangledynamics ] summarizes the key characteristics of classical triangular billiards and their pos , and treats their implications for dielectric resonators .
section [ sec : exptechnique ] explains the fabrication of our microlasers and the experimental setup for their characterization .
the experimental results for different triangular microlasers are presented in , starting with highly symmetric and well - understood cases and going all the way to completely asymmetric cavities .
section [ sec : conclusions ] concludes with a summary of the results .
triangles are the simplest type of polygons and hence can serve in many respects as paradigms for more general polygons .
the dynamics of a classical polygonal billiard with @xmath1 vertices depends on the internal angles @xmath2 at the vertices , where @xmath3 is the index of the vertex .
a polygonal billiard is called rational if all of its angles can be written as a rational multiple of @xmath4 , that is , @xmath5 , where the @xmath6 and @xmath7 are coprime integer numbers .
all other polygons are called irrational .
the topology of the phase space of a rational polygon is determined by its genus @xmath8 , which is a function of the @xmath6 and @xmath7 @xcite . for @xmath9
the phase space has the topology of a torus and the billiard is integrable .
the three cases of integrable triangle billiards are the equilateral triangle , the right isosceles triangle , and the triangle with angles @xmath10 , @xmath11 , and @xmath12 , which is an equilateral triangle cut in half @xcite . for @xmath13 ,
the phase space resembles a @xmath8-handled sphere and its dynamics are said to be pseudointegrable @xcite .
irrational polygons have ergodic dynamics @xcite .
it should be noted that the subset of rational polygons is dense in the set of all polygons , that is , each irrational polygon can be approximated by a rational triangle to arbitrary precision .
one type of trajectory that plays an important role in the dynamics of both classical and wave - dynamical billiards are periodic orbits , that is , orbits that retrace themselves after a finite number of reflections .
one example is shown in .
while it has been proven that pos exist in any billiard with a smooth ( @xmath14 ) contour @xcite , there is no such theorem for polygonal billiards : it is not known whether any po exists in an arbitrary polygonal billiard .
there are , however , several results for specific cases , one of the oldest concerning acute triangles .
it was proven back in the 18th century by fagnano that the po connecting the feet of the three altitudes ( thick solid red line in ) is the shortest of all possible pos in such triangles @xcite .
it is hence called fagnano s orbit .
the existence of other pos in irrational acute triangles is , however , not evident .
the pos in rational and several other types of triangles can be constructed with the so - called unfolding technique as demonstrated in ( a ) for an isosceles triangle with a top angle of @xmath15 .
we follow a trajectory [ red ( dark gray ) line ] that starts perpendicularly to the height of the triangle by reflecting the triangle each time that the boundary is encountered so the trajectory unfolds into a straight line .
the trajectory returns to its starting point after a finite number of reflections .
the actual po is obtained by folding the red ( dark gray ) line back into the triangle as shown in ( b ) .
furthermore , there are estimates for the number of pos up to a given length in rational triangles , and it has been shown that the pos are dense in their phase space @xcite . one particular class of rational triangles for which even stronger theorems concerning the pos are known
are the so - called veech triangles @xcite .
much less is known for irrational triangles .
the pos in right and isosceles triangles can also be constructed by unfolding @xcite .
it should be noted , though , that even in cases where pos can be constructed with the unfolding technique , the shortest po can be long and complicated compared to the simple examples shown in figs .
[ fig : fagnano ] and [ fig : dbowtiepo ] .
furthermore , the existence of pos in triangles with all angles smaller than @xmath0 has been proven @xcite .
however , no general theorems are known concerning the existence of pos in obtuse triangles with one angle larger than @xmath0 .
all trajectories in polygonal billiards are marginally stable with respect to perturbations of their initial conditions .
pos with an odd number of reflections like fagnano s orbit or the double bow - tie orbit in are isolated , while pos with an even number of reflections are part of a family of pos with parallel trajectories @xcite .
hence , repeating an isolated po an even number of times yields a nonisolated po .
for example , the po indicated as dotted red line in belongs to the family of the twice - repeated fagnano s orbit .
the complete family of a po can be found by unfolding since its other members cover a strip parallel to the po .
this strip is called the po channel .
the po channel of the double bow - tie orbit in the isosceles @xmath15 billiard is indicated as a gray strip in .
the po channel is bounded by two trajectories ( thin black lines ) that touch the corners of the triangle .
these lines are the optical boundaries of the po channel .
parallel trajectories beyond these line are not members of that family as can be verified by unfolding them .
depending on the po and the triangle , the po channel can cover either a part of the billiard like in or the complete billiard .
for example , the pos of the equilateral and right isosceles triangles cover them completely since these triangles tessellate the plane when unfolding .
the flat organic microlasers studied in this article are treated as 2d passive open dielectric resonators since the lasing modes close to threshold can usually be well described by the modes of the passive cavity .
a passive resonator is governed by the scalar helmholtz equation @xmath16 \psi(x , y ) = 0\ ] ] where @xmath17 is the free - space wave number and @xmath18 is the effective refractive index @xmath19 for @xmath20 inside the resonator and the refractive index of the surrounding medium ( air with @xmath21 ) on the outside .
the wave function @xmath22 corresponds either to the @xmath23 component of the electric field , @xmath24 , for transverse magnetic ( tm ) modes or to that of the magnetic field , @xmath25 , for transverse electric ( te ) modes .
the wave functions inside and outside of the resonator are connected by the usual boundary conditions for dielectric interfaces @xcite .    since the typical size of the cavities considered here is in the range of several hundred wavelengths , the resonators are in the so - called semiclassical regime which is the transition regime from classical physics ( @xmath26 ray optics ) to quantum mechanics ( @xmath26 wave optics ) .
semiclassical methods permit to explain various properties of the resonators using concepts and quantities from the dynamics of the corresponding classical billiard systems @xcite , the pos playing an important role in such approximations .
two well - known cases are trace formulas that connect the density of state with the pos @xcite and resonant states localized on pos , so - called scars @xcite and superscars @xcite .
trace formulas and scars were mainly investigated in the context of closed resonators with dirichlet boundary conditions , but the underlying principles can be extended to dielectric resonators .
in fact , a trace formula for dielectric resonators has been developed @xcite , and modes of dielectric resonators localized on classical trajectories are often observed .
this includes gaussian modes that are localized on stable pos @xcite , scar states localized on unstable pos @xcite , and superscar states localized on families of marginally stable pos @xcite or classical tori @xcite .
it has been shown in ref .
@xcite for resonators with dirichlet boundary conditions that superscar states are localized inside the po channel ( see ) due to repeated diffraction at the corners that define the optical boundaries .
a similar effect has been proposed in ref .
@xcite for dielectric resonators even though the diffraction at dielectric corners is not understood @xcite . therefore we expect to find superscarred lasing modes in triangular dielectric resonators with pseudointegrable classical dynamics .
the influence of a po on the properties of a resonator depends on several factors , among them its length , its stability or , in the case of nonisolated orbits , the area covered by its family and its refractive losses @xcite .
the losses depend on the refractive index and the angles of incidence of a po since a ray traveling in a dielectric resonator is reflected and refracted at the side walls according to the fresnel formulas
. the emission directions of the refracted rays are determined by snell s law .
the most long - lived modes of passive dielectric cavities hence mainly exhibit the influence of the shortest and best - confined pos @xcite .
on the other hand , modes based on pos that are not confined by total internal reflection can be observed for laser cavities . the threshold condition for a ray traveling along a po with length @xmath27 in an active medium with linear gain @xmath28
is given by @xmath29 where @xmath30 is the fresnel reflection coefficient for the reflection at the @xmath31th vertex and the product runs over all vertices of the po .
the threshold gain @xmath32 is hence @xmath33 we use as a simple estimate for the threshold of a mode localized on a po . in practice , however , also other parameters can be of importance like the overlap between the gain region and the mode profiles @xcite or the coupling between the molecules of the gain medium and the electric field of a mode @xcite .
therefore , we do not expect quantitative agreement of the measured lasing thresholds with .
more sophisticated approaches are necessary for a quantitative understanding of the lasing thresholds .    it should be noted that even though the pos of a classical billiard can explain many properties of the corresponding resonators , wave effects such as tunneling and diffraction can also have a significant influence .
one example is the existence of so - called diffractive orbits that have one vertex at a diffractive corner of the billiard .
diffractive corners are corners with an angle that is not equal to @xmath34 , where @xmath35 is an integer .
the reflection of a ray at such a corner is not defined in classical mechanics .
in contrast , a wave impinging on it is diffracted and hence scattered into various directions .
this enables to close ray trajectories that impinge on a diffractive corner and thus diffractive orbits can appear in wave - dynamical billiards ( i.e. , resonators ) .
for example , they contribute to the trace formula for resonators with diffractive corners or point scatterers @xcite .
another important point is that the dynamics of wave systems is less sensitive to geometric perturbations than that of classical systems .
therefore the influence of a po on a wave system can survive geometric perturbations even though the perturbation completely eliminates that po in the classical dynamics @xcite .
the organic microlasers consisted of poly(methyl methacrylate ) ( pmma ) doped with @xmath36 wt% of the laser dye dcm .
a solution of pmma and dcm was spin coated on a silicon wafer with a @xmath37-@xmath38m - thick layer of silica .
the thickness of the pmma layer was about @xmath39  nm .
the desired cavity shapes were written by @xmath40-kv electron - beam lithography .
this process allowed us to define the cavity boundaries with nanometric precision and achieve vertical side walls and sharp corners and edges @xcite .
scanning electron microscope ( sem ) images of an equilateral triangle cavity are presented in .
the triangular microlasers considered here had typical side lengths @xmath41 in the range of @xmath42 to @xmath43 m , i.e. , several hundred times larger than the wavelength @xmath44  nm . they are considered two - dimensional ( 2d ) systems with an effective refractive index of @xmath45 since they are only about one wavelength thick and support only a single vertical excitation for each polarization @xcite .
the experimental setup was similar to the one described in ref .
the microlasers were pumped by a pulsed frequency - doubled nd : yag laser ( @xmath46  nm , @xmath47  ps , @xmath48  hz , teem photonics png-002025 - 140 ) that impinged perpendicularly to the cavity plane .
the intensity and the polarization of the pump beam were controlled independently using half- and quarter - wavelength plates and polarizers .
a circularly polarized pump beam was used , and the pump intensity was normally chosen only slightly above threshold .
the pump beam had an approximately gaussian intensity profile and its diameter was adjusted to cover the complete area of a single microlaser .
the lasing emission in the plane of the microlasers was collected in the far field by a lens @xmath49  cm away from the sample and transferred by a fiber to a spectrometer ( spectra pro 2500i , acton research ) and a cooled ccd camera ( pixis 110b , princeton instruments ) .
the spectra were integrated over @xmath48 pump pulses .
the samples could be rotated to record the spectra in all possible directions in the plane of the cavity and thus measure the azimuthal far - field distributions .
the polarization of the lasing emission was determined using a linear polarization filter @xcite .
all of the microlasers presented in the following emitted transverse electrically ( te ) polarized light , i.e. , the lasing emission had an electric field parallel to the plane of the cavities .
furthermore , a complementary metal - oxide semiconductor sensor camera ( , ids imaging ) with a zoom lens ( zoom 6000 , navitar ) was used to take photographs of the lasing cavities .
the observation angle of the camera was chosen at a @xmath50 tilt angle above the plane so the whole cavities could be surveyed .
the pos that the lasing modes might be localized on can be deduced from the various experimental observables
. it should be noted , however , that not all lasing modes are in fact localized on specific classical trajectories .
therefore , a careful analysis of all available data is needed to determine the nature of the observed resonant states .
the lasing thresholds are related to the lifetime of the cavity modes and hence to the losses of a possible underlying po [ cf .  ] .
the far - field distributions are often concentrated around a few specific directions . from these directions
one can infer the possible trajectories within the resonator via snell s law .
the photographs indicate from which parts of the cavities the light is emitted .
the spectra typically exhibit multimode lasing with several tens of resonances .
they are often organized in sequences of equidistant resonances .
if a set of lasing modes is localized on a certain po , their resonance wave numbers are given by @xmath51 where @xmath35 is an integer and @xmath52 a constant phase shift .
the resonance spacing @xmath53 is hence inversely proportional to the optical length @xmath54 of this po , which can be conveniently obtained from the fourier transform ( ft ) of the spectrum that exhibits peaks at @xmath54 and its multiples @xcite .
it can also be deduced from the free spectral range ( fsr ) , @xmath55 , via the relation @xmath56 where @xmath57 is the wavelength of the lasing emission .
the geometric length of the po , @xmath27 , is related to the optical length by @xmath58 , where @xmath59 is the group refractive index .
the latter differs from the effective refractive index since it also takes into account dispersion .
it has a value in the range of @xmath60 to @xmath61 depending on the sample @xcite .
the precise value for each sample can be determined from calibration measurements with ribbon - shaped fabry - prot cavities since they sustain only a single type of po , the well - known bouncing ball or fabry - prot orbits .
seven different triangles with varying degrees of symmetry and different types of classical dynamics were investigated .
we start with the simplest and most symmetric triangles and gradually pass to less symmetric and accordingly more complicated ones .
the first examples are the equilateral and the right isosceles triangle that have both integrable classical dynamics .
the next two triangles are isosceles triangles with top angles @xmath0 and @xmath15 .
they are rational triangles with pseudointegrable classical dynamics . since these four ( pseudo- ) integrable triangles feature short and simple pos it is expected to find lasing modes localized on some of these orbits .
the fifth triangle is an irrational right triangle .
it also exhibits a short and simple po due to its right angle but has ergodic classical dynamics in contrast to the previous examples . since the rational triangles are a dense subset of all triangles , it is interesting to see if the properties of rational and irrational triangle microlasers significantly differ .
the final two triangles are perturbations of the equilateral triangle and the @xmath0 isosceles triangle , respectively .
they are irrational and hence have ergodic dynamics .
they were chosen because no simple pos are known for them except for fagnano s orbit in the quasiequilateral triangle .
in addition , they permit us to study the influence of geometric perturbations on the microlasers properties and in particular the effect of breaking their mirror symmetry .
the equilateral triangle ( et ) is the triangle with the highest degree of symmetry , and microlasers with equilateral triangular shape have been intensely studied @xcite .
the equilateral triangle has integrable classical dynamics , and the corresponding cavity problem with dirichlet or neumann boundary conditions can be solved analytically @xcite .
this problem was already investigated in the context of vibrating membranes by lam in the 19th century @xcite .
there is , however , no analytical solution in the case of the dielectric boundary conditions considered here .
all pos of the et are known , but none of them is confined by total internal reflection .
this is due to the relatively low value of @xmath62 that corresponds to a critical angle of @xmath63 .
the same is true for the other triangular microlasers .
the lack of good confinement required comparably large cavities to provide sufficient gain .
the side length of the et microlaser was @xmath64 m .
figure [ fig : tequispectrum](a ) shows the lasing spectrum @xmath65 of the et microlaser for a pump energy just above the threshold in the direction @xmath66 [ see ( a ) for the definition of the azimuthal angle ] .
the spectrum exhibits a clear structure of equidistant peaks .
the fourier transform of the spectrum , @xmath67 , is plotted with respect to the optical length @xmath54 in ( b ) .
it features several equidistant peaks with decreasing amplitude as expected for a series of equidistant resonances .
the first peak at @xmath68 m corresponds to the fsr , @xmath69  nm , of the lasing spectrum , and the further peaks are harmonics .
the two shortest types of pos in the et billiard are fagnano s orbit with @xmath70 m [ see ( a ) ] and the so - called quasi - fabry - prot ( qfp ) orbits shown in ( b ) .
a qfp orbit has two reflections with perpendicular incidence and two reflections with an angle of incidence of @xmath71 with respect to the surface normal .
it has a length of @xmath72 m .
the family of the qfp orbits in the et is constructed as follows : first , one qfp orbit [ e.g. , the solid red line in ( b ) ] is shifted perpendicularly to its trajectory ( yielding , for example , the dashed blue lines ) and , second , the qfp orbits can be rotated by @xmath73 ( yielding , for example , the dotted green lines ) .
the family of qfp orbits covers the whole surface of the et .
the calculated optical lengths corresponding to fagnano s orbit and the qfp orbit are @xmath74 m and @xmath75 m , respectively .
they are indicated by the arrows in ( b ) .
obviously , the observed optical length corresponds to the qfp orbit while fagnano s orbit is too short .
all other pos are significantly too long .
further evidence is gained from the azimuthal far - field distribution in ( c ) .
the lasing emission is concentrated in the three directions perpendicular to the cavity side walls .
this is the behavior expected from modes localized on the qfp orbit as shown in ( b ) .
note that there is no emission from the reflections with angle of incidence @xmath71 since this angle is larger than the critical angle . on the other hand
, fagnano s orbit would correspond to six emission directions with an angle of @xmath76 with respect to the surface normals as indicated in ( a ) , but no emission was found in these directions .
it should be furthermore noted that the lasing threshold of the et microlaser is only about @xmath77 higher than that of a fabry - prot ( fp ) cavity of corresponding width .
altogether this proves that the observed lasing modes are localized on the qfp orbit .
finally , a photo taken from the direction @xmath66 and presented in shows that the whole side wall of the cavity is lasing .
this was expected since the family of the qfp orbit covers the whole triangle , though also all other po families do so .
photos taken from @xmath78 and @xmath79 show the same behavior , whereas no lasing light was observed with the camera in all other directions . in summary ,
the qfp orbit was identified unambiguously from the experimental data as the orbit supporting the lasing modes .
it is at first surprising that the dominant lasing modes are localized on the qfp orbit and not on fagnano s orbit like in refs .
first , however , the refractive index of the semiconductor materials used in refs .
@xcite was significantly higher so fagnano s orbit was confined by total internal reflection , which is not the case here .
second , the lasing modes that we observed were te polarized , i.e. , their electric field was parallel to the plane of the resonator .
in fact , te polarized modes are favored by the properties of the lasing dye and the pumping scheme that is used here @xcite . since the angle of incidence of fagnano s orbit , @xmath80 , is close to the brewster angle @xmath81 , a te mode localized on this po would suffer from very high losses . from we
calculate @xmath82 as the threshold of fagnano s orbit and @xmath83 for the qfp orbit .
thus , the dominance of the qfp modes can be well explained by taking into account the peculiarities of the organic microlasers .
it should be noted that the modes of the other triangular microlasers considered in the following were all te polarized as well .
the second triangle that was investigated is the right isosceles triangle ( rit ) , which is essentially a square cut in half along a diagonal .
the rit billiard is integrable like the square and equilateral triangle billiards .
both classical and quantum right triangle billiards have been studied and their pos investigated @xcite .
a property well known by opticians is the fact that a corner with a right angle sends a ray back parallel to its initial direction regardless of the angle of incidence . from this
follows directly the existence of a family of pos that impinge perpendicularly on the hypotenuse as indicated by the thick solid red line in .
we call this orbit quasidiamond orbit in analogy to the diamond po in the square billiard , indicated by the dashed red line in .
another important po is indicated by the dash - dotted blue line .
it is reflected perpendicularly at the two short sides of the rit and with an angle of incidence of @xmath84 at the hypotenuse .
since it corresponds to the fabry - prot orbit of the square billiard ( indicated by the dotted blue lines in ) , it is also called the qfp orbit .
these are the two shortest pos of the rit , and both po families cover the whole area of the billiard .
the qfp orbit exists in all isosceles triangles , with the angle of incidence on the long side and the area covered by its family depending on the top angle of the triangle .    two spectra measured in the directions perpendicular to the hypotenuse and one of the short sides , respectively ,
are shown in . both spectra consist of families of equidistant resonances that have , however , a different fsr each .
the largest peak in the ft of the spectrum for @xmath85 is at @xmath86 m .
as expected for this direction , the corresponding po is the quasidiamond orbit with a geometric length of @xmath87 m , where @xmath88 m is the length of the hypotenuse , and a calculated optical length of @xmath89 m .
the ft also features a smaller peak at half this optical length .
it stems from the slight modulation of the resonance amplitudes [ see ( a ) ] , i.e. , the fact that every second resonance has a somewhat smaller amplitude than its neighbors .
it is interesting to note that there is an isolated po along the height of the triangle , indicated as thin solid red line in , that has half the length of the quasidiamond orbit .
the physical origin of the modulation of the resonance amplitudes and whether it is connected to this po remains , however , unclear .
in contrast , the ft of the spectrum at @xmath90 shows a peak at @xmath91 m .
this corresponds to the qfp orbit with a geometric length of @xmath92 m and an optical length of @xmath93 m .
an unexpected find is a small peak at @xmath86 m [ indicated by the black arrow in the inset of ( b ) ] that corresponds to the quasidiamond orbit .
it originates from another family of barely visible resonances in the spectrum .    the far - field distribution shown in features several emission lobes with differing amplitudes .
the strongest emission lobe is in the direction of @xmath85 and is due to the quasidiamond orbit .
it was cut off in the inset of since its amplitude of @xmath94 counts far exceeds that of the other emission lobes .
it should be noted that the emission lobe at @xmath95 is about @xmath36 times larger than that at @xmath84 , whereas the two lobes are expected to have equal amplitudes due to the mirror symmetry of the triangle . a significant asymmetry of the microcavity itself was excluded .
further experiments demonstrated that the ratio between the amplitudes at @xmath96 and @xmath95 depended sensitively on how precisely the cavity was pumped , that is , for example , on the size and position of the pump beam .    according to
, modes localized on the qfp orbit should emit in the directions of @xmath84 and @xmath95 , and indeed the spectra measured in these directions exhibit a fsr corresponding to its optical length .
however , small contributions also of the quasidiamond orbit were found in the spectra at @xmath84 and @xmath95 as demonstrated in ( b ) .
in addition , two small lobes at @xmath97 and @xmath98 were found that are also related to the quasidiamond orbit .
the origin of these can be elucidated by the photos shown in .
the photos taken from @xmath96 show a strong emission from the side wall perpendicular to the camera perspective as predicted for the qfp orbit .
classically , it is expected that the whole side wall emits like in since the qfp orbit family covers the whole triangle . why this is not observed experimentally
remains unclear .
it could also be related to the strong sensitivity of the emission to the pumping conditions .
in addition , a weak emission from the side wall parallel to the camera perspective was observed .
this is best seen in the top right panel of .
a similar grazing emission was also observed at @xmath90 , @xmath99 , and @xmath98 .
it is not expected classically since the quasidiamond orbit is totally reflected at the two short sides .
the same kind of grazing emission is also observed for the diamond orbit modes of square organic microlasers and will be discussed elsewhere @xcite .
the photo taken at @xmath85 finally shows emission from the whole hypotenuse as expected for modes localized on the quasidiamond orbit .
the key characteristics of the rit microlaser can be well explained by simple pos of the corresponding billiard as in the case of the equilateral triangle .
but in contrast to the equilateral triangle , the rit microlaser features two families of modes localized on different pos that coexist .
a calculation of the thresholds according to reveals that they are very close to each other since both pos have the same losses and nearly the same lengths .
this prediction agrees qualitatively with the measured thresholds .
another interesting observation is the grazing emission of the modes localized on the quasidiamond orbit that is not expected from the ray dynamics .
it demonstrates that some properties of the triangular microlasers need a more careful treatment taking into account wave - dynamical effects .
the most general class of triangles with a symmetry are isosceles triangles .
microlasers with an isosceles triangle shape have been investigated , for example , in ref .
the pos of isosceles triangles can be constructed by the unfolding technique @xcite .
the simplest po that exists in all isosceles triangles is the qfp orbit already known from the right isosceles triangle .
another , more complicated example is the double bow - tie orbit shown in .
it exists for top angles between @xmath100 and @xmath101 .
the first of two obtuse isosceles triangles that are discussed here is the one with top angle @xmath102 ( abbreviated it100 in the following ) .
its classical dynamics is pseudointegrable .
the spectrum of the it100 microlaser measured in the direction @xmath103 , i.e. , perpendicular to one of the short side walls , is shown in .
the fsr of the equidistant resonance family corresponds to an optical length of @xmath104 m .
the underlying po is hence the qfp orbit shown in ( a ) that has a length of @xmath105 m , where the length of the long side is @xmath106 m , which corresponds to an optical length of @xmath107 m .
the far - field distribution presented in ( b ) shows four major emission lobes .
their directions , @xmath103 , @xmath108 , @xmath109 , and @xmath110 , are precisely the ones expected classically for the qfp orbit .
there are also some smaller emission lobes at , e.g. , @xmath79 and @xmath111 , hardly visible in ( b ) .
the spectra in these directions also have a fsr corresponding to the qfp orbit .
these directions can not , however , be explained by the ray dynamics of the qfp orbit .
figure [ fig : tiso100photos ] presents photographs of the it100 microlaser taken from two of the major emission directions .
the photographs taken from @xmath112 in ( a ) show that the whole side wall is emitting with the exception of a small part next to the corner with the long side to the right .
this is in fact the part of the boundary that is expected classically to emit , indicated by the thick red ( dark gray ) line in the sketch in the right panel of ( a ) , because the family of the qfp orbit covers only a part of the billiard ( indicated as the gray area ) in contrast to the right isosceles triangle .
the photographs taken from @xmath113 in ( b ) demonstrate that the most intense emission into that direction originates from the middle part of the long side ( to the left in the photo ) , again in good agreement with the classical prediction shown in the right panel .
however , also the smaller side wall ( to the right in the photo ) emits light , though with lesser intensity .
this emission can not be explained classically with the properties of the qfp orbit family . in summary ,
most of the observed lasing characteristics of the it100 microlaser are in very good agreement with the classical predictions for the qfp orbit , while only some details are beyond a simple ray - dynamical analysis as in the case of the right isosceles triangle microlaser .
the second isosceles triangle that was investigated is the one with a top angle of @xmath15 ( abbreviated it110 in the following ) .
it is also pseudointegrable .
the spectrum of the it110 microlaser at @xmath114 is shown in ( a ) .
it shows a single family of resonances , the fsr of which corresponds to @xmath115 m .
the corresponding po is again the qfp orbit with geometric length @xmath116 m and optical length @xmath117 m , where the length of the long side is @xmath106 m .
the spectrum measured at @xmath118 is presented in ( b ) .
its structure is less clean than that of the one at @xmath119 , but its ft ( see inset ) shows clear peaks at @xmath120 m and multiples of this length .
the corresponding fsr of @xmath121  nm is the fsr of the dominant family of modes [ see ( b ) ] .
the spectrum also exhibits a second family of modes with smaller amplitude and the same fsr .
this optical length as well as the emission direction correspond well to the double bow - tie orbit ( shown as inset ) with a geometric length of @xmath122 = 706.4~\mu$]m and an optical length of @xmath123 m .
it should be noted that the measured threshold of the modes localized on the qfp orbit is almost @xmath124 times higher than that of the modes localized on the double bow - tie orbit .
this agrees qualitatively with , which predicts a @xmath125 times higher threshold .
the far - field distribution of the it110 microlaser is presented in .
the six principal emission directions agree very well with those calculated for the double bow - tie orbit that are indicated by the black arrows [ see also inset of ( b ) ] .
the amplitudes of the emission lobes lack , however , the expected symmetry as in the case of the right isosceles triangle microlaser .
no emissions lobes corresponding to the qfp orbit were observed since the microlaser was pumped slightly above the threshold of the double bow - tie orbit modes but well below the threshold of the qfp orbit modes .
the photos shown in were taken with the same pump intensity as for the measurement of the far - field distribution .
the directions of observation correspond to the major emission directions .
all photos show that the most intense part of the laser emission originates from those parts of the side walls that are covered by the family of the double bow - tie orbit [ indicated by the red ( dark gray ) lines in the drawing of the it110 ] .
this confirms that the lasing modes with the lowest threshold are localized on the double bow - tie orbit .
however , the photos also show weak emission from other parts of the side walls that are not covered by the family of the double bow - tie orbit .
a possible explanation is that even though the field distributions of superscarred resonant states are strongly concentrated inside the po channel , they also have a nonvanishing field outside of the po channel due to coupling to nonscarred states @xcite .    in conclusion ,
the it110 microlaser is another example like the right isosceles triangle where two families of modes localized on different pos coexist .
in contrast to the case of the right isosceles triangle , however , there is a significant difference between the thresholds of the two mode families so one of them is easily selected by keeping the pump intensity sufficiently low .
furthermore , it should be noted that the dominant po , the double bow - tie orbit , is selected because it has the lowest losses even though it is longer and more complicated than the qfp orbit .
the double bow - tie orbit also exists in the it100 billiard but has a higher threshold than the qfp orbit in that case .
lasing modes based on the double bow - tie orbit were not found experimentally for the it100 microlasers even at considerably higher pump energies .
this demonstrates that the change of a single geometric parameter , the top angle in this case , can significantly modify the lasing characteristics .
the fifth triangle that was investigated is a pythagorean triangle ( pt ) with side length ratio @xmath126 as shown in .
its hypotenuse is @xmath127 m long .
in contrast to the other triangles studied so far , it is irrational and exhibits no symmetry .
it exhibits a qfp orbit that impinges perpendicularly on the hypotenuse as shown in in analogy to the right isosceles triangle .
the pos of right triangles can also be constructed by unfolding @xcite even if they are irrational like the pt .
the lasing spectrum of the corresponding pt microlaser observed in the direction perpendicular to the hypotenuse is presented in ( a ) .
it exhibits a single family of equidistant resonances .
its ft , shown in ( b ) , features a dominant peak at @xmath128 m that corresponds to the fsr of the spectrum .
the length of the qfp orbit is @xmath129 m , which yields an optical length of @xmath130 m that is in good agreement with the optical length observed in the ft . in complete analogy to the case of the right isosceles triangle , there is an additional , smaller peak at half this optical length that stems from a modulation of the resonance amplitudes , and there is a po along the height of the triangle with half the length of qfp orbit ( indicated as dashed blue line in ) .
it should be noted that the relative amplitude of the peaks at @xmath131 m and @xmath132 m depends sensitively on the pump beam position and other details of the excitation scheme .
the far - field distribution presented in has three major emission directions , @xmath133 , @xmath134 , and @xmath135 .
the lasing spectra in these directions exhibit the same modes as the one shown in ( a ) .
the emission directions predicted for the qfp orbit according to snell s law for @xmath62 are @xmath133 , @xmath136 , and @xmath137 , respectively ( see ) .
they are indicated in by the black arrows and agree quite well with the observed ones .
it should be noted that no emission is expected from the left ( smallest ) side of the pt since the angle of incidence of the qfp orbit on it is larger than the critical angle .
the amplitude of the emission lobe at @xmath135 is considerably smaller than that of the lobe at @xmath134 while an equal amplitude is expected classically .
in fact , the relative amplitude of the lobes strongly depended on the position and size of the pump beam as in the case of the right isosceles triangle microlaser .    furthermore , the directions of these two lobes are not symmetric as well . according to the geometry of the qfp orbit , their angles with respect to the surface
normal should be equal , but in reality they are @xmath138 and @xmath139 , respectively , and thus differ by @xmath140 .
this is a significant deviation that is within the resolution of the setup .
the third emission lobe on the contrary has precisely the expected direction perpendicular to the hypotenuse .
in contrast to the previously considered triangles , however , the pt cavity itself exhibits no symmetry , and symmetric emission directions are only expected due to the properties of the underlying qfp orbit .    a completely unexpected experimental result are the three smaller and broader emission lobes around @xmath141 , @xmath142 , and @xmath143 .
the spectra in these directions feature the same structure and fsr as those in the three major emission lobes ; however , these three emission directions can not be related to the dynamics of the qfp orbit .
photographs of the pt microlaser taken from the three major emission directions are shown in the left panels of .
the sketches in the right panels indicate the parts of the side walls that are expected to emit in these directions according to the geometry of the qfp orbits . indeed
, the brightest areas of emission in the photos correspond well to these classical predictions .
in contrast , also some weak emissions are observed from the side wall on the left side at @xmath144 and from the side wall at the right side at @xmath145 , which are not expected classically .    in summary , the spectra , the far - field distribution , and the photos clearly evidence that the observed lasing modes are localized on the qfp orbit family . in detail , however , there are deviations from the ray - dynamical predictions . first , the directions of the two emission lobes at @xmath134 and @xmath135 do not exhibit the expected symmetry ; second , there are three additional emission lobes the directions of which have no apparent connection to the qfp orbits ; and , third , parts of the lasing emission stem from sections of the cavity boundary that are not expected to emit classically .
so even though the underlying po could be identified , the properties of the pt triangle microlaser can not be explained with the same precision and completeness as in the previous cases .
this might be related to the fact that the pt billiard is neither rational nor symmetric , and in fact its pos can be constructed easily only due to its right angle .      while all the triangles considered so far had symmetries or other properties that enabled an easy construction of their pos
, the case considered in the following has none of these .
it is a deformation of the equilateral triangle and hence called quasiequilateral triangle ( qet ) in the following .
its side lengths are @xmath146 m , @xmath147 m , and @xmath148 m , respectively , as indicated in .
it is an irrational triangle and we hence can not easily construct any po besides fagnano s orbit .
in particular , the qfp orbits no longer exist due to the lack of symmetry .
a typical spectrum , observed at @xmath149 , is shown in ( a ) .
the structure of the spectrum is not as clear as for the previously shown ones ; nonetheless , a small sequence of equidistant resonances can be identified .
the ft of the spectrum shows a peak at @xmath150 m and its approximate multiples .
the corresponding fsr of @xmath69  nm matches that of the resonance sequence in the spectrum .
this optical length is close to that of the qfp orbit in an equilateral triangle with side length @xmath41 , but the qfp orbit no longer exists in the qet .
there are , however , three diffractive pos along the heights of the triangle that have nearly the same lengths .
they are depicted in and called @xmath151 , @xmath152 , and @xmath153 in the following .
they are called diffractive orbits because one of their vertices is at a diffractive corner of the billiard ( cf .  ) .
it should be noted that for an equilateral triangle , the orbits along the heights belong to the family of the qfp orbit .
thus , these three diffractive pos can be considered as the remnants of the qfp orbit family .
their lengths are @xmath154 m , @xmath155 m , and @xmath156 m .
the corresponding optical lengths are indicated by arrows in ( b ) and are indeed close to the peaks observed in the ft of the spectrum .
the far - field distribution of the qet microlaser is plotted in ( a ) .
it features three broad bundles of emission lobes the centers of which are roughly perpendicular to the cavity side walls .
the ft of the spectra measured at different azimuthal angles is shown in ( b ) .
the dominant optical lengths vary somewhat with @xmath157 , but they always stay close to the lengths of the three diffractive pos , the optical lengths of which are indicated by the vertical white lines .
this corresponds to the fact that the spectra always feature a similar fsr even though their structure and quality varies significantly , leading to a relatively high noise level in the ft .
no evidence of fagnano s orbit was found like in the case of the equilateral triangle since its angles of incidence are close to the brewster angle .
so while the fsrs approximately match those corresponding to the diffractive orbits , other observations do not indicate that the modes are localized on them .
for example , the directions of maximal emission are not exactly perpendicular to the side walls as one would naively expect .    the photos shown in enable a better understanding of the nature of the resonant modes .
they show the lasing qet from the directions approximately perpendicular to the side walls .
all of them show that the points of origin of the lasing emission are more or less broadly distributed over the side walls .
the expectation for modes localized on the diffractive pos , in contrast , would be that the origin of emission is strongly concentrated around the feet of the heights that are roughly in the middle of the side walls .
another expectation for this kind of mode would be strong emission coming from the corners of the triangle .
this is , however , not observed in any direction .
in fact , the images shown in are typical also for other directions in which the qet microlaser emits .    in conclusion , the lasing characteristics of the qet microlaser can not be explained by any po of the qet billiard . the observed fsrs correspond to an optical length similar to that of the qfp orbits of an equilateral triangle .
a possible explanation is that the observed modes are the perturbed modes of the equilateral triangle microlaser that were localized on the qfp orbit . due to the small perturbation ,
their fsr stays approximately the same , but the far - field distribution broadens around the emission directions of the qfp orbit that are perpendicular to the side walls .
the same effect was observed in ref .
@xcite for triangular resonators with dirichlet boundary conditions .
some modes were shown to be localized on so - called ghost pos , i.e. , the pos of a geometrically different , but similar , triangle .
an analogous case is the persistence of the influence of the bouncing ball orbits in a quantum stadium billiard when the originally parallel side walls of the stadium are slightly tilted @xcite .
the reason for these effects is that while the classical dynamics can exhibit singular behavior with respect to perturbations of the billiard geometry , e.g. , pos suddenly vanishing completely , wave - dynamical systems react in a continuous manner to geometric perturbations , essentially smoothing out the singularities of classical mechanics . the data presented here
lead us to believe that the modes of the quasiequilateral triangle microlasers are localized on the ghost qfp orbit , but numerical investigations of the wave functions will be necessary to confirm this notion .
the last triangle is again a deformation of one of the previous triangles .
the triangle was constructed by moving the right vertex of the isosceles @xmath0 triangle by @xmath158 m to the left while keeping the other two vertices fixed as demonstrated in ( a ) .
hence it is called quasi - isosceles triangle ( qit ) .
it is irrational and therefore no simple pos are known since the qfp and the double bow - tie orbit of the isosceles @xmath0 triangle are destroyed by the geometric perturbation .
there are , however , three diffractive pos [ see ( b ) ] similar to the case of the quasiequilateral triangle .
one is along the height of the triangle ( called the height orbit in the following ) , and the other two have one perpendicular reflection at a short side wall and one at the top vertex .
they can be considered as the remnants of the qfp orbits of the isosceles @xmath0 triangle and are called diffractive quasi - fabry - prot orbits in the following .
two spectra of the qit microlaser observed in different directions are plotted in .
both spectra show a relatively clear structure of equidistant modes .
however , their fsrs and the corresponding optical lengths differ as can be seen in the fts shown as insets .
the optical length for the spectrum at @xmath159 is @xmath160 m .
this can only correspond to the height orbit with @xmath161 m and @xmath162 m .
the optical length for the spectrum at @xmath163 is @xmath164 m and corresponds approximately to the lengths of the two diffractive qfp orbits .
the geometric and optical length of the left diffractive qfp orbit [ solid red line in ( b ) ] are @xmath165 m and @xmath166 m , respectively , and are identical to those of the qfp orbit in the isosceles @xmath0 triangle , while the geometric and optical length of the right diffractive qfp orbit [ dashed green line in ( b ) ] are @xmath167 m and @xmath168 m , respectively . so the qit microlaser exhibits ( at least ) two different families of modes that have similar thresholds but different fsrs .
next , we investigated what family of modes emitted in which directions .
the far - field distribution in shows a large number of emission lobes with varying amplitudes , though the emission is not as broadly distributed as in the case of the quasiequilateral triangle .
the directions in which the spectra exhibit the optical length of the diffractive qfp orbits and the height orbit are indicated by the solid red and dotted blue arrows , respectively . note that for some of the smaller emission lobes the corresponding optical length could not be clearly determined due to indistinct spectra .
in fact , some of the most prominent emission lobes point approximately in the directions expected for a diffractive qfp orbit which are indicated by the black double arrows [ see also ( b ) ] .
there remain , however , significant deviations between the emission directions of the presumed diffractive qfp modes and the directions predicted by ray optics . these deviations can not be consistently explained by a different refractive index either . regarding the presumed height orbit modes ,
there are several emission lobes around , though not precisely in , the direction of @xmath169 predicted classically .
furthermore , the two emission lobes close to @xmath100 seem reasonable for such modes , too .
the strong emission in the direction of @xmath170 from the presumed height orbit modes and in the direction of @xmath171 from both families of modes , however , defy any simple ray - dynamical explanation .
thus , the dynamics of the diffractive pos can explain the observed emission directions at best on a qualitative level .
the photos taken in the main emission directions of the presumed diffractive qfp modes are presented in .
they show that the emission originates from almost the whole side walls and not just small portions of it as would be expected for the isolated diffractive qfp orbits .
the photos rather resemble those of the equilateral and isosceles triangle microlasers .
in addition , a part of the emission originates from those side walls that are not directly facing the camera .
emission from these side walls is not at all expected from any fabry - prot - like modes .
the photos in were analogously taken in the main emission directions of the presumed height orbit modes .
the lasing emission is not as broadly distributed along the side walls as in the cases presented in but shows nonetheless no single points of concentration . in particular
, the photos at @xmath172 and @xmath173 demonstrate that the emission in these directions is not predominantly originating from the vertex of the triangle as expected for a mode localized on the diffractive height orbit .    in summary , the photos evidence that both families of lasing modes cover more or less the whole microlaser and are not strongly concentrated along the isolated diffractive orbits shown in .
we are hence led to believe that one family of modes is localized on the ghost qfp orbits , i.e. , they are the perturbed qfp modes of the isosceles @xmath0 triangle , analogously to the case of the quasiequilateral triangle .
the situation for the other family of modes with an fsr corresponding to the height orbit is less clear since a ( nondiffractive ) height po does not exist in the isosceles @xmath0 triangle , and no family of modes with a similar fsr was observed for the isosceles @xmath0 triangle microlaser . in any case ,
further investigations are necessary for a full understanding of the qit microlaser s behavior .
it is interesting to note that even though a po is known to exist in the qit since its angles are all less than @xmath0 , the dominant lasing modes are clearly not localized on any classical po .
it can be presumed that even the shortest po in the qit is too long or has too high losses to support lasing modes with a reasonably low threshold and that therefore other types of lasing modes are predominant .
while there are still open questions regarding the existence and properties of pos in classical triangular billiards and this remains a domain of active research , it is also very interesting to explore the influence of the pos on the properties of the corresponding wave - dynamical billiards . of particular interest
are billiards for which the pos can not be easily constructed or are not even known because they are irrational and lack symmetries .
one of the objectives of the experiments presented here was therefore to see whether the resonant states in these cases are localized on pos , and , if yes , on which ones or if , on the contrary , they have no relation to specific pos .
we investigated this ray - wave correspondence in experiments with organic microlasers of triangular shape .
the shapes that were chosen correspond to different types of classical dynamics featuring diverse types of pos and were hence expected to exhibit very different lasing characteristics . even though each triangle microlaser had properties distinct from the others and had to be treated separately , it is possible to draw some general conclusions from the experiments .
several examples of triangles with well - understood classical dynamics and relatively simple pos were studied , namely the equilateral triangle and several rational isosceles triangles .
they all exhibit lasing modes that are clearly localized on pos .
the most important observables were the fsr of the spectra and the directions of the emission lobes in the far field that could be measured with high precision and showed excellent agreement with ray - optical predictions .
in contrast , the amplitudes of the emission lobes proved to be very sensitive to small changes in the experimental setup .
photographs of the lasing cavities allowed further insight into the characteristics of the lasing modes and showed good qualitative agreement with the classical predictions .
it seems to be a general rule for polygonal billiards having simple , short pos with not too high losses that the modes of the corresponding microlasers are localized on these pos @xcite .
nonetheless , even for these simple and symmetric triangles some details of the observed spectra and far - field distributions were beyond simple ray - dynamical explanations .
this underlines the need to further refine the different models of modal localization in dielectric resonators to better account for wave - dynamical effects like diffraction .
the agreement with the ray - optical calculations was less good for the case of an irrational right triangle , namely the pythagorean triangle . due to its right angle
, it still features a simple po , and the fsr of the spectra corresponded to its optical length . the directions of the emission lobes , however , showed deviations from the classical expectations , and several emission lobes in addition to the predicted ones were found . also the photographs of the pt microlaser
only partly agreed with the classical expectations .
finally , two examples of irrational triangles lacking any symmetry and hence any simple pos ( with the exception of fagnano s orbit for the quasiequilateral triangle ) were investigated .
both cases were deformations of previously investigated triangles .
their modes seemed to be localized on ghost pos @xcite , i.e. , they resembled the modes localized on pos of the undeformed triangles .
but even though their modes retained some of the features of the unperturbed triangles , their properties like the major emission directions could only be explained qualitatively , if at all .
this means that , contrary to expectations , the lasing modes of triangular microlasers are not based on a po if no simple po exists in the corresponding billiard .
it also demonstrates that there seems to be a significant , qualitative difference between microlasers with the shape of rational and irrational triangles which is surprising since any irrational triangle can be approximated by a rational one .
the examples of irrational triangles demonstrate that the predictive power of ray optics decreases for irrational triangles like the pythagorean triangle , and quantitative predictions are no longer possible for irrational triangles without symmetries and simple pos .
there is apparently a transition from symmetric triangles with simple classical dynamics for which the lasing properties can be very well explained by ray optics to less symmetric triangles with increasingly complicated classical dynamics for which ray - optical methods no longer yield a good description of the lasing properties .
so the complexity of the classical dynamics is directly reflected in the complexity of the corresponding microlasers .
it remains an interesting challenge to better understand the latter cases by further experimental and numerical studies .
the authors thank e.  bogomolny and a.  nosich for illuminating and fruitful discussions .
s.  b.  gratefully acknowledges funding from the european union seventh framework programme ( fp7/2007 - 2013 ) under grant
this work was supported by a public grant from the laboratoire dexcellence physics atom light matter ( labex palm ) overseen by the french national research agency ( anr ) as part of the investissements davenir program ( reference no .
anr-10-labx-0039 ) .
81ifxundefined [ 1 ] ifx#1 ifnum [ 1 ] # 1firstoftwo secondoftwo ifx [ 1 ] # 1firstoftwo secondoftwo `` `` # 1''''@noop [ 0]secondoftwosanitize@url [ 0 ]
 + 12$12  & 12#1212_12%12@startlink[1]@endlink[0]@bib@innerbibempty link:\doibase 10.1007/bf01388890 [ * * , ( ) ] link:\doibase 10.1103/physreve.52.2066 [ * * ,   ( ) ] @noop *
* ,   ( ) link:\doibase 10.1090/s0002 - 9947 - 98 - 02089 - 3 [ * * ,   ( ) ] link:\doibase 10.1007/s000140050113 [ * * ,   ( ) ] link:\doibase 10.1080/10586458.2006.10128961 [ * * ,   ( ) ] link:\doibase 10.1007/s10711 - 007 - 9129 - 9 [ * * ,   ( ) ] @noop _ _  ( ,  ,  ) @noop _ _  ( ,  ,  ) link:\doibase 10.1103/physrevlett.67.785 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.64.2215 [ * * ,   ( ) ] in  @noop _ _ ,  , vol .  ,
( ,  ,  )  pp .
link:\doibase 10.1103/physrevlett.83.2171 [ * * ,   ( ) ] @noop `` , ''   ( ) ,  ed . , @noop _ _  ( ,  ,  ) link:\doibase 10.1364/oe.10.000752 [ * * ,   ( ) ] link:\doibase 10.1364/ol.27.000824 [ * * ,   ( ) ] link:\doibase    10.1103/physreve.67.015207 [ * * ,   ( ) ]
link:\doibase 10.1103/physreva.76.023830 [ * * ,   ( ) ] link:\doibase 10.1103/physreva.62.013816 [ * * ,   ( ) ] link:\doibase    10.1364/ao.46.002969 [ * * ,   ( ) ] link:\doibase 10.1103/physreve.83.036208 [ * * ,   ( ) ] link:\doibase 10.1007/s003400000210 [ * * ,   ( ) ] link:\doibase 10.1364/ol.29.000005 [ * * ,   ( ) ] link:\doibase    10.1109/jstqe.2006.883150 [ * * ,   ( ) ] \doibase doi:10.1364/oe.16.016516 [ * * ,   ( ) ] link:\doibase 10.1364/ol.33.000509 [ * * ,   ( ) ] link:\doibase    10.1364/ol.34.001810 [ * * ,   ( ) ] link:\doibase    10.1364/ol.35.002723 [ * * ,   ( ) ] link:\doibase    10.1103/physreve.83.016208 [ * * ,   ( ) ] link:\doibase 10.1103/physreva.72.063806 [ * * ,   ( ) ] link:\doibase 10.1063/1.2779094 [ * * ,   ( ) ] link:\doibase    10.1063/1.3276069 [ * * ,   ( ) ] link:\doibase    10.1364/ol.29.002758 [ * * ,   ( ) ] link:\doibase 10.1140/epjp / i2012 - 12008 - 0 [ * * ,   ( ) ] link:\doibase 10.1109/aps.1991.174903 [ * * ,   ( ) ] link:\doibase 10.1364/josaa.24.002822 [ * * ,   ( ) ] link:\doibase 10.1109/tap.2009.2038181 [ * * , ( ) ] link:\doibase 10.1016/0167 - 2789(81)90024 - 5 [ * * ,   ( ) ] link:\doibase 10.2307/1971280 [ * * ,   ( ) ] @noop _ _  ( ,  ,  ) link:\doibase 10.1090/s0002 - 9939 - 2012 - 11379 - 6 [ * * ,   ( ) ] link:\doibase 10.1063/1.532762 [ * * , ( ) ] @noop _ _  ( ,  ,  ) @noop _ _  ( , ,  ) link:\doibase 10.1103/physreva.77.013804 [ * * ,   ( ) ] link:\doibase 10.1063/1.1665328 [ * * , ( ) ] link:\doibase 10.1063/1.1665596 [ * * , ( ) ] link:\doibase 10.1103/physrevlett.53.1515 [ * * ,   ( ) ] @noop `` , ''   ( ) ,   link:\doibase 10.1103/physrevlett.97.254102 [ * * ,   ( ) ] link:\doibase 10.1103/physreve.78.056202 [ * * ,   ( ) ] link:\doibase 10.1103/physreve.86.026202 [ * * ,   ( ) ] link:\doibase    10.1103/physreve.81.066215 [ * * ,   ( ) ] link:\doibase 10.1088/1751 - 8113/44/15/155305 [ * * ,   ( ) ] link:\doibase 10.1364/ol.36.001023 [ * * ,   ( ) ] link:\doibase    10.1103/physreva.88.023834 [ * * ,   ( ) ] link:\doibase    10.1103/physreve.88.062906 [ * * ,   ( ) ] link:\doibase 10.1088/0951 - 7715/16/6/310 [ * * ,   ( ) ] link:\doibase 10.1109/tap.2010.2103031 [ * * ,   ( ) ] link:\doibase 10.1016/0167 - 2789(88)90075 - 9 [ * * ,   ( ) ] link:\doibase 10.1109/jqe.2010.2055836 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.109.033903 [ * * ,   ( ) ] link:\doibase 10.1103/physreva.86.043817 [ * * ,   ( ) ] link:\doibase 10.1088/0951 - 7715/9/4/012 [ * * , ( ) ] link:\doibase 10.1103/physreve.61.3689 [ * * ,   ( ) ] link:\doibase 10.1103/physreve.50.1886 [ * * ,   ( ) ] link:\doibase 10.1088/0305 - 4470/27/13/018 [ * * ,   ( ) ] link:\doibase 10.1063/1.4720474 [ * * ,   ( ) ] link:\doibase 10.1364/oe.22.012316 [ * * ,   ( ) ] link:\doibase 10.1103/physreve.68.066207 [ * * , ( ) ] link:\doibase 10.1088/1464 - 4258/7/9/009 [ * * ,   ( ) ] link:\doibase 10.1364/josab.23.001586 [ * * ,   ( ) ] link:\doibase 10.1103/physreva.76.023822 [ * * , ( ) ] link:\doibase 10.1109/jqe.2009.2024006 [ * * , ( ) ] link:\doibase 10.1364/josab.26.002449 [ * * ,   ( ) ] @noop _ _ ( ,  ,  ) link:\doibase 10.1088/0305 - 4470/34/40/306 [ * * ,   ( ) ] @noop link:\doibase 10.1103/physrevlett.100.204101 [ * * ,   ( ) ] link:\doibase    10.1103/physreve.78.045201 [ * * ,   ( ) ]