Infrared upconversion spectrometer for the mid-IR range

10753795 · 2020-08-25

Assignee

Inventors

Cpc classification

International classification

Abstract

The invention provides an infrared upconversion spectrometer for determining a mid-IR spectrum of received infrared light with a high resolution. The spectrometer applies upconversion to transform light in the mid-IR to the near-IR range where efficient detectors are available. The up-conversion causes divergence of the light, and in addition, the invention applies an extra dispersive element to record a spectrum.

Claims

1. An infrared upconversion spectrometer for determining a mid-IR spectrum of received infrared light, wherein the spectrometer comprises: a nonlinear material and first optical element configured to couple the received infrared light into the nonlinear material; a laser configured to couple into the nonlinear material in at least substantially the same direction of propagation as, and in spatial overlap with, in-coupled infrared light; where the nonlinear material, the first optical element, and the laser are configured for non-collinear phase-matching in-coupled mid-IR light and upconverted light to be transmitted by the nonlinear material; and a spatially resolved detector with multiple pixels being sensitive to upconverted light; the infrared upconversion spectrometer further comprising: imaging optics arranged between the nonlinear material and the detector to image the spatial overlap between the laser light and the in-coupled infrared light inside the nonlinear material onto the detector; and a dispersive element configured to disperse upconverted light so that different wavelength sub-ranges of the upconverted light are imaged onto different pixels of the detector; and in that a spectral resolution of the spectrometer is determined by the laser, the dispersive element and a virtual slit being a spatial overlap between the laser light and the in-coupled infrared light inside the nonlinear material.

2. The infrared upconversion spectrometer according to claim 1, wherein the spatial overlap between the laser light and the in-coupled infrared light inside the nonlinear material is at least 1 mm long in the direction of propagation.

3. The infrared upconversion spectrometer according to claim 1, wherein the imaging optics comprises a diverging lens followed by a converging lens for expanding and collimating the upconverted light, respectively; wherein the dispersive element is positioned to receive the collimated upconverted light, and wherein the detector is a line-detector.

4. The infrared upconversion spectrometer according to claim 1, wherein the imaging optics comprises a cylindrical lens or a concave cylindrical mirror for focusing the dispersed upconverted radiation on the detector, and wherein the detector is a 2-dimensinal detector.

5. The infrared upconversion spectrometer according to claim 1, wherein the laser is an asymmetric laser field having an asymmetric emitter area with an aspect ratio of more than 3.

6. The infrared upconversion spectrometer according to claim 5, wherein the nonlinear material is a nonlinear crystal with a length (z-axis) to thickness (x-axis) ratio of more than 20.

7. The infrared upconversion spectrometer according to claim 1, wherein the nonlinear material is poled at a two different superposed frequencies.

8. The infrared upconversion spectrometer according to claim 1, further comprising a means for preventing the upconverted light having the largest angles in the x-axis direction from being dispersed by the dispersive element.

9. A method for determining a mid-IR spectrum of received infrared light, comprising: providing an infrared upconversion spectrometer comprising: a nonlinear material and first optical element configured to couple received infrared light into the nonlinear material; a laser configured to couple into the nonlinear material in at least substantially the same direction of propagation as, and in spatial overlap with, in-coupled infrared light; wherein the nonlinear material, the first optical element, and the laser are configured for non-collinear phase-matching in-coupled mid-IR light and upconverted light to be transmitted by the nonlinear material; and a spatially resolved detector with multiple pixels being sensitive to upconverted light; receiving infrared light and coupling received infrared light and the laser into the nonlinear material whereby in-coupled infrared light and the laser light generates upconverted light via a nonlinear optical process in the nonlinear material, the upconverted light forming a concentric pattern with radii as a function of wavelength; collimating the upconverted light transmitted by the nonlinear material; dispersing the collimated light; and focusing dispersed light onto the detector; wherein a spectral resolution of the spectrometer is determined by the laser, the dispersive element and a virtual slit being a spatial overlap between the laser light and the in-coupled infrared light inside the nonlinear material.

Description

BRIEF DESCRIPTION OF THE FIGURES

(1) The infrared upconversion spectrometer according to the invention will now be described in more detail with regard to the accompanying figures, wherein same reference numerals are used to designate identical or corresponding features. The figures show one way of implementing the present invention and is not to be construed as being limiting to other possible embodiments falling within the scope of the attached claim set.

