Design, simulation and characterization of integrated photonic spectrographs for Astronomy II: Low-aberration Generation-II AWG devices with three stigmatic points

In the second part of our series on integrated photonic spectrographs for astronomy, we present theoretical and experimental results on the design, simulation and characterization of custom-manufactured silica-on-silicon Arrayed Waveguide Gratings (AWGs) constructed using the three-stigmatic-point method. We derive several mid-to-high resolution field-flattened AWG designs, targeting resolving powers of 11,000 - 35,000 in the astronomical H-band, by iterative computation of differential coefficients of the optical path function. We use numerical simulations to study the imaging properties of the designs in a wide wavelength range between 1500 nm and 1680 nm. We theoretically discuss the design-specific degradation of spectral resolving power at far-off-centre wavelengths and suggest possible solutions. In the experimental section, we provide characterization results of seven manufactured AWG devices of varying free spectral range and resolution. We obtain estimates on spectral resolving powers of up to 27,600 for polarized input at 1550 nm from measurements of the channel transmission bandwidth. Furthermore, we numerically predict expected resolving powers of up to 36,000 in polarized mode and up to 24,000 in unpolarized mode for direct continuous imaging of the spectrum.


Introduction
Integrated photonic spectrographs on the basis of Arrayed Waveguide Gratings are a subject of active development driven mainly by astronomy and astrophotonics [1][2][3][4], whose special requirements call for customized AWG designs for e.g. high-resolution spectroscopy in the near infrared. Integrated photonics plays an increasingly important role in the miniaturization and mass reduction of spectroscopic and sensing instruments, which is especially beneficial for space-exploration oriented applications. Arrayed Waveguide Gratings, originally introduced as wavelength demultiplexers in telecommunications, are modified for spectroscopic use by removal of the output waveguide array and polishing of the image plane facet to optical quality. The spectral image is captured by a near infrared (NIR) camera on a flat image detector. Canonical AWG designs based on the Rowland geometry known from telecommunications were studied in a previous work [5], hereafter paper I, with the intent to determine the performance of various AWG designs in the resolving power range between 15, 000 and 60, 000 on a low-index silica platform. The Rowland-type AWG suffers from a nonuniformity of the spectral resolution across the flat image plane due to defocus aberration, causing the spectral resolving power to decrease substantially towards the centre of the spectral image. Defocus-corrected AWG-designs using a field-flattening lens in the free propagation region (FPR) have been proposed as a solution [6], which, however, adds complexity to the device fabrication. An elegant alternative is provided through the three-stigmatic-point geometry, which achieves field-flattening by modifying the grating shape and array waveguide lengths [7][8][9]. Working field-flattened devices using three stigmatic points have been demonstrated on polymer and SOI platforms [10,11]. This work is the second paper in a series covering the design of integrated diffractive elements for the "Potsdam Arrayed Waveguide Spectrograph" (PAWS) [12], a semi-integrated photonic spectrograph using various AWG designs as well as photonic echelle gratings (PEGs) with a modified focal plane using two stigmatic points.
We present the implementation of the world's first field-flattened AWG devices on a lowrefractive-index-contrast silica platform (Δ = 0.02) designed specifically for mid-to highresolution spectroscopy in astronomy. The new anastigmatic AWG designs are a modification of the canonical Rowland-type Generation-I designs with spectral resolving powers up to = /Δ = 34, 800 at 1550 nm. In paper I, higher resolving powers were found impractical to achieve experimentally without phase error trimming and were subsequently avoided in the following generation of devices. Stigmatic versions of the Gen-I designs were obtained by defining three stigmatic points on a straight line perpendicular to the principal optical axis, located a focal distance above the grating centre. The general properties, such as waveguide core width = 3.4 μm, grating pitch = 10.57 μm and free spectral range (FSR) were retained in order to minimize the required geometry changes. Like their Rowland-counterparts, the anastigmatic AWGs were designed to operate in a region of the astronomical H-band between 1500 nm and 1700 nm. Fabrication of the designs was accomplished in collaboration with the photonics foundry Enablence USA Components Inc. The devices were fabricated using UV-photolithography and atmospheric pressure chemical vapour deposition (APCVD) of 2 on a substrate (silica-on-silicon). In this work, we present the construction and simulation of anastigmatic AWG designs operating at a central wavelength of 1550 nm, as well as experimental characterization results of fabricated anastigmatic AWG designs. We show that defocus of the spectral image occurs at off-centre wavelengths due to a difference in dispersion between the waveguide array and free propagation region of the anastigmatic AWG. We provide a theoretical condition for defocus cancellation in the wavelength range of operation.

