19 December 2013 Whispering-gallery modes of μ-optic silicon bottle resonator examined using a Fourier–Bessel eigen-state approach
Author Affiliations +
J. of Micro/Nanolithography, MEMS, and MOEMS, 12(4), 043007 (2013). doi:10.1117/1.JMM.12.4.043007
The whispering-gallery mode space of a cylindrically symmetric silicon micro-optic bottle resonator with an axial Gaussian dielectric profile is theoretically examined. Maxwell’s wave equation for the magnetic field, expressed in cylindrical coordinates, is converted into an eigenvalue problem using a Fourier–Bessel basis space to expand the inverse dielectric and field components. The order of the eigen-matrix is reduced using symmetry arguments and azimuthal order segmentation such that the matrix can be solved on a desktop PC. The eigen-space reveals that several states can be confined in the bottle region, whereas another set of states extends the axial length of the structure similar to ridge waveguide modes. The hollow core bottle is also presented as a design modification. The micro-optic bottle resonator structure size makes it directly compatible with silicon-based photonics increasing its potential as a key component in sensing, communication, and other applications.
Gauthier: Whispering-gallery modes of μ-optic silicon bottle resonator examined using a Fourier–Bessel eigen-state approach



The micro-optical resonator is a structure that plays a key role in many photonic devices and can take on several geometries, such as defect sites in photonic crystals,1,2 ring waveguides,3 laser cavities,4,5 microspheres,6 capillaries,7 disks,8 and of particular interest here, the bottle resonator.9 The most common bottle resonator configuration consists of a radially enlarged region on a solid or hollow glass fiber. The radial and axial profile of the enlarged region determines the optical properties of the modes that may circulate in the azimuthal direction.10,11 The resulting bottles are usually quite large compared to the wavelength at which they are designed to operate as they are constructed in glass. We examine theoretically the modal properties of a silicon-based bottle resonator configuration in air. The high-index contrast that exists permits a smaller bottle radius and is on a size scale compatible with silicon photonic devices. For this article, we have chosen the silicon bottle enlargement to have a Gaussian profile in the axial direction. The structure is analyzed using the cylindrical coordinate representation of the full vector wave equation cast as an eigenvalue problem by making use of Fourier–Bessel (FB) basis functions to represent the fields and inverse dielectric profile. The resonator frequencies (wavelengths) and field profiles are computed and presented for various azimuthal mode orders.

There are a number of numerical techniques available for the calculation of the frequencies and field profiles of resonator states. In this article, we are concerned with whispering-gallery modes in a cylindrically symmetric dielectric structure. If the dielectric profile is such that the separation of variable techniques can be applied, Maxwell’s wave equation may be directly solved and field continuity applied at the various interfaces.12 In the case involving more general dielectric profiles, one usually resorts to either the finite-difference-time-domain (FDTD) or finite-element-method (FEM) techniques.13,14 The dielectric structure is decomposed into a three-dimensional grid and the modes are obtained. The computer resources and computation time using these techniques scale with the complexity of structure under analysis. In a number of instances, the dielectric profile examined is such that simplifications in the analysis equations can be performed resulting in a numerical approach that provides the states in a more time efficient manner. The technique presented here is general and can be placed on the same level of complexity as FDTD and FEM. As will be discussed, the present cylindrical symmetry permits a significant reduction in computer resources and computation time.

In the next section, the key mathematical steps leading to the eigen-matrix for the resonators expressed in cylindrical coordinates and using a FB basis function set are provided. Section 3 presents the geometry of the micro-optic silicon bottle resonator and examines the expansion coefficient space related to FB basis. In Sec. 4, the eigenvalue state space and corresponding field profiles, determined through the eigenvectors, are represented. The dominant features of the modes are provided. In Sec. 5, the hollow core silicon bottle resonator is examined as a design modification to the solid core bottle. Appendix A provides a set of equations required to populate the eigen-matrix.


General Mode Solver Development

The mode solver suitable for obtaining the frequencies and field profiles of the states supported by a cylindrically symmetric dielectric structure, “bottle configuration of interest here,” is based on solving Maxwell’s wave equations using an eigenvalue formulation. When the curl equations are combined for a charge-free, current-free, and nonmagnetic medium with a time variation expressed as ejωt, the equations are



In Eq. (1), c is the free-space speed of light and εr is the relative dielectric constant. As for the familiar plane-wave analysis for photonic crystals,15 a set of basic functions is selected from which a series expansion of the inverse dielectric and field components are proposed. In cylindrical space, the basic functions selected are



