Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Seeking celestial positronium with an OH-suppressed diffraction-limited spectrograph

Open Access Open Access

Abstract

Celestially, positronium (Ps) has been observed only through gamma-ray emission produced by its annihilation. However, in its triplet state, a Ps atom has a mean lifetime long enough for electronic transitions to occur between quantum states. This produces a recombination spectrum observable in principle at near IR wavelengths, where angular resolution greatly exceeding that of the gamma-ray observations is possible. However, the background in the near IR is dominated by extremely bright atmospheric hydroxyl (OH) emission lines. In this paper, we present the design of a diffraction-limited spectroscopic system using novel photonic components—a photonic lantern, OH fiber Bragg grating filters, and a photonic TIGER 2D pseudo-slit—to observe the Ps Balmer alpha line at 1.3122 µm for the first time, to our knowledge.

© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. INTRODUCTION

It has long been known that positrons are produced in the Milky Way galaxy, through the detection of the 511 keV $\gamma $-rays produced as the positrons annihilate with ambient electrons [13]. The emission is concentrated towards the galactic center, but has wide extension in the bulge and (at lower intensity) along the plane of the galaxy [3].

However, the precise nature and location of the source(s) of the positrons are unclear. The $\gamma $-ray distribution is consistent with a population of discrete sources in the galactic bulge. Microquasars consist of an accretion disk surrounding a stellar-mass black hole and are a possible source because annihilation $\gamma $-rays have been detected from one of them [4]. Alternatively, the positrons could arise from objects such as exploding stars (supernovae) or processes occurring near the galactic center black hole itself (with subsequent diffusion through the interstellar medium before annihilation). More exotic sources could include annihilating dark matter [3]. The main reason for this uncertainty is the coarse angular resolution (${\sim}{2}\deg)$ of the best available $\gamma $-ray observations. It is this problem that we aim to overcome with the present project.

Before annihilation, the positrons and electrons briefly combine to form positronium (Ps) atoms, in one of two spin states. In the singlet state (para-Ps), the electron and positron have anti-parallel spins, and the system undergoes rapid decay into two oppositely directed 511 keV $\gamma $-rays. But the triplet state (ortho-Ps), with parallel spins, has a much greater mean life of 142 ns in the ground state and 3.8 µs in the ${n} = {3}$ quantum state. It decays into three $\gamma $-rays, which results in a wide range of photon energies rather than the identifiable sharp 511 keV line [5]. However, the longer lifetime allows for the possibility of optical transitions from excited states before annihilation [2]. The energy levels for Ps are almost exactly half those of hydrogen, resulting in spectral lines with twice the wavelength of the hydrogen spectrum.

The best candidate for observation using ground-based telescopes is the analog of hydrogen Balmer $\alpha $ (656 nm), with Ps Balmer $\alpha $ ($n = 3 \to 2$) having a wavelength of 1.3122 µm. However, this wavelength lies in the band subject to strong hydroxyl (OH) emission from the Earth’s atmosphere. The emission is variable in time and across the sky, making sky subtraction after detection unsuccessful. Furthermore, the OH lines are so strong that any desired spectral features can be swamped by a nearby OH line in the wings of the instrumental response function (line spread function, LSF). Figure 1 shows the atmospheric emission in the region of interest [6].

 figure: Fig. 1.

Fig. 1. Atmospheric OH emission spectrum [6]. The dashed vertical line shows the wavelength of positronium Balmer $\alpha $ emission for a source with zero radial velocity relative to the Earth. The OH lines are shown with width $\Delta \lambda$ equal to the 0.1 nm expected width of the FBG notches (i.e., $\lambda/\Delta\lambda = 13{,}100$).

Download Full Size | PDF

We propose to overcome this problem using fiber Bragg grating (FBG) technology to suppress the OH lines at high resolution, and then to observe the relevant spectral region with a small spectrograph. Our proposal is based on the photonic integrated multimode spectrograph (PIMMS) concept [7,8], which enables a small inexpensive spectrograph to receive the light from any telescope, while at the same time gaining the advantages from photonic technology, such as FBGs.

Observing in the near IR (NIR) does have a crucial advantage over the visible band, namely, the greater transmission through dust obscuration, which is prominent towards the galactic center. It also benefits from the availability of adaptive optics systems at major telescopes. In any case, there are no Ps lines in the visible range, with its Lyman series lying in the UV.