Theory
Aberration theory of the arrayed waveguide grating [7] provides a simple and versatile method for designing various types of AWGs and PEGs with suppressed optical aberrations. Designs known as field-flattened AWGs and two-stigmatic-point PEGs [13] are of special interest for integrated spectroscopy applications due to their uniformity of the point spread function (PSF) across the entire image plane. Low aberration geometries can be derived by defining a set of at most three stigmatic points (i.e. points of perfect constructive interference) and their respective wavelengths in the image plane. Figure 1 schematically illustrates the structure of a field-flattened AWG star coupler. The stigmatic points are arranged in a straight line, minimizing defocus in the entire region between 1 and 3 . Note the non-circular shape of the array-FPR interface on the left side.
A stigmatic-point device geometry is derived by forcing all aberration coefficients in the stigmatic points to zero. The optical path length function (OPF) and its differential coefficients are given by the set of equations [7] Fig. 1. Schematic layout of a field-flattened star coupler with three stigmatic points 1 , 2 and 3 corresponding to shorter and longer wavelengths. Light emerges from the waveguide array on the left side; optimal focus of the spectral image at = .
where and are coordinates of the th stigmatic point, and are the slab waveguide and array waveguide effective index, respectively, and are geometrical lengths of the paths in the input/output FPR (denoted A and B) connecting the stigmatic points with the waveguide array, ( ) is the array waveguide length function, ( ) is the grating period chirp function, ( ) (0) and ( ) (0) are the differential coefficients of the grating curve function ( ) and the free propagation lengths ( ) at = 0, respectively, and is the grating order of the waveguide array. The optical path function ( ) is required to be constant in the range − grat /2 < < grat /2, where grat is the width of the grating (see Figure 1), which implies that its derivatives with respect to must be zero at = 0. We introduce an additional degree of freedom into the aberration theory by Wang et al. by allowing distribution of the three stigmatic points over different spectral orders. Introducing the reduced effective index = ( )/ ( ) and the reduced stigmatic wavelength Λ = / ( ), equation (1b) can be written as The geometry of the optical paths in the AWG is uniquely defined by three stigmatic points, which are required to exactly satisfy (2). Each stigmatic point is defined by its coordinates ( , ), the corresponding wavelength and the spectral order of the stigmatic point. Inserting the stigmatic point properties into (2), together with (1e), yields the differential coefficients of the grating curve function ( ) Once the coefficients ( ) (0) are known, the coefficients of ( ) and ( ) follow immediately from equation (2). The differential coefficients are recursively computed from equations (1f)-(3) using the condition (0) (0) = (1) (0) = 0 and equations (1c)-(1e). The classical Rowland-geometry emerges as a special case when the stigmatic points are located on a circle. The minor generalization of variable spectral orders introduced in this work enables an alternative way of defining the AWG geometry by specifying three stigmatic points in three different spectral orders at one single wavelength as opposed to a definition by three different wavelengths in one single order. In this work, however, the latter approach was implemented.

Design of field-flattened AWG devices
We begin the design process by defining the target material platform to be a low-refractive-index (RI) silica-on-silicon platform with an RI contrast Δ = ( 2 core − 2 clad )/2 2 core = 0.02, cladding refractive index clad = 1.444 and core refractive index core = 1.4738 at = 1550 nm. The waveguide structure used in this work is identical to the structure described in paper I, with a guiding layer thickness of 3.4 μm and a waveguide core width of 3.4 μm -9.07 μm. Waveguide mode fields and wavelength-dependent effective indices of the waveguide core were determined by numerical simulations using a commercial implementation (RSoft) of a full-vectorial 3D-beam propagation method (BPM). The simulated optical properties of the basic building blocks, such as channel waveguides, bends and tapers, have been covered in detail in paper I. In summary, we have defined the fundamental design parameters as listed in Table 1.