For the radial coordinate, only the lowest order Bessel function is needed as a series expansion in Jo() can represent any function over a finite interval 0rR.16 The radial expansion index, p, selects the p’th zero of the Bessel function, ρp. Within the interval for r, the lowest order Bessel functions are orthogonal. In the azimuthal direction, the complex exponential function is employed as it is orthogonal over the 2π interval with m as an integer, m+. In the axial direction, a complex exponential is also selected. Here, it is assumed that the dielectric is periodic in z with a period T. The Gn in the exponential is equivalent to reciprocal lattice “vectors” (Ref. 17) for the axial direction (thus scalar) and take on the form n(2π/T) with n+. The axial basis function contributions are orthogonal over the periodic interval T. As the development to follow will involve rather lengthy expressions, the basis function notation is simplified to Fo(), which requires three indices (p,m,n) and within the bracket, a parameter indicating if the basic function is related to the inverse dielectric [Ω=(1/εr)] or field component (rHr,φHφ,zHz). The lower subscript indicates the order of the Bessel function. Radial derivatives give rise to the first-order Bessels and the associated basis function notation is F1(). With these basis functions, the inverse dielectric expansion is


where κΩ is the expansion coefficients and the summation is performed over the indices pΩmΩnΩ now implicitly written as belonging to the inverse dielectric. From a mathematical perspective, it is more practical to develop the eigen-matrix expression using the H wave equation. Once the matrix element generating expressions are obtained, the E matrix expression can be directly obtained by setting all derivatives of the inverse dielectric to zero and substituting the symbol E for H. As such, the eigen-matrix development will proceed using the top of Eq. (1). The H field in cylindrical coordinates is expressed as H=Hrr^+Hφφ^+Hzz^, where the series expansion for each component is of the form:



Here, i is used to represent the field component and the summation is taken over the indices related to that particular field component. The next step is to expand the double curl of Eq. (1) and form three separate equations, one for each of the field components specified by the RHS field. The derivative expressions related to (Hr,Hφ,Hz) are







When the derivatives of the series expansions for the field components and inverse dielectric are performed, the following equations are obtained:







The set of Eqs. (8)–(10) can be cast into an eigenvalue problem by making use of the orthogonality properties of the basic functions and integrating over the cylindrical volume defining the dielectric profile. The resulting eigen-matrix is cast into the following general form:



The column matrix [RΦZ] represents the expansion coefficients collected by field component and structured based on the summation index sequence. The square matrix is segmented into three row blocks and three column blocks. The first row is generated from Eq. (8), the second row from Eq. (9), and the third row from Eq. (10). The expressions required to generate the matrix elements within each block are provided in Appendix A. A fair bit of structuring is required in order to keep the square matrix elements linked to the expansion coefficients of the column matrix. I is the identity matrix. As will be shown in the next sections, symmetry arguments can be used to greatly reduce the order of the matrix in Eq. (11) and can be solved on a desktop PC. Note that expressions provided in Appendix A used to produce the matrix elements in Eq. (11) have no simplifications or assumption introduced other than those indicated prior to Eq. (1) in their derivation; therefore, they can be applied to the determination of the steady states of any nonmagnetic, charge, and current free-dielectric structure with overall cylindrical symmetry. The derivation can be slightly modified to include an axial propagation contribution, kz, by replacing ejGniz with ej(Gni+kz)z in the field expressions [Eq. (4)] and replacing Gni by (Gni+kz) in the matrix element generating expression in Appendix A.

The eigenvalues obtained may be real or complex depending on which type of state they are related to, (ωc)2=Wreal+jWimag. The real and imaginary parts of the state frequencies can be obtained through the following:



The corresponding free space wavelength is determined from



The eigenvectors represent the expansion coefficients of the field profile and may also contain real and complex contributions. These are obtained at the same time as the eigenvalues using, for instance, the eig() function in MATLAB©. Examination of the eigenvector coefficients and determination of the dominant field contribution provide insight into the state’s properties.


Silicon Bottle Dielectric Profile

The silicon bottle resonator configuration is shown in Fig. 1. The structure can be considered as the contact of a uniform dielectric silicon cylinder of height T=3μm and radius D=1.75μm and a silicon Gaussian profile extension centered on the z axis with an amplitude of 0.5 μm (L-D) and σ of 0.33 μm. The structure presents a cylindrical symmetry and can be efficiently characterized using the (r,φ,z) coordinate system. For computation purposes, the structure is taken as periodic in the z direction with a height sufficiently large such that the resonator bottle modes are well isolated from the axial plane periodic interfaces at +T/2 and T/2. The external free-space region extends to a radius of R=3μm and is sufficient to isolate the bottle resonator modes from the radial computation domain edge. In a previous publication dealing with photonic crystals decomposed using a FB basis functions for the (r,θ) polar coordinates, we have shown that the rotational symmetry of the dielectric imposes restrictions on the nonzero expansion coefficients for the inverse dielectric.18 Only the base symmetry of the dielectric and its integer multiples produce nonzero expansion coefficients. The dielectric structure in Fig. 1 is uniform with respect the azimuthal coordinate φ and can thus be regarded as having infinite rotational symmetry. This property implies that the decomposition of the dielectric should be independent of the azimuthal coordinate and imposes the condition that the only nonzero inverse dielectric expansion coefficients in Eq. (3) are obtained when mΩ=0·=0. The next index pair would be mΩ=±1·, which is not attained in a finite series expansion.