2. INSTRUMENT CONCEPT

Figure 2 shows a schematic of the proposed system. Objects will be imaged by a telescope on to the multimode fiber (MMF) end of a photonic lantern, which transitions the light into 19 single-mode fibers (SMFs) [911]. In each SMF, a custom FBG will be imprinted, with sharp and deep rejection notches at the wavelengths of the bright OH lines near 1.312 µm. This technology has been demonstrated to successfully suppress the OH emission [1214].

 figure: Fig. 2.

Fig. 2. Schematic of the proposed system. The target is imaged by a telescope to a multimode fiber, which transitions in a photonic lantern to 19 single-mode fibers, each of which has a custom fiber Bragg grating to reject specific wavelengths of OH lines. The SMFs are arranged in a hexagonal pattern at the entrance plane of a small diffraction-limited spectrograph, with spectra imaged on a C-RED-One detector.

Download Full Size | PDF

The advantage of the FBG approach is that the atmospheric OH lines are very narrow and also stable in wavelength—so rejection notches at high resolving power can remove them, while losing a small fraction of the flux of the desired object, which is observed at lower resolving power. The FBGs are temperature stabilized to ensure that the notch wavelengths remain aligned with the atmospheric emission lines. Importantly, the notch profile can be designed with steep sides and a flat bottom, which for this purpose is superior to the typical LSF of a high-resolution spectrograph. For example, Trinh et al. [15] demonstrated notches of 28 dB depth, resolving power 10,400 and square sides. It will be necessary to suppress only some six to 10 OH lines for the present project, which will enable the FBGs to be configured to give narrower notches than would be the case for FBGs seeking to suppress OH over a wider band.

In an OH suppression spectrograph, the resolving power is generally determined by the expected linewidths for the target, and not by the requirement to resolve or separate the OH lines themselves, since these have been suppressed at higher resolution by the FBGs. For Ps forming through radiative combination of positrons with free thermalized electrons, the linewidth Δλ is given by [16]

$$\frac{{\lambda}}{{{\Delta \lambda}}} = 73774.4{T^{- 0.44}}.$$
So for typical interstellar medium temperatures ${T} = {{10}^3} {-} {{10}^5}\;{\rm K}$, the resolving power $\lambda /\Delta \lambda$ required is $R \approx {500 {-} 3500}$. However, in the present case, the Ps line is extremely faint compared to the OH lines, and therefore we wish to ensure the Ps signal is well separated from the residual OH lines. The nearest OH line is at 1.31278 µm, separated by 0.58 nm from Ps Balmer $\alpha $. A resolving power of ${R}\approx {4500}$ ensures a separation of two resolution elements, and will adequately resolve any Ps line.
Tables Icon

Table 1. Components of the Spectrograph

Since the light emerges in SMFs, we can obtain good resolving power with a small diffraction-limited custom spectrograph. A diffraction-limited spectrograph with an SMF input offers two important advantages for this experiment. First, it allows a very stable and compact spectrograph at moderate resolving power (${R}\sim{4000}$) to be built for a relatively low cost [17]. Second, it provides a spectrograph LSF with very low scattering [18], which allows clean sky subtraction—still an important consideration for faint signals even with OH suppression [13]. This latter advantage outweighs the increase in detector noise due to the multiple SMF inputs spreading the light over more pixels than the equivalent MMF injection (as used in previous OH suppression instruments), especially given the low readout noise of the detector.

The detector is a C-RED-One avalanche photodiode (APD) array (www.first-light-imaging.com) [19]. This has the crucial advantage of sub-electron readout noise, meaning it is effectively a photon-counting array. It also has high detection quantum efficiency and low dark noise. The array consists of ${320} \times {256}$ pixels, at a pitch of 24 µm.

