Design of cycloidal rays in optical waveguides in analogy to the fastest descending problem

In this work, we present the design of cycloidal waveguides from a gradient refractive index (GRIN) medium in analogy to the fastest descending problem in classical mechanics. Light rays propagate along cycloids in this medium, of which the refractive index can be determined from relating to the descending speed under gravity force. It can be used as GRIN lenses or waveguides, and the frequency specific focusing and imaging properties have been discussed. The results suggest that the waveguide can be viewed as an optical filter. Its frequency response characteristics change with the refractive index profile and the device geometries.


Introduction
The connection between classical mechanics and optics is a long lasting and still intriguing problem.The movement of an object under the action of an arbitrary force, as dictated by the Euler-Lagrange equation, appears to be deeply connected with the propagation of light in a special class of medium with gradient refractive index (GRIN), as dictated by the famous Fermat's principle.Both follow the least action principle, which can be described by the variational method.Guided by this beautiful analogy, people formulated the peer law of Newtonian mechanics in optics [1] , developed general methods of ray tracing in a symmetric GRIN profile from Fermat's ray invariants [2,3] and optical Binet equation [4] , and solved the optical Lagrange equation using the first integral method to control light trace in GRIN media [5] .
In recent years, progress in material science has made it possible to fabricate optical medium of GRIN type by a variety of techniques, such as ion-exchange process [6,7] , photo-thermal reaction in glass [8] , nano-scale porous silicon [9] , gradient-structured metamaterials [10] , and additive manufacturing with nanocomposite materials [11,12] .These artificial materials enable people to apply the well-known methods and conclusions in mechanics in novel optical design, which has potential application in invisibility cloaks [13,14] , perfect imaging lenses [15,16] , and transformation optical media [17,18] .
In this work, we propose the design of a GRIN optical device based on the analogy to the famous problem of fastest descent (also known as the Brachistochrone problem).The light paths follow cycloidal curves, just like the descending object, if the refractive index profile is devised to imitate the effect of gravitational force acting on the object.In this medium, we can construct the refractive index profile to be cylindrically symmetric, such that light is confined transversally and directed along the cylinder axis, which can be applied as a waveguide.Furthermore, we provide field simulation results to support our design.The light field distribution reveals the physical nature behind the geometric optics that light rays obey in the medium.

Theoretical model
To restate the fastest descending problem, we need to seek a track for an object to descend between two fixed points within the shortest time.In Cartesian coordinates, the object descends at a speed under the action of gravitation  = √2 0  , where  0 is the gravitational acceleration.Therefore, the duration time for the object to descend from point A to point B is: where the actual path is determined by the function  = ().
Through variational method, the path that minimizes  is a cycloid, which satisfies the parametric equations: in which  is a parameter determining the spatial extension of the curve.The schematics of a cycloid is shown in Fig. 1(a).When it comes to optics, the speed of light propagating in a dielectric medium is  = / , where  is the refractive index of the medium.It is found that if we prescribe a refractive index profile of the form: where  0 is a constant, the medium plays a similar role as gravitation does in the fastest descending problem.Consequently, light rays will simply follow cycloids described by Eq. ( 2), which minimize the length of optical path for the light rays.This fact agrees with Fermat's principle.On the other hand, light rays can only propagate along cycloids when light feels a GRIN profile of the medium, as dictated by Eq. ( 3).It represents a specific class of GRIN media, and we call them the cycloidal media in following discussion for convenience.If we rotate the cycloidal medium around the  = 0 axis, we can construct a cylindrical waveguide.Fig. 2 shows its cross-section, in which light propagation is restricted in-plane due to symmetry in the azimuthal direction.We can also insert a uniform layer with constant refractive index in the center, serving as the core layer of the waveguide.