(2) FIG. 1A illustrates an embodiment of the invention, and 1B is a close-up of the nonlinear material.

(3) FIGS. 2-6 illustrates different layouts of the spectroscopy part of the spectrometer according to the invention.

(4) FIGS. 7 and 8 illustrates different layouts of embodiments where the upconversion occurs intracavity in the laser.

(5) FIGS. 9 and 10 illustrates different embodiments of a multichannel spectrometer according to embodiments of the invention.

(6) FIGS. 11A and B are calculated curves showing the phase matched upconverted wavelength as a function of angle of the upconverted light relative to the upconversion laser for angles external to the nonlinear material (11A) and angles internal in the nonlinear material (11B).

(7) FIGS. 12A and B are plots of the intensity distribution of upconverted light when propagated (back or forth) to the center of the non-linear crystal illustrated in FIG. 12C.

(8) FIGS. 13A-C are contour plots showing the phase match curves as they can appear in the plane of the grating. The grating considered here is a 25 mm grating placed at an angle of 60 to the Z-axis. In this example the grating acts as a spatial aperture removing larger angles in the x-direction, whereas larger angles in the y-direction

DETAILED DESCRIPTION OF AN EMBODIMENT

(9) FIG. 1A illustrates part of an embodiment 10 of the infrared upconversion spectrometer according to the invention used in the following for explanatory purposes. Here, solid and dashed lines represent different wavelengths .sub.1 and .sub.2 of received polychromatic infrared light, and are used to show and explain the different radiation paths for different wavelengths.

(10) The first optical element for coupling the received infrared light 1 into the nonlinear material 3 is here embodied by a spherical lens 2. Preferred choices for the nonlinear material are various nonlinear crystal, such as periodically-poled Lithium Niobate (PP:LN), aperiodically-poled Lithium Niobate (AP:LN), AgGaS.sub.2, AgGaSe.sub.2, OP:GaAs, OP:GaN, BNA, and DAST. Inside the nonlinear material, the incoupled light is overlapped with a beam 9 (see FIG. 1B) from a laser 11 coupled in via mirror 12, and with proper orientation or poling of the nonlinear material, a nonlinear optical process such as sum frequency generation (SFG) converts the incoupled light to upconverted light 13 having a shorter wavelength.

(11) The angular dependent (non-collinear) phase-matching relied upon in the present invention means that the phase matching condition will be fulfilled for different angles of incidence, , for different wavelengths. Thereby, for say wavelength .sub.1, only light having a large angle of incidence will be phase matched, whereas for .sub.2 (<.sub.1), the phase matching condition will be fulfilled at smaller angles of incidence. Light with wavelength .sub.1 but a small angle of incidence will not fulfil the phase matching condition and will thus not be upconverted. The spatial distribution of the resulting upconverted light will therefore form a concentric pattern with radii as a function of wavelength, i.e. each constituent frequency will emerge with a different angle, similar to dispersion in a prism or a grating.

(12) Adding to this is the conservation of momentum which means that the angle of the generated upconverted non-collinear photons scales with the relative shift in wavelength compared to the incoming mid-IR light. This will narrow the cone of the concentric pattern.

(13) The upconversion process can take place through the entirety of the overlap inside the nonlinear material, here equal to the length L of the nonlinear material. This means that a photon converted near the input can travel along with its new angle as compared to a photon with the same wavelength converted near the exit of the nonlinear material.

(14) The upconverted light 13 is transmitted by the nonlinear material, and the overlap between incoupled light and laser inside the nonlinear material should be imaged onto the detector 8. In the embodiment shown in FIG. 1A, this imaging is performed by imaging optics comprising lenses 4 and 7. In the infinity corrected plane (or Fourier plane) between the lenses, the dispersive element in the form of grating 6 with the grating lines in the Y direction is positioned. The dispersion causes the rings corresponding to different wavelengths to be displaced differently in the X-direction, and the focusing by lens 7 results in a line of points, each point corresponding to a ring of different wavelength. The referral to rings and points corresponding to a distinct wavelength is used for the purpose of illustration only. Typically, the received infrared light will contain a continuous wavelength spectrum and the ring pattern and line of points will be continuous features as well.