μm
Core width (taper) 9.07 μm Bend radius ≥ 1.5 mm In paper I, we have studied canonical Rowland-type AWG designs with target resolving powers of up to 60, 000. This work focuses on AWG designs with modest resolving powers of up to 35,000 to avoid excessive degradation of the PSF due to phase errors. In this work, we have studied eight different custom-designed and fabricated AWGs with FSRs of 11.9 nm -47.2 nm, covering a range of spectral resolving powers of = 11, 500 -34, 800. Since the spot size of the input waveguide and the angular dispersion of the waveguide grating are fixed by the design, the spectral resolution is determined only by the FSR and focal length . Stigmatic points define the locations and associated wavelengths of zero aberration in the focal plane of the AWG and are represented by 1 = (− /2, ) at = 0 − /2, 2 = (0, ) at = 0 , 3 = ( /2, ) at = 0 + /2, where is the width of the dispersed spectral image in the focal plane and 0 is the central wavelength. The three stigmatic points define a focal plane perpendicular to the AWG's optical axis at a distance from the centre of the grating. We have determined for each design by iterative tuning of the parameter and computing of Equations (1) until the grating pitch condition (1) (0) = 1/ was satisfied.
We have generated AWG designs with focal lengths of 5 mm, 7 mm, 9 mm, 12.728 mm and 13.5 mm, enumerated 1 -5, and FSRs 11.9 nm, 16.1 nm, 22.6 nm, 32.2 nm and 47.2 nm, enumerated A -E. In this work, we have implemented the FSR-combinations A-1, A-2, B-3, B-5, C-4, D-3 and E-5. Each AWG was designed for a central wavelength 0 = 1550 nm. The design parameters obtained for these eight combinations are listed in Table 2.

Construction of the field-flattened AWGs
Equations (1)-(3) were solved numerically for each configuration of stigmatic points, yielding the differential coefficients of the Taylor series expansions of the grating curve ( ), array waveguide length ( ) and grating chirp function ( ). Grating point coordinates, i.e. locations of waveguides at the array-FPR boundary, and corresponding array waveguide lengths were obtained by evaluating the Taylor polynomials of ( ) and ( ) on a discrete set of points { | ∈ R, ( ) ∈ Z} in the range − grat /2 < < grat /2. The AWG devices in this work have been designed with grating aperture angles of 25.6 • -27.5 • , resulting in grat = 2.191 mm ( = 5 mm)grat = 5.794 mm ( = 13.5 mm). The optical delays of the AWGs were realized as conventional three-segmented circular waveguide arrays consisting of two straight segments and one circular segment. Narrow-FSR AWGs were implemented in a folded configuration of two overlapping FPRs crossing at an angle ≥ 90 • . This configuration provides sufficient room for a consistent waveguide array structure at arbitrarily large path length increments. Wide-FSR AWGs, since they require less space for path length accumulation, were implemented in the traditional horseshoe shape known from commercial telecom devices. Arrays of input/output waveguides with a wavelength spacing of 1 nm were added to each device. The I/O waveguides were not tapered in order to minimize the mode field diameter to achieve the maximum possible spectral resolution. Lithographic layouts of three selected AWG designs C-4, E-5 and B-5, which represent the main device families, are shown in Figure 2. The remaining designs are structurally similar. Design D-3 is a modification