Fig. 1

Three-dimensional representation of the high dielectric region of the silicon bottle resonator configuration. The bottle region has a Gaussian profile with design parameters as indicated in the text. For computation purposes, the structure is periodic in the axial direction with sufficient height such that the bottle states are well isolated from the periodic boundaries. The structure domain is extended to a radius of 3 μm, black line external cylinder, with air filling the region between silicon and radial boundary.


As indicated above, the expansion coefficients of the inverse dielectric are required in the eigen-matrix. Figure 2 shows a plot of the absolute value for the real part of the expansion coefficients for up to 100 Bessel Functions (1pΩ100) and 101 axial exponential (50GnΩ50) with mΩ=0. As expected, the dominant expansion coefficients are observed for the zero axial order and lowest Bessel order since the structure is primarily a right cylinder. All nonzero expansion coefficients with axial order other than zero define the bottle dielectric region. The range of the indices and number of basis functions needed to accurately represent the inverse dielectric is determined by examining how well the inverse dielectric profile is reproduced through Eq. (3). The use of 100 Bessel and 101 azimuthal functions reproduces the inverse dielectric bottle region with <1% difference when compared to the original inverse dielectric profile.

Fig. 2

Plot of the absolute value of the real part of the expansion coefficients for the inverse dielectric given in Eq. (3). The Bessel order, pΩ, ranges from 1 to 100 and the axial order, nΩ, ranges from 50 to 50. Rotational symmetry of the structure imposes that the only nonzero expansion coefficients occur when mΩ=0. The imaginary parts to the coefficients (not shown) are 2 to 3 orders of magnitude smaller.



Eigen-State Determination


Matrix Properties/Order Reduction

Determination of the eigen-state requires that the elements of the matrix be computed using the expressions given in Appendix A. The dielectric decomposition has nonzero expansion coefficients only when mΩ=0 for the structure considered with infinite rotational symmetry. This plays a significant role in simplifying the eigen-matrix and reducing the order of the matrix that needs to be solved. The expressions in the eigen-matrix produce a zero matrix element unless mΩ=mLHSmRHS. If all possible azimuthal mode orders are considered and the resulting matrix is examined, the column and rows for the nonzero elements of one azimuthal mode order are padded with zeros produced from the other azimuthal mode orders. As a result, the matrices for each azimuthal order are independent, which can be extracted into separate submatrices and solved independently. This significantly reduces the order of the matrix to which the eigenvalues and eigenvectors are to be determined. In addition, these submatrices may have the order further reduced by one-half, recognizing that the matrix contribution produced for (mLHS,mRHS) both positives is separable from the matrix contribution with (mLHS,mRHS) both negatives. Both of these submatrices (+,+), (,) produce the same eigenvalue space but the eigenvector space is the complex conjugate of each other. The full eigenvector is produced by combining both. As an example of matrix order reduction, using 100 Bessel functions and 101 axial exponentials and requiring the properties, the 50 lowest azimuthal order modes result in a full matrix of order 3,060,300 while through matrix order reduction, the individual azimuthal mode properties can be determined from a matrix with order 30,300, which can be accommodated on a desktop PC. Note: The plane-wave analysis technique applied to this type of dielectric structure would require a unit cell well separated such that the steady states are isolated between cells. A square unit cell would be defined in Fig. 1 rather than a right cylinder. If 100 exponentials were to be used in the expansions for each of the coordinate directions, the matrix would have an order of 3,000,000. Clearly, the FB approach in cylindrical coordinates provides a much more resource efficient and time-managed approach to the steady-state determination process for the cylindrical symmetric dielectric profiles. Using a basis function set related to the symmetry of the dielectric to reduce the order of the eigenmatrix is suggested in Refs. 19 and 20.


Mode Solutions by Azimuthal Mode Number

