# Patent

## Title
Room-temperature quantum noise limited spectrometry and methods of the same

## Abstract

In one embodiment, a heterodyne detection system for detecting light includes a first input aperture configured to receive first light from a scene input, a second input aperture configured to receive second light from a local oscillator input, a broadband local oscillator configured to provide the second light to the second input aperture, a dispersive element configured to disperse the first light and the second light, and a final condensing lens coupled to an infrared detector. The final condensing lens is configured to concentrate incident light from a primary condensing lens onto the infrared detector, and the infrared detector is a square-law detector capable of sensing the frequency difference between the first light and the second light. More systems and methods for detecting light are described according to other embodiments.

## Background
RELATED APPLICATIONS
This application is a continuation of U.S. application Ser. No. 13/076,162, filed Mar. 30, 2011; which claims priority to a U.S. Provisional Patent Application filed Mar. 30, 2010, under Appl. No. 61/319,130, from each of which priority is claimed and which are all herein incorporated by reference.


The United States Government has rights in this invention pursuant to Contract No. DE-AC52-07NA27344 between the United States Department of Energy and Lawrence Livermore National Security, LLC for the operation of Lawrence Livermore National Laboratory.


FIELD OF THE INVENTION
The present invention relates to spectrometry, and more particularly, to quantum noise limited performance in spectrometers.
BACKGROUND
Infrared spectral detection systems are useful in a variety of applications, especially for ground and space based persistent surveillance. Impediments to field deployment of hyper-spectral detection systems include the detector's overall size, weight, and energy usage, which are driven in large part by the requirements of the cooling system for the spectrometer and detector. Therefore, room temperature operation of hyper-spectral detection systems has long been sought to eliminate the vast cooling requirements of conventional hyper-spectral detection systems. However, intrinsic semiconductor detector noise from dark current and radiation from the spectrometer walls presents a significant limitation on system signal-to-noise performance.
One conventional solution to this problem for high-spectral resolution is the use of signal multiplexing, such as Fourier transform spectrometry. This approach detects the entire signal spectral band simultaneously, thereby increasing the competitive position of the signal current with respect to the noise sources. However, even signal multiplexing is overwhelmed by large background flux at elevated temperatures, and at high spectral resolution, this approach fails to achieve quantum noise limited (QNL) performance. QNL performance, where detection sensitivity is limited only by the quantum shot noise of the signal, is a widely used measure with which to determine the performance of signal detection in spectrometry applications.
Another approach uses thermal bolometer detectors for room temperature operation. However, Johnson-Nyquist noise imposes a requirement on these systems for broadband operation, with bandwidths in a range of about 10-90 cm−1. In addition, the thermal response times for bolometers impose limits on the operational speed of a spectrometer.
Another approach has been used to achieve QNL performance by heterodyne detection by using a single frequency laser source. In heterodyne detection, a known light source is combined with an incoming signal on a non-linear detector to produce beat frequencies that are amplified and detected. If the known light source generates more detector signal than the sources of noise, a QNL results. Heterodyne detection has historically depended on sufficiently bright light sources provided by narrow beam lasers. These heterodyne approaches rely on tuning the narrow band laser over the spectral region of interest. This tuning of the laser over the spectral band of interest can considerably reduce the efficiency in gathering spectral information over an extended spectral region, however. Conventional incandescent light sources, while potentially bright enough, are extremely inefficient since a substantial fraction of their output is outside the spectral range of interest.
Another approach illustrated in the prior art is to use an immersion condensing lens to reduce the size of the detector. Because dark current is proportional to detector size, this discriminates against background and dark current. A factor of 16 in focal area reduction has been achieved by VIGO for HgCdTe detectors in the long-wave region through the use of high index immersion lens technology. However this reduction is not enough to overcome dark currents sufficiently to allow QNL performance.
Therefore, all these approaches have failed to overcome the imposed limitations on size, weight, and energy usage imposed by cooling requirements. Even if QNL performance is achievable with hyper-spectral detection systems at elevated temperatures, the size, weight, and energy usage limitations placed on implementation of these devices makes their use impractical. Therefore, a hyper-spectral detection system which could reach QNL performance while achieving practical implementation regarding size, weight, and energy usage would be very beneficial, particularly in military and law enforcement efforts to deploy highly sensitive instruments capable of ground-based persistent surveillance and micro-power space surveillance.
SUMMARY
In one embodiment, a heterodyne detection system for detecting light includes a first input aperture configured to receive first light from a scene input, a second input aperture configured to receive second light from a local oscillator input, a broadband local oscillator configured to provide the second light to the second input aperture, a dispersive element configured to disperse the first light and the second light, and a final condensing lens coupled to an infrared detector. The final condensing lens is configured to concentrate incident light from a primary condensing lens onto the infrared detector, and the infrared detector is a square-law detector configured to sense the frequency difference between the first light and the second light.
In another embodiment, a method for detecting light includes receiving first light from a desired scene in a first input aperture, introducing second light produced by a broadband local oscillator to a second input aperture, passing the first light from the desired scene and the second light produced by the broadband local oscillator through a dispersive element, concentrating the light from the dispersive element onto a detector pixel of an infrared detector using a condensing lens coupled to the infrared detector, and detecting simultaneously the first light and the second light using the infrared detector. The concentrating reduces an intrinsic detector noise resulting from dark current by a factor of at least 20.
Other aspects and embodiments of the present invention will become apparent from the following detailed description, which, when taken in conjunction with the drawings, illustrate by way of example the principles of the invention.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 shows a schematic of a heterodyne spectrometer, according to one embodiment.
FIG. 2A shows a plot of current (electrons/(pixel·sec)) versus wavelength (microns) expected from various sources, according to one embodiment.
FIG. 2B shows a plot of current (electrons/(pixel·sec)) versus wavelength (microns) expected from various sources, according to one embodiment using a condensing superlens.
FIG. 3 shows a schematic of a system including a condensing superlens capable of condensing electromagnetic radiation in a cavity mode, according to one embodiment.
FIG. 4 shows a plot of streamlines for a cavity mode that propagates to the smallest radii, according to one embodiment.
FIG. 5 shows a plot of the change in energy density from a mouth of a cavity mode, according to one embodiment.
FIG. 6 is a flow diagram of a method for detecting light according to one embodiment.