Numerical simulation
The AWG designs were numerically studied for their transmission characteristics using a combination of 3D-BPM and a model of the AWG utilizing scalar Fraunhofer diffraction in the FPR slab waveguide, which we assume to be homogeneous and isotropic in the plane of propagation. In the first stage of the simulation, we have used 3D-BPM to obtain the mode field distributions of the untapered input/output waveguides, the tapered array waveguides at the FPR-interface and the fundamental propagating mode of the FPR slab. The wavelengthdependent effective indices of the channel waveguides and the slab waveguide were calculated in the wavelength range 1500 nm -1700 nm. Material dispersion was included in the model by using refractive index data for 2 from the RSoft dielectric material library. The calculated mode profiles of the I/O waveguides and tapered array waveguides were transformed into their corresponding far fields in the FPR slab and subsequently used in the scalar diffraction model. The power transmission from input waveguide to output waveguide was simulated by numerically computing [14] where Ω is the signal phase after propagation, , and ℎ are angular far field functions of the input waveguides, array waveguides and output waveguides, respectively; and are emission angles relative to the axis of the input waveguide and array waveguide , respectively; and are effective indices of the waveguide array and the FPR slab, respectively; and are the distances between the waveguide array and input/output waveguides; and is the length of the th array waveguide.
Transmission spectra were calculated for each AWG design between 1500 nm and 1680 nm. The 3-dB bandwidths Δ of the transmission curves were of special interest, as they provide a rough lower estimate of the spectral resolution. Figure 3 shows simulation results for designs A-1, A-2, D-3 and B-5. Figures 3(a), (b) and (c) show the transmission curves of output channels on the short-wavelength edge, the centre and the long-wavelength edge of the output region, respectively. The transmission curves are remarkably similar for edge channels and central channels, to the point of being almost indistinguishable, which indicates a very high uniformity of spectral resolution across the flat focal plane. The sidelobe intensity is -47 dB in all designs.
We have evaluated the spatial uniformity of the transmission peak bandwidth in the spectral image region of the AWG by simulating the transmission for each of the output waveguides distributed equidistantly across the focal plane and measuring its 3-dB bandwidth in the wavelength range 1500 nm -1680 nm. Uniformity of the transmission bandwidths of each channel across the output waveguide array indicates uniformity of the focused beam profile in the focal plane. Figure 4(a) shows the spectral resolving power measured across all output waveguides in the wavelength interval ( 0 − /2, 0 + /2), whereby each waveguide transmits a particular wavelength depending on its location in the focal plane (see top of Figure 3). The resolving power is nearly constant across one FSR in all designs. In the main spectral order, simulated variations of = /Δ , defined as Δ = ( max − min )/ max were found to be on the order of 10 −4 -10 −2 . Designs with shorter focal lengths exhibit less variation of in comparison with larger AWGs with longer focal lengths. The stigmatic point method ensures a low variation of spectral resolving power in the main spectral order spanning one FSR. In practice, however, the AWG devices are expected to operate in a wide wavelength range of a few 100 nm, covering multiple spectral orders. Figure 4(b) shows the variation of spectral resolving power between 1500 nm and 1680 nm. A degradation of spectral resolution is observed at off-centre wavelengths. Again, the severity of non-uniformity increases with focal length, being strongest in design B-5 with = 13.5 mm. Table 3 presents the simulation results for all seven AWG designs. We have observed simulated 3-dB transmission bandwidths of 44.5 pm -133.9 pm, corresponding to resolving powers of 11,500 -34,800 at 1550 nm. The analysis of the simulated transmission presented in the previous section shows a significant increase of transmission peak bandwidth, respectively spectral resolving power at off-centre wavelengths, which becomes more pronounced in AWG designs with larger star couplers. This type of behaviour is very critical for high-resolution astronomical applications due to the requirement of operation in a wide wavelength range. Strong non-uniformity of the diffraction peak width across the spectral orders severely limits the usable wavelength range at high resolutions. In this section, it is shown that the observed degradation of the PSF in anastigmatic AWG designs is caused by the difference of waveguide array and star coupler slab waveguide dispersion and a condition for zero-chromatic-defocus is obtained. The condition for constructive interference at the central point on the focal plane is defined as where 0 and are the central wavelength and grating order of the AWG, respectively; and are the array waveguide effective index and the slab waveguide effective index, respectively; is the length of the ℎ array waveguide; is the propagation length between the ℎ array waveguide and the centre of the focal plane; 0 = 0 + 2 0 is the total optical length of path = 0. The waveguide effective index of silica channel waveguides and slab waveguides in the studied wavelength range can be approximated by In equation (5), the geometrical lengths of the array and slab sections in the anastigmatic case can be expressed as = 0 + Δ + and = 0 + , where 0 + Δ and 0 represent the array waveguide length and free propagation length in the Rowland-geometry, respectively, while and represent the length deviations from the Rowland geometry. With the approximation (6) and ( 0 )Δ = 0 , the constructive interference condition (5) becomes Evaluation of equations (7) at = 0 yields the relation between the modulations of the geometrical path lengths as ( 0 ) + 2 ( 0 ) = 0. Finally, equation (7a) becomes It is immediately clear that the above equation only holds in the cases (i) = 0 and (ii) / = / = 0, as the right hand side is constant with respect to . At off-centre wavelengths, waveguide dispersion gives rise to aberrations, which are represented by the two terms in square brackets on the left hand side of equation (8). The first term depends only on the variation of the waveguide array effective index and adds a constant excess path length increment to the array waveguide lengths, which causes a shift of the central wavelength. The second term depends on the variation of both the waveguide array and slab waveguide effective index and is proportional to the modulation of the array waveguide lengths. In an anastigmatic AWG design, is non-linear and symmetric with respect to the centre of the grating. The aberrations resulting from the second term at off-centre wavelengths tend to defocus the image. The condition for zero chromatic defocus is obtained from equation (8) as The effect of condition (9) on the uniformity of spectral resolving power is shown in Figure 5. = −5.5 × 10 −3 μm −1 . This analysis shows several ways to mitigate or eliminate chromatic defocus. It follows from Equation (8) that small values of result in smaller aberrations. Therefore, designs with smaller foot-print are advantageous in terms of resolving power uniformity. Alternatively, waveguide cores with a large aspect ratio, such as strip waveguides, can be used to achieve a smaller value. A promising candidate for this purpose is the SiN based TriPleX technology [15,16]. Another possibility is the fabrication of AWGs with different refractive indices for the array waveguides and slab waveguides, such that condition (9) is satisfied.