We first consider obtaining the eigen-states for a mode profile that has a rotational symmetry in the azimuthal direction equal to 22. For this mode order, the only possible nonzero field expansion coefficients occur when mi=±22. The eigen-matrix of the H field for this mode order can be generated using Eq. (11) and the expressions are provided in Appendix A. Sufficiently converged eigenvalues and field profiles are obtained using 75 Bessel terms and 39 axial exponentials for the bottle structure examined here. (Note: The convergence is based on examining the variation in the eigenvalues and eigenvectors of the desired bottle states.) The matrix solution returns 8775 eigenvalues that extend from very small to very large values, real and complex. Within these are contained the desired bottle states. It is expected that the bottle states should have a wavelength located within the scale of the dielectric structure. The eigenvalues [converted to free-space wavelength using Eq. (13)] are plotted in Fig. 3 for the 1 to 3 μm range. The eigenvalues for which the eigenvector expansion coefficient space is dominated by the Hz field component are plotted in the upper trace. The lower trace is plotted for the remaining eigenvalues dominated by the (Hr,Hφ) field components and are related to the other set of hybrid states and in fact are decoupled from Hz for a uniform axial dielectric. Solving the E field equation in Eq. (1) and examining the eigenvectors dominated by the Ez field component expansion coefficients provide the second set of hybrid states. The remainder of the presentation focusses on states dominated by Hz. For the wavelength range depicted in Fig. 3(a), the wavelength states of 2.23 μm and larger (not shown) have eigen frequencies with a large imaginary contribution. These form the set of unphysical states and are similar to physically unrealistic states of slabs and other waveguides.21 The imaginary frequency of the states with wavelengths of 1.71 μm and lower are of the order of 109 and are effectively zero (numerically). These form the eigen-state of the resonator with several of the states localized in the bottle region and other states extending the full axial length of the dielectric structure. The field profiles, reconstructed using the eigenvectors, are examined below and facilitate the distinction between the two types. It is worth noting that the field expression provided in Eq. (4) makes use of the lowest order Bessel function for the entire radial extent. The radial series expansion is suitable for modeling the field in the high dielectric region and in the low dielectric region even if the low dielectric field region were to be of an “evanescent” nature.

Fig. 3

Wavelengths computed for modes with azimuthal order 22 which lie in the 1 to 3 μm range. Top trace is plotted for the eigenvectors dominated by the Hz field expansion coefficients. Within these, three types of states are observed, unphysical, bottle, and dielectric confined states. Lower trace is for the Hr,Hφ dominated states.


Figure 4 shows the intensity of the Hz field for the states at wavelengths 1.707, 1.557, 1.507, and 1.426, plotted in the (r,z) plane, top profiles, and in the (r,φ) plane, lower profiles. The solid line represents the interface between the high and low dielectric regions. Here, the modes are labeled as R0Z0, R0Z1, R1Z0, and R1Z1, which are based on the number of zero crossing of the field in the radial and axial directions. As the RZ null number increases, the fields are less confined to the bottle region typical of optical waveguide solutions in slab and ridge geometries. For the traces in the azimuthal plane, the states are plotted using the axial value that corresponds to the maximum in the field value. Azimuthal mode order of 22 is confirmed as there are 44 high intensity regions per 2π rotation. In general, the excitation of bottle states is accomplished using an external waveguide (fiber) or a focused laser beam. Knowledge of the mode profile and polarization of the bottle states is crucial in designing a large overlap integral between the incident and supported fields. The FB eigen-approach presented here provides a means of accurately calculating the modes frequency, wavelength, and field profile. A parameter highly desired in resonator design is the quality factor, Q. The dielectric structure under analysis has no loss mechanisms presented. Thus, the eigen-frequencies of the bottle states are real and indicate an infinite Q. However, in a real application of the bottle resonator, input and output power coupling mechanisms must be included and the Q of the resonator would be finite. An estimation of the Q factor should be possible by taking the field profile of the bottle state of interest and calculating its overlap with the power bleeding environment (waveguide, lossy medium, …).

Fig. 4

Field profiles in the (r,z)top and (r,φ)lower planes for the bottle confined modes with azimuthal order of 22. Counting the null regions in the radial and axial direction permits the labeling of the states as R0Z0, R0Z1, R1Z0 and R1Z1. The ghosting that is observed in the azimuthal plots is related to limiting the expansion series to 75 Bessel and 39 axial exponentials. It constitutes <1% of the field value.


Six of the modes with wavelengths below 1.426 μm are plotted in Fig. 5. They extend the full axial length and present either an even symmetry (top) or odd symmetry (lower) with respect to the structures’ center and upper/lower axial periodic boundaries. These modes resemble, to a large extent, the guided modes of a ridge waveguides resting on a low dielectric substrate. To a first approximation, whispering-gallery states have been analyzed using the Helmholtz equation with periodic azimuthal boundary conditions.22 It is observed that the modes inner caustic, Rc, decreases in relation to a decrease in mode wavelength. This is a manifestation of the conservation of angular momentum for these modes, all with the same azimuthal symmetry order mi and angular momentum L=mih=Rck. Although these states are available and would produce field profiles contained primarily in the high dielectric region, they would not normally be significantly excited using a tapered fiber. A more efficient way to couple to these states would be, for instance, to adopt an end-fire coupling approach by cutting out a portion of the high dielectric region.