Results and discussion
Imagine that light enters the waveguide from left.Without losing generality, we assume light enters in the lower half side of the waveguide at  = −ℎ 1 .Starting from the entrance point, the light rays follow cycloidal trajectories, propagating down from inner to peripheral layer of the waveguide until it reaches  = −ℎ 2 , at which the light rays become parallel to the horizontal axis.Then it is deflected back to the inner layer, and goes straight through the core layer where () is a constant.
Due to the rotation symmetry of the refractive index profile, the light path in the upper half side for  > ℎ 1 is symmetric to that for  < −ℎ 1 .In this way, the propagation of light in the waveguide is periodic, where the period length  depends on the entrance angle of each light ray, denoted as .We plot the first few periods of light trajectories with different  in Fig. 2. If absorption loss is neglected, the light just continues to propagate in the waveguide and can be guided to the far end of the device.
Among all the light rays, there are a group of rays with periods that satisfy rational ratios between each other.In this case, we can find a least common multiple of periods for these rays, so they constitute a commensurate group.After propagating for multiple periods, these light rays will arrive at the same point, similar to the focusing behavior of an ordinary lens.However, unlike the fact that there are infinitely many light rays that are focused for the ordinary lens, the number of light rays in the cycloidal waveguide that arrive at the same focus point is limited.
For example, let us examine all the light rays with the entrance angle  from about -70 o to -20 o .The variation of  versus  in this range is plotted in Fig. 3(a).As the data shows, the ratio between the maximum and minimum of  for these rays is about 12, and therefore the possible rational ratios in this condition are: which includes periods: , 1) ×   (5)   It contains 12 elements in total, each corresponding to a certain  , where   and   are the maximum and minimum of period in the group with the largest and smallest || respectively.
in which (, ) is the light intensity distribution on the image plane, and (, ) is the discrete frequency spectrum of the transmission function on the object plane.From the Kirchhoff diffraction formula of a thin lens [19] , the spatial frequency of an object under coherent illumination equals to the tangent of the entrance angle  multiplied by the wave number  0 of light: where the index ,  run over the items in the commensurate group, and ∆ , in Eq. ( 8) are intervals of discrete frequencies related to sampling precision, which can be treated as constant coefficients.In this sense, the waveguide acts as a spatial frequency selector along its axis, as certain positions on the axis correspond to certain combinations of spatial frequency components.
On the Fourier plane of the lens, these components overlap to form a filtered image of the light source, which may be distinct from the original appearance.To shed light on the filtering effect, we assume that light is emitted from a coherent light source of square shape, as shown in Fig. 3(b).It is worth noticing that, in order to observe the frequency specific image in experiments, an iris has to be used and placed at the focus point, blocking out all other rays that do not belong to the commensurate group., (10)   For instance, assuming g() is a quadratic function of : In this case, () takes the form: The denominator is the square root of a polynomial of variable  with parameters  1 ,  2 and  3 .We can change  1 ,  2 and  3 to modify the refractive index profile, and hence changing the light trajectories.It is found that for certain values of  1 ,  2 and  3 , the difference of period  between light rays emitting at different  becomes larger, such that the range of   /  becomes larger.
For example, when  1 = 1.2 × 10 −4  0 /mm 3 ,  2 = 1.7 × 10 −2  0 /mm 2 and  3 = 0.82 0 /mm ,   /  can reach 24.As Fig. 4 shows, light follows trajectories that are quasi-cycloidal, and in this case the commensurate group becomes: Compared to the case when   /  = 12, the number of light rays with rationally harmonic frequencies increases.The result illustrates that the waveguiding properties of the cycloidal medium can be manipulated through changing the functional form of () , in analogy to changing the gravitational force in the fastest descent problem.
The partial imaging properties of the cycloidal waveguide is reminiscent of an important self-imaging phenomenon in wave optics, the fractional Talbot effect in Fresnel diffraction [20] .The Talbot image is also space dependent, consisting of a series of superposed images of a periodic object.Each image represents a different subset of spatial frequency components of the object.It is found in Fig. 6 that the maximum of the electric field magnitude follows a deflected trace as the refractive index changes gradually, which coincides largely with the corresponding cycloids.

Field simulation result
Moreover, it is observed that as the wavelength decreases to 200m, the propagation of light in the waveguide shows less scattering and dispersion loss, forming a farther-reaching trace along the light path obtained in rays tracing.This feature is in accordance with the fact that geometrical optics becomes a better approximation when the wavelength decreases, wherein the eigen-modes of the waveguide reduce to light rays, as described by cycloidal curves.
In general, the concept of geometrical optics is well established

Conclusion
In this work, we transplant the classical problem of fastest descent into optics, which enables the design of a GRIN medium in which light rays propagate along cycloids.The refractive index of this medium can be prescribed in analogy to the descending speed under gravitation.With rotational symmetry, the medium can function as a waveguide, and its light transmitting, focusing and imaging properties have been discussed.
If we further change the functional form of the refractive index, the trajectories of light rays can be controlled, hence altering the ratio between the distance light rays propagate in one period.In this way, the number of rational ratios in the commensurate group changes, leading to different combinations of light rays at different distances, so the waveguide functions as an optical filter.Specific spatial frequency components of the light source can be guided and selected out for further utilization.Field simulation results support the geometrical optical properties obtained in our method.

3. 1
Constant gravitation model In Fig. 1(b), we depict the trajectories of light rays in the cycloidal medium with  0 = 1 0 /mm.Light emits at different angles, and each path is a cycloid from the same starting point at  = 0, where  denotes the distance light goes in the -direction.If we reverse the light propagation, the paths end with convergence, showing that parallel light rays are focused to the starting point, i.e. the image of an infinitely distant object.Therefore the medium itself can be used as a lens, with its focus length depending on the incident position of light rays.