Fabrication and testing
Like their predecessor generation, the field-flattened AWG designs have been fabricated by an external foundry on a Δ = 2%-contrast silica-on-silicon material platform using UV-photolithography and APCVD technology. A 3.4 μm thick Ge-doped silica guiding layer was deposited on top of a 20 μm thick thermal oxide buffer, followed by application of the contact mask and reactive ion etching (RIE) of the waveguide structure. Finally, a 15 μm thick top cladding layer was deposited. Figure 6 shows photographs of three fabricated chips of type C-4, E-5 and A-2. The fabricated AWG chips were tested in our laboratory using a T100S-HP CL-band tunable laser source (TLS)    causes, which can be (i) a transversal periodic modulation of the waveguide array amplitude profile and (ii) a transversal periodic modulation of the waveguide array phases due to systematic fabrication errors. The exact origin of the sidelobe anomaly in the folded-type design is as yet unknown and will be investigated in a future study.
The third group of AWGs consists of devices A-1 and A-2, which are modified B-x type AWGs with a reduced FSR of 11.9 nm. The central-channel insertion loss of the fabricated designs is 1.45 dB in device A-1 and 2 dB in device A-2. The devices have transmission peak 3-dB bandwidths of 96 pm and 70 pm, or, = 16, 100 and = 22, 000, respectively, corresponding to 92.5% and 89.4% of the theoretical resolving power. The polarization sensitivity is lowest in these two devices, with = 28 pm (A-1) and = 27 pm (A-2). Figure 9 shows the measured transmission of AWG A-2. Similarly to the case B-5, this device exhibits symmetrical sidelobes with a crosstalk level of −20 dB relative to the transmission maximum. Similar spectral anomalies can be observed in folded AWG designs of the Rowland-type, presented in paper I, which excludes the field-flattening geometry as a possible cause.
In order to test the uniformity of spectral resolving power across the flat image plane of the AWGs, we have measured the power transmission curves of multiple output waveguides distributed equidistantly in the flat focal plane of the output FPR. The array of output waveguides was designed to cover the entire region of the diffraction image, thereby allowing the characterizaton of the PSF uniformity in the entire focal plane by measuring the transmission curves on each output waveguide in the wavelength range 1500 nm -1680 nm. After obtaining the measurements for all AWG designs, we have calculated estimates of the spectral resolving power in the entire spectral image region (see Figure 4 for comparison), as shown in Figure 10 (left). The results show a low variation of the resolving power across the focal plane, with relative variations of 1.7% (A-1), 1.8% (A-2), 2% (B-3), 2.1% (B-5), 2.4% (C-4), 4.1% (D-3) and 5.8% (E-5). The long-range variation of is shown on the right side of Figure 10. Devices A-2 and B-3 exhibit the smallest long-range variation of with a total difference of 6% between maximum and minimum in both AWGs, closely followed by devices D-3, A-1 and E-5 with variations of 7.5%, 8% and 8.4%, respectively. The strongest variation of was observed in devices B-5 and C-4 with 14% and 15.6%, respectively. The transmission peak bandwidths were measured with an accuracy of ±2 pm. Correspondingly, the resolving power was measured with an accuracy ≤ 1000, which is reflected by the measured short-range variation of .
The AWG devices in this work have been characterized by measuring the transmission from an input waveguide to an array of output waveguides, which only shows the properties of the AWGs as multiplexing/demultiplexing optical filters, but does not reveal the focused diffraction image directly. For the purpose of direct imaging of the spectrum by a camera, the output waveguide array will be removed by dicing of the output facet and the output facet will be polished to optical quality in a pending post-processing step. The procedure was performed on two selected AWGs of the previous generation and the results were presented in paper I. In order to predict the spectral resolving power in the direct imaging regime after removal of the output waveguides for the devices in this work, we have conducted a numerical analysis, presented in the next section.