Fig. 5

Mode profiles for six eigen-states with azimuthal order 22. Although the mode profiles are complicated in appearance, they all extend the full axial length and penetrate below the bottle region into the uniform dielectric cylinder. The inner caustic, Rc, is directly related to the mode wavelength and is a manifestation of angular momentum conservation. Top three profiles display a symmetric symmetry while the lower three states are anti-symmetric with respect to the center line.


Applying the matrix order reduction technique, the matrices for azimuthal orders 5 to 30 were produced and the state space was determined. From the eigenvectors, the field profiles were computed and those states which are highly contained in the bottle region have their wavelengths plotted in Fig. 6. The lines between data points link states with the same general mode profiles. As the azimuthal order is reduced, the wavelength associated within a mode class increases as well as the field’s extent into the cylindrical dielectric region. The “cut-off” of the modes results when the state is converted into a mode that extends the full axial height of the dielectric structure. For example, Fig. 7 shows the mode profile for the fundamental state for several azimuthal orders. As the azimuthal order is decreased (wavelength increased), the state extends further into the cylindrical dielectric region. At azimuthal order of 4, the fundamental state extends the full axial range and is no longer considered as confined to the bottle dielectric region.

Fig. 6

Plot the bottle confined mode wavelengths versus the azimuthal mode order. Cut-off occurs for modes as the azimuthal order is reduced. The solid light blue line (top) represents a wavelength of fundamental whispering-gallery modes determined using geometrical optics considerations. Good agreement is observed in the high rotational order range as geometrical and wave optics are in better agreement when the wavelength is small.


Fig. 7

Intensity of the fundamental mode at various azimuthal orders. As the azimuthal order is reduced, the field extends further into the dielectric. At cut-off, the mode extends the full axial direction and is no longer considered as a bottle resonator state.


The solid light blue line (top most line in figure) in Fig. 6 is produced by considering the fundamental mode in the geometrical regime, R0Z0 geometrical. Based on the circumference, C, of the bottle resonator and azimuthal mode order, the free-space wavelengths that are plotted where calculated from λ=CnSi/mi. Good agreement between geometrical and wave optics is observed when the wavelength is small as is the general trend in the traces shown in Fig. 6. The ratio of the R0Z0 FB-computed wavelength to the R0Z0 geometrical wavelength multiplied by the index of refraction of silicon gives the R0Z0 mode’s effective index.

The effect of reducing the amplitude of the Gaussian profile making up the bottle resonator region is shown in Fig. 8 for the mode with azimuthal order of 22. As the amplitude decreases, the modes become less and less confined and extend further into the solid dielectric region. At sufficiently low amplitude, the mode may be cut-off and extend the full axial direction. The figure further shows that by adjusting the amplitude, the wavelength of the fundamental mode can be tuned to that of 1.55 μm. A similar type of curve can be generated for the other azimuthal mode orders and through a combination of selecting the amplitude and σ of the Gaussian profile, the wavelength and number of bottle modes can be adjusted.

Fig. 8

The effect of reducing the amplitude of the Gaussian function that defines the bottle resonator region for the modes with azimuthal order of 22 is plotted. The amplitude of 0.5 corresponds to the amplitude used in previous computations. As the amplitude is decreased, states are forced to extend further into the solid dielectric region and eventually are cut-off as the bottle can no longer confine the mode. Below cut-off, the state extends the full axial direction.


The multimode nature of the bottle confined states permits the specification of two types of free-spectral ranges;2324.25.26 wavelength spacing of similar states with adjacent azimuthal order; and wavelength spacing of adjacent modes within an azimuthal order. Figure 9 shows the free-spectral ranges computed for the fundamental mode, R0Z0, plotted on a log–log trace. The solid black line is the free-spectral range computed using geometrical optics for the R0Z0 mode. For the Gaussian bottle parameters chosen, the FB analysis indicates that both free-spectral ranges can be several hundred nanometers.

Fig. 9

Free-spectral range for the fundamental mode with azimuthal order 22. Red (line with squares) indicates FSR of the fundamental mode relative to adjacent azimuthal order. Green (line with triangles) represents FSR for an azimuthal order between fundamental and next higher order mode (R0Z1). Black (line with diamonds) represents geometrically computed FSR of the fundamental mode relative to adjacent azimuthal order.



Hollow Silicon Bottle Resonator