The limited number of pixels does influence the design of our instrument: if the fibers were arranged along a linear pseudo-slit, at the minimum feasible separation of ${\sim}{60}\;\unicode{x00B5}{\rm m}$ (due to engineering limitations in the fabrication; assuming etched 80 µm fiber such as the Thorlabs SM1250G80), the slit length becomes 1.14 mm, and fitting the full slit length on the detector places an upper limit on the magnification provided by the collimator/camera combination. The result is that in the spectral direction, the profiles would be severely undersampled, at ${\sim}{1.35}$ pixels/full width at half maximum (FWHM). We avoid this by arranging the fibers in a hexagonal pattern at the spectrograph “slit” plane, in what has become known as the photonic TIGER pattern [20,21]. The result is an effective slit length of 0.3 mm, allowing a higher magnification and optimal sampling of the LSF. This technology has been demonstrated in a compact high-resolution diffraction-limited Raman spectrograph [22] using multicore fiber, but the principle is also valid for a bundle of individual fibers as proposed here. Furthermore, the sensor size limits this approach to a small number of SMF inputs; hence, we use a 19-SMF photonic lantern, although with larger sensors, the number of SMFs could be increased. The 2D arrangement instead of a linear pseudo-slit results in different wavelength ranges on the detector for different fibers, and a consequent reduction in the wavelength range that is common to all fibers. The common range is still sufficient for our purpose, however. Having the individual spectra staggered in this way limits the use of on-chip binning along the spatial axis, but binning is in any case unnecessary with a nearly noise-free detector.

 figure: Fig. 3.

Fig. 3. Front end of the spectrograph. The injection from the fibers is at the left-hand side. The beam encounters the three-element collimator, VPH grating, and the camera lens. The detector is off the diagram to the right. The rays shown are for ${\rm NA} = {0.11}$, from the on-axis fiber.

Download Full Size | PDF

A. Spectrograph Design

In the first instance, we have sought a spectrograph design using only commercial off-the-shelf components. Table 1 lists the selected components. All lenses are anti-reflection coated for 1.05–1.7 µm. The transmission diffraction grating is of volume phase holographic (VPH) type. Figure 3 shows the optical layout. The main design problem is the collimator, which is quite fast: $f/{4.5}$ to the ${{1/e}^2}$ points, but about $f/{3.4}$ to avoid excessive truncation of the numerical aperture $({\rm NA}) = {0.11}$ Gaussian beams [23]. The three-lens combination for the collimator has an effective focal length of 61 mm, giving a nominal beam diameter of 13.5 mm to the ${{1/e}^2}$ points.

Figure 4 shows the detector plane as well as the window and four filters (including cold-stop) of the C-RED-One.

 figure: Fig. 4.

Fig. 4. Back end of the spectrograph, showing the warm sapphire Dewar window and four cooled Hereaeus Infrasil 302 filters of the C-RED-One detector. The detector plane is at the right-hand side. The same rays as in Fig. 3 are shown, with wavelengths 1.299 µm (blue), 1.312 µm (green), and 1.325 µm (red). The cold stop is at filter 2 and has a diameter of 9.5 mm.

Download Full Size | PDF

 figure: Fig. 5.

Fig. 5. Pseudo-color representation of the spectra, showing the full detector area. The white squares show the imaged location of each fiber at the center wavelength of 1.312 µm, in the inclined hexagonal photonic TIGER pattern; numbered squares refer to the rows in Fig. 6 and Table 2. The individual spectra are shown with colors ranging from magenta (1.283 µm) to red (1.342 µm), and with the calculated Gaussian profile in the spatial direction. Detector pixelization is not shown here. The figure shows that the layout does not use the full extent of the detector’s spatial axis.

Download Full Size | PDF

3. SPECTROGRAPH PERFORMANCE

A. Format

Figure 5 shows a simulation of the format of the spectra on the detector. The hexagonal pattern of the fibers at the injection plane has a pitch of 80 µm, and is inclined by 10.9° so that spectra do not overlap [21]. The pattern is magnified by the ratio of the camera and collimator focal lengths, which is 12.3.

B. Image Quality

Figure 6 shows ray-trace spot diagrams at the detector plane, for the center and four extremes of the hexagonal pattern, at the center wavelength plus the extremes of the common range. The actual beams are Gaussian, and will be shown below. But for assessing aberration performance, it is useful to first examine the ray-trace results using a beam with uniform apodization, and ${\rm NA} = {0.11}$. (The actual Gaussian beams will be tapered to 13.5% of peak intensity at the same angle as the edge of the uniform beam, but will then extend somewhat beyond.)

 figure: Fig. 6.

Fig. 6. Matrix spot diagram showing ray-trace results at the center wavelength (middle column) and the two extreme wavelengths of the common range. The five rows show the central fiber and the four with the greatest displacements in both spatial and wavelength directions. The slit plane $xy$ positions of the inputs for each row are shown in Table 2 (and annotated in Fig. 5 at the center wavelength). The circles show the Airy disk radius; the box side is 200 µm. These plots assume rays are uniform within ${\rm NA} = {0.11}$. The individual traced rays have merged into a single circular patch in each case, because the aberrations are minimal.