(15) An example of a line-detector could be e.g. the Xlin detector series from Xenics includes multiplexed InGaAs line arrays with 1024 or 2048 pixels and a 12 m pitch. The detectors are equipped with one-stage Peltier cooling; three-stage cooling is also available on an optional basis. With the 1024-pixel detector, a line rate of 40 kHz can be achieved; whereas with the 2048-pixel detector, this value is 10 kHz.

(16) The infrared upconversion spectrometer of the invention consists overall of an infrared signal input, represented by a set of mixed wavelengths, an upconversion module that converts the infrared signal to higher energy levels, and a spectroscopic part that disperse the upconverted signal further and measure the specific amount of each spectral component. In the following, a number of different layouts for spectroscopic part will be described in relation to FIGS. 2 through 6. In these Figures, the input set of wavelengths are upconverted into a set of upconverted wavelengths represented by full, dashed and dotted lines respectively. The illustrations track the trajectories of each wavelength.

(17) In the layout 20 illustrated in FIGS. 2, 1 and 5 are optional filter elements to separate the signal components from noise and background. 2 is a focusing element to match the incoming infrared signal with the upconversion process in the nonlinear material 3. The upconverted light is collimated by the focusing element 4 and diffracted by the grating 6. The diffracted light is collected by the focusing element 7 and detected by detector 8.

(18) The layout 30 illustrated in FIG. 3 is similar to 20 described above, with a different configuration of the focusing element for the upconverted signal. In this configuration a diverging element, 14, together with a converging element 4, expands and collimate the upconverted signal before it reaches the grating. This results in a more compact construction of the spectrometer

(19) The layout 40 illustrated in FIG. 4 is similar to 30 described above, with a general dispersive element, 6. This could be a reflection grating as in 30, or other such as a transmission grating, a prism or a holographic element.

(20) The layout 50 illustrated in FIG. 5 is similar to 20 described above, where the focusing elements for the upconverted signal are included in the dispersive element, 6.

(21) The layout 60 illustrated in FIG. 6 is similar to 50 described above, with a 4f-scaling of the upconverted signal by focusing elements 4 and 7. The dispersive element, 6, includes further a design that allow for optimal flat field detection at detector 8.

(22) In the following, different layouts for intracavity upconversion are described in relation to FIGS. 7 and 8.

(23) FIG. 7 illustrates a layout having three high reflective mirrors, 16, 17, and 18, that together create a laser cavity 70. Inside the laser cavity a gain medium, 19, is placed such that the gain support lasing inside the cavity. The laser crystal 19 is pumped by a diode laser through the dichroic mirror, 18. A focusing element 21 scales the field inside the laser cavity the nonlinear crystal 3 converts the incoming infrared signal to a higher photon energy output through a nonlinear sum frequency process between the intracavity laser field and the incoming infrared light. The infrared light is coupled in through the dichroic mirror 16 and out through the dichroic mirror 17.

(24) The layout of the laser cavity 80 illustrated in FIG. 8 is similar to 70 described above, but the infrared light is coupled to the crystal with an intracavity parabolic mirror 22, where the laser beam pass through a hole in the parabolic mirror. The upconversion output is done through the dichroic mirror 16.

(25) FIG. 9 shows an embodiment of a multi-channel infrared upconversion spectrometer 90 using a technique which may be useful in combination with other embodiments of the invention. Herein, an angular dependence of the phase-matching conditions within the nonlinear material 3 is used to form different up-conversion channels. This is achieved by propagating the in-coupled infrared light along multiple paths within the nonlinear material 3. In this way, different spectral ranges of the infrared light are matched by the different phase-matching conditions along the different paths of the in-coupled light. Then, the nonlinear material 3 is rotatably mounted, so that a rotation of the nonlinear material according to arrow 91 directly changes the angle. By rotating the nonlinear material 3 between e.g. two rotational positions, different up-conversion channels may be formed sequentially in time. The number of channels may easily be increased by using more rotational positions. In this embodiment, the demultiplexer is the rotational mount of the nonlinear material 3, which may selectively couple the infrared light into a first up-conversion channel, a second up-conversion channel, and any other up-conversion channels. Detection of light from the respective up-conversion channels may then be performed by time-gated detection while the nonlinear material 3 dwells in the respective rotational position. Thus, light incident on the detector (not shown in this figure) at any given time will only be the up-converted light corresponding to one channel. The skilled person will realize that a continuum of channels may be achieved in this way by scanning the rotational position of the nonlinear material 3.