DETAILED DESCRIPTION
The following description is made for the purpose of illustrating the general principles of the present invention and is not meant to limit the inventive concepts claimed herein. Further, particular features described herein can be used in combination with other described features in each of the various possible combinations and permutations.
Unless otherwise specifically defined herein, all terms are to be given their broadest possible interpretation including meanings implied from the specification as well as meanings understood by those skilled in the art and/or as defined in dictionaries, treatises, etc.
It must also be noted that, as used in the specification and the appended claims, the singular forms “a,” “an” and “the” include plural referents unless otherwise specified.
In one general embodiment, a heterodyne detection system for detecting light includes a first input aperture adapted for receiving first light from a scene input, a second input aperture adapted for receiving second light from a local oscillator input, a broadband local oscillator adapted for providing the second light to the second input aperture, a dispersive element adapted for dispersing the first light and the second light, and a final condensing lens coupled to an infrared detector. The final condensing lens is adapted for concentrating incident light from a primary condensing lens onto the infrared detector, and the infrared detector is a square-law detector capable of sensing the frequency difference between the first light and the second light.
In another general embodiment, a method for detecting light includes receiving first light from a desired scene in a first input aperture, introducing second light produced by a broadband local oscillator to a second input aperture, passing the first light from the desired scene and the second light produced by the broadband local oscillator through a dispersive element, concentrating the light from the dispersive element onto a detector pixel using a condensing lens coupled to the infrared detector, and detecting simultaneously the first light and the second light using the infrared detector.
According to one embodiment, broadband infrared spectrometry based on the principle of heterodyne detection that improves spectral light gathering efficiency while maintaining sensitivity may be achieved. To accomplish this, light from a broadband local oscillator source, such as a broadband quantum cascade laser or photonic bandgap light source, may be introduced into a dispersive spectrometer at an entrance aperture (window) that is physically or spatially displaced from an input aperture (window) for receiving the desired scene light. The physical displacement between signal and local oscillator translates, after dispersion, into a physical displacement of their spectra on a square-law detector array. Consequently, each pixel of the array responds to two frequencies creating an intermediate beat frequency in the detector output. In this way, hundreds of individual spectral channels may simultaneously record a high resolution infrared spectrum, substantially improving the throughput of heterodyne spectrometry. Additionally, this approach allows for operation at a significantly higher detector temperature than conventional broadband spectrometers.
A broadband heterodyne approach with emerging bright light sources is promising for facilitating low-noise, elevated-temperature infrared detection. One advantage of using a heterodyne detection approach with a bright light source sufficiently powerful for use in hyper-spectral detection systems is evident from the signal-to-noise expression for heterodyne detection, presented below as Equation 1.

S/N=(Ns1/2NLO1/2)/[(Ns+NLO+Nspec+Ndark+Nother)1/2]  (Eq. 1)