Download Full Size | PDF

Figure 6 shows all spots are well within the Airy disk, i.e., all images are close to diffraction limited. They show good correction of spherical aberration, but do have a little longitudinal chromatic aberration, as expected from the singlet meniscus lens and camera lens.

Zemax Physical Optics Propagation (www.zemax.com) has been used to find the image properties when using Gaussian beams with NA $= {0.11}$ from the fiber cores. The 2D images in the focal plane were found, then exported to Matlab (www.mathworks.com) where they were projected to the wavelength axis, giving the simulated LSF. Matlab was also used to find the accurate FWHM values for the LSFs.

Figure 7 shows the projected LSFs for two cases—the on-axis fiber at the center wavelength, and fiber no. 3 in Fig. 5 and Table 2, observed near the detector edge at 1.325 µm. As expected from the spot diagrams in Fig. 6, there is very slight broadening for the off-axis case. Both have a clean Gaussian profile.

 figure: Fig. 7.

Fig. 7. Monochromatic profiles of the line spread function as calculated using Gaussian beam inputs, for the center of the detector (green curve, location 2 in Fig. 5) and off-center (blue curve, fiber 3, at $\lambda $ 1.325 µm).

Download Full Size | PDF

Tables Icon

Table 2. Slit Plane $XY$ Positions of Spot Diagram Fibers

The FWHMs are 57.9 µm on axis and 60.2 µm for the other. Thus there is a 4% broadening away from the detector center. The resolving power (λ/FWHM) is 4410 for the on-axis case, and 4240 for the off-axis case. With 24 µm pixels, the detector sampling is 2.4–2.5 pixels/FWHM, a near-ideal balance between accurate characterization of spectral features and wavelength coverage.

C. Wavelength Range

Figure 8 shows the number of fiber spectra available at any given wavelength. Coverage begins at 1.2827 µm, with one fiber, and ends at 1.3420 µm, again with one fiber. Within the common range of 1.303–1.322 µm, all 19 spectra are on the detector. The common range is equivalent to a radial velocity range of ${-}{2070}$ to $+2270\; {\rm km/s}$ for a spectral line with rest wavelength 1.312 µm. This will be adequate for any source within the Milky Way galaxy, but for extragalactic sources such as active galactic nuclei, it extends only to redshift 0.0076. (However, coverage with fewer spectra extends to 6700 km/s or redshift of 0.023.) If desired, the Ps line wavelength of 1.312 µm could be moved towards the shorter wavelength end of the detector by changing the grating and camera tilt angles, so increasing the available range of redshifts, at the expense of coverage shortward of the line.

 figure: Fig. 8.

Fig. 8. Number of individual fiber spectra available at each wavelength. Only 1.303–1.322 µm is common to all 19 fibers.

Download Full Size | PDF

4. INSTRUMENT PERFORMANCE

A. On-Sky Angular Capture

Due to atmospheric turbulence (“seeing”) a ground-based telescope is not in general diffraction limited, and hence cannot efficiently couple a star image into an SMF. For this reason, our system uses an MMF, which transitions to 19 SMFs. We now evaluate the angular size of the telescope’s seeing disk that can be captured by our 19-mode system.

For a Gaussian beam, the far-field beam divergence half-angle (to the ${{1/e}^2}$ point) is

$$\theta = \frac{\lambda}{{\pi {r_0}}}$$
(e.g., [24]), where ${r_0}$ is the mode-field radius. For moderate NAs, our value of 0.11, $\theta $ (in radians) is approximately equal to the NA. For individual SMFs, we therefore have
$${r_0}\,{{\rm NA}_{\text{fib}}} = \frac{\lambda}{\pi}.$$

We will make use of the principle of etendue (${\rm A}\Omega$, area × solid angle) conservation. The ${\rm A}\Omega$ product will be increased over that for a single fiber by a factor equal to the number $N$ of fibers (modes) included. Since the above equations relate to radii and plane angles rather than areas and solid angles, the increase will be a factor of $\surd{N}$.

Consider a telescope with primary mirror diameter ${2}R$, forming star images with angular FWHM $\Gamma $ arscec (the usual measure of astronomical seeing). Seeing disks have approximately Gaussian profiles, so a meaningful comparison of ${\rm A}\Omega$ products can be made.