Numerical prediction of the spectral resolving power after dicing
In order to obtain the spatial profile of the focused beam at the interface between the FPR and the output waveguides in the fabricated AWGs, we have used the model outlined in section 4 to numerically recreate the measured transmission of the devices. In a Monte-Carlo process, which minimized the residual sum of squares between the output of the model and the measured transmission curves by random variation of the optical path lengths [5,17], we have fitted the AWG model with the experimental results, obtaining numerical models which closely match the fabricated AWGs in transmission behaviour. Next, we have used the fitted models to calculate the spatial intensity distribution in the focal plane at several wavelengths for each of the AWG devices. The wavelength was varied in 0.5 nm steps around the central wavelength 1550 nm. The fitting process was repeated multiple times to ensure stability and the obtained resolving power values were averaged over the multiple runs. A variation of the resolving power estimate in a range 1% -4% between fitting runs was observed. Figure 11 shows the results of model fitting and numerical simulation of the diffraction image of AWG A-1. The wavelength axis is shown relative to the output waveguide transmission maximum and the fitted data region contains transmission peaks of two consecutive spectral orders, ensuring a correct FSR in the fitted model. The measured power transmission on the central output waveguide of the device is shown along with the fitted model output in Figure 11(a). The measured data matches the fit well down to a level of -33 dB, below which discrepancies appear. The quality of the fit is sufficient for an estimate of the full-width-half-maximum (FWHM) of the output beam in the focal plane, as the region of interest of the measured curve is approximated with an accuracy < 0.1 dB, which is in the range of the experimental measurement accuracy. Figure 11(b) shows the simulated diffraction images in the focal plane of the AWG at incremental wavelengths from 1549.5 nm to 1550.5 nm in 0.5 nm steps. We have obtained estimates of the spectral resolving power in the direct-imaging regime as where is the wavelength shift and is the corresponding spatial shift of the diffraction image. The estimated resolving power of device A-1 was calculated as 20,540 at 1550 nm, which exceeds the value estimated from the channel transmission measurement by a factor of 1.28. The analysis was performed for all remaining AWGs, each case yielding a higher resolving power than previously estimated from measurement. The lowest predicted direct was calculated for device E-5 as 12,880, while the highest values were achieved by devices A-2 and B-5 with 29,490 and 36,000, respectively.
The above results assume strictly polarized input signals. As mentioned earlier, the resolving capability of the AWG devices suffers from birefringence when operating with unpolarized light. We have evaluated the degradation of resolving power for mixed TE/TM input by using the values obtained in section 5 to calculate the spatial polarization splitting of the PSF. Assuming fully unpolarized light, we added the intensities of the two polarization components and measured the FWHM of the composite peak. In the case of device B-5, we found that operation in the unpolarized regime is not possible, since the TE/TM components appear as two distinct lines due to the large in this AWG. For the remaining devices, we calculated resolving powers between 9,600 (E-5) and 24,000 (A-2) by inserting the composite FWHM into equation (10). We found that the resolving power decreased by 9.4% -35.2% in comparison with the polarized regime. Figure 12 shows the cases of device A-1 (a), which exhibits the lowest ratio /Δ , and device B-5 (b) with the largest /Δ . The summarized characterization results are presented in Table 4.