The demonstration of hollow core rolled up semiconductor structures can be traced to the early work reported by Prinz et al.27,28 The initial structures were designed in the nanometer scale diameter regime and produced by allowing the strain mismatch in the freed layers to provide a means of rolling the sheet into a hollow structure. The technique has been extended into the micron-scale optical domain by increasing the diameter of the hollow tubes and through the utilization of semiconductor materials common employed in integrated optic devices.29,30 A great deal of effort is directed toward the determination of the modal properties that may circulate the circumference of the hollow tube.31,32 In this section of the article, the dependence of the wavelength of the fundamental mode, azimuthal order 22, with respect to the inner radius of the air-filled hollow silicon bottle resonator is examined. The hollow internal region provides an access port for fluid flow33,34 and for mechanical positioning.35 The Gaussian shape bottle resonator region is retained in the analysis and the presence of an inner radius equal to 1.75 μm completely removes the supporting cylindrical region. Figure 10 is produced by modifying the dielectric profile, recomputing the inverse dielectric expansion coefficients and then computing the eigenvalues for the fundamental mode of azimuthal order 22. The wavelength of the mode is plotted versus hollow core radius. The inner air region starts to influence the mode’s wavelength (and profile) when the radius of the hollow region approaches and exceeds that of the original mode’s inner caustic. The trace was produced using 75 Bessel and 39 axial indices. The fluctuations on the trace (amplified here due to the small y axis range) represent <1% of the wavelength value.

Fig. 10

Wavelength of the fundamental mode with azimuthal order 22 plotted as a function of the inner hollow radius. Fluctuation (<1%) in data represents the numerical error in the eigen-mode solver when using 75 Bessel and 39 axial indices. Results indicate that the mode’s wavelength is insensitive to hollow regions with radius <1.2μm. When the radius starts to affect the mode’s profile, the field adjusts to the reduced index value in the central region and shifts to lower wavelength. At an inner radius of 1.75 μm, the entire central cylindrical support has been replaced with low dielectric and only the Gaussian bottle remains.




In this work, we have examined the steady states that are available in a micro-optic silicon optical bottle resonator. The structures are examined using a novel FB basis function set to expand the inverse dielectric and field components in order to cast Maxwell’s wave equations into an eigenvalue formulation. The expressions are provided in a general form and may be applied to solving for steady states in a large number of cylindrically symmetric dielectric configurations. Examination of the eigenvalue space and dominant field component of the associated eigenvector facilitates the determination of bottle confined states and additional extended states. For the parameters of the Gaussian bottle resonator presented, four bottle confined eigen-states are observed when the azimuthal mode order is 22. As the azimuthal order is reduced, the modal wavelengths increase and the extention of the field profile into the high dielectric region increases. A cut-off is defined when the bottle state transforms from one confined to the Gaussian region to one that extends the full axial direction. A hollow variation of the bottle resonator is examined and it is shown that the fundamental mode is unaltered by the air region provided the modal field remains confined to the high dielectric region. The analysis technique developed here along with the micro-optic bottle resonator presented in silicon indicates a rich variety of confined and extended states. Variations on the structure parameters promise interesting micro-optical components and devices on the Silicon platform.


Appendix A:

Generating Expressions for the Matrix Elements of Eq. (11)

When the orthogonality integration is applied to Eqs. (8)–(10), the resulting expressions written below are obtained. The indices (pi*,mi*,ni*) relate to the application of the complex conjugate of a basis function on the right and left hand side prior to integration of the cylindrical domain.



















For each of the nine matrix element expressions, the conditions on the indices to produce nonzero matrix elements are given within the bracket. These conditions arise from the orthogonality properties of the exponential portions of the basis functions. The |pimini indicates the matrix element based on the LHS field indices. The orthogonality condition of the Bessel function integrals on the LHS is not met and produces finite values. The resulting integrals involve the product of three Bessel functions, indicated as (S,T,U,V) in Eqs. (14) to (22) and are given below:















The integrals are evaluated over the normalized range x=rR and involve various combinations of the zero and first-order Bessel functions as well as different powers of x. These integrals, written in the form provided, can be evaluated numerically and the results stored in look-up tables. Since the integrals are independent of the dielectric profile and particular field components required, they only need to be calculated once and can be used in the steady state determination using the expressions presented here for any dielectric profile where the radial direction is expanded in the lowest order Bessel series.

Appendix B:

Eigenvector Space for Fundamental Mode of Azimuthal Order 22

Figure 11 shows a plot for the real and imaginary parts of the expansion coefficients for the Hz field component for the fundamental azimuthal mode of order 22 and wavelength at 1.707 μm. The expansion coefficients are collected in terms of axial index nz as labeled along the x axis of the plot. Within each of these are collected the Bessel indices in increasing order up to pz=75. Such a figure shows that the eigenvector space is well converged.