For a Gaussian, the ${{1/e}^2}$ points are separated by ${1.699} \times {\rm FWHM}$ [23]. Expressing the seeing angle as an effective NA of the telescope,

$${{\rm NA}_{\text{tel}}} = \frac{{1.699{\Gamma}}}{2} \times \frac{1}{{3600}} \times \frac{\pi}{{180}}.$$

Etendue conservation is expressed as

$$\sqrt N {r_0}\;{{\rm NA}_{\text{fib}}} = {{\rm RNA}_{\text{tel}}},$$
which leads to
$${\Gamma} = 7.73 \times {10^4}\frac{{\sqrt N \lambda}}{R},$$
where $\Gamma $ is in arcseconds.

This is the seeing disk diameter that could be efficiently coupled into our system. For ${\rm R} = {2}$ (i.e., a 4 m telescope) and $N = {19}$, we find $\Gamma = 0.22$ arcsec. This is smaller (better) than the natural seeing at most observatories. For an 8 m telescope, good coupling requires seeing of 0.11 arcsec. In both cases, an adaptive optics system would in general be needed to achieve good coupling, but without it, we can still observe using the core region of a natural seeing disk. A larger detector would enable us to use a larger photonic lantern, e.g., 37 SMFs, which would increase the angular capture (field of view) of our instrument.

B. Sensitivity Estimate

We have made preliminary estimates of the signal/noise $({\rm S}/{\rm N}$) ratio obtainable from microquasar point sources of Ps Balmer $\alpha $ radiation, using estimated positron fluxes from [2].

The instrument throughput depends strongly on the coupling efficiency achieved at the injection to the MMF. This in turn depends on the seeing and can be greatly improved through the use of adaptive optics. At the present early design stage, we use indicative figures to estimate the overall system efficiency. Experience with the GNOSIS system [15] shows that throughput of the lantern/FBG is ${\sim}{70}\% ;$ lenses are all anti-reflection coated, with ${\sim}{0.5}\%$ loss, the VPH grating has efficiency ${\sim}{90}\% ,$ and the detector quantum efficiency is at least 70%. Combining with telescope throughput, we have an indicative system efficiency of 0.38 (excluding coupling losses).

The OH suppression improves the S/N by a factor of 1.5–2 for typical sources, although this is sensitive to the relative brightness of the Ps line versus the source continuum. The reason for the improvement is the reduction in the night sky background—the fainter the source, the greater the relative improvement. If a target source had a very bright continuum, the OH suppression would be no advantage, due to the extra losses from the FBG unit, although the diffraction-limited LSF would still provide a major benefit by making sky subtraction more accurate. However, our likely target sources are faint, and will benefit from OH suppression.

Note that we seek to maximize S/N—not simply flux from the source. The optimum field of view therefore depends on the site conditions (i.e., seeing), and the level of adaptive optics correction, as well as the throughput of the photonic lanterns, which is a function of injected focal ratio [25]. The best S/N is achieved by injecting the maximum target light into the smallest possible field of view, so minimizing the night sky contamination of the signal.

With a 4 m telescope, assuming seeing of 0.5 arcsec, and instrument and telescope throughput of 0.38, we obtain S/N in a 10 h observation equal to seven for the brightest example (XTE ${\rm J1118} + {480}$). Figure 9 shows a simulated observation for a similar source and demonstrates the value of OH suppression. For other sources towards the galactic center region, S/N of five to six could be expected in 10 h. The OH suppression is crucial for such objects. As noted, our system would be ideally suited for use behind an adaptive optics facility at an 8 m telescope, which would greatly increase the coupling by matching the delivered image size to our system’s angular capture, without increasing the amount of night sky light. Fainter objects would then be detectable.

 figure: Fig. 9.

Fig. 9. Simulation of a source with an incident Ps Balmer $\alpha $ flux of ${50}\;{{\rm ph/s/m}^2}/\unicode{x00B5}{\rm m}$, with continuum of ${\rm J} = {18}\;{\rm mag}$, in 10 h on a 4 m telescope. The OH suppression notches are assumed to have depth 30 dB and width 100 pm. The result has ${\rm S}/{\rm N} = {8}$ when OH suppressed (black), which is ${2.1} \times$ better than the non-suppressed case (gray). The positronium spectral line is the peak at 1.312 µm, and the imperfect subtraction of the OH lines in the non-suppressed case is clearly detrimental to detection.