(26) FIG. 10 illustrates another embodiment of a multi-channel infrared upconversion spectrometer 100 that may be useful in combination with embodiments of the infrared spectrometer according to the invention. In this case, the nonlinear material 101 comprises multiple regions (here three are illustrated, in the form of periodically poled regions 102, 103, 104), each having different phase-matching conditions. Each region may correspond to a single up-conversion channel and be selectable by transverse translation of the nonlinear material 101 according to arrow 105. In this way, each channel may be formed one at a time. Thus, the demultiplexer is the translational mount of the nonlinear material 101 in this embodiment.

(27) The following is an example of how to calculate the expected spectral resolution in an infrared upconversion spectrometer according to the invention, and also contains a lot of considerations relating to important design parameters which may be generalized to other embodiments of the invention.

(28) The upconverted wavelength will vary according to phasematch relations. In particular, non-collinear phasematch at varying angles to the laser will lead to different wavelengths being phasematched at different angles. An example calculation of a periodically poled Lithium Niobate crystal gives the relation shown in FIG. 11A expressing the crystal external angles of the upconverted infrared light.

(29) It is important to note that the different output angles corresponds to different input angles as well (conservation of transverse momentum). The above example is performed for PPLN crystal at 60 C. with a poling period of 22 m, mixing with a 1.064 m laser source.

(30) We note that most of the wavelength coverage is for angles between 2 and 3 degrees, and remember that it is (almost) rotationally symmetric (neglecting a small variation from change in refractive index).

(31) Virtual Slit Size Estimation

(32) In the following, we do an example calculation on how to determine the effective virtual slit size from theoretical considerations and later how a system can be modified to minimize the virtual slit size. The calculation is based on an example relating to FIGS. 11A, 11B, 12A, 12B and 12C.

(33) Firstly, FIG. 12C illustrates a nonlinear crystal 3 of length L, in which a mixing laser field 9 and a plane mid-infrared signal are overlapped. The mid-infrared signal is a plane wave extending over the full crystal and is not illustrated here. The virtual slit 15, is computed from the shape of the mixing laser, the incoming field and the collective shape of the upconverted light 13 at a specific angle as seen from the center part of the crystal in the z-direction. I.e. a convolution of a top hat and a Gaussian distribution, that is the width of the virtual slit is a combination of the width of the mixing laser and the angular spread of the upconverted signal in the x-direction.

(34) The angle of the upconverted light inside the crystal plays a crucial role in determining the virtual slit size in the crystal. In the photon efficiency calculation we assume that the infrared light being upconverted is incoherent light coming from e.g. a large thermal source such as a hot filament, which is able to fill the non-linear crystal with infrared light from all angles in the entire volume of the laser. Since the angles of the infrared light are quite large, it is common to have total internal reflections on the facets in the x and/or y directions. This helps significantly in filling the whole length of the crystal with infrared light.

(35) An external angle at around 2 degrees correspond to an internal crystal angle of about 1 degree (since the refractive index of PPLN is about 2). An accurate calculation gives the curves in FIGS. 11A and B showing the phase matched upconverted wavelength as a function of angle of the upconverted light relative to the upconversion laser for angles external to the nonlinear material (11A) and angles internal in the nonlinear material (11B).

(36) The calculation is performed for a PPLN crystal with a poling period of 22 m and a temperature of 60 C. The plot shows the angle of optimal conversion efficiency, although each wavelength will be upconverted in a range of angles near this maximum. The angles external to the crystal (11A) are relevant for designing size and focal length of lenses/mirrors and gratings. The internal crystal angle (11B) plays an important role in determining the size of the virtual slit in the center of the crystal.