Fig. 11

Plots of the real and imaginary parts of the eigenvector expansion coefficients for the Hz field component of the fundamental mode with azimuthal order 22. The x axis labeling corresponds to the z axis index nz and within each axial index is plotted the Bessel index in increasing order up to pz=75. The real parts are an order of magnitude larger than the imaginary parts.



This work is supported by a grant provided by Natural Sciences and Engineering Research Council of Canada (NSERC).



S. RobinsonR. Nakkeeran, “Couple mode theory analysis for circular photonic crystal ring resonator-based add-drop filter,” Opt. Eng. 51(11), 114001 (2012).OPEGAR0091-3286http://dx.doi.org/10.1117/1.OE.51.11.114001Google Scholar


M.-K. Seoet al., “Controlled sub-nanometer tuning of photonic crystal resonator by carbonaceous nano-dots,” Opt. Exp. 16(13), 9829–9837 (2008).OPEGAR0091-3286http://dx.doi.org/10.1364/OE.16.009829Google Scholar


C. KusalajeerungS. ChlanggaS. Putikwongsaporn, “Nonlinear switching in silicon-based ring resonators,” Opt. Eng. 50(2), 024601 (2011).OPEGAR0091-3286http://dx.doi.org/10.1117/1.3533033Google Scholar


A. E. Siegman, “Laser beams and resonators: the 1960s,” IEEE J. Sel. Top. Quant. Electron. 6(6), 1380–1388 (2000).IJSQEN1077-260Xhttp://dx.doi.org/10.1109/2944.902192Google Scholar


A. E. Siegman, “Laser beams and resonators: Beyond the 1960s,” IEEE J. Sel. Top. Quant. Electron. 6(6), 1389–1399 (2000).IJSQEN1077-260Xhttp://dx.doi.org/10.1109/2944.902193Google Scholar


N. M. Hanumegowdaet al., “Refractometric sensors based on microsphere resonators,” Appl. Phys. Lett. 87(20), 201107 (2005).APPLAB0003-6951http://dx.doi.org/10.1063/1.2132076Google Scholar


I. M. Whiteet. al., “Versatile waveguide-coupled optofluidic devices based on liquid core optical ring resonators,” Appl. Phys. Lett. 91(24), 241104 (2007).APPLAB0003-6951http://dx.doi.org/10.1063/1.2824843Google Scholar


B. Minet al., “High-Q surface-plasmon-polariton whispering-gallery microcavity,” Nature 457, 455–459 (2009).NATUAS0028-0836http://dx.doi.org/10.1038/nature07627Google Scholar


G. Kakarantzaset al., “Miniature all-fiber devices based on CO2 laser microstructuring of tapered fibers,” Opt. Lett. 26(15), 1137–1139 (2001).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.26.001137Google Scholar


M. Sumetsky, “Whispering-gallery-bottle microcavities: the three-dimensional etalon,” Opt. Lett. 29(1), 8–10 (2004).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.29.000008Google Scholar


G. S. Muruganet. al., “Hollow-bottle optical microresonators,” Opt. Exp. 19(21), 20773–20784 (2011).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.19.020773Google Scholar


G. Arfken, Mathematical Methods for Physicists, pp. 86–90, Academic Press, New York (1970).Google Scholar


N. KanedaB. HousmandT. Itoh, “FDTD analysis of dielectric resonators with curved surfaces,” IEEE Trans. Microwave Theory Technol. 45(9), 1645–1649 (1997).IETMAB0018-9480http://dx.doi.org/10.1109/22.622937Google Scholar


M. Oxborrow, “Traceable 2-D finite-element simulation of the whispering-gallery modes of axisymmetric electromagnetic resonators,” IEEE Trans. Microwave Theory Technol. 55(6), 1209–11218 (2007).IETMAB0018-9480http://dx.doi.org/10.1109/TMTT.2007.897850Google Scholar


K. Buschet al., Eds., Photonic Crystals, pp. 267–270, Wiley-VCH Weinheim (2004).Google Scholar


F. Bowman, Introduction to Bessel Functions, Dover Publications New York (1958).Google Scholar


N. W. AshcroftN. D. Mermin, Solid State Physics, p. 88, Saunders College Philadelphia (1976).Google Scholar


S. R. NewmanR. C. Gauthier, “Representation of photonic crystals and their localized modes through the use of Fourier-Bessel expansions,” IEEE J. Photonics 3(6), 1133–1141 (2011).1943-0655http://dx.doi.org/10.1109/JPHOT.2011.2175909Google Scholar