Download Full Size | PDF

5. CONCLUSION

We have presented a feasible concept for a small instrument to be attached to a telescope to search for NIR recombination line emission from Ps. This paper outlines the photonic technologies employed, concentrating on the optical design of the proposed spectrograph.

The limited angular capture means that our instrument is suited to detection of recombination radiation from point sources but not truly diffuse regions. However, an upper limit from non-detection can still provide discrimination between point sources and diffuse emission, which is an important test, for both truly diffuse emitters of positrons (such as dark matter) and of positron transport—do positrons annihilate close to their sources or not?

We are working towards completion of the design, to be followed by construction and on-sky use initially at the Anglo-Australian Telescope, which is situated at an ideal latitude for observations of the galactic center region.

Funding

University of Sydney Faculty of Science Dean’s International Postgraduate Research Scholarship; German Federal Ministry for Education and Science (03Z22AI1).

Acknowledgment

We thank the reviewers for comments that have improved the paper.

Disclosures

The authors declare no conflicts of interest.

Data Availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

REFERENCES

1. N. Prantzos, C. Boehm, A. M. Bykov, R. Diehl, K. Ferrière, N. Guessoum, P. Jean, J. Knoedlseder, A. Marcowith, I. V. Moskalenko, A. Strong, and G. Weidenspointner, “The 511 keV emission from positron annihilation in the galaxy,” Rev. Mod. Phys. 83, 1001 (2011). [CrossRef]  

2. S. C. Ellis and J. Bland-Hawthorn, “The search for celestial positronium via the recombination spectrum,” Astrophys. J. 707, 457 (2009). [CrossRef]  

3. S. C. Ellis and J. Bland-Hawthorn, “Astrophysical signatures of leptonium,” Eur. Phys. J. D 72, 18 (2018). [CrossRef]  

4. T. Siegert, R. Diehl, J. Greiner, M. G. H. Krause, A. M. Beloborodov, M. C. Bel, F. Guglielmetti, J. Rodriguez, A. W. Strong, and X. Zhang, “Positron annihilation signatures associated with the outburst of the microquasar V404 Cygni,” Nature 531, 341–343 (2016). [CrossRef]  

5. S. Berko and H. N. Pendleton, “Positronium,” Annu. Rev. Nucl. Part. Sci. 30, 543 (1980). [CrossRef]  

6. S. C. Ellis and J. Bland-Hawthorn, “The case for OH suppression at near-infrared wavelengths,” Mon. Not. R. Astron. Soc. 386, 47 (2008). [CrossRef]  

7. J. Bland-Hawthorn and A. Horton, “Instruments without optics: an integrated photonic spectrograph,” Proc. SPIE 6269, 62690N (2006). [CrossRef]  

8. J. Bland-Hawthorn, J. Lawrence, G. Robertson, S. Campbell, B. Pope, C. Betters, S. Leon-Saval, T. Birks, R. Haynes, N. Cvetojevic, and N. Jovanovic, “PIMMS: photonic integrated multimode microspectrograph,” Proc. SPIE 7735, 77350N (2010). [CrossRef]  

9. S. Leon-Saval, T. A. Birks, J. Bland-Hawthorn, and M. Englund, “Multi-mode fiber devices with single-mode performance,” Opt. Lett. 30, 2545–2547 (2005). [CrossRef]  

10. S. Leon-Saval, A. Argyros, and J. Bland-Hawthorn, “Photonic lanterns: a study of light propagation in multimode to single-mode converters,” Opt. Express 18, 8430–8439 (2010). [CrossRef]  

11. T. A. Birks, I. Gris-Sánchez, S. Yerolatsitis, S. G. Leon-Saval, and R. R. Thomson, “The photonic lantern,” Adv. Opt. Photon. 7, 107–167 (2015). [CrossRef]  

12. S. C. Ellis, S. Bauer, C. Bacigalupo, J. Bland-Hawthorn, J. J. Bryant, S. Case, R. Content, T. Fechner, D. Giannone, R. Haynes, E. Hernandez, A. J. Horton, U. Klauser, J. S. Lawrence, S. G. Leon-Saval, E. Lindley, H.-G. Loehmannsroeben, S.-S. Min, N. Pai, M. Roth, K. Shortridge, L. Waller, P. Xavier, and R. Zhelem, “PRAXIS: an OH suppression optimised near infrared spectrograph,” Proc. SPIE 10702, 107020P (2018). [CrossRef]  