In Equation 1, Ns, NLO, Nspec, Ndark are net signal, local oscillator, spectrometer wall, and dark current photoelectrons, respectively, as collected over the data integration time. This expression exhibits achieving quantum noise limit (QNL) performance when the photoelectrons from the local oscillator exceed other noise sources. NLO then dominates the denominator and cancels itself in the numerator. The expression predicts that elevated-temperature QNL performance is achievable with efficient and bright broadband light sources. This is especially important in the infrared spectral region.
Now referring to FIG. 1, one embodiment of a spectrometer 100 utilizing heterodyne detection is shown, where an input window or aperture 102 for the scene input 104 (source of light to be detected) may be spatially/physically displaced with respect to the input aperture 106 for a bright local broadband oscillator 108. A dispersive element 110, such as a grating, may be used to disperse the signals from the scene input 104 and the bright local broadband oscillator 108 such that the signals are dispersed according to wavelengths. This physical displacement between signal and local oscillator translates, after dispersion, into shifts in the wavelength in the focal plane of the spectrometer 100, e.g., a physical displacement of their spectra. Consequently, each pixel of the array responds to two frequencies creating an intermediate beat frequency in the detector output. An infrared detector 112 registers the displacement shifts in the wavelength. The detector 112 may be an array of photodetectors as known in the relevant art, according to one embodiment, and more preferably is a square-law detector capable of sensing a frequency difference between the light from the aperture 102 and the light from the bright local broadband oscillator 108. The photodetectors may be mixing elements of a type known in the art.
In one approach, where the two frequencies meet spatial and spectral coherency requirements, they coincide at each detector pixel and combine to create a beat frequency with a bandwidth determined by the spectral resolution of the spectrometer 100.
According to various approaches, either, neither, or both apertures 102, 104 may have a lens placed therein for affecting light from the respective input. The lens may be for light shaping, filtering, focusing, concentrating, condensing, etc., and may be of any type known in the art.
In a preferred embodiment, the optical configuration meets an antenna theorem condition such that the product of the aperture area and the collection solid angle is approximately equal to the square of the light source wavelength. Furthermore, spectral coherency is established by the bandwidth of the light intercepted by each pixel, ranging from about 3.0 GHz for a 0.1 cm−1 spectrometer to about 30 GHz for a 1.0 cm−1 spectrometer. Also, the physical displacement at the entrance may preferably be at least twice the spectral bandwidth, in some approaches.
Recently, the development of photonic bandgap light sources and broadband quantum cascade laser light sources makes it practical to consider using broadband light sources in the spectrometer configuration shown in FIG. 1. Sources that emit light in an adjustable band centered on the wavelength region of interest can be bright and efficient. Commercial photonic bandgap (PBG) sources are equivalent to 1000 K blackbodies in bands of several hundred wavenumbers, from at least about 500 cm−1 to about 5000 cm−1. Some virtues of the photonic crystal emitters include targeted spectral ranges and output efficiencies at about 5%. The intensity of these sources is sufficient to outpace spectrometer wall flux at ambient temperature for operation at wavelengths shorter than 5 microns. This relaxes the requirement for cooling the entire spectrometer and requires only detector cooling. The intensities of wall radiation and coherent PBG sources are presented in FIG. 2A. This shows that the PBG source has sufficient intensity to serve as a LO for a room temperature spectrometer in the mid-infrared region.
The above analysis assumes that IR detectors have high frequency operation and sufficiently low dark current noise. Recent developments in band-structure engineered materials such as quantum well infrared photodetectors (QWIPs) have the requisite operating frequencies, up to 75 Ghz, but they have high dark currents that require operation at cryogenic temperatures. More recently, new options have opened up in “quantum engineering” of materials, particularly with the development of quantum dot detectors, quantum dot IR photodetectors (QDIPs), and atomic clusters including “metallic semiconductors,” that are opening up whole new vistas for searching for combinations of high speed and low dark currents. Limiting the location and density of states inhibits one of the major sources of dark current, Auger processes. Quantum clusters that assign angular momentum to electron trajectories further restrict the electron coupling options by conservation rules.
Developments along these lines are being pursued with Type II InAs/GaSb superlattice structures. The dark current of these detector types can be significantly lower than that of QWIPs. One example is shown in FIG. 2A for both 200K operation and 300K operation. While the PBG source is bright enough to operate the spectrometer at ambient temperature, a brighter source is needed to overcome detector current.
The development of broadband quantum cascade infrared lasers, offers a potential robust local oscillator option to implement this approach. A 4μ to 5μ cw laser operating at 285K with 1 Watt output and 5% wall plug efficiency has been described. The 400 cm−1 span provides about 2.5 mWatts per cm−1 wave-number interval. This meets the local oscillator requirement by a wide margin over the detector dark current at ambient temperature. Lowering the overall power consumption to 1 Watt yields about 2×1015 photons/pixel·sec, more than 2000 times higher than detector dark current. This allows for more than 200 spatial pixels creating the opportunity for an imaging spectrometer. Adjusting the well width of the QCL may provide a LO in the more important 3 to 4 micron spectral region. Comparable QCL sources in the long-wave infrared region provide sufficient power to allow room temperature spectrometer operation. However, detector dark current in this spectral region is currently much too high to operate using uncooled detectors.
FIG. 2B illustrates the relative currents for signal, noise sources, and local oscillators expected with the addition of a plasmonic condensing lens (see below) that reduces the detector volume. This reduction directly translates into lower background radiation and lower detector dark current. Using a photonic bandgap broadband source as a local oscillator, a room temperature spectrometer coupled with a detector cooled to 200K is possible. In the long-wave infrared, a room temperature spectrometer with a cryo-cooled detector is possible. In the mid-infrared shorter than 5 microns, the combination of a QLC coupled with a plasmonic condensing lens makes possible QNL operation with a room temperature spectrometer and uncooled detector.
According to one exemplary embodiment, a spectrometer includes a detector, two entrance apertures, and a dispersive element. It may have a resolution of about 1 cm−1 that may provide about a 30 GHz band pass system. The two entrance apertures may have a spatial separation of about 30 microns that translates into an overlap of two colors separated by about 30 GHz at the detector. This matches a frequency response of fast quantum well detectors (about 32 GHz) that are available. These detectors may be fabricated as single elements or arrays exceeding 250×250 pixels and are multiquantum well structures of composition InAs (9)/GaAs (10) with a lattice period of 58 A. Custom design filter/amplifier following circuits may be included, such as a single detector element with external signal processing or an array with on-board pixel circuit elements. Of course, this is just one exemplary embodiment, and many other approaches may be used in designing dimensions and materials of construction for a spectrometer according to various embodiments.
A spectrometer, according to any embodiment described herein, may include diagnostic/performance measuring and analysis capability, either on-board or accessible through an interface with an external device. To assess the performance of broadband spectrometer designs, a model may be developed that includes key physics elements as comparative values. Inputs to the model may include spatial and temporal characteristics of the broadband source and the signal, filter function, spectral dispersion, spatial (frequency) offset of the local oscillator and signal light, and general detector characteristics, such as efficiency, bandwidth, dark current, etc. The model may be embodied as computer readable code, which may be capable of calculating the heterodyne efficiency due to matching of the local oscillator and signal field distributions at each pixel, the base band detector photocurrent in the presence of Poisson shot noise, detector dark current, wall radiation, and other noise sources (as would be apparent to one of skill in the relevant art), and ultimately the signal-to-noise ratio as a function of the fundamental inputs. Performance may be assessed versus a filter spectral bandwidth and shape, as well as the heterodyne offset (intermediate frequency).
Additionally, a model may be produced to predict the improvement expected from using immersion lenses to focus light onto each pixel. The performance model may leverage a time-domain heterodyne detection code for vibrometry applications, to which may be added the finite local oscillator bandwidth associated with the filtered broadband source.
Now referring to FIG. 3, in another embodiment, heterodyne detection may be further improved by placing a condensing superlens or immersion micro-lens 302 in contact with detector pixels 304 such that the superlens 302 receives the scene input 310 and condenses the light, thereby increasing the signal-to-noise ratio (SNR) by a factor of n2, where n is the refractive index of the detector pixels 304 (n can be as large as four). In some embodiments, the superlens 302 may increase the SNR and allow QNL in the long-wave infrared range.
In yet another embodiment, a condensing superlens 302 may be used to reduce dark current (Ndark in Equation 1), which results from detection events in the complete absence of any light source.
In a preferred embodiment, the condensing superlens 302 may be fabricated from silicon (Si) by reactive ion etching, but any formation technique may be used with any material suitable for transmission of an input signal in a desired wavelength through a superlens. According to one embodiment, fields of Si cones 306 may be formed on SiO2 and coated with highly reflective metal 308 as known in the art, such as gold, copper, silver, platinum, any combination thereof, etc. After coating the Si cone 306 with the reflective metal 308, properly formed cones 306 may be selected for focused ion beam (FIB) modification in order to permit transmittance of light therethrough. In a preferred embodiment, the FIB modification may produce a superlens 302 with an aperture diameter approximately one-seventh ( 1/7) that of the wavelength of light that is to be detected by the detector pixels 304. Furthermore, preferred embodiments may place the superlens 302 in direct contact with the detector pixels 304 to improve light absorption therein.
Some exemplary dimensions are indicated in FIG. 3, with a length of the superlens being about 10 μm, a diameter at a receiving end being about 3.5 μm, and a diameter at a transmitting end being about 0.47 μm. Of course, any dimensions may be used, as would be apparent to one of skill in the art, particularly when taking into account a SNR increase desired, a size of individual pixels in the photodetector, the material used, the wavelength of light to be detected, etc.
General derivatives of the preferred embodiment of the condensing superlens 302 may be fabricated from any material known in the art, with the only requirement being that the material be transparent to the wavelength of light to be detected. For example, IR light is transmittable through Si, which makes Si a possible material for a superlens used to condense IR light for detection.
An electromagnetic wave can propagate in a conical cavity until the diameter of the cavity approaches one-half wavelength of light in the propagating medium. The focal spot of a lens would be substantially larger, by twice the F-number, for a wave in a perfectly conducting cavity.
This discussion is based on treatment of cylindrical waveguides and fibers, which extends to a conical geometry. If the field cannot be free of an axial component, the coefficient of this component satisfies a wave-like equation, and from it the remainder of the field can be determined. The resulting solutions can be expressed analytically in terms of Bessel functions. Electric and magnetic fields are assumed to have the time dependence ε=Ee−ωt and β=Be−ωt, respectively.
Maxwell's equations in a uniform nondissipative medium of permittivity ε, permeability μ, wave number k=(με)1/2ω are shown below in Equation 8.

