Laser Transverse Modes with Ray-Wave Duality: A Review

: We present a systematic overview on laser transverse modes with ray-wave duality. We start from the spectrum of eigenfrequencies in ideal spherical cavities to display the critical role of degeneracy for unifying the Hermite–Gaussian eigenmodes and planar geometric modes. We subsequently review the wave representation for the elliptical modes that generally carry the orbital angular momentum. Next, we manifest the ﬁne structures of eigenfrequencies in a spherical cavity with astigmatism to derive the wave-packet representation for Lissajous geometric modes. Finally, the damping effect on the formation of transverse modes is generally reviewed. The present overview is believed to provide important insights into the ray-wave correspondence in mesoscopic optics and laser physics.


Introduction
In paraxial approximation, the wave equation for spherical resonators was verified to be analogous to the Schrödinger equation for two-dimensional (2D) harmonic oscillators [1]. This analogy has been fruitfully employed to visualize the quantum wave functions from generating various laser transverse modes in spherical cavities [2]. The eigenmodes of the spherical cavities are analytically found to be Hermite-Gaussian (HG) modes in rectangular coordinates and Laguerre-Gaussian (LG) modes in polar coordinates [3]. The development of diode-pumped solid-state lasers enables us to experimentally generate both high-order HG and LG modes with the selectively focused excitation [4][5][6][7][8][9][10][11]. Besides the HG and LG eigenmodes, the so-called geometric modes with ray-wave duality can be systematically generated in the degenerate cavities [12][13][14][15][16]. The observations of eigenmodes or geometric modes are confirmed to be dependent on the transverse order N and the frequency ratio ∆ f T /∆ f L , where ∆ f L is the longitudinal mode spacing and ∆ f T is the transverse mode spacing. The laser resonator that satisfies the condition of ∆ f T /∆ f L = P/Q is called the degenerate cavity, where P and Q are coprime positive integers. The optical ray tracing in a degenerate cavity with ∆ f T /∆ f L = P/Q can be found to be back to the initial point after the Q round tracings. To be brief, a periodic orbit comprises the ray tracing with the Q round trips in a degenerate cavity with ∆ f T /∆ f L = P/Q.
Experimental observations displayed that if the transverse order N was less than the number Q, the observed modes were often dominated by the HG or LG eigenmodes [4,5]. However, the lasing modes were found to be governed by the geometric modes under the circumstance of N >> Q [13][14][15][16]. One representation for the geometric modes was theoretically derived from the inhomogeneous wave equation [14]. The other representation proposed for expressing the geometric modes was the Fox-Li method by using an initial field matching the gain distribution [17]. Nevertheless, both methods cannot simultaneously represent the eigenmodes and geometrical modes. Recently, we have successfully exploited the formulism of quantum harmonic oscillators to develop an analytical theory to unify the representation for the eigenmodes and geometric modes [18].
The marriage between the resonant modes and the eigenfrequency distributions is of interest for exploring novel mesoscopic optics [19][20][21][22]. Due to the analogy between quantum mechanics and optics, large-area vertical-cavity surface emitting lasers have been developed to explore the subtle relationship between level degeneracies and periodic orbits in mesoscopic quantum billiards [23][24][25][26][27]. Furthermore, two dimensional (2D) photonic systems analogous to quantum Hall [28] and quantum spin Hall [29] systems have been proposed and experimentally demonstrated in optical regime. Quantum level degeneracies have been found to link the conductance fluctuation to the prevalence of classical periodic orbits [30]. Analogous to mesoscopic quantum phenomena, the emergence of geometric modes has been evidenced to be related to the degeneracies of eigenfrequencies in spherical cavities [31][32][33][34]. Laser geometric modes have been paid much attention because of various applications including ray-wave connection, light-matter interaction, optical communication, optical microscopy, and material processing [13,14,17,[35][36][37][38].
One of the most striking geometric modes is the transverse intensities exhibiting the Lissajous figures [31][32][33][34][39][40][41]. In 1964, Herriot et al. [42] proved that under foldback conditions, an incident laser beam can be reflected forth and back between two spherical mirrors to exhibit an elliptical pattern. Herriott and Schulte further [43] proved that the use of a pair of astigmatism lenses can lead to elliptical precession, forming a Lissajous dot pattern under reentrant conditions. Herriott-type multiple-pass cells have been used in various experiments, such as optical delay lines [44], absorption spectroscopy [45], Raman conversion [46][47][48], and high-power laser systems [49]. Moreover, an analysis was developed by McManus et al. [50] to illustrate the Lissajous dot pattern on a mirror [51].
We previously used the SU(2) quantum theory to derive the Lissajous geometric mode to be a superposition of different order HG modes combined with different longitudinalorder modes subject to the degeneracy [39,41]. Recently, we have exploited the generalized wave-packet state representation to construct the ray-wave connection for the multiple-pass Lissajous dotted patterns [52]. Another important aspect of the wave-packet representation is to explore the damping effect that causes the asymmetric structure of the laser mode [53]. Asymmetric laser modes [54] have been employed in manipulating biological cells [55] and generating entangled states of quantum communication [56].
In this article, we systematically summarize the laser transverse modes with raywave duality in spherical resonators with and without astigmatism. We start from the eigenfrequency spectrum in an ideal spherical cavity to comprehend the important role of degeneracy in the formation of geometric modes. We then review the Gaussian wavepacket representation for unifying the HG eigenmodes and the planar geometric modes. We further extend the wave-packet representations to the case related to elliptical modes carrying orbital angular momentum. Next, we present an overview the fine structure of the characteristic frequency in a spherical cavity subject to astigmatism to obtain the parametric ray equations for Lissajous geometric modes. Based on the parametric ray equations, we clearly review the modification of the wave packet representation for considering the astigmatism effect. Finally, we review the damping effect on the generation of asymmetric geometric mode from the Gaussian wave-packet representation. This review is believed to provide important insights into the ray-wave duality in mesoscopic optics and laser physics. For example, the astigmatism of the laser cavity can lead the eigenfrequency spectrum to be topologically identical to the 2D commensurate harmonic oscillators. Consequently, the astigmatic degenerate cavities can be exploited to generate the geometric modes with transverse patterns to manifest the quantum wave functions.

Eigenfrequency Spectrum
The round-trip propagation matrix for a stable spherical resonator is generalized as with det(M) = 1. Two eigenvalues of the matrix M can be found to be Conveniently, two eigenvalues are expressed as λ ± = e ± i θ , where From Sylvester's theorem, M Q with a positive integer Q can be shown as [57] If the parameter θ/2π is equal to a rational number P/Q with P and Q to be coprime positive integers, Equation (4) can be used to confirm that M Q = I, where I is the unit matrix. The feature of M Q = I means that an optical ray inside the cavity is back to the initial position after the Q round trips. The fraction θ/2π represents the ratio of the phase variations between the transverse and longitudinal parts. Therefore, the fraction P/Q is just the mode-spacing ratio ∆ f T /∆ f L , where ∆ f T and ∆ f L are the frequency spacings for the transverse and longitudinal modes, respectively. Note that the value of Q is infinite for the ratio ∆ f T /∆ f L to be an irrational number. From Equation (3), the mode-spacing ratio can be explicitly given by Considering a concave-flat cavity with a mirror with radius of curvature R c at z = −L c and a flat output coupler at z = 0, the round-trip propagation matrix starting from z = 0 is Using Equations (5) and (6), the frequency ratio ∆ f T /∆ f L can be obtained as Consequently, the eigenfrequencies of the cavity can be expressed as where n = (n 1 , n 2 , n 3 ), the integers n 1 and n 2 are the transverse mode indices, and the integer n 3 is the longitudinal mode index. Let n c be the index for the central mode, the relative eigenfrequency distribution can be expressed as where δn = n − n c = (δn 1 , δn 2 , δn 3 ). The connection between geometric modes and eigenfrequency degeneracies can be manifested from the relative spectrum F(δn; L c ). Figure 1 shows the spectrum F(δn; L c ) as a function of L c /R c for |δn i | ≤ 10 with i = 1, 2, 3. Numerous degeneracies and gaps with ∆ f T /∆ f L = P/Q can be seen to appear at the effective cavity lengths of L c = L P/Q = R c sin 2 (πP/Q).

Parametric Equations for Geometric Rays
By using the propagation matrix for a degenerate cavity with , the characteristic equations for periodic ray trajectories can be verified as [18] ( ) is the running index for the different rays, In the spherical cavity, the properties for the x-and y-direction are the same. Therefore, the trajectory equation for the y-direction can be expressed as

Parametric Equations for Geometric Rays
By using the propagation matrix for a degenerate cavity with ∆ f T /∆ f L = P/Q, the characteristic equations for periodic ray trajectories can be verified as [18] x where θ s = (P/Q)2πs, s = 0, 1, 2 · · · Q − 1 is the running index for the different rays, is the Rayleigh range, W is the amplitude of the transverse waist, φ x is associated with the initial condition, and − and + in the symbol of ∓ represent the leftward and rightward tracings, respectively. Let the amplitude W = √ N x w o , Equation (10) can be expressed as where w(z) = w 0 1 + (z/z R ) 2 is the Gaussian beam parameter with the waist w o , θ G (z) = tan −1 (z/z R ) is the Gouy phase, and N x is the transverse order. Defining the dimensionless variable x = √ 2x/w(z), the ray equation in Equation (11) can be given by In the spherical cavity, the properties for the xand y-direction are the same. Therefore, the trajectory equation for the y-direction can be expressed as where y = √ 2y/w(z), φ y is associated with the initial position, and N y is the transverse order. If one of the N x and N y is zero, the ray traces form a planar periodic orbit. For convenience, we set N x = 0 to discuss the feature of planar periodic orbits. Figure 2a shows the planar periodic orbits in the yz plane for P/Q = 1/4, 2/7, and 3/11 with φ y = 0. The planar periodic orbits can be seen to be the reciprocating paths. It can be confirmed that if the phase factor φ y is the integer multiple of π/Q, the planar periodic tracing is a reciprocating type. On the other hand, when φ y is not the integer multiple of π/Q, the planar tracing is a generic type, as an example seen in Figure 2b for φ y = π/(2Q). The total numbers of rays for the reciprocating and generic periodic orbits can be easily found to be Q and 2Q, respectively. In addition to planar orbits, when both of N x and N y are not zero and φ x = φ y , the ray orbits are generally nonplanar, as shown in Figure 3 for P/Q = 1/4, 2/7, and 3/11 with N x = N y and φ x − φ y = π/2. The total number of rays in the nonplanar orbits can be seen to be Q.
. The planar periodic orbits can be seen to be the reciprocating paths. It can be confirmed that if the phase factor y φ is the integer multiple of π/Q, the planar periodic tracing is a reciprocating type. On the other hand, when y φ is not the integer multiple of π/Q, the planar tracing is a generic type, as an example seen in Figure 2b for The total numbers of rays for the reciprocating and generic periodic orbits can be easily found to be Q and 2Q, respectively. In addition to planar orbits, when both of x N and y N are not zero and x y φ φ ≠ , the ray orbits are generally nonplanar, as shown in Figure 3 for / P Q = 1/4, 2/7, and 3/11 with x y N N = and / 2 x y φ φ π − = . The total number of rays in the nonplanar orbits can be seen to be Q.

Wave Representation for Unifying the Eigenmodes and Geometric Modes
In this section, we review that the Gaussian wave-packet state can be developed from the coherent state of quantum harmonic oscillator to unify the representation for the HG eigenmodes and geometric modes. From quantum theory, the coherent state for the onedimensional (1D) harmonic oscillator is given by [58]  y . The planar periodic orbits can be seen to be the reciprocating paths. It can be confirmed that if the phase factor y φ is the integer multiple of π/Q, the planar periodic tracing is a reciprocating type. On the other hand, when y φ is not the integer multiple of π/Q, the planar tracing is a generic type, as an example seen in Figure 2b for The total numbers of rays for the reciprocating and generic periodic orbits can be easily found to be Q and 2Q, respectively. In addition to planar orbits, when both of x N and y N are not zero and x y φ φ ≠ , the ray orbits are generally nonplanar, as shown in Figure 3 for / P Q = 1/4, 2/7, and 3/11 with x The total number of rays in the nonplanar orbits can be seen to be Q.

Wave Representation for Unifying the Eigenmodes and Geometric Modes
In this section, we review that the Gaussian wave-packet state can be developed from the coherent state of quantum harmonic oscillator to unify the representation for the HG eigenmodes and geometric modes. From quantum theory, the coherent state for the onedimensional (1D) harmonic oscillator is given by [58]

Wave Representation for Unifying the Eigenmodes and Geometric Modes
In this section, we review that the Gaussian wave-packet state can be developed from the coherent state of quantum harmonic oscillator to unify the representation for the HG eigenmodes and geometric modes. From quantum theory, the coherent state for the one-dimensional (1D) harmonic oscillator is given by [58] g( where ψ n (x) is the HG eigenstates given by H n (·) are the Hermite polynomials of order n, and u is a complex parameter related to the moving path. The maximal point of the Gaussian coherent state g(x; u) in Equation (14) is just at x = Re( √ 2u). The 2D HG eigenmodes in a spherical cavity with the transverse orders of n x and n y are given by [1] where Θ ± (r) comprises the quadratic phase and Gouy phase, r = ( x, y, z), and the symbol of ± represents the propagating directions. From Equation (16), the transverse part of the HG mode can be given by ψ (HG) n x ,n y ( x, y) = ψ n x ( x) ψ n y ( y). Using Equations (12)- (17), the 2D Gaussian wave-packet state with the maximal point at (12) and (13) can be generalized as The Gaussian wave-packet state G(r; u ± s , v ± s ) can be seen to be the product of two 1D quantum coherent states except for the phase term exp[−i Θ ± (r)]. In terms of G(r; u ± s , v ± s ), the wave representation for the geometric mode with the leftward and rightward parts of a complete tracing can be expressed as [18] where the phase term exp[i(N x + N y ) θ s ] indicates the phase arising from the transverse component. From Equation (19), the intensity of the geometric mode within −L c ≤ z ≤ 0 can be calculated by where the intensity emitted from the output coupler only considers the forward part Ψ + N x ,N y (r) 2 .

Experimental Results and Theoretical Confirmations
The experimental setup was a diode-end-pumped Nd:YVO 4 laser with a concave-flat cavity. The gain medium was an a-cut 2.0 at.% Nd:YVO 4 crystal with a length of 2 mm and an aperture of 10 × 10 mm 2 . The aperture size of the used Nd:YVO 4 crystal was much larger than the conventional one for generating the transverse modes as high as possible. Both end surfaces of the Nd:YVO 4 crystal were coated as antireflective at 1064 nm (R < 0.2%). The laser crystal was wrapped with indium foil and mounted into a copper holder with water cooling. A concave mirror with a radius-of-curvature of 20 mm was employed as the input mirror whose entrance face was coated as antireflective at the pump wavelength of 808 nm and the second surface was coated as highly reflective at 1064 nm (R > 99.8%) and highly transmissive at 808 nm (T > 95%). A flat mirror with a partial reflection of 97% at 1064 nm was used as the output coupler. The pump source was a 3.0 W 808 nm fiber-coupled laser diode. The core diameter and the numerical aperture of the coupled fiber were 100 µm and 0.16, respectively. The laser crystal was placed behind the input mirror 1-2 mm. The pump beam from the fiber output was focused into the laser crystal with an average pump radius of approximately w p = 100 µm.
The wave-packet state Ψ ± N x ,N y (r) in Equation (19) can reveal the spatial structure with ray-wave duality. Moreover, the state Ψ ± N x ,N y (r) can unify the representation for HG eigenmodes and geometric modes. Let us set N x = 0 and N y = N to demonstrate the geometric modes associated with the planar rays. From Equation (12), N x = 0 indicates u ± s = 0 and g(x; 0) = ψ 0 (x). As a result, Equation (19) becomes Figure 4 depicts the calculated patterns I(r) = Ψ + 0,N (r) by using Equation (20) with N = 100 to correspond to the trajectories shown in Figure 2. The calculated distributions agree very well with the geometric rays for all cases. On the other hand, if the transverse order N is not high enough, the distribution cannot exhibit the ray feature. When N < Q/2, the mathematical form of Ψ ± 0,N (r) can be verified to approach to the HG mode ψ (HG) 0,N ( x, y), as discussed later. = Ψ + Ψ r r r by using Equation (20) with 100 N = to correspond to the trajectories shown in Figure 2. The calculated distributions agree very well with the geometric rays for all cases. On the other hand, if the transverse order N is not high enough, the distribution cannot exhibit the ray feature. When  By using Equation (20), the calculated far-field patterns 2 0, from Equation (20) with N = 100, corresponding to the results shown in Figure 2.
By using Equation (20), the calculated far-field patterns I( x, y) = Ψ + 0,N ( x, y, z = 20) 2 with P/Q = 2/7, φ y = 0, and z = z/z R are shown in Figure 5a. The values of N in Figure 5a cover from 0 to 370 to display the variation from the HG modes to the geometric modes. The experimental results are shown in Figure 5b for comparison. The experimental results were obtained by using a diode-pumped solid-state laser with the off-axis pumping [4,[14][15][16].
The experimental results agree very well with the calculated patterns for the variation from the HG modes to the geometric modes. Note that the transverse order N can be determined from ∆y by finding the integer closest to (∆y/ω c ) 2 , where ω c is the fundamental mode size in the gain medium.
N I x y x y z  (20) can be rewritten as the integration: For a nondegenerate cavity, corresponding to Q → ∞ , we can confirm that the state Ψ ± 0,N (r) is directly related to the HG mode ψ (HG) 0,N ( x, y). For Q → ∞ , the summation in Equation (20) can be rewritten as the integration: with Applying the orthogonal property to Equation (14), we can obtain Equation (24) implies that the Gaussian wave packet g( y, v ± ) can be used to express the eigenmodes ψ N ( y) with an integral transformation. Substituting Equation (24) into Equation (21), we can obtain where a N = N N/2 e −N/2 / √ N!. Equation (25) means that the origin of HG eigenmodes is related to a nondegenerate cavity under the off-axis local pumping. Figure 6 shows the calculated patterns for Ψ ± 0,N (r)  Besides the planar rays, the representation in Equation (19) can be extended to the nonplanar tracings. The calculated results by using Equation (19) with  Figure 7, all the cases one-to-one corresponding to the nonplanar geometric modes shown in Figure 3. The calculated results are in good agreement with the geometric spots for all cases. Similar to the planar modes, if the transverse orders x N and y N are less than Q, the wave intensity cannot be localized on the geometric rays. For can be found to exhibit the spatial pattern very close to the LG mode. Figure 8 shows  Besides the planar rays, the representation in Equation (19) can be extended to the nonplanar tracings. The calculated results by using Equation (19) with N x = N y = 50 for the transverse patterns on the flat mirror are shown in Figure 7, all the cases one-to-one corresponding to the nonplanar geometric modes shown in Figure 3. The calculated results are in good agreement with the geometric spots for all cases. Similar to the planar modes, if the transverse orders N x and N y are less than Q, the wave intensity cannot be localized on the geometric rays. For N < Q/2, the state Ψ ± N,N (r) with φ x − φ y = π/2 can be found to exhibit the spatial pattern very close to the LG mode. Figure 8 shows the calculated transverse patterns of Ψ ± N,N (r) on the flat mirror with P/Q = 2/7 and φ x − φ y = π/2 for N = 1, 3, 14 to reveal the variation from the LG modes to the geometric modes. The present demonstration can be extended to the elliptical mode by considering φ x − φ y = π/2. sults are in good agreement with the geometric spots for all cases. Similar to the planar modes, if the transverse orders x N and y N are less than Q, the wave intensity cannot be localized on the geometric rays. For can be found to exhibit the spatial pattern very close to the LG mode. Figure 8 shows the calculated transverse patterns of , ( )

Eigenfrequency Spectrum and Lissajous Geometric Modes
Next, we review the eigenfrequency spectrum subject to the influence of a gain medium with birefringence. In the vicinity of In terms of the Rayleigh range can be approximated as Nd-doped vanadate crystals are frequently utilized as gain materials for implementing diode-pumped solid-state lasers. Due to the large birefringence, the Nd-doped vanadate crystal unavoidably leads to significant astigmatism inside the resonator without

Eigenfrequency Spectrum and Lissajous Geometric Modes
Next, we review the eigenfrequency spectrum subject to the influence of a gain medium with birefringence. In the vicinity of L P/Q , let the cavity length be given by L c = L P/Q + δL c . The variation of the ratio ∆ f T /∆ f L due to an infinitesimal change δL c can be verified as In terms of the Rayleigh range z R = (R c /2) sin(2πP/Q), the frequency ratio ∆ f T /∆ f L for the cavity length near L P/Q can be approximated as Nd-doped vanadate crystals are frequently utilized as gain materials for implementing diode-pumped solid-state lasers. Due to the large birefringence, the Nd-doped vanadate crystal unavoidably leads to significant astigmatism inside the resonator without compensation. Consequently, there is a considerable difference between the effective cavity lengths L c,x and L c,y , where L c,x and L c,y are the effective cavity lengths for the transverse parts in the directions of x//b and y//c, respectively. By using the new variables of d = L c,y − L c,x and L c = (L c,y + L c,x )/2, the frequency ratios for x and y directions can be individually given by The difference between ∆ f T,x and ∆ f T,y leads the relative eigenfrequency spectrum to be The calculated results for the relative spectrum F(δn; L c , d) as a function of L c /R c for |δn i | ≤ 10 with i = 1, 2, 3 are shown in Figure 9. The values of d/R c used in Figure 9a,b are 3.3 × 10 −3 and 6.6 × 10 −3 , respectively. It can be found that the astigmatism causes the degeneracies at L c = L P/Q to split. Furthermore, the width of the frequency splitting can be seen to increase with increasing the value of d that measures the degree of astigmatism. Even though the astigmatism destroys the original degeneracies at L c = L P/Q , a variety of new fine degeneracies can be generated in the vicinity, as shown in Figure 9c for the cavity length near P/Q = 1/4.
Equations (27)- (29) can be exploited to obtain the cavity lengths for the new fine degeneracies. From Equations (27)-(29), the frequency spacings ∆ f T,x and ∆ f T,y for the cavity length near L P/Q can be explicitly approximated as Substituting Equations (31) and (32) into Equation (30), the cavity lengths for the new fine degeneracies can be derived from the condition of F(δn; L c , d) = 0. Accordingly, the requirements for the new degeneracies are where p and q are integers. Equations (33) and (34) demand the integer pair (p, q) to satisfy p + q = QK with K to be an integer. Furthermore, the cavity length δL c with respect to L P/Q for the new degeneracy of (p, q) can be obtained from Equation (33) as Figure 10a-c shows the relative spectrum F(δn; L c , d) as a function of δL c /d for |δn i | ≤ 10 with i = 1, 2, 3 for three cases of P/Q = 1/3, 2/5, and 3/8, respectively. 9a,b are 3.3 × 10 and 6.6 × 10 , respectively. It can be found that the astigmatism causes the degeneracies at , a variety of new fine degeneracies can be generated in the vicinity, as shown in Figure 9c for the cavity length near / 1/4 P Q = . can be explicitly approximated as

Parametric Equations for Lissajous Geometric Rays
From Equations (31), (32), and (35), the frequency spacings ∆ f T,x and ∆ f T,y for the new degeneracy (p, q) can be rewritten as where β = d/(2πz R KP). The transverse patterns of the high-order modes at the degeneracy (p, q) have been observed to be localized on the Lissajous figures. Consequently, these observed modes are called the Lissajous geometric modes. The Rayleigh ranges for x and y directions can be given by (39) spect to / P Q L for the new degeneracy of ( , ) p q can be obtained from Equation (33) as   with z R,x and z R,y , the beam waists are w o,x = λz R,x /π and w o,y = λz R,y /π, where λ is the laser wavelength. In terms of ∆ f T,x , ∆ f T,y , z R,x and z R,y , the parametric equations for an astigmatic cavity can be derived as [18,41] where s = 0, 1, 2, · · · is the running index for different paths, w x (z) = w o, w y (z) = w o,y 1 + (z/z R,y ) 2 , θ G,x (z)= tan −1 (z/z R,x ), and θ G,y (z)= tan −1 z/z R,y , φ x and φ y are associated with the initial condition, N x and N y are the transverse orders in the xand y-directions, and − and + in the symbol of ± correspond to the backward and forward rays, respectively. The parameter β for the geometric mode needs to be a rational number m/M for forming a complete orbit with finite rays, where m and M are coprime positive integers. By using β = m/M, the mode spacings ∆ f T,x and ∆ f T,y for the degeneracy (p, q) are given by ∆ f T,y /∆ f L = (PM y )/(QM) (43) where M x = M − mq and M y = M + mp. Substituting Equations (42) and (43) into Equations (40) and (41), the parametric equations for the periodic rays of the Lissajous geometric mode related to the degeneracy (p, q) in the vicinity of L P/Q could be derived as y s (z) = N y w y (z) cos θ p,s + φ y ± θ G,y (z) (45) where s = 0, 1, 2, · · · , MQ − 1 is the running index, θ s = 2πs(P/Q), θ q,s = θ s (M − mq)/M, , and θ G,y (z) = tan −1 (z/z R,y ). The total ray number for forming a geometric period is MQ. Typically, the value of M ranges from 10 2 to 10 4 .

Wave Representation for Lissajous Geometric Modes
Next, we review that the parametric equations in Equations (44) and (45) for geometric rays can be implanted into the coherent state in Equation (14) to derive the Gaussian wave-packet representation for Lissajous geometric modes. By using the dimensionless variables x = √ 2x/w x (z) and y = √ 2y/w y (z), the parametric equations in Equations (44) and (45) can be written as Considering the effect of astigmatism in Equations (16)- (19), these equations for the HG eigenmode and geometric modes are revised as Φ ± n x ,n y (r) = e −i Θ ± (r) ψ n x ( x) e ∓i n x θ G,x (z) ψ n y ( y) e ∓i n y θ G,y (z) (48) Equations (49)-(51) stand for the Gaussian wave-packet representation for the Lissajous geometric modes.
The spatial characteristics of Lissajous geometric modes can be deeply comprehended from the parametric equations in Equations (44) and (45) and the wave representation in Equations (49)- (51). We can use Equations (44) and (45) to verify that the total number of output rays is N = MQ/C d , where C d is the common divisor among (M − qm), (M + pm), and MQ. The total path length on the transverse plane for Lissajous geometric modes can be expressed as y(θ; z) = N y w y (z) cos[p(θ + φ y )] Numerical calculations revealed that the path length L p,q (N x , N y ; z) relies primarily on the parameters (p, q) and N x , N y and is nearly independent of the phase factors (φ x , φ y ). Moreover, the path length L p,q (N x , N y ; z) can be accurately approximated with Since both values of w x (z) and w y (z) are rather close, we can use √ 2w x (z) ≈ √ 2w y (z) as the effective width of the Gaussian wave packets. Based on the number of Gaussian wave packets N = MQ/C d and the overall length L p,q (N x , N y ; z) in Equation (55), the degree of the average overlapping between two contiguous Gaussian spots can be estimated as where n is the greatest common divisor of p and q. The value of n means the number of Lissajous curves appearing in the transverse plane of the corresponding geometric mode. When D 1, the Gaussian wave packets in Equation (51) substantially overlap with each other to display a feature of continuous curves, not the original discrete spots. In contrast, when the value of D is considerably less than 1.0, the transverse pattern can be seen to reveal a feature of separate spots corresponding to geometric rays.

Experimental and Numerical Results
We present an overview of the comparison between experimental and theoretical results for Lissajous geometric modes. It has been demonstrated that Lissajous geometric modes can be well generated from an off-axis diode-end-pumped Nd:YVO 4 laser. The detailed description for the experiment can be obtained from previous papers [14,40,41]. When a short laser crystal with length of 1~2 mm was used in experiment, the value of M was found to be greater than 300 by far. Since the transverse orders of N x and N y for Lissajous lasing modes can be approximately 10 2 ∼ 10 3 by varying the off-axis displacement of the pumping, the degree of the overlapping D in Equation (56) is then significantly larger than 1.0. The first experimental case is shown in Figure 11a for the Lissajous geometric mode corresponding to the parameters of (P, Q) = (1, 3), (p, q) = (4, −1), and φ x = φ y = π/3 for the transverse patterns with propagation dependence. The theoretical results are depicted in Figure 11b by using Equations (14), (46), (47), (50) and (51) with the parameters and M = 400, N x , N y = (400, 200), and m = 1. Theoretical patterns can be seen to agree very well with experimental results. From Equation (56), the value of D for the case in Figure 11 can be obtained to be approximately 7.07. As expected, the transverse pattern displays a characteristic of continuous curve due to large overlapping between Gaussian wave packets.
The condition for the critical overlapping between Gaussian wave packets can be specified as D = 1. Figure 12 shows the influence of the degree of the overlapping D on the transverse far-field patterns with the parameters of (P, Q) = (1, 4), (p, q) = (1, 7), φ x = φ y = 0, M = 100, and m = 1. The values of N x , N y in Figure 12a-d are (200, 200), (400, 400), (600, 600), and (800, 800), respectively. Accordingly, the values of D for Figure 12a-d are 1.414, 1.000, 0816, and 0.707, respectively. Figure 12 clearly reveals that the transverse pattern varies from the continuous curves to the discrete spots for the value of D from greater to less than 1.0. = = /3 for the transverse patterns with propagation dependence. The theoretical results are depicted in Figure 11b by using Equations (14), (46), (47), (50) and (51) with the parameters and = 400, , = (400, 200), and = 1. Theoretical patterns can be seen to agree very well with experimental results. From Equation (56), the value of D for the case in Figure 11 can be obtained to be approximately 7.07. As expected, the transverse pattern displays a characteristic of continuous curve due to large overlapping between Gaussian wave packets.  The condition for the critical overlapping between Gaussian wave packets can be specified as = 1. Figure 12 Figure 12 clearly reveals that the transverse pattern varies from the continuous curves to the discrete spots for the value of D from greater to less than 1.0. Experimentally, the astigmatism needs to be increased to reduce the value of M for realizing the condition of 1. One practical way of increasing the astigmatism is to use a gain medium with longer length. By using a length of 10 mm, the value of M can be reduced to the range between 40 and 120. In addition to reducing the value of M, Equation (56) indicates that the degree of the overlapping D can be varied by the factor which Experimentally, the astigmatism needs to be increased to reduce the value of M for realizing the condition of D < 1. One practical way of increasing the astigmatism is to use a gain medium with longer length. By using a length of 10 mm, the value of M can be reduced to the range between 40 and 120. In addition to reducing the value of M, Equation (56) indicates that the degree of the overlapping D can be varied by the factor C d which is the common divisor among (M − qm), (M + pm), and MQ. The influence of the factor C d is shown in Figure 13 for two cases of M = 100 and M = 101 for the far-field transverse patterns with the parameters of (P, Q) = (1, 4), (p, q) = (1, 7), N x , N y = (500, 300), φ x = φ y = 0, and m = 1. All far-field patterns were measured at the position of z = 100 L c . The results in Figure 13a with M = 100 can be found to correspond to C d = 1 and D = 0.898. On the other hand, the case shown in Figure 13b with M = 101 corresponds to C d = 2 and D = 0.454. It can be found that even though the difference of the astigmatism is very small for M = 100 and M = 101, the degree of the overlapping D varies significantly due to the influence of the factor C d . Figure 14a   The results shown in Figures 13 and 14 are for the cases with n = 1. The value of n arising from the greatest common divisor of p and q determine the number of transverse Lissajous curves in the geometric mode. An example for the far-field transverse pattern with n = 2 is depicted in Figure 15a with the parameters of (P, Q) = (1, 4), (p, q) = (−2, 6), N x , N y = (300, 200), φ x = φ y = π/2, M = 102, and m = 1. The result in Figure 15a corresponds to C d = 4 and D = 0.67. Although the value of n determines the number of Lissajous curves in the geometric mode, some phase values of φ x , φ y may lead the corresponding curves to coincide together, as depicted in Figure 15b. All parameters for Figure 15a,b are the same except for φ x = φ y = π/4. Once again, all theoretical results are in good agreement with experimental patterns.   The results shown in Figures 13 and 14 are for the cases with = 1. The value of n arising from the greatest common divisor of p and q determine the number of transverse Lissajous curves in the geometric mode. An example for the far-field transverse pattern with = 2 is depicted in Figure 15a

Damping Effect
Experimental patterns shown in Figures 11-15 were generated by using an output coupler with a transmission as low as 2% at the lasing wavelength. It has been experimentally observed that the spatial structures of lasing modes can be significantly affected by the transmission of the output coupler. It has been demonstrated [53] that an output coupler with high transmission can be used to deliberately generate the asymmetric lasing

Damping Effect
Experimental patterns shown in Figures 11-15 were generated by using an output coupler with a transmission as low as 2% at the lasing wavelength. It has been experimentally observed that the spatial structures of lasing modes can be significantly affected by the transmission of the output coupler. It has been demonstrated [53] that an output coupler with high transmission can be used to deliberately generate the asymmetric lasing modes. In theory, the influence of output transmission on the spatial structures of Lissajous lasing modes can be directly taken into account from Equation (51). Let the effective damping coefficient be considered as a loss factor from the output transmission T oc . The loss factor on the sth ray in Equation (51) can be given by exp(−γθ s ) with γ = T oc /2π. Taking the damping effect into account, the Gaussian wave-packet representation for Lissajous geometric modes can be expressed as [38,53] Figures 16a and 17a show experimental results obtained, respectively, with T oc = 2% and 15% for the far-field patterns of Lissajous geometric modes with the cavity length near (P, Q) = (1, 4) for four different (p, q). The first and second row in Figures 16a  and 17a are the cases for φ x = φ y = 0 and φ x = φ y = π/4, respectively. Note that the far-field patterns were measured at the position of z = 50 L c . Making a comparison between the patterns of Figures 16a and 17a, a high output transmission can lead the Lissajous transverse pattern to exhibit a feature of broken line. The asymmetric structure arising from a high output transmission can be found to depend on not only (p, q) but also φ x , φ y . By using Equations (51) and (57), the calculated patterns corresponding to Figures 16a and 17a are shown in Figures 16b and 17b with experimental parameters and M = 200, N x , N y = (150, 100), and m = 1. All the calculated patterns are in excellent agreement with experimental results to manifest the damping effect on the spatial structures of Lissajous geometric modes. Finally, experimental results for the far-field patterns of Lissajous geometric modes generated by using an output coupler with a higher transmission of T oc = 25% are shown in Figure 18a. By using Equation (57) with the parameters of M = 200, N x , N y = (75, 50), and m = 1, theoretical patterns for T oc = 25% are calculated and shown in Figure 18b. The superiority of the Gaussian wave-packet representation can be confirmed from the good agreement between experimental and theoretical patterns in manifesting the damping effect arising from cavity losses.
imental results to manifest the damping effect on the spatial structures of Lissajous geometric modes. Finally, experimental results for the far-field patterns of Lissajous geometric modes generated by using an output coupler with a higher transmission of = 25% are shown in Figure 18a. By using Equation (57) with the parameters of = 200 , , = (75, 50), and = 1, theoretical patterns for = 25% are calculated and shown in Figure 18b. The superiority of the Gaussian wave-packet representation can be confirmed from the good agreement between experimental and theoretical patterns in manifesting the damping effect arising from cavity losses.   (57) with the parameters used in Figure 16 (b). Figure 17. (a) Experimental results for the far-field patterns of Lissajous geometric modes generated by using an output coupler with T oc = 15%. (b) Numerical patterns obtained from Equation (57) with the parameters used in Figure 16b.  (57) with the parameters used in Figure 16 (b).

Conclusions
We presented a thorough overview of the laser transverse modes with ray-wave duality in spherical cavities with and without astigmatism. We first reviewed the eigenfrequency spectrum in an ideal spherical cavity to manifest the critical role of degeneracy in the ray-wave correspondence. The derivation of Gaussian wave-packet states was given to understand the unification of the HG eigenmodes and the planar geometric modes. We extended the wave representation of the planar geometric modes to the case for the elliptical modes with orbital angular momentum. Furthermore, we reviewed the fine structure of the characteristic frequency in a spherical cavity with astigmatism. Taking the astigmatism into account, the wave packet representation for Lissajous geometric modes has been completely introduced. Finally, we reviewed the damping effect on the formation of transverse modes to reveal the generation of asymmetric geometric modes. It is believed that this review can offer useful insights into the ray-wave duality in mesoscopic optics and laser physics.  Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.
Data Availability Statement: All of the data reported in the paper are presented in the main text. Any other data will be provided on request.