13. S. C. Ellis, J. Bland-Hawthorn, J. S. Lawrence, A. J. Horton, R. Content, M. M. Roth, N. Pai, R. Zhelem, S. Case, E. Hernandez, S. G. Leon-Saval, R. Haynes, S. S. Min, D. Giannone, K. Madhav, A. Rahman, C. Betters, D. Haynes, W. Couch, L. J. Kewley, R. McDermid, L. Spitler, R. G. Sharp, and S. Veilleux, “First demonstration of OH suppression in a high efficiency near-infrared spectrograph,” Mon. Not. R. Astron. Soc. 492, 2796 (2020). [CrossRef]  

14. J. Bland-Hawthorn, S. Ellis, S. Leon-Saval, R. Haynes, M. M. Roth, H.-G. Löhmannsröben, A. J. Horton, J.-G. Cuby, T. A. Birks, J. S. Lawrence, P. Gillingham, S. D. Ryder, and C. Trinh, “A complex multi-notch astronomical filter to suppress the bright infrared sky,” Nat. Commun. 2, 581 (2011). [CrossRef]  

15. C. Q. Trinh, S. C. Ellis, J. Bland-Hawthorn, J. S. Lawrence, A. J. Horton, S. G. Leon-Saval, K. Shortridge, J. Bryant, S. Case, M. Colless, W. Couch, K. Freeman, H.-G. Löhmannsröben, L. Gers, K. Glazebrook, R. Haynes, S. Lee, J. O’Byrne, S. Miziarski, and M. Roth, “GNOSIS: the first instrument to use fiber Bragg gratings for OH suppression,” Astron. J. 145, 51 (2013). [CrossRef]  

16. P. Wallyn, W. A. Mahoney, P. Durouchoux, and C. Chapuis, “The positronium radiative combination spectrum: calculation in the limit of thermal positrons and low densities,” Astrophys. J. 465, 473 (1996). [CrossRef]  

17. C. H. Betters, A. Murray, J. Bland-Hawthorn, and S. G. Leon-Saval, “Precision radial velocities with inexpensive compact spectrographs,” Proc. SPIE 9908, 99081B (2016). [CrossRef]  

18. C. H. Betters, S. G. Leon-Saval, J. G. Robertson, and J. Bland-Hawthorn, “Beating the classical limit: a diffraction-limited spectrograph for an arbitrary input beam,” Opt. Express 21, 26103–26112 (2013). [CrossRef]  

19. P. Feautrier, J.-L. Gach, T. Greffe, F. Clop, S. Lemarchand, T. Carmignani, E. Stadler, C. Doucouré, and D. Boutolleau, “C-RED one and C-RED 2: SWIR advanced cameras using Saphira e-APD and snake InGaAs detectors,” Proc. SPIE 10209, 102090G (2017). [CrossRef]  

20. R. Bacon, G. Adam, A. Baranne, G. Courtes, D. Dubet, J. P. Dubois, E. Emsellem, P. Ferruit, Y. Georgelin, G. Monet, E. Pecontal, A. Rousset, and F. Sayede, “3D spectrography at high spatial resolution. I. Concept and realization of the integral field spectrograph TIGER,” Astron. Astrophys. Suppl. Ser. 113, 347 (1995).

21. S. G. Leon-Saval, C. H. Betters, and J. Bland-Hawthorn, “The photonic TIGER: a multicore fiber-fed spectrograph,” Proc. SPIE 8450, 84501K (2012). [CrossRef]  

22. C. H. Betters, J. Bland-Hawthorn, S. Sukkarieh, I. Gris-Sanchez, and S. G. Leon-Saval, “A multi-core fiber photonic lantern-based spectrograph for Raman spectroscopy,” IEEE Photon. Technol. Lett. 32, 395 (2020). [CrossRef]  

23. J. G. Robertson and J. Bland-Hawthorn, “Compact high-resolution spectrographs for large and extremely large telescopes: using the diffraction limit,” Proc. SPIE 8446, 844623 (2012). [CrossRef]  

24. C. A. Bennett, Principles of Physical Optics (Wiley, 2008).

25. A. Horton, R. Content, S. Ellis, and J. Lawrence, “Photonic lantern behaviour and implications for instrument design,” Proc. SPIE 9151, 915122 (2014). [CrossRef]  