(37) In this example we assume a crystal length of 20 mm. Light generated in the beginning of the crystal at an angle of 1 degree will move sideways to the center of the crystal with: 10 mm*Tan(1)=174 m. If we for example have a radius of the laser of 90 m, we will then achieve an effective slit width dominated by the sideways translation of the upconverted light. However, to the untrained eye, this relation may seem trivial, but there are complexities arising from phase match conditions. If we consider the exactly phase matched process, the field generated in the beginning of the crystal will be in phase with the field generated everywhere else in the crystal. However, much of the contribution of upconverted light is for non-perfect phasematching. This means that the field in the virtual slit will add with different phase, dependent on how far it has been (back or forth) propagated within the crystal. This will lead to non-trivial virtual field distribution in the center of the crystal (virtual since some of the field hasn't been generated yet, and only exists in the center of the crystal as a mathematical abstraction).

(38) The virtual field (converted to intensity) can look like the profiles shown in FIGS. 12A and B for a phasematched (12A) and somewhat non-phasematched (12B) case (k*L=7.6), for a laser with w.sub.0=90 m.

(39) Here, a single wavelength is modelled as a plane wave coming from a single non-collinear angle that may be phase matched (FIG. 12A), or not perfectly phase matched (12B). The plots show example calculation for a laser spot radius of 90 m, and an internal angle of 1, a crystal length of 20 mm. The longer the upconverted crystal is, and the larger the internal angle of propagation is, the larger the virtual slit will appear to be. It is important to note that only the field propagated forward to the center of the crystal has a physical presence in the virtual slit, as opposed to the field generated by SFG later in the crystal. When calculating the total virtual slit size for a given wavelength, the intensity for all angles and corresponding phase match conditions should be added. If the IR source is incoherent, such as light from a hot filament, it is the intensities that should be added. If the IR source is a spatially coherent source, the fields should be added before calculating the generated field. FIG. 12A, is calculated as the convolution of a tophat function and a Gaussian function, whereas FIG. 12B is the same convolution but added with different phase contribution depending on the z-coordinate from which it stems, as outlined in FIG. 12C. The tophat function originates from the fact that the light is generated at different angles and z-positions in the crystal and the Gaussian distribution is the product of the incoming plane infrared wave and the mixing laser.

(40) The total amount of light is an integration of all the individual intensity distribution over all incident angles with correspondingly varying phasematch condition. This integration should be done separate for each wavelength under consideration.

(41) When all these intensity distributions are added together for a given wavelength, it is seen that for the example case we have an effective slit radius of about 200 m.

(42) Note that the effective slit width is reduced considerably by using a shorter crystal. E.g. the sideways translation will be halved by using a half as long crystal. However, using a shorter crystal is detrimental to the amount of photons being upconverted. The expert in the field will know that for large beam diameters and collinear mixing the quantum efficiency will generally decrease by the crystal length squared. For the non-collinear case with angles larger than w.sub.0/L, the situation is quite different and a detailed analysis reveals that for the case of thermal light upconversion half as long crystal leads only to half as many upconverted photons, where the significant decrease in quantum efficiency is partly alleviated by a corresponding increase in angular acceptance. Surprisingly, this holds true for even very large incidence angles of infrared light to a narrow laser beam in a relatively long crystal.

(43) Now, to calculate the spectral resolution achievable in a setup like the one illustrated in FIG. 2, where we use an f=200 mm as the first lens after the crystal and the grating is tilted 60 relative to the z-axis. the radius of the illuminated area on the grating will be: 200 mm*sin(2.9)/cos(60)=20 mm. Hence, it will fit nicely within a 50 mm grating (although the grating needs only be 20 mm high).

(44) For optimal resolution, it is better to use a grating with more lines per mm (up to a point). For this calculation we consider a blazed grating with 1200 lines/mm using the first order diffraction. The blazed grating is particular efficient in making sure that most of the light is generated in the desired spectral order, and are available at low cost from many different vendors.

(45) The equation for diffraction on a grating can be expressed as:
a[sin .sub.m+sin .sub.i]=m

(46) Where a is the grating constant (1/a=grating lines/mm), m is the diffraction order, is the wavelength of the light after upconversion and .sub.m and .sub.i are the output and incidence angles to the grating respectively.

(47) The grating with 1200 lines/mm has a wavelength dispersion of around 0.74 nm/mrad. The light we send to the grating appears to come from an area of up to 400 m in diameter (the effective slit width), giving an uncertainty of the angle that is hitting the grating of 400 m/200 mm=2 mrad (i.e. +/1 mrad). This uncertainty in angle of incident light is a fundamental limiting factor of how the achievable resolution in the spectrometer.

(48) Let's assume we intend to image the whole upconverted spectrum in the above example from 3.2 m to 5.5 m, corresponding to 798 nm to 891 nm in upconverted wavelengths. These wavelengths will be distributed by the grating between .sub.m(798 nm)=0.0917 rad to .sub.m(891 nm)=0.2046 rad. It can be shown that the uncertainty in the angle in .sub.i scales with Cos(.sub.i)/Cos(.sub.m) to .sub.m, and hence the uncertainty on .sub.m in this case is about half of what it is in .sub.i. Hence .sub.m uncertainty is about +/0.5 mrad. The bandwidth (891 nm798 nm=93 nm) is distributed over 0.2046 rad0.0917 rad=0.1129 rad. With an uncertainty of em of +/0.0005 rad, this gives a resolution on the upconverted side of 93 nm*0.0005 rad/0.1129 rad=+/0.4 nm.+/0.4 nm uncertainty on the upconverted wavelength corresponds to a quite different uncertainty (or rather resolution) for the infrared light. From energy conservation we have:

(49) 1 up = 1 laser + 1 IR

(50) Which by differentiation leads us to the bandwidth in upconverted wavelength .sub.up corresponds to a bandwidth much larger at the infrared wavelength .sub.IR:

(51) IR = ( IR up ) 2 up

(52) Hence, if the laser is 1064 nm, and the IR wavelength is 3500 nm, then .sub.up is 816 nm. This scales the +/0.4 nm @ 816 nm to +/7.4 nm @ 3500 nm. Sometimes it is easier to consider bandwidth in units that are energy-equivalent such as MHz or cm.sup.1, where the bandwidth is conserved in the wavelength transformation process (fact of energy conservation).

(53) Now, to image the light diffracted from the grating to a detector we need an optical component such as a lens. The lens should be at a position to catch all the light diffracted in the desired order by the grating, and focus it to a line on the detector. In this case, the detector should be placed in the focal plane of the lens, since it should convert the spectral information contained in the angles of the light to positions on the detector. If we want to distribute the light over e.g. 2000 pixels of 14 m width, we can choose a lens with focal length of 2000*14 m/0.1129 rad=248 mm. So a lens with a focal length of 200 mm would be a reasonable choice to image the light to the line-detector. The uncertainty in output angles of +/0.5 mrad through a 200 mm lens, corresponds to position uncertainty on the camera of 0.5 mrad*200 mm=+/100 m. This is well larger than the pixel size of the example line-detector, so one might design the system with a shorter focal length without significant decrease in spectral resolution.

(54) One solution to achieve better spectral resolution is to avoid using the larger upconverted angles in the plane where the grating is active. One simple solution for this purpose is to use a 25 mm grating in the above described setup. The 25 mm grating is high enough to capture all the light in the vertical direction, but since it is tilted 60 degrees it is not wide enough to capture the light at large angles in the horizontal plane. In the example above a 25 mm grating would limit the horizontal angles to about 25 mm/2*Cos(60)/200 mm=31.25 mrad=1.8. This would limit the largest horizontal angle to be 1.8 degree rather than 2.9 degree (external angle). Again this would limit the width of the apparent horizontal aperture to be a factor 1.8/2.9 times smaller, i.e. 174 m*1.8/2.9=108 m, thus leading to significant improvement in the spectral resolution, at the cost of lower signal at the wavelengths phase matched above 1.8 degree in FIG. 9A. It is important to realize that the broadening of the virtual slit from the incident angles is relevant only for the angle projected to the x-axis. In particular, upconverted light with only an angle in the y-direction will not give rise to transverse broadening of the effective slit width. Hence, by making sure that preventing light upconverted at larger angles in the x-direction from reaching the dispersive element, we thus limit the virtual slit dimension and hence improve the spectral resolution obtained, at the cost of some degradation of the signal for the wavelengths cut off.

(55) The effect of the grating acting as an angle limiter, would remove some of the upconverted light. To be more precise, the horizontal edges of the phase match curve would be removed. as indicated in the contour plots shown in FIGS. 13A-C. Contour plot 13A and 13B shows the lower and upper half of the phasematch curve in FIG. 11A, as it would appear on a grating that is intentionally too narrow to accept the full circular pattern. FIG. 13C shows the combination of FIGS. 13A and 13B. In fact, a 2525 mm grating tilted 60 relative to the z-axis, has an apparent size of 12.525 mm when viewed along the z-axis.