Fig. 1 (
Fig. 1 (a) Cycloidal track in the fastest descending problem, which is the trace of a point on a rotating circle.(b) Light trajectories in white curves in the cycloidal medium.The notations of ℎ 1 , ℎ 2 and  mentioned in the main text are labeled for one of the curves.On the background mesh colors are used to indicate the value of refractive index.

Fig. 2
Fig. 2 Cycloidal trajectories of light rays (white curves) within the first two periods in the waveguide under the constant gravitation model.The entrance angle  ranges from -66.7 o to -33.9 o , and the distance  that light travels in one period differs with .

1 (
7 o and   = 8.12 m m,  = −21.2o .Given those ratios in Eq. (5), the least common multiple of the numerators is 2310, which means that at  = 2310  all these light rays are converged together after propagating for a multiple number of periods.In addition, prior to this point there are other positions at which only part of the light rays are converged.For instance, at  = 30  , the light rays with the periods: and at  = 6  , the combination becomes: that the position of the focus point depends on the selection of a certain combination of rational ratios.If we place an ordinary lens behind the waveguide, it will form the image of those light rays from a selected focus point.From the viewpoint of Fourier optics, each light ray represents a certain spatial frequency component.The refraction and imaging of light rays through a lens corresponds to a superposition of different spatial frequency components.It can be described by the inverse discrete Fourier transform with finite items: (, ) = 2) 2 ∑ (, ) (  +  ) ∆  ∆  , ,

Fig. 3 (
Fig. 3 (a) Entrance angle  versus period length  in constant gravitation model.(b) Geometry of the light source at  = 0. (c) and (d) Normalized images reconstructed from sum of spatial frequency components in the commensurate group at  = 30  and 2310  , respectively.

Fig. 4
Fig. 4 Cycloidal trajectories of light rays (white curves) within the first two periods in the waveguide under the variable gravitation model.The entrance angle  ranges from -68.8 o to -29.2 o .
3(c) and 3(d), which tell us about the evolution of imaging properties in frequency domain.When  is increased, there are more spatial frequency components from the light source making a contribution, and the image in Fig. 3(d) resembles the source profile to a greater extent than that in Fig. 3(c), as reflected in the ratio    /   taken between local intensity maximum in the central area and near the edge on the image plane.The result is 0.989 in Fig.
3(c) and 1.542 in Fig. 3(d), which shows that light intensity is increasingly distributed around the central area as  increases.Since the selected light rays emit at relatively large angles, they represent high frequency components of spatial structure.In Fourier transform, spectrum in high frequency range comes from rapid change in space domain.Consequently, it is the edge of the light source that is highlighted in the image, while the clarity varies at different distances, as the number of spatial frequency components included in the commensurate group changes.From another perspective, the image can be understood as the interference pattern between light rays arriving at the image plane with different wavefronts.

3. 2
Functional gravitation model In the problem of fastest descent, gravitational force is treated as a constant.If the gravity field changes, the extremal track of the descending object does change its path.This reminds us that for optical design, we can extend the model to introduce space dependent gravitation generated by either real or artificial mass distribution, under which () is an arbitrary function of .As a more general result, the trajectories of light rays follow the equation:  = sign(ℎ 2 − ℎ 1 ) The focus point shifts approximately to  = 2.2 × 10 8   for all these light rays to meet together.In fact, we need not to go as far as the full focus length, since most rays already encounter each other at a smaller distance.The corresponding imaging results at  = 6  , 30  and 2310  are shown in Fig. 5.

Fig. 5 (
Fig. 5 (a) Entrance angle  versus period length  in functional gravitation model.(b)-(d) Normalized images reconstructed from summing up finite items of spatial frequency components in the commensurate group at  = 6  , 30  and 2310  , respectively.

Fig. 6
Fig. 6 shows the field simulation results under the constant and the functional gravitation model at the wavelength of 500m and 200m respectively.It reveals how light propagates in the cycloidal medium from wave optics perspective.A plane wave of width 1.5mm and initial field magnitude of 1V/m is incident from the left side.The wavefront then follows the eikonal equation and gets distorted by the GRIN profile on the light path.Finite element method was used to calculate the static distribution of transverse electric field when light passes through the cycloidal medium.Perfectly matching layers were assigned to the boundaries.Note that compared with Fig. 2 and Fig. 4, the figure has been rotated by 45 o for the convenience of setting up boundary conditions.

Fig. 6
Fig. 6 Electric field of plane waves propagating in the cycloidal medium under the constant gravitation model (a, c) and the variable gravitation model (b, d).The wavelength of input light is 500m (a, b) and 200m (c, d) respectively.To distinguish electric field from the background substance, we use thermal color scale to represent the electric field magnitude, and gray scale to represent the refractive index.
for GRIN media if the wavelength of light satisfies |∇λ| ≪ 1, which holds true in our model.It asserts that the spatial variation of wavelength should be small enough.The results imply that Fermat's principle together with its analogy to classical mechanics is an instructive clue to modeling and designing the proper refractive index profile given the path of light rays.It provides us with an ancient but powerful tool, bringing to optical design more flexibilities and possibilities.With respects to fabrication techniques of the waveguide proposed in this work, the broad index range the GRIN profile covers is challenging for the fabrication process.Future research is called for achieving a large refractive index span of the GRIN profile in optical materials, which allows for expanding the design space of GRIN optical elements and devices.