∇×E=iωB ∇·E=0 ∇×B=(k2/iω)E ∇·B=0  (Eq. 8)

However, because of the similarity of the E and B fields, Ψ may represent either value, thus:

∇×∇×Ψ=k2′Ψ  (Eq. 9a)

∇·Ψ=0  (Eq. 9b)

In spherical coordinates (r, ≧, φ) orthonormal unit vectors {e1, e2, e3} are used in the {r, θ, φ} directions, for which Ψ=e1Ψ1+e2Ψe2+e3Ψ3, which is referred to as {“radial,” “axial,” “circular,”}, respectively. Similarly, “transverse” means non-radial on the e2, e3 surface. With s=sin(θ),

∇·Ψ=(1/r2)(r2Ψ1),r+(1/rs)(sΨ2),θ+(1/rs)Ψ3,φ=0  (Eq. 10a)

∇×Ψ=(e1/rs)[(sΨ3),θ−Ψ2,φ]+(e2/rs)[Ψ1,φ−(rsΨ3),r]+(e3/r)[(rΨ2),r−Ψ1,θ]  (Eq. 10b)

Thus, according to Equation 10a which allows Ψ2 and Ψ3 to be replaced with Ψ1, the wave equation for Ψ1 is:

(1/r2)[(r2ψ1),rr+(1/s)(sψ1,θ),θ+(1/s2)ψ1,φφ]+k2ψ1=0  (Eq. 11)