Conclusion
In our second paper in a series of works on integrated photonic spectrograph design, we have presented a new generation of optimized low-aberration AWG devices for applications in spectroscopy utilizing the three-stigmatic-point construction method. In the theoretical part of this work, we have presented a modified version of AWG aberration theory by Wang et al., introducing the possibility of distributing stigmatic points over multiple spectral orders. We have applied aberration theory of arrayed waveguide gratings on existing conventional AWG designs on a 2%-contrast silica platform to obtain AWG designs with uniform spectral resolving power on a flat focal plane in the near infrared spectral region between 1500 nm and 1680 nm, overcoming the defocus aberration of conventional AWG designs. The AWG designs reported in this work have theoretical resolving powers between = 11, 800 and = 34, 800, determined numerically. Simulations of the AWG designs have shown low variation of spectral resolving power within the main spectral order of the AWG on the order of 10 −4 -10 −2 , as well as variations of the resolving power across the entire wavelength range of operation of 13.9% in the smallest AWG design A-1 to 30% in design B-5. We have theoretically shown that the observed degradation of resolving power is caused by the difference in dispersion between the waveguide array and the free propagation region which causes a chromatic defocus effect and is a property unique to three-stigmatic-point AWG devices. From the analysis, we have obtained a condition for zero chromatic defocus and proposed possible solutions to the issue. The AWG designs constructed in our work have been fabricated by an external foundry on a 6-inch silicon wafer using APCVD technique. The fabricated devices were tested for their spectral transmission characteristics using a T100S-HP tunable laser and a CT440 component tester in the wavelength range 1500 nm -1680 nm.
The FSRs of the fabricated devices were in very good agreement with the design values, indicating that material and geometrical properties of the fabricated waveguide structures are close to the requirements. All devices exhibited insertion losses relative to a reference waveguide at or below 3 dB, the lowest being observed in the smallest high-order AWG A-1 at 1.45 dB.
In terms of spectral performance, the AWGs showed good transmission characteristics with 3-dB transmission bandwidths of 56 pm -152 pm for TE polarized input signals, corresponding to estimates of spectral resolving powers between 10, 200 and 27, 600. Designs A-1 and D-3 performed close to theoretical specifications, reaching 92.5% and 92% of their theoretical resolving power, respectively. Designs A-2, C-4, B-3 and E-5 performed slightly worse at 89.4%, 87%, 86.6% and 86%, respectively. Design B-5, while exhibiting the highest resolution, was most severely affected by fabrication errors, reaching 79.3% of its theoretical performance. The measured variation of the wavelength-dependent resolving power was found to be 1.7% -5.8% in the main spectral order and 6% -15.6% in the entire wavelength range of the measurements. The lowest variation of was observed in the designs A-1, A-2, B-3, D-3 and E-5 with 6% -8.4%. In the designs B-5 and C-4, the variation of was strongest at 14% and 15.6%, respectively.
Measurements of the transmission with TM-polarized input signals showed a polarization dependent wavelength shift between 27 pm in design A-2 and 65 pm in design E-5 despite a square waveguide core profile. In a numerical study of the experimental results, we have used a Monte-Carlo algorithm to fit the models of the AWG devices with measured data in order to calculate the expected spectral image after removal of the output waveguides. Simulations of the fitted AWG models have shown resolving powers of 12,880 (E-5) -36,000 (B-5) to be expected in the direct-imaging regime. Furthermore, the obtained model fits were used to study the impact of polarization sensitivity on the resolving power of the AWGs. The polarization dependence resulted in a degradation of spectral resolving power by 9.4% in design A-1 -35.2% in design C-4. Design A-2 achieved the highest resolving power in unpolarized mode with = 24, 000.
We have successfully demonstrated functional custom-tailored AWG designs with an anastigmatic geometry, achieving spectral resolving powers in the near infrared of up to 36, 000 in the polarized case and up to 24, 000 in the unpolarized case.
The new AWG prototypes presented in this work mark an important milestone on the way towards highly integrated, miniaturized, medium to high resolution spectrographs for astronomy, that may also find applications in other fields. In view of the saving of volume and mass in comparison to bulk optics spectrographs, they are particularly promising for airborne and space applications. The anastigmatic spectroscopic AWG devices will be tested on-sky in the PAWS spectrograph [12] in the near future, supplanting the first generation of Rowland-type AWG designs presented in paper I.

Design, simulation and characterization of integrated photonic spectrographs for
Astronomy II: Low-aberration Generation-II AWG devices with three stigmatic points: supplemental document

LITHOGRAPHIC AWG DEVICE LAYOUTS
This section contains the lithographic AWG layouts not included in the main document.

CHANNEL TRANSMISSION CURVES
Here, we show the channel power transmission curves omitted in the main document.