T. M. Monroet al., “Holey optical fibers: an efficient modal model,” J. Lightwave Technol. 17(6), 1093–1102 (1999).JLTEDG0733-8724http://dx.doi.org/10.1109/50.769313Google Scholar


T. M. Monroet al., “Modeling large air fraction holey fibers,” J. Lightwave Technol. 18(1), 50–56 (2000).JLTEDG0733-8724http://dx.doi.org/10.1109/50.818906Google Scholar


R. SymsJ. Cozens, Optical Guided Waves and Devices, p. 130, McGraw Hill, Amsterdam (1992).Google Scholar


F. LiZ. Mi, “Optically pumped rolled-up InGaAs/GaAs quantum dot microtube lasers,” Opt. Exp. 17(22), 19933–19939 (2009).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.17.019933Google Scholar


H. QuanZ. Guo, “Analysis of whispering-gallery modes in small resonators,” J. Micro/Nanolith. MEMS MOEMS 8(3), 033060 (2009).JMMMGF1932-5134http://dx.doi.org/10.1117/1.3213247Google Scholar


Z. Tianet al., “Selective polarization mode excitation in InGaAs/GaAs microtubes,” Opt. Lett. 36(17), 3506–3508 (2011).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.36.003506Google Scholar


M. SumetskyY. DulashkoR. S. Windeler, “Super free spectral range tunable optical microbuble resonator,” Opt. Lett. 35(11), 1866–1868 (2010).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.35.001866Google Scholar


M. Pöllinger, “Bottle microresonators for applications in quantum optics and all-optical signal processing,” Dissertation Thesis (2010).Google Scholar


V. Y. Prinzet al., “Free-standing and overgrown InGaAs/GaAs nanotubes, nanohelices and their arrays,” Physica E 6(1–4), 828–831 (2000).PELNFM1386-9477http://dx.doi.org/10.1016/S1386-9477(99)00249-0Google Scholar


V. Y. Prinz, “A new concept in fabricating building blocks for nanoelectronic and nanomechanic devices,” Microelectron. Eng. 69(2–4), 466–475 (2003).MIENEF0167-9317http://dx.doi.org/10.1016/S0167-9317(03)00336-8Google Scholar


R. Songmuanget al., “SiOx/Si radial superlattices and microtube optical ring resonators,” Appl. Phys. Lett. 90(9), 091905 (2007).APPLAB0003-6951http://dx.doi.org/10.1063/1.2472546Google Scholar


A. Bernardiet al., “On-chip Si/SiOx microtube refractometer,” Appl. Phys. Lett. 93(9), 094106 (2008).APPLAB0003-6951http://dx.doi.org/10.1063/1.2978239Google Scholar


G. S. Huanget al., “Optical properties of rolled-up tubular microcavities from shaped nanomenbranes,” Appl. Phys. Lett. 94(14), 141901 (2009).APPLAB0003-6951http://dx.doi.org/10.1063/1.3111813Google Scholar


C. Strelowet al., “Light confinement and mode plitting in rolled-up semiconductor microtube bottle resonators,” Phys. Rev. B 85(15), 155329 (2012).PRBMDO1098-0121http://dx.doi.org/10.1103/PhysRevB.85.155329Google Scholar


G. Huanget al., “Rolled-up optical microcavities with subwavelength wall thicknesses for enhanced liquid sensing applications,” ACS Nano 4(6), 3123–3130 (2010).1936-0851http://dx.doi.org/10.1021/nn100456rGoogle Scholar


K. Dietrichet al., “optical modes excited by evanescent-wave-coupled PbS nanocrystals in semiconductor microtube bottle resonators,” Nano Lett. 10(2), 627–631 (2010).NALEFD1530-6984http://dx.doi.org/10.1021/nl9036966Google Scholar


Z. Tianet al., “Single rolled-up InGaAs/GaAs quantum dot microtubes integrated with silicon-on-insulator waveguides,” Opt. Exp. 19(13), 1216–12170 (2011).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.19.012164Google Scholar


Robert C. Gauthier obtained his PhD in physics in 1988, Dalhousie University. He is a registered member of the professional engineers of Ontario and is currently a faculty member in the Department of Electronics at Carleton University. His research interests include theoretical studies photonic crystals, quasi-crystals, and resonators. Theoretical efforts have resulted in the development of 3-D eigen-mode solvers using special basis functions. Past research interests include laser trapping, fiber optic sensors, and integrated optics.

Robert C. Gauthier, "Whispering-gallery modes of μ-optic silicon bottle resonator examined using a Fourier–Bessel eigen-state approach," Journal of Micro/Nanolithography, MEMS, and MOEMS 12(4), 043007 (19 December 2013). http://dx.doi.org/10.1117/1.JMM.12.4.043007




Bessel functions


Micro optics


Back to Top