Next a separable solution is sought to the wave Equation 11 having the form:

ψ1∝R1(r)Θ1(θ)Φ1(φ)  (Eq. 12)

This leads to the following eigenvalue equations:

(r2R1),rr+(k2r2−l(l+1)R1=0  (Eq. 13c)

(1/s)(sΘ1,θ),θ+(l(l+1)−m2/s2)Θ1=0  (Eq. 13b)

Φ1,φφ+m2Φ1  (Eq. 13a)

Substituting x=kr, X=x3/2R1, or R1=x−3/2X into Equation 13a yields Bessel's equation:

x2X,xx+xX,x+(x2−l(l+1)−1⁄4)X=0  (Eq. 14)

However, the large value of t in present applications renders the 1⁄4 term negligible, and solutions are the Bessel and Hankel functions, and with ν=[l(l+1)]1/2, the functions have asymptotic forms as x approaches infinity:




  
    
      
        
          
            
              
                
                  X
                  υ
                
                ⁡
                
                  (
                  x
                  )
                
              
              =
              
              ⁢
              
                
                  
                    
                      J
                      υ
                    
                    ⁡
                    
                      (
                      x
                      )
                    
                  
                  →
                  
                    
                      
                        (
                        
                          
                            2
                            /
                            π
                          
                          ⁢
                          
                            
                          
                          ⁢
                          x
                        
                        )
                      
                      
                        1
                        /
                        2
                      
                    
                    ⁢
                    
                      cos
                      ⁡
                      
                        (
                        
                          x
                          -
                          
                            
                              1
                              2
                            
                            ⁢
                            πυ
                          
                          -
                          
                            
                              1
                              4
                            
                            ⁢
                            π
                          
                        
                        )
                      
                    
                    ⁢
                    
                      
                    
                    ⁢
                    
                      υ
                      2
                    
                  
                
                =
                
                  l
                  ⁡
                  
                    (
                    
                      l
                      +
                      1
                    
                    )
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    Y
                    υ
                  
                  ⁡
                  
                    (
                    x
                    )
                  
                
                →
                
                  
                    
                      (
                      
                        
                          2
                          /
                          π
                        
                        ⁢
                        
                          
                        
                        ⁢
                        x
                      
                      )
                    
                    
                      1
                      /
                      2
                    
                  
                  ⁢
                  
                    sin
                    ⁡
                    
                      (
                      
                        x
                        -
                        
                          
                            1
                            2
                          
                          ⁢
                          πυ
                        
                        -
                        
                          
                            1
                            4
                          
                          ⁢
                          π
                        
                      
                      )
                    
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  H
                  υ
                  
                    (
                    1
                    )
                  
                
                ⁡
                
                  (
                  x
                  )
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    
                      J
                      υ
                    
                    ⁡
                    
                      (
                      x
                      )
                    
                  
                  +
                  
                    
                      iY
                      υ
                    
                    ⁡
                    
                      (
                      x
                      )
                    
                  
                
                →
                
                  
                    
                      (
                      
                        
                          2
                          /
                          π
                        
                        ⁢
                        
                          
                        
                        ⁢
                        x
                      
                      )
                    
                    
                      1
                      /
                      2
                    
                  
                  ⁢
                  
                    e
                    
                      i
                      ⁡
                      
                        (
                        
                          x
                          -
                          
                            
                              1
                              2
                            
                            ⁢
                            πυ
                          
                          -
                          
                            
                              1
                              4
                            
                            ⁢
                            π
                          
                        
                        )
                      
                    
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  H
                  
                    (
                    2
                    )
                  
                
                ⁢
                
                  υ
                  ⁡
                  
                    (
                    x
                    )
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    
                      J
                      υ
                    
                    ⁡
                    
                      (
                      x
                      )
                    
                  
                  -
                  
                    
                      iY
                      υ
                    
                    ⁡
                    
                      (
                      x
                      )
                    
                  
                
                →
                
                  
                    
                      (
                      
                        
                          2
                          /
                          π
                        
                        ⁢
                        
                          
                        
                        ⁢
                        x
                      
                      )
                    
                    
                      1
                      /
                      2
                    
                  
                  ⁢
                  
                    e
                    
                      -
                      
                        i
                        ⁡
                        
                          (
                          
                            x
                            -
                            
                              
                                1
                                2
                              
                              ⁢
                              πυ
                            
                            -
                            
                              
                                1
                                4
                              
                              ⁢
                              π
                            
                          
                          )
                        
                      
                    
                  
                
              
            
          
        
      
    
    
      
        (
        
          
            Eq
            .
            
              
            
            ⁢
            15
          
          ⁢
          
            
          
          ⁢
          a
        
        )
      
    
  




Wherein the functions are related by the expression:

Xν′=Xν−1−(ν/x)Xν  (Eq. 15b)

With the implicit eiωt time dependence, the Hankel functions H(1), H(2), represent outgoing, ingoing traveling waves, respectively, and the Bessel functions represent standing waves with equal components of outgoing and ingoing waves.

Since Yν and the Hankel functions diverge at small radii, they cannot be physical solutions to the propagation into an empty loss-free cavity because “loss-free” implies perfect reflection. The Bessel function of the first kind, however, corresponds to a well-behaved standing wave with perfect reflection; so, in a loss-free cavity that reaches a sufficiently small diameter Xν=Jν.
The solutions to Equation 13b are the Legendre polynomials (m=O) and associated Legendre functions (m≠0)ν. The solutions to Equation 13c are trigonometry functions.

Θ1=Plm(c) c=sin θ  (Eq. 16a)

Φ1={ cos(mφ), sin(mφ)}  (Eq. 16b)

Hence, ψ1 is composed of:

(ψ1)lm=ψ1lm{x−3/2Xν(x)}{Plm(c)}{ cos(mφ), sin(mφ)};  (Eq. 17)

    
    
        ψ1lm is a constant.

By analogy, in a waveguide, distinguish modes with transverse electric “TE” fields (E1=0) and transverse magnetic “TM” fields (M1=0). By this construction, Ψ1≠0, so Ψ is the “non-transverse” field, and ∇xΨ is the “transverse” field.

    
    


This choice of modes allows construction of the transverse component of P.

(∇×Ψ)1=(1/rs)[(sψ),θ−ψθ,φ]=0  (Eq. 18)

Which will be satisfied for a scalar field ξ with:

ψ2=ξ,θ, ψ3=(1/s)ξ,φ  (Eq. 19)

∇·Ψ=0 implies

(1/s)(sξ,θ),θ+1/s2)ξ,φφ=−(1/r)(r2ψ1),r  (Eq. 20)