Data Availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (9)

Fig. 1.
Fig. 1. Atmospheric OH emission spectrum [6]. The dashed vertical line shows the wavelength of positronium Balmer $\alpha $ emission for a source with zero radial velocity relative to the Earth. The OH lines are shown with width $\Delta \lambda$ equal to the 0.1 nm expected width of the FBG notches (i.e., $\lambda/\Delta\lambda = 13{,}100$ ).
Fig. 2.
Fig. 2. Schematic of the proposed system. The target is imaged by a telescope to a multimode fiber, which transitions in a photonic lantern to 19 single-mode fibers, each of which has a custom fiber Bragg grating to reject specific wavelengths of OH lines. The SMFs are arranged in a hexagonal pattern at the entrance plane of a small diffraction-limited spectrograph, with spectra imaged on a C-RED-One detector.
Fig. 3.
Fig. 3. Front end of the spectrograph. The injection from the fibers is at the left-hand side. The beam encounters the three-element collimator, VPH grating, and the camera lens. The detector is off the diagram to the right. The rays shown are for ${\rm NA} = {0.11}$ , from the on-axis fiber.
Fig. 4.
Fig. 4. Back end of the spectrograph, showing the warm sapphire Dewar window and four cooled Hereaeus Infrasil 302 filters of the C-RED-One detector. The detector plane is at the right-hand side. The same rays as in Fig. 3 are shown, with wavelengths 1.299 µm (blue), 1.312 µm (green), and 1.325 µm (red). The cold stop is at filter 2 and has a diameter of 9.5 mm.
Fig. 5.
Fig. 5. Pseudo-color representation of the spectra, showing the full detector area. The white squares show the imaged location of each fiber at the center wavelength of 1.312 µm, in the inclined hexagonal photonic TIGER pattern; numbered squares refer to the rows in Fig. 6 and Table 2. The individual spectra are shown with colors ranging from magenta (1.283 µm) to red (1.342 µm), and with the calculated Gaussian profile in the spatial direction. Detector pixelization is not shown here. The figure shows that the layout does not use the full extent of the detector’s spatial axis.
Fig. 6.
Fig. 6. Matrix spot diagram showing ray-trace results at the center wavelength (middle column) and the two extreme wavelengths of the common range. The five rows show the central fiber and the four with the greatest displacements in both spatial and wavelength directions. The slit plane $xy$ positions of the inputs for each row are shown in Table 2 (and annotated in Fig. 5 at the center wavelength). The circles show the Airy disk radius; the box side is 200 µm. These plots assume rays are uniform within ${\rm NA} = {0.11}$ . The individual traced rays have merged into a single circular patch in each case, because the aberrations are minimal.
Fig. 7.
Fig. 7. Monochromatic profiles of the line spread function as calculated using Gaussian beam inputs, for the center of the detector (green curve, location 2 in Fig. 5) and off-center (blue curve, fiber 3, at $\lambda $ 1.325 µm).
Fig. 8.
Fig. 8. Number of individual fiber spectra available at each wavelength. Only 1.303–1.322 µm is common to all 19 fibers.
Fig. 9.
Fig. 9. Simulation of a source with an incident Ps Balmer $\alpha $ flux of ${50}\;{{\rm ph/s/m}^2}/\unicode{x00B5}{\rm m}$ , with continuum of ${\rm J} = {18}\;{\rm mag}$ , in 10 h on a 4 m telescope. The OH suppression notches are assumed to have depth 30 dB and width 100 pm. The result has ${\rm S}/{\rm N} = {8}$ when OH suppressed (black), which is ${2.1} \times$ better than the non-suppressed case (gray). The positronium spectral line is the peak at 1.312 µm, and the imperfect subtraction of the OH lines in the non-suppressed case is clearly detrimental to detection.

Tables (2)

Tables Icon

Table 1. Components of the Spectrograph

Tables Icon

Table 2. Slit Plane X Y Positions of Spot Diagram Fibers

Equations (6)

Equations on this page are rendered with MathJax. Learn more.

λ Δ λ = 73774.4 T 0.44 .
θ = λ π r 0
r 0 N A fib = λ π .
N A tel = 1.699 Γ 2 × 1 3600 × π 180 .
N r 0 N A fib = R N A tel ,
Γ = 7.73 × 10 4 N λ R ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.