This invention relates to a stationary waveguide broadband spectrum analyser, i.e. an analyser that in operation requires no relatively moving parts (although this does not rule out parts that can be adjusted prior to operation).
The invention is applicable especially, but not necessarily, to an optical waveguide spectrum analyser for the spectral analysis of light. In the present specification, light includes any wavelength which can be transmitted by internal reflection in a transparent medium, and therefore includes UV and IR wavelengths as well as visible wavelengths.
Spectrometers are widely used in a variety of disciplines for measurement of a diverse quantity of measurands. Interferometry-based spectrometer systems consist of either temporally scanned configurations or spatially scanned configurations that are typically stationary.
Many commercially available spectrometers are based on diffraction grating technology. These devices are not inherently capable of characterization of wide spectral bandwidths without scanning, which again introduces moving parts to the system.
A fibre device is known which act as a Young's slits device and is used to produce a diffraction pattern at the detector (see “Electronically scanned optical-fiber Young's white-light interferometer”, S. Chen et al., Optics Letters, Vol. 16, No. 10, 15 May 1991, Page 761). This device uses the interference fringe visibility profile to detect optical path length changes and is therefore not inherently capable of the spectral analysis of sources.
A second fibre device has been reported which acts as a wavemeter (see ‘James J. Snyder and Stephen L. Kwiatkowski, “Wavelength measurement with a Young's interferometer, Opt. Eng. 44, 083602 (Aug. 23, 2005); doi:10.1117/1.2030947”’). This device is based on an unbalanced interferometric configuration. It is, therefore, incapable of characterizing low coherence sources as it does not scan through a zero delay point and cannot produce an interferogram when a broadband source is used to illuminate the interferometer. This wavemeter embodiment also uses fringe counting to determine wavelength and is therefore not capable of interrogating more than a single wavelength in any one scan.
U.S. Pat. No. 6,016,197 (Krivoshlykov) describes a spectral analyser which uses a lens for performing a Fourier transformation. The use of a bulk optical component is undesirable.
Marcenac, D. D. et al: “Maximum-entropy optical spectrum analyser”, Optics Letters, 20, 1995, pp. 1074-1076 concerns the examination of laser lines and monochromatic light in a narrow spectral range. It is therefore not directed to a broadband solution.
To date, as far as we are aware, stationary interferometric Fourier transform devices for broadband applications are bulk optic and generally use a form of a Michelson interferometer.
According to the present invention there is provided a waveguide spectrum analyser comprising:
Preferably the data processing means transforms the captured image data to the frequency domain using at least one of HTP and DFTS processing.
The invention permits the analysis of broadband sources at least over the range 200 nm to 2500 nm, depending on the materials selected.
In a preferred embodiment the electromagnetic radiation is light and the input and output waveguides are optical waveguides, most preferably optical fibre waveguides but alternatively optical waveguide chips.
The input waveguide may be a single-mode waveguide. Alternatively it may be a multi-mode waveguide and the input radiation is distributed by the splitter between more than two single-mode waveguides.
The spectrum analyser further includes means to provide a reference image data to the processing means, the reference image being formed on the detector array by radiation of known wavelength.
The invention also provides a method for the spectral analysis of a beam of electromagnetic radiation, the method comprising:
Although the embodiments of the invention described herein are designed for the spectral analysis of light, the invention is applicable to the spectral analysis of electromagnetic radiation in general, it being understood that the nature of the waveguides and detector array will be selected according to the wavelengths concerned.
Embodiments of the invention will now be described, by way of example, with reference to the accompanying drawings, in which:
In use, a beam of light to be interrogated (i.e. spectrally analysed) is coupled into the input arm 10 and split, independently of wavelength, equally between the two output arms 14, 16 by the splitter 18. The input arm 12, when present, is used to receive a monochromatic reference beam, as will be described. This too is split equally between the two output arms 14, 16 by the splitter 18 and is used to calibrate the spectrum analyser. Each output arm 14, 16 has a respective exit port 20, 22 facing a CCD array 24 of individual photosensitive detector elements (pixels) 30,
The interference fringe frequency is wavelength dependent. Thus the total interferogram generated across the detector array can be Fourier processed to provide a spectral characterization of the light illuminating the interferometer.
This is done by dedicated processing electronics (or a suitably programmed PC) 26 which captures an image of the fringes from the array by sampling the array at above the Nyquist rate, and transforms the captured image data to the frequency domain by Fourier analysis.
Preferably, but not necessarily, the two single-mode optical fibre arms 14, 16 are substantially identical in material, waveguide dimension, index contrast and length so that the light propagating along each arm experiences the same optical path delay and dispersion in each arm. In the embodiment shown in
As the light exits the port 20 or 22 it will diffract into the surrounding medium with a divergence angle governed by the wavelength dependent numerical aperture (NA) of the fibre.
NA=√{square root over (nco2−ncl2)} (1)
where nco is the wavelength dependent refractive index of the single-mode optical waveguide core and ncl is the wavelength dependent refractive index of the single-mode optical waveguide cladding.
In the present embodiment the array 24 is linear and the exit ports 20, 22 are spaced apart in the longitudinal direction of the array. At the ports 20, 22 the optical axes of the arms 14, 16 are mutually parallel and normal to the array 24. However, in other embodiments the relative delay between the beams from the two exit ports can be varied by changing the angles at which the light is launched from the exit ports 20, 22 (which need not be parallel) and/or by changing the separation of the two exit ports. The relative delay should be such that the zero delay position is scanned to enable spectral analysis of low-coherence sources.
The divergent output beams from the two exit ports 20, 22 may be collimated so that the beams do not diverge beyond the limits of the array 24.
Referring to
A(Q, t)=K1A(P1,t−t1)+K2A(P2,t−t2) (2)
where P1 and P2 are the centres of secondary disturbances and t1 and t2 are the times taken for the light to travel from P1 to Q and from P2 to Q respectively. The intensity I(Q) at Q is then given by the general interference law for partially coherent light as
I(Q)=I1(Q)+I2(Q)+2√{square root over (I1(Q))}√{square root over (I2(Q))}[{tilde over (γ)}12(τ) (3)
where I1(Q) and I2(Q) are the irradiances from the respective optical fibre arms, denotes the real part of a complex number, {tilde over (γ)}12(τ) is the complex mutual degree of coherence function and τ represents delay.
In the case where I1≅I2 and illumination is by a narrow line-width source of mean optical frequency
I(τ)=2I0[1+V(τ)cos φ(τ)] (4).
where I0=I1=I2. Here, V(τ) is a slowly varying irradiance envelope function that is the interferogram visibility, and the temporal evolution of phase is given by φ(τ)=
Considering
I(x)=2I0[1+V(x)cos φ(x)] (5).
In this equation, the visibility V(x) varies slowly with x due to the quasi-monochromatic nature of the source, and the spatial evolution of phase φ(x) is given by
φ(x)={k(x)φ(τ)+φC} (6).
In a manner analogous to the stationary Michelson interferometer, the term k(x) is constant in the case of negligible separation d between the fibre output arms and in the case of perfectly plane surfaces on the fibre ends. In practice, however, ideal surfaces cannot be guaranteed and a non-negligible separation d is deliberately applied in order to achieve a scan of the optical path delay; the x dependence of k(x) is therefore explicitly included. A phase correction term φC is introduced to allow an arbitrary choice of delay origin. An expression for the intensity distribution across the pixel array I(x) in terms of phase evolution φ(τ) is then
I(x)=2I0[1+V(x)cos{k(x)φ(τ)+φC}] (7).
When a reference laser beam, of optical frequency
I
R(x)=2I0R[1+VR(x)cos{kR(x)φR(τ)+φCR}] (8)
where subscript ‘R’, is used to denote terms relating to the reference laser.
As the numerical aperture (NA) of the optical fibre is wavelength dependent, the phase evolution φ(x) of the interference pattern for each source wavelength will be unique and will give rise to a distinct interference pattern.
For a fixed source frequency or wavelength, fixed NA fibres and coplanar fibre ports, the evolution of interferogram phase φ(τ) across the detector array plane is governed by the separation between the ends of the fibre ports, d and the distance, z between the fibre ends and the detector array plane, as shown in
For a cylindrical optical fibre, the condition for single-mode operation is
where V is the waveguide parameter, λ is the vacuum wavelength of light, a is the core radius, nco is the core refractive index and ncl is cladding refractive index.
Silica has a transmission window that spans the wavelength range from ˜200 nm to ˜2000 nm. As such, silica is the preferred material for optical fibres in this wavelength region. For extended wavelengths, there are a range of soft glass materials such as germanate and chalcoginide that provide transmission windows out into the mid IR.
In one embodiment a visible/near IR stationary optical fibre waveguide spectrum analyser based on
In another embodiment a telecoms wavelength version of the stationary optical fibre waveguide spectrum analyser uses smf-28 optical fibre that operates single-mode for wavelengths above 1260 nm and an InGaAs detector array with spectral bandwidth (900-1700) nm.
It is possible to get enhanced InGaAs detector arrays that can operate up to 2200 nm and 2500 nm.
There are a number of methods of referencing the interferometer. When the input arm 12 is present the device allows a reference beam of known wavelength to be launched to co-illuminate the detector array 24 with the light to be interrogated, as described above. In such a case the wavelength of the reference beam should be separable from the wavelength range of the light to be interrogated in the frequency domain. Alternatively, the device can be precalibrated by introducing the reference beam into the arm 10 or, if present, the arm 12, with reference images produced thereby stored in a database for subsequent use by the data processing electronics 26.
As shown in
The path difference, ΔPD, between the beams exiting from port 1 and port 2 is determined by the separation of the two ports, d, the distance to the detector array, z, and the spatial extent of the detector array, x.
In the observation plane the fringe intensity varies sinusoidally with x as
where λ is the wavelength of the light.
As mentioned above, the sinusoidal variation I(x) is sampled above the Nyquist rate at the highest frequency of light to be interrogated to prevent aliasing.
The material of the detector array 26 depends on the wavelengths (or frequencies) of the light to be interrogated. For example, a silicon CCD or CMOS detector array can be used to interrogate light at visible and NIR wavelengths. An InGaAs detector array can be used to interrogate light from ˜1000 nm through to ˜2000 nm. (The determination of the detector array type depends on the responsivity to the wavelength of the light to be interrogated and the optical power. It has been shown for example, that where high optical powers are used InGaAs detector arrays can also be used for detection of visible wavelengths.)
In one example, the fibre has a numerical aperture of 0.13@1550 nm, the fibres exit port ends are separated by ˜5 mm, the distance to the detector array is ˜200 mm, the detector array has 512 pixels with a pitch p (
The foregoing assumes that the input arm 10 is a single mode fibre. Multi-mode light that is incident on single-mode fibre is spatially filtered so that only a single mode of the incident light is propagated in the single-mode fibre. This leads to a loss of optical input power to ensure high optical resolution of the signal.
Recently, the photonic lantern has been used to capture multi-mode light in a multi-mode fibre and then distribute it with minimum loss into a plurality of single-mode fibres (see “Multimode fiber devices with single mode performance”, S. G. Leon-Saval et al., Optics Letters, Vol. 30, No. 19, 1 Oct. 2005, Page 2545). Using this approach, it is possible to extend the stationary waveguide interferometer to interrogate multi-mode light. Possible options to resolve the spectral content of this multi-mode light are shown in
In
In
The interferograms are captured and Fourier transform processed to retrieve information, such as the spectral content, of the light illuminating the detector array.
DFTS
Having detected an interferogram or interferogram ensemble at the detector array 26, it is captured using a frame grabber and stored for processing on the processing electronics or computer 26. Dispersive Fourier Transform Spectrometry (DFTS) can be used to recover optical frequency and phase information. From Equation (3), the oscillatory part of the interferogram, Ios is
I
os=2√{square root over (I1(Q))}√{square root over (I2(Q))}[{tilde over (γ)}12(τ)] (12)
As a consequence of the optical equivalent of the Wiener-Khintchine theorem, the complex degree of coherence function {tilde over (γ)}12(τ) and the normalised complex interferogram spectrum {tilde over (G)}(ω)exp[−i φ12(ω)] are Fourier transform pairs. This may be expressed mathematically as
{tilde over (γ)}12(τ)=∫−∞∞{{tilde over (G)}(ω)exp[−i φ12(ω)]}exp[−i ωτ]dω (13)
Defining the complex spectrum {tilde over (S)}(ω) of the interferogram as
{tilde over (S)}(ω)={tilde over (G)}(ω)exp[−i φ12(ω)] (14)
it follows that
{tilde over (γ)}12(τ)=−1{{tilde over (S)}(ω)} (15)
Then, using the linearity of the Fourier transform { }
where −1{ } denotes the inverse Fourier transform, * denotes the complex conjugate and {tilde over (S)}(ω) is Hermitian. For negative optical frequencies {tilde over (G)}(ω) and therefore {tilde over (S)}(ω) are zero-valued so that from Equation (16)
{tilde over (S)}(ω)={2[{tilde over (γ)}12(τ)]}, for all ω>0 (17)
Combining Equation (12) with Equation (17):
{tilde over (G)}(ω)exp[−i φ12(ω)]∝{Ios(τ)};
and {tilde over (G)}(ω)∝|{Ios(τ)}| (18)
The spectral phase is then found from the argument of the Fourier transform of the interferogram
φ12(ω)=arg {Ios(τ)} (19).
Equation (19) forms the basis of DFTS processing.
HTP
The Hilbert transform processing (HTP) approach is used for measurements of mean wavelengths or frequencies. The technique is based on the comparison between recovered phase values from the respective analytic signals of the signal and reference sources that illuminate the interferometer. In the case of an interferogram generated using quasi-monochromatic light, the spectral amplitudes will only have appreciable values in a frequency range Δω that is small compared to the mean frequency
A
os(τ)=x(τ)+ih(τ) (20)
where x(τ) is the real part of the signal and ih(τ) is its associated imaginary part. The real and imaginary parts are Hilbert transforms of each other, as expressed by
x(τ)=∫0∞a(ω)cos[φ(ω)−
ih(τ)=i∫0∞a(ω)sin[φ(ω)−
Recovery of the analytic signal from a captured interferogram is therefore achieved as follows. A Fourier Transform {tilde over (X)}(ω) of the real signal x(τ) contains both positive and negative frequency components. A Fourier transform of the analytic signal is one-sided with no negative frequency components, and equal to twice the single-sided Fourier transform of the real signal
The analytic signal is found then from twice the inverse transform of the positive frequency components of the Fourier transform of the real signal
A
os(τ)=2−1[{tilde over (X)}POS(ω)] (23).
Taking the argument of this complex temporal representation of the interferogram yields the desired values of temporal phase φ(τ) for the interferogram in question
The invention is not limited to the embodiments described herein which may be modified or varied without departing from the scope of the invention.
Number | Date | Country | Kind |
---|---|---|---|
1107227.9 | May 2011 | GB | national |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/EP2012/057826 | 4/27/2012 | WO | 00 | 4/10/2014 |