For a separable solution ξ=ξlmRξ(r) Θ,(θ)Φξ(φ)

ξlm=ψ1lm/l(l+1)  (Eq. 21a)

Rξ=(x2R1),x/x  (Eq. 21b)

Θξ=Θ1, Φξ=Φ1  (Eq. 21c)

Thus, the entire field for an eigenmode is

(ψ1)lm=l(l+1)ξlmR1{Plm(c)}{ cos(mφ), sin(mφ)}

(ψ2)lm=ξ,θ=ξlmRξ{Plm(c),θ}{ cos(mφ), sin(mφ)}

(ψ3)lm=ξ,φ/s=ξlmRξ{Plm(c)/s}{−sin(mφ), cos(mφ)}

(∇×Ψ)1lm=0  (Eq. 22d)

(∇×Ψ)2lm=(1/rs)[ψ1,φ−(rsψ3),r]=ξlm{kxR1}{mPlm(c)/s}{−sin(mφ), cos(mφ)}  (Eq. 22e)

(∇×Ψ)3lm=(1/r)[(rψ2)−ψ1,θ]=ξlm{kxR1}{Plm(c),θ}{−cos(mφ),−sin(mφ)}  (Eq. 22f)

R1=x−3/2Xν(x)

Rξ=(x2R1),x/x=x−1/2Xν−1(x)−(ν−1⁄2)R1 

Plm(c),θ=(l+m)(l−m+l)Plm−1(c)−(mc/s)Plm(c)  (vi)

The fields satisfy boundary conditions appropriate for a conductor at the surface of the cone, θ=θ0:

E1(θ0)=E3(θ0)=0; B2(θ0)=0  (Eq. 23)

Which is satisfied where

Plm(cos θ0)=0(TM); Plm(cos θ0),θ=0(TE)  (Eq. 24)

For a narrow cone with θ0<<1, l will be large. In this case,

Plm(cos θ0)=lmJm(lθ)  (Eq. 25)

However, (Equation 13b) may be approximated using s=sin θ≈θ

α2Plm,αα+αPlm,θ+(y2−m2)Plm≈0 y=νθν2=l(l+1)  (Eq. 26)

Which yields a more accurate approximation in the form of Bessel's equation:

Plm(c)≈νmJm(y)

mPlm(c)/s≈νm+1mJm(y)/y 

Plm(c),θ≈νm+1mJm′(y)=νm+1[Jm−1(y)−mJm(y)/y]  (Eq. 27a)

In this small θ approximation, for TE

E1=0

E2=ExR1 mJm(y)/y{−sin(mφ), cos(mφ)}

E3=ExR1Jm′(y){−cos(mφ),−sin(mφ)}

B1=BνR1Jm(y){ cos(mφ), sin(mφ)}

B2=BRξ Jm′(y){ cos(mφ), sin(mφ)}

B3=BRξ mJm(y)/y{−sin(mφ), cos(mφ)}

R1=x−3/2Xν(x)

Rξ=(x2R1),x/x=x−1/2Xν−1(x)−(ν−1⁄2)R1 

B=νm+1ξlm E=(iω/k)B  (Eq. 27b)

For TM, B replaces E in Equation 27b, and the eigenvalue ν2=l(l+1) satisfies the boundary conditions at θ0:

νnm=ynm/θ0  (Eq. 28)

where ynm is the nth zero of Jm′(TE) or Jm(TM). This simple relation between ν and θ0 motivates application of the small-angle limit.








TABLE 1







The Zeroes of Jm' and Jm








ynm (TE)
ynm (TB)














n/m
0
1
2 
n/m
0
1
2

















0
3.832
1.841
3.054
0
2.045
3.832
5.136


1
7.016
5.331
6.076
1
5.520
7.016
8.417


2
10.173
8.536
9.969
2
8.864
10.173 
11.620









In one embodiment, the mode TE11, for which ν=1.842/θ0, and TM10, for which ν=2.405/θ0, since such modes propagate to the smallest radii. Where the cavity is driven by a linearly polarized external field, only the m=1 mode(s) will be excited. Streamlines for the lowest mode TE01 will overlap with the external field and comprises the cavity mode that propagates to the smallest radii, as depicted in FIG. 4.
As an indicator of propagation into the cavity, use the energy density in the electromagnetic field. The volume energy density is:

u=1⁄2(ε|E|2+|B|2/μ)  (Eq. 29)

Wherein the linear density may be expressed by:

dU/dr=r2∫dθ sin θ∫dφu≈(x/kv)2∫dyu∫dφu  (Eq. 30)

From Equation 20, the density component for single modes may be evaluated as:

dU/dR≈u*/(νk)2{2π:m=0,π:m≠0}Inm(2/π)ν(x)  (Eq. 31a)

u*=1⁄2ξE2=1⁄2 B2/μ  (Eq. 31b)

Inm=∫Ynm[Jm′2+(mJm/y)2]ydy  (Eq. 31c)

ν(x)=(π/2)[(x2+fν2)|xR1|2+x2|Rξ|2]  (Eq. 31d)

F=∫YnmJm2ydy/Inm=1  (Eq. 31e)

Where energy density of a given mode at the mouth of cavity is known, the change in energy density from the mouth of the cavity is contained within ν(x). Thus for a standing wave X=J with perfect reflection from the cavity, the energy density function is as shown in FIG. 5.
In any of the embodiments described above, the broadband oscillator may produce a range of wavelengths (e.g., broadband quantum cascade laser), as opposed to a single-frequency laser. In preferred embodiments, the broadband oscillator may further reduce energy usage by including the ability to limit the range of frequencies produced that correspond to a desired detection range.
Now referring to FIG. 6, a method 600 for detecting light is shown according to one embodiment. The method 600 may be carried out in any desired environment, including those shown in FIGS. 1 and 3, among others. Of course, the method 600 may include more or less operations than those shown in FIG. 6, as would be apparent to one of skill in the art upon reading the present descriptions.
In operation 602, first light is received from a desired scene in a first input aperture.
In operation 604, second light produced by a broadband local oscillator is introduced to a second input aperture. In one embodiment, the local broadband oscillator may include a broadband quantum cascade laser having broad spectral bands between about 500 cm−1 and about 5000 cm−.
In one embodiment, the second light from the broadband local oscillator may be limited to a frequency that corresponds to a desired detection range, which may reduce power usage of the broadband local oscillator.
In one preferred approach, the first input aperture may be spatially offset from the second input aperture. The spatial offset may be based on a desired frequency and/or wavelength detection range, in some approaches.
In operation 606, the first light from the desired scene and the second light produced by the broadband local oscillator is passed through a dispersive element.
In operation 608, the light from the dispersive element is concentrated onto a detector pixel using a condensing lens coupled to the detector. In one embodiment, the condensing lens may include a plasmonic IR condensing lens, as previously described.
According to some approaches, the plasmonic IR condensing lens may have a conical shape tapering inward toward the sensor, and may include a Si core (or any other material as known in the art that is compatible with the light being concentrated) having a conductive metal cladding on radial surfaces thereof. The conductive metal cladding may be any material which operates to concentrate the light within the plasmonic lens, as known in the art, such as Au, Ag, Cu, etc. The conductive metal cladding may also be selected based on a gain provided by the conductive metal cladding, in more approaches.
In one embodiment, an angle measured from perpendicular to a planar surface of the sensor and a sidewall of the conical shape may be between about 5° and about 12°, for example about 8.5°.
According to another embodiment, the conical shape may have a first diameter at a distal end adjacent the detector of less than about one-half wavelength of incident light in Si. Also, the conical shape may have a second diameter at an end receiving the incident light from the dispersive element of approximately a diffraction limit set by the primary condensing lens. Of course, the dimensions may be chosen to increase the likelihood of light being concentrated by the plasmonic lens, according to techniques that would be apparent to one of skill in the art upon reading the present descriptions.
In operation 610, the first light and the second light are simultaneously detected using the detector. In one approach, the detector may include an array of photodetectors, as previously described.
In one preferred approach, the array of photodetectors may be sufficient in number to detect a desired frequency range with a desired resolution, in order to conserve energy usage of the array of photodetectors, among other advantages.
According to one embodiment, a frequency of the first light and a frequency of the second light may coincide at each pixel of the detector with a bandwidth determined by a spectral resolution of a spectrometer.
In an optional operation, a coherent beat frequency may be created using the two frequencies, wherein each area of the first and second input apertures times a collection solid angle of the spectrometer is about a wavelength squared (about λ2).
While various embodiments have been described above, it should be understood that they have been presented by way of example only, and not limitation. Thus, the breadth and scope of a preferred embodiment should not be limited by any of the above-described exemplary embodiments, but should be defined only in accordance with the following claims and their equivalents.

## Claims

1. A heterodyne detection system for detecting light, comprising:
a first input aperture configured to receive first light from a scene input;
a second input aperture configured to receive second light from a local oscillator input;
a broadband local oscillator configured to provide the second light to the second input aperture;
a dispersive element configured to disperse the first light and the second light; and
a final condensing lens coupled to an infrared detector,
wherein the final condensing lens is configured to concentrate incident light from a primary condensing lens onto the infrared detector,
wherein the infrared detector is a square-law detector configured to sense a frequency difference between the first light and the second light; and
wherein each area of the first and second input apertures times an input collection solid angle of the system is about a wavelength squared (about λ2).

2. The heterodyne detection system as recited in claim 1, further comprising a radio frequency (RF) waveguide, wherein the RF waveguide is configured to facilitate optical upconversion of RF energy generated by the infrared detector.

3. The heterodyne detection system as recited in claim 1, wherein at each detector pixel, a frequency of light from the scene input and a frequency of light from the oscillator input coincide with a bandwidth determined by a spectral resolution of the system, and
wherein the detection system is configured such that the two frequencies create a coherent beat frequency.

4. The heterodyne detection system as recited in claim 1, further comprising a radio frequency (RF) waveguide, wherein the RF waveguide is coplanar to the infrared detector.

5. The heterodyne detection system as recited in claim 1, wherein the infrared detector comprises:
a photodetector element; and
an RF detector element decoupled from the photodetector element.

6. The heterodyne detection system as recited in claim 5, wherein the photodetector element comprises a quantum well infrared photodetector (QWIP), and wherein the RF detector element comprises a carbon nanotube antenna.

7. The heterodyne detection system as recited in claim 1, wherein the broadband local oscillator comprises a broadband quantum cascade laser having broad spectral bands between about 500 cm−1 and about 5000 cm−1.

8. The heterodyne detection system as recited in claim 1, wherein the first input aperture is spatially offset from the second input aperture, wherein the first and second light pass through the first and second apertures in a same direction relative to each other.

9. A heterodyne detection system for detecting light, comprising:
a first input aperture configured to receive first light from a scene input;
a second input aperture configured to receive second light from a local oscillator input;
a broadband local oscillator configured to provide the second light to the second input aperture;
a dispersive element configured to disperse the first light and the second light; and
a final condensing lens coupled to an infrared detector;
wherein the final condensing lens is configured to concentrate incident light from a primary condensing lens onto the infrared detector;
wherein the infrared detector is a square-law detector configured to sense a frequency difference between the first light and the second light;
wherein the final condensing lens comprises a plasmonic infrared (IR) condensing lens;
wherein the plasmonic IR condensing lens has a conical shape tapering inward toward the infrared detector, and comprises a Si core having a conductive metal cladding on radial surfaces thereof; and
wherein the conical shape has a first diameter at an end adjacent the infrared detector of less than about one-half wavelength of the incident light in Si, and wherein the conical shape has a second diameter at an end receiving the incident light from the dispersive element of approximately a diffraction limit set by the primary condensing lens.

10. A method for detecting light, the method comprising:
receiving first light from a desired scene in a first input aperture of a spectrometer;
introducing second light produced by a broadband local oscillator to a second input aperture;
passing the first light from the desired scene and the second light produced by the broadband local oscillator through a dispersive element;
concentrating the light from the dispersive element onto a detector pixel of an infrared detector using a condensing lens coupled to the infrared detector; and
detecting simultaneously the first light and the second light using the infrared detector,
wherein the concentrating reduces an intrinsic detector noise resulting from dark current by a factor of at least 20.

11. The method as recited in claim 10, further comprising:
generating radio-frequency (RF) energy using the infrared detector in response to detecting the first light and the second light; and
optically upconverting the generated RF energy.

12. The method as recited in claim 11, further comprising creating a coherent beat frequency using the two frequencies, wherein each area of the first and second input apertures times a collection solid angle of the spectrometer is about a wavelength squared (about λ2).

13. The method as recited in claim 11, wherein the upconverting comprises:
channeling the generated RF energy to a close-coupled RF waveguide; and
interacting the channeled RF energy with a laser beam, wherein the interacting generates an optical electric field.

14. The method as recited in claim 13, wherein the interacting is performed at an electro-optic modulator coupled to an RF waveguide.

15. The method as recited in claim 10, wherein the broadband local oscillator comprises a broadband quantum cascade laser having broad spectral bands between about 500 cm−1 and about 5000 cm−1.

16. The method as recited in claim 10, wherein the first input aperture is spatially offset from the second input aperture, and
wherein the first and second light pass through the first and second apertures in a same direction relative to each other.

17. The method as recited in claim 10, wherein the condensing lens comprises a plasmonic infrared (IR) condensing lens.