This invention relates to a system and method to measure chromatic dispersion in short lengths of waveguides using a three wave interference pattern and a single arm interferometer.
One of the main components in a photonic device is a waveguide or an optical fiber which serves as a low-loss medium for light transmission. An important characteristic of waveguides such as optical fiber is the dispersion that light experiences as it travels inside the waveguide. Dispersion is the phenomenon that causes different frequencies of light to travel at different velocities. The phenomenon of dispersion is commonly observed through the spreading of light by a prism. When white light, which contains a broad spectrum of frequencies, enters a prism the different wavelengths are bent at different angles since each frequency sees a different index of refraction, a phenomenon first quantified by Newton in the 17th century. Inside a fiber, this variation in the index of refraction with frequency is what causes the frequency dependence of the velocity.
A more modern example of the phenomenon of dispersion is the affect it has on the performance of photonic devices used in communication systems. In these systems, dispersion, or more specifically second order dispersion, leads to a broadening of the pulses used to represent 1 or 0 in a digital communication system. Pulse broadening causes adjacent bits to overlap and leads to intersymbol interference. Intersymbol interference occurs when a pulse is broadened beyond its allocated bit slot to such an extent that it begins to overlap with adjacent bits and it is no longer possible to determine whether or not a specific bit contains a 1 or a 0.
As a result of intersymbol interference the allocated bit slots must be widened and this effectively lowers the number of bits that can be transmitted in a given period of time and reduces the system bandwidth. As a result modern communication systems have evolved methods to mitigate the effects of dispersion.
Current methods of countering the effects of dispersion in an optical fiber use dispersion compensating devices such as chirped fiber Bragg gratings and dispersion compensating fiber (DCF). In order to effectively use these techniques it is important to know the exact magnitude of the dispersion that is being compensated for. As a result, knowledge of the dispersion in both the transmission system and the dispersion compensation system is important to the design of the overall communication system.
Knowledge of dispersion in a waveguide is also significant for the study of fiber based nonlinear wave interaction phenomena. An optical soliton is a pulse that maintains a constant shape (width) as it propagates along a fiber (first order soliton) or has a shape that is periodic with propagation (higher order soliton). This is due to the fact that the effects of dispersion and self phase modulation (SPM) are in balance. SPM is the effect whereby the phase of a given pulse is modified by its own intensity profile. Knowledge of the dispersion in an optical fiber allows for the determination of the required intensity for the formation of an optical soliton. This effect has also been used in the area of soliton effect pulse compression where the combination of the chirping effect of SPM and subsequent distributed compression effect of negative dispersion is used to compress an optical pulse. Knowledge of dispersion is also important for the study of nonlinear effects such as second harmonic generation, three-wave mixing and four-wave mixing since it determines the interaction lengths between the various wavelengths. Dispersion is particularly important in techniques that aim to extend this interaction length such as in Quasi Phase Matching (QPM) devices.
Theory on Chromatic Dispersion of a Waveguide
Dispersion is the phenomenon whereby the index of refraction of a material varies with the frequency or wavelength of the radiation being transmitted through it. The term ‘Chromatic Dispersion’ is often used to emphasize this wavelength dependence. The total dispersion in a waveguide or an optical fiber is a function of both the material composition (material dispersion) and the geometry of the waveguide (waveguide dispersion). This section outlines the contributions of both material and waveguide dispersion, identifies their physical source and develops the mathematical terminology for their description.
Dispersion in a Waveguide
When light is confined in a waveguide or an optical fiber the index is a property of both the material and the geometry of the waveguide. The waveguide geometry changes the refractive index via optical confinement by the waveguide structure. The refractive index is therefore a function of both the material and waveguide contributions. For this reason in a fiber or a waveguide the index is known as an effective index.
The dispersion parameter, D, which represents second order dispersion since it describes how the second derivative of the effective index varies with respect to wavelength:
The dispersion parameter is important since it is related to pulse broadening which greatly limits the bit rate of a communication system.
The dispersion parameter of a waveguide such as an optical fiber is given by the total dispersion due to both the material and waveguide contributions. The total dispersion is the combination of the material dispersion and the waveguide dispersion and thus the dispersion parameter of a waveguide is given by:
The next two sections discuss the contributions that both material and waveguide dispersion make individually to the total dispersion.
Material Dispersion
Material dispersion originates from the frequency or wavelength dependent response of the atoms/molecules of a material to electromagnetic waves. All media are dispersive and the only non-dispersive medium is a vacuum. The source of material dispersion can be examined from an understanding of the atomic nature of matter and the frequency dependent aspect of that nature. Material dispersion occurs because atoms absorb and re-radiate electromagnetic radiation more efficiently as the frequency approaches a certain characteristic frequency for that particular atom called the resonance frequency.
When an applied electric field impinges on an atom it distorts the charge cloud surrounding that atom and induces a polarization that is inversely proportional to the relative difference between the frequency of the field and the resonance frequency of the atom. Thus the closer the frequency of the electromagnetic radiation is to the atoms resonance frequency the larger the induced polarization and the larger the displacement between the negative charge cloud and the positive nucleus.
The material dispersion is then determined by taking the derivative of the group index of the material with respect to wavelength or equivalently the second derivative of the absolute index with respect to wavelength:
Waveguide Dispersion
Waveguide dispersion occurs because waveguide geometry variably affects the velocity of different frequencies of light. More technically, waveguide dispersion is caused by the variation in the index of refraction due to the confinement of light in an optical mode. Waveguide dispersion is a function of the material parameters of the waveguide such as the normalized core-cladding index difference, Δ=(ncore−ncladding)/ncore and geometrical parameters such as the core size, a. The index in a waveguide is known as an effective index, neff, because of the portion of the index change caused by propagation in a confined medium.
Confinement is best described by a quantity known as the V parameter, which is a function of both the material and geometry of the waveguide. The V parameter is given by Eq. 4:
Propagation in a waveguide is described by a quantity known as the normalized propagation constant, b, which is also a function of the material and geometry of the waveguide. This quantity is given in Eq. 5:
The contribution of the waveguide to the dispersion parameter depends on the confinement and propagation of the light in a waveguide and hence it is a function of both the V parameter and the normalized propagation constant, b. The waveguide dispersion can be calculated via knowledge of V and b via Eq. 6:
In most cases the main effect of the waveguide dispersion in standard single mode fibers is a reduction in dispersion compared to dispersion in bulk. In comparison to material dispersion the contribution of waveguide dispersion is quite small and in most standard single mode fibers it only shifts the zero dispersion wavelength from 1276 nm to 1310 nm.
In summary, the dispersion in a waveguide or an optical fiber is caused not only by the material but also by the effect of confinement and propagation in the waveguide. Thus accurate knowledge of the dispersion in a waveguide cannot be made by simple knowledge of the material dispersion but must include the effect of the waveguide. As a result either the dimensions of the waveguide must be known to a high degree of accuracy so that the waveguide dispersion can be calculated (which is not easy since fabrication processes are hardly perfect) or the dispersion must be measured empirically. Accurate measurement of the (total) dispersion parameter, D, is important to the design of photonic systems.
Polarization Mode Dispersion
In addition to the above, optical waveguides may suffer from polarization mode dispersion (PMD). PMD may exist in fibers with asymmetrical cores. In optical fibers, the light that travels along one of the two polarization axis travels at a right angle to light traveling along the other axis. In asymmetrical optical fibers, the light travels along the two axes at different speeds. This causes pulses to spread, which can cause them to become undetectable at the detector.
Conventional Measurement Techniques
There are 3 categories of dispersion measurement techniques: Time of flight (TOF), Modulation phase shift (MPS) and Interferometric. TOF and MPS are the most widely used commercial dispersion measurement techniques. Interferometric techniques are not widely used commercially but have been used in laboratories for dispersion measurements. Interferometric techniques come in two forms; temporal and spectral. The existing techniques differ in measurement precision and fiber length requirements.
Time of Flight Technique
In the TOF technique the second order dispersion parameter, D, hereafter referred to simply as the dispersion parameter, can be determined either by measuring the relative temporal delay between pulses at different wavelengths or by measuring the pulse broadening itself. The relative temporal delay between pulses at different wavelengths is measured to determine the group velocity which can then be used to determine the dispersion parameter using Eq. 7:
The above equation can also be used to determine the dispersion parameter from the pulse broadening itself if Δt is the measured pulse broadening and Δλ is the bandwidth of the wavelengths in the pulse. The measurement precision achievable by the TOF technique is on the order of 1 ps/nm.
One of the main problems with the TOF technique is that it generally requires several kilometers of fiber to accumulate an appreciable difference in time for different wavelengths. Another issue with the TOF technique when the pulse broadening is measured directly is that the pulse width is affected by changes in the pulse shape which leads to errors in the measurement of the dispersion parameter. As a result, in order to measure the dispersion parameter with a precision near 1 ps/nm-km several kilometers of fiber are required.
Modulation Phase Shift Technique
The MPS technique is another dispersion characterization technique that requires long lengths of fiber. In the MPS technique, a continuous-wave optical signal is amplitude modulated by an RF signal, and the dispersion parameter is determined by measuring the RF phase delay experienced by the optical carriers at the different wavelengths.
The RF phase delay information is extracted by this technique, and by taking the second derivative of the phase information, the dispersion parameter can be determined. Measurement precision achievable by the MPS technique is on the order of 0.07 ps/nm. Due to its higher precision, MPS has become the industry standard for measuring dispersion in optical fibers. However, MPS has several disadvantages. The first is that it is expensive to implement since the components required such as an RF analyzer and a tunable laser, are costly. The second is that its precision is limited by both the stability and jitter of the RF signal.
MPS has several limitations on the minimum device length that it is capable of characterizing. In the MPS method the width of the modulated signal limits the minimum characterizable device length. This method also typically requires fiber lengths in excess of tens of meters to obtain a precision to better than 1 ps/nm-km. Therefore it is not desirable for the characterization of specialty fibers or gain fibers, of which long fiber lengths are expensive to acquire or not available. Also, when fiber uniformity changes significantly along its length, the dispersion of a long span of fiber cannot be used to accurately represent that of a short section of fiber. In such cases, dispersion measurement performed directly on short fiber samples is desirable. As a result a technique for measuring the dispersion of short lengths of fiber is desired.
Dispersion Measurements on Short Length
Interferometric techniques are capable of characterizing the dispersion on fiber lengths below 1 m. There are two categories of interferometric techniques for making dispersion measurements on fiber of short length: temporal and spectral.
Temporal Interferometry (Fourier Transform Spectroscopy)
Dual Arm temporal interferometry employs a broadband source and a variable optical path to produce a temporal interferogram between a fixed path through the test fiber and variable air path. It involves moving one arm of the interferometer at a constant speed and plotting the interference pattern as a function of delay length (time). The spectral amplitude and phase are then determined from the Fourier transform of the temporal interferogram.
A temporal interferogram gives the phase variation as a function of time. The spectral phase variation can be extracted from the temporal interferogram if a Fourier Transform is applied to it. The spectral phase contains the dispersion information which can be indirectly obtained by taking the second derivative of the spectral phase. A precision of 0.0015 ps/nm measured on a 0.814-m-long photonic crystal fiber was recently reported using temporal interferometry. The main disadvantage of temporal interferometry is that it is susceptible to noise resulting from both translation inaccuracy and vibration of the optics in the variable path. A tracking laser is typically required to calibrate the delay path length. Another problem with this technique is that a second derivative of the phase information must be taken to obtain the dispersion parameter which means that it is less accurate than a method that can obtain the dispersion parameter directly.
Spectral Interferometry
Spectral interferometry, like temporal interferometry, is capable of characterizing the dispersion in short length fiber (<1 m). In spectral interferometry, instead of stepping the length of one of the arms, a scan of the wavelength domain performed to produce a spectral interferogram. Spectral interferometry is generally more stable than temporal interferometry since the arms of the interferometer are kept stationary. Thus it is simpler than temporal interferometry since no tracking laser is necessary.
There are two types of spectral interferometry, one is general and does not require balancing, and another, the special case, is ‘balanced’. In the balanced case it is possible to directly measure the dispersion parameter from the interferogram. This makes it more accurate than temporal interferometry and it is for this reason that spectral interferometry is discussed as a dispersion measurement technique.
In general spectral interferometry the dispersion parameter is obtained from the interference spectrum produced by two time delayed light pulses/beams in an unbalanced dual arm interferometer. Two pulses/beams from the two arms of the interferometer are set up to interfere in a spectrometer and a spectral interferogram is produced. The interference pattern produced for a given time or phase delay is given by:
The last two terms in Eq. 8 result in spectral interference pattern via a cos(Δφ(ω)+ωτ) term.
There are several ways to extract the phase information from the cosine term but the most prevalent way to do so is to take the Inverse Fourier transform of the spectral interference pattern. Note that f(ω)=F.T.f(t)=|E*o(ω)E(ω)|exp(iΔφ(ω)) contains all the phase information on the spectral phase difference Δφ(ω). Therefore, if f(ω) can be extracted from the interference pattern then the phase difference information can be known. If an Inverse Fourier Transform of the spectral interference is performed on the interference pattern the following is obtained:
F.T.−1(I(ω))=E*o(−t){circumflex over (x)}Eo(t)+E*(−t){circumflex over (x)}E(t)+f(t−τ)+f(−t−τ) Eq. 9
If all terms except the f(t−τ) term get filtered out via a band pass filter then the phase information can be extracted from a Fourier Transform on f(t−τ).
The phase information can then be extracted by applying a Fourier Transform to the filtered component f(t−τ) thereby transferring it back to the spectral domain. The complex amplitude therefore becomes f(ω)=|Eo(ω)∥E(ω)|exp(iΔφ(ω)+ωr). The phase of this complex amplitude minus the linear part (ωτ) that is due to the delay, yields the spectral phase difference between the two pulses as a function of ω and is independent of the delay between the two pulses. In this way the phase difference between the two pulses can be obtained.
If one of the pulses travels through a non-dispersive medium such as air and the other pulse travels through a dispersive medium such as an optical fiber then the phase difference spectrum will be directly related to the dispersion in the fiber. Thus the dispersion parameter plot can be determined by taking the second derivative of the phase difference spectrum with respect to wavelength.
The main issue with this form of spectral interferometry, however, is that the dispersion parameter is not determined directly but rather via a second order derivative of the phase information with respect to wavelength. Therefore, like temporal interferometry, this general unbalanced method of spectral interferometry is not as accurate as the balanced method capable of measuring the dispersion parameter directly of the present invention.
In balanced spectral interferometry the arm lengths of an interferometer are kept constant and they are balanced for a given wavelength called the central wavelength such that the group delay in both arms is the same. This allows for the removal of the effect of the large linear dispersion term in the interferogram. Balanced interferometry measures the dispersion parameter D at the wavelength at which the group delay is the same in both arms. This wavelength is henceforth referred to as the central wavelength. The dispersion parameter in this case can be directly determined from the phase information in the spectral interferogram without differentiation of the phase. For this reason it is more accurate than both unbalanced general spectral interferometry and temporal interferometry. As a result balanced spectral interferometry is often used to obtain accurate dispersion measurements in short length waveguides and fibers.
Both forms of spectral interferometry are considered to be less susceptible to noise since the arms of the interferometer are kept still and there are no moving parts. It is for this reason that spectral interferometry in general is considered to be more accurate than temporal interferometry. Spectral interferometry is therefore considered to be the technique of choice for measuring the dispersion of photonic components and spectral depth resolved optical imaging. One well known and important class of spectral interferometry is optical coherence tomography (OCT).
Balanced dual arm spectral interferometers are typically found in a Michelson or a Mach-Zehnder configuration in which the path lengths are equalized at the given wavelength in which the dispersion is to be measured. The most often used configuration, however, is the Michelson and the discussion that follows considers the Michelson interferometer. In a balanced Michelson interferometer the dispersion is measured from the interference between two waves: one that passes through the test fiber and another that passes through an air path. Balancing the air path length with the fiber eliminates the effect of the group index of the fiber in the interference pattern. This allows for the measurement of the second derivative of the effective index with respect to wavelength directly from the interference pattern.
The main disadvantage of this Michelson configuration is that two types of path balancing must occur simultaneously. The first type of path balancing is coupler arm balancing wherein the coupler arms need to be balanced exactly or an extra set of interference fringes will be created from the reflections at the two end facets of the coupler arms as shown in
The second type of balancing is test fiber-air path balancing to ensure that the optical path length in the air path exactly equals that of the fiber for a given central wavelength. This ensures that the central wavelength in the interference pattern is within the viewable bandwidth of the OSA.
The main problem in implementing a Michelson interferometer is that the arms of the coupler cannot be balanced exactly and as a result the effect of the extra set of reflections produced at the coupler facets cannot be removed.
Comparison of Dispersion Measurement Techniques
There have been several techniques developed for the measurement of chromatic dispersion in fiber. Especially important are those developed for the measurement of short lengths of fiber. One reason that short length characterization techniques are important stems from recent developments in the design and fabrication of specialty fiber.
Specialty fiber such as Twin Hole Fiber (THF) (
Based on the above discussion, the technique of choice for dispersion measurement is balanced spectral interferometry since it will provide the most accurate measurements. As a result the new technique will employ balanced spectral interferometry.
The two important parameters in comparing dispersion measurement techniques are the minimum device length that each is capable of characterizing and the precision to which the characterization is achieved. It is generally desirable to characterize as short a fiber as possible with as high a precision as possible. It is also desirable to perform the measurement in the simplest way possible.
Therefore, what is needed is a new method for the measurement of dispersion that does not require the cancellation of any extra fringes. What is also needed is a method to measure the dispersion parameter in short lengths of optical fiber. The initial need for a short length characterization scheme came from the need to measure the dispersion of a specialty fiber such as THF, PCF, or gain fibre. This requirement is based on the expense of fibre, nonlinear wave interaction phenomena in fibre, and non-uniform dispersion along the length of a fibre.
According to the present invention, a system and method to determine chromatic dispersion in short lengths of waveguides using a three wave interference pattern and a single-arm interferometer has been developed.
In a first aspect, the present invention relates to an interferometer system for obtaining a measure of the chromatic dispersion of a waveguide comprising a radiation source operable to emit radiation connected to a means for separating incident and reflected waves; the means for separating incident and reflected waves having an output arm adjacent to a first end of the waveguide; the means for separating incident and reflected waves further connected to a detector; a collimating means positioned at a second end of the waveguide; and a reflecting means positioned at a balanced distance from the collimating means operable to reflect a test emission from the radiation source back through the collimating means, the waveguide, and the means for separating incident and reflected waves thereby generating an interference pattern that is recorded by the detector;
In a second aspect of the present invention said interference pattern consists of three waves wherein a first wave is a reflection of the test emission from one facet of the waveguide, a second wave is a reflection of the test emission from a second facet of the waveguide and a third wave is a reflection of the test emission from the reflecting means.
In another aspect, the present invention relates to an interferometric method for obtaining a measure of the chromatic dispersion of a waveguide comprising the steps of:
In a still further aspect of the present invention relates to a method for increasing a maximum length of a waveguide for which chromatic dispersion may be measured using an interferometer comprising the steps of:
The invention will be better understood with reference to the drawings.
The accompanying drawings, which are incorporated in and constitute a part of the specification, illustrate various example systems and methods and are not intended to limit the scope of the specification.
System Operation
1: Development of Single Arm Interferometry
A Single Arm Interferometer (SAI) can be produced by folding the two arms of a Michelson interferometer (
1 1 A New Concept
A balanced Single-Arm Interferometer (SAI) can directly measure dispersion in short fibers. A balanced SAI is depicted in
The SAI system may operate in several different embodiments, the source (11)/detector (23) pairing may be a broadband source operable with an Optical Spectrum Analyzer (OSA) or the source (11)/detector (23) pairing may be a tunable laser (31) operable with a detector (23) system.
The test fiber (17) is a waveguide. As previously mentioned waveguides encompass a large class of materials and include but are not limited to the following: optical fiber, photonic crystal fiber, nanowire, nanofiber, twin-hole fiber and etched waveguides.
This configuration is not only much simpler than a dual arm interferometer but as further evidenced by the experimental setup detailed in
The SAI is a balanced interferometer since the group delay in the test fiber (17) is the same as the group delay in the air path (27). It will be shown mathematically that this balancing of the group delay in each path (17, 27) allows the dispersion parameter to be measured directly from the interference pattern. The conceptual difference between SAI and Dual Arm interferometers is that, in SAI, the interference pattern is produced by three waves: two from the reflections (Uo, U1) at the facets of the test fiber (17) and one (U2) from a mirror (21) placed behind it (as shown by Uo, U1, and U2 in
From the phase information in this envelope the dispersion parameter can be extracted. Both dual and single arm balanced interferometers have in common this ability to directly measure the dispersion parameter from the interference pattern.
Optionally, the launch fiber (15) may be configured as a polarization controller. This may be accomplished by a number of means known to those skilled in the art, such as looping the launch fiber (15) into three coils. Alternately, the polarization controller may be implemented by placing a linear polarizer in the air path (27). By providing a polarization controller, the SAI may be operable to measure PMD.
The SAI configuration appears similar to common path interferometers, often used for depth imaging as in Common-Path Optical Coherence Tomography (CP-OCT). The SAI, however, is fundamentally different from CP-OCT since it utilizes 3 reflections (Uo, U1, and U2), and extracts the dispersion parameter directly from the envelope of the interference pattern. The differences between the Michelson Interferometer, CP-OCT and balanced Single Arm Interferometry are outlined in Table 1.
1.2 Mathematical Description
1.2.1.1 Equal Amplitude Case
Dispersion measurements can be made using a single-arm interferometer by extracting the second derivative of the effective index with respect to wavelength from the envelope of the interference pattern generated by three waves Uo, U1 and U2 depicted in
The extra reflection from the launch fiber (15) is minimized using angle polished connector (25) as shown in
U1=U0e−j2βL
U2=U0e−j2βL
In Eq. 10, Lf and Lair are the lengths of the test fiber (17) and the air path (27), respectively. β and ko are the propagation constant of the fundamental mode in the fiber and the propagation constant in free space. The interference pattern is produced by the interference of the three reflections (Uo, U1, and U2) is given by Eq. 11:
Eq. 11 contains two fast terms, with a phase φ1=+(βLf+koLair) and φ2=2(βLf+koLair). Since φ1 is slower than φ2 it will amplitude modulate the faster term. As a result the period of the ‘carrier’ will be that of the slowest of the fast terms, φcarrier=φ1. This carrier is then itself amplitude modulated by the slower term φenvelope=(βLf−koLair) to produce the ‘envelope’ of the interference pattern. This envelope is equivalent to the interference pattern produced by Michelson interferometer and it can be written as:
Uo2(5+4|cos(φenvelope)|) Eq. 12
The calculated interference pattern generated by the setup for a 39.5 cm SMF28™ test fiber is illustrated in
Applying a Taylor expansion to the phase of the slow envelope and replacing β with 2πneff/λ, where neff is the effective index of the fiber, gives the phase relation in Eq. 13:
The first term in Eq. 13 (in the square brackets) disappears when Lair is adjusted to balance out the group delay of the test fiber (17) at λo, the balanced wavelength. Taking the difference between the phases at two separate wavelengths; λ1 and λ2 results in:
Note that m is the number of fringes between the two wavelengths. If this phase difference is taken using a different pair of peaks/troughs (i.e. λ3 & λ4) the result is a system of equations in which d2neff/dλ2|λ
Note that, if the third-order dispersion is ignored, then only two wavelengths (e.g. λ1 and λ2) are required to calculate the second-order dispersion. This, however, would be less accurate. The dispersion parameter D can then be found as follows:
1.2.1.2 Unequal Amplitude Cases
In reality the reflections (Uo, U1, and U2) from the three facets of the interferometer as illustrated in
In general the reflections from each facet (Uo, U1) of the test fiber (17) and the reflection from the mirror (U2) shown in
U1=aU0e−j2βL
U2=bU0e−j2βL
In Eq. 16 Lf and Lair are the lengths of the test fiber (17) and the air path (27), respectively. β and ko are the propagation constant of the fundamental mode in the fiber and the propagation constant in free space. ‘a’ is the fraction of the amplitude reflected from the second facet in terms of the first and ‘b’ is the fraction of the amplitude reflected from the mirror (21) in terms of the fraction reflected from the first facet. The interference pattern of the spectral interferogram can be expressed as:
The expression in Eq. 17 can be treated as a fast-varying “carrier” (with respect to frequency or wavelength) modified by an upper and a lower slow-varying envelope, as shown in
Uo2(1+a2+b2+2a(b−1)+2b+4a|cos(φenvelope)|) Eq. 18
It will now be shown that although the magnitude of the interference pattern is not the same as the envelope for cases in which b≠1, the peak and trough locations of the two match exactly. As a result the phase information of the interferogram is preserved and the dispersion information can be extracted from the interferogram. Note that a=b=1 is a special case of this more general analysis and was presented in the previous section.
Since the phase of the upper envelope, φenvelope (and therefore the dispersion information) is unaffected by the magnitude of the reflections from the facets (Uo, U1) and the mirror (U2), the method for determining the dispersion parameter as presented in Eqs. 13 to 15 is valid even in the general case. The dispersion parameter, therefore, can always be obtained from an SAI.
As mentioned earlier, the main difference between the fringes produced in the SAI and those produced by dual arm interferometers is the presence of a fast carrier (Eq. 17) slowly modulated by the desired envelope. The presence of this carrier sets extra operational constraints that will be discussed in the next section.
1.3 System Parameters
There may be four factors of interest with regard to the dispersion measurement of the present invention. These factors may determine the quality and range of the output of the dispersion measurements. The first factor of interest is the wavelength resolution of the measurement, the second is the minimum required bandwidth of the source, the third is the measurable bandwidth of the dispersion curve, and the fourth is the test fiber length. The sections that follow discuss how each of these factors affect the output of the dispersion measurement.
1.3.1 Wavelength Resolution of the Dispersion Measurement
The wavelength resolution of the points in the plot of the dispersion parameter is determined by the minimum step size of the translation stage. With smaller step increments in the translation stage there are smaller step increments in the plot of the dispersion parameter vs. wavelength. This is because variation of the air path (27) changes the wavelength where the air path (27) and test fiber (17) are balanced and produces a new interferogram from which the dispersion parameter can be determined. Examination of Eq. 13 shows that the first term can be removed if the group delay in the air path (27) is equal to that in the test fiber path for the central wavelength, λo (central wavelength at which the group delay in test fiber (17) and air paths (27) are balanced). The relationship between the air path length (Lair) and the fiber length (Lf) at the wavelength λo is given by Eq. 19:
Taking the derivative of Lair with respect to λo and using the definition given by Eq. 15:
Therefore the change of λo with respect to the change of Lair can be written as
Thus, the relationship between a change in the central (balanced) wavelength and the change in the air path length (Lair) is given by:
The minimum amount by which the air path (27) can change sets the minimum increment of the central wavelength in the interferogram. This amount must be several times smaller than the bandwidth of the source. Thus the minimum step size of the air path (27) sets the wavelength resolution of the measured dispersion curve. Note the wavelength resolution is also inversely proportional to the dispersion-length product of the test fiber (17).
The dependence of the wavelength resolution on the dispersion length product will now be shown. As a numerical example, for a step size of 0.1 μm, assuming a 50-cm-long SMF28™ test fiber, the wavelength resolution is 0.1 nm, which is sufficient for most applications. As a graphical example the wavelength resolution is plotted against the dispersion-length product of standard SMF28™ fiber.
1.3.2 Minimum Required Source Bandwidth
A minimum number of envelope fringes are required for accurate measurements of dispersion. As long as the balanced wavelength, λ0, and four other wavelengths corresponding to the peaks (or troughs) of the envelope fringes are captured within the source bandwidth, Bsource, (
For a given test fiber, the dispersion-length product is fixed. Therefore, the only factor that limits the number of envelope fringes is the source bandwidth, Bsource. The longer the test fiber (17), or the larger the dispersion, the more closely spaced the envelope fringes, and hence the smaller the required bandwidth. In order to determine Bmin quantitatively, it is required to determine the maximum value for the wavelength spacing (λ2−λ0), as shown in
Applying the definition of dispersion in Eq. 15, the upper bound of the wavelength spacing (λ1−λ0) is:
Next, the wavelength spacing between λ1 and λ2 is examined. From 4-5, ignoring the 3rd-order term and applying Eq. 15 gives:
Combining Eqs. 24 and 25, the upper bound for the wavelength spacing λ2−λ0 is:
The minimum required source bandwidth Bmin should be not less than the upper bound of 2(λ2−λ0), therefore,
It is clear that the dispersion-length product of the test fiber (17) also affects the minimum required bandwidth. Using a similar numerical example, assuming a 50-cm-long SMF test fiber and 1550 nm as the balanced wavelength, the minimum required bandwidth is 85 nm.
1. 3. 3 Measurable Bandwidth of the Dispersion Curve Bmea
Since each spectral interferogram produces one dispersion value at the balanced wavelength, λ0, to obtain dispersion versus wavelength, a number of interferograms are recorded at various balanced wavelengths by setting the appropriate air path lengths (27). Since each interferogram should be taken over a bandwidth of at least Bmin, from
Alternatively, if two of the troughs are not required to be on each side of 4, then the measurable bandwidth Bmea can be larger. In order to accurately determine λ0, the central fringe (from λ−1 to λ1 in
Either equation Eq. 28 or 29 give the lower bound for the measurable bandwidth, which assumes the widest possible central fringe. In practice, since φenvelope(λ0) cannot be controlled, the width of the central fringe can be anywhere between zero and twice the limit of Eq. 29. Therefore, Bmea can be as large as Bsource in certain cases.
Examination of Eq. 28 or 29 shows that increasing the dispersion-length product of the test fiber (17) increases Bmea. Note that for a given measurement system, Bsource is fixed, so the only parameter that can be used to extend Bmea is Lf. The dispersion length product is, in fact, the main independent variable in determining the system parameters.
The dispersion length-product has been shown to be the main independent variable in determining the measurable bandwidth and the minimum bandwidth. But the range of this parameter is itself affected by the source used. The bandwidth of the source determines the minimum fiber length that can be characterized using this technique and the minimum wavelength step of the source leads to a maximum characterizable fiber length. The next section discusses how the source bandwidth and minimum wavelength step size affect the range of fiber lengths that can be measured using the SAI technique.
1. 3. 4 Minimum Fiber Length
The bandwidth of the source determines the minimum fiber length that can be characterized using SAI. A smaller fiber length produces a wider spectral interferogram as determined by Eq. 27. Thus in order for a certain fiber length to be characterizable using SAI the interferogram produced must fit inside the source bandwidth. Therefore the requirement is that,
Bmin≦Bsource Eq. 30
Using Eq. 27:
Note that for a longer fiber there will be a greater measurement bandwidth (according to Eq. 28 or 29) and a higher wavelength resolution (Eq. 22). As a numerical example, for a source bandwidth of 130 nm, the minimum length for a SMF28 fiber is 0.23 m. The minimum fiber length is plotted as a function of the source bandwidth in
1. 3. 5 Maximum Fiber Length
The SAI method uses the slow-varying envelope function to obtain dispersion. Though the “carrier” fringes are not of interest, they still need to be resolved during measurement otherwise the envelope shape cannot be preserved. The carrier fringe spacing is directly affected by the length of the fiber under test, Lf. A longer fiber will lead to narrower carrier fringes.
The minimum step size of the tunable laser, however, sets a limit on the minimum carrier fringe period that can be detected due to aliasing. Since a longer fiber length has a higher frequency carrier this minimum detectable fringe period results in a limit on the maximum fiber length. The carrier fringe period is the wavelength difference that causes the fast varying phase to shift by 2π. The Fast phase term in Eq. 11 for a balanced air path, Lair=Ng(λo)Lf, can be written as:
φ=(koneffLf+koNg(λo)Lf) Eq. 32
Using a first order approximation of neff and Ng
Ng(λo)≈neff≈n Eq. 33
Where n is the core index, the phase term is written as
The fringe period, Δλ, corresponds to a 2π phase shift
Hence,
In order to detect one fringe accurately, the Nyquist criterion is applied where at least 2 sample points have to be included in one fringe. This sets the following limit over the fiber length:
where Δλ is the minimum wavelength step size of the tunable laser.
If the fiber length limit is exceeded aliasing occurs.
The preceding analysis assumes that it is necessary to avoid aliasing to ensure that all of the peaks of the interferogram are sampled in order to accurately plot the envelope of the interferogram. It is this assumption that leads to the upper limit in the fiber length given in Eq. 37. This upper limit however can be exceeded by dividing the interferogram into small window sections and selecting a single point in each window to plot the envelope. The theory behind this technique, called wavelength windowing, will be explained in detail in the next section.
1. 4 The Effect of Wavelength Windowing
The problem with trying to measure a fiber longer than Eq. 37 allows is that the period of the carrier gets shorter with increasing fiber length. According to Nyquist theory the sampling period, determined by the average step size of the tunable laser, must be at least 2 times smaller than the period of the carrier in order to avoid aliasing. This ensures that all the sampled peaks of the carrier match the true envelope of the interference pattern.
Aliasing is a phenomenon that prevents every peak of the carrier from being sampled but it does not mean that some of the peaks in a given wavelength window range will not be sampled. The interferogram can be divided into small window sections, as shown in
Detailed statistical analysis (developed in the next section) shows how the probability that at least one of the peaks will be sampled within a wavelength window is determined. This technique shows that the upper limit in Eq. 37 can be exceeded by many folds by wavelength windowing.
1. 5 Model Development
This technique uses a tunable laser (31) system to sample the peaks of an interferogram. A real world tunable laser (31) system, however, does not step the wavelength with equal step sizes but has a certain standard deviation in its step size. In order to produce an accurate modeling of a real world process this variation in the step size of the tunable laser (31) must be taken into account by the model. The tunable laser (31) system used in the experiments was the Agilent 8164A™ which has an average step size of 1 pm and a standard deviation of 0.17 pm as determined from the histogram and the Gaussian probability density function (PDF) in
In order for the model to accurately determine the probability of a sampled point matching at least one peak of the carrier wave within a certain wavelength window, certain parameters must be determined. The model that will be developed requires knowledge of the fiber length, the width of wavelength window, the average step size of the tunable laser, the standard deviation of this step size and the tolerance in detecting the peak as a percentage of the carrier period.
In this model the fiber length as Lf, the wavelength window within which it is desired to detect a peak as W, the average step size of the tunable laser as μ, the standard deviation of the step size of the tunable laser as σ and the tolerance in detecting the peak as a percentage of the carrier period as ε can be designated. If λo is the separation between the first carrier peak and the maximum sampling probability density of the first step, as shown in
Assuming the location of the first carrier peak to be at λ1, as shown in
Therefore the probability that the first peak is sampled by the first step is:
Here ε, shown in
If N is the number of steps of the tunable laser in a given window size then the probability of not sampling the first peak with any of the N steps is given by:
If M is the number of peaks of the carrier in a given window size then the probability of not sampling any of the M peaks with any of the N steps is given by:
Where λm is the location of the mth peak in the wavelength window and is given by m ki and Λ+ and Λ− are the normalized wavelength parameters given by:
Since the model assumes a fixed relationship between the first carrier peak and the maximum of the probability density function this probability should be averaged for λo varying over one carrier wave period. This gives the probability that no carrier peak is sampled in a given window for a random alignment between the carrier peaks and the maximum of the probability density function. The result is given as:
Thus the probability that at least one of the peaks is sampled for a given window size is determined as:
1.6 Simulation Results
1.6.1 Probability vs. Fiber Length
The probability that at least one peak is sampled in a given window size, W, is shown in
2: Experimental Results
2.1 Experimental Process
The first step in the experiment is to set up the Single Arm Interferometer and to assemble the control and data acquisition hardware (29). The second step in the experiment is to test the technique by using it to measure the dispersion of test fibers (17) for which the dispersion curves are known or that can easily be measured using conventional techniques. To do this, the dispersion curves of Single Mode Fiber (SMF28™) and Dispersion Compensating Fiber (DCF) were measured. After careful analysis of the results for the experiments on SMF28™ and DCF the new technique was then used to measure the dispersion of a test fiber (17) that has never before been characterized (twin-hole fiber). The entire experimental process for this project is outlined in
2. 2 Experimental Instrumentation & Specific Limits
The experimental set up is shown in
Optionally, the launch fiber (15) may be configured as a polarization controller. This may be accomplished by a number of means known to those skilled in the art, such as looping the launch fiber (15) into three coils. Alternately, the polarization controller may be implemented by placing a linear polarizer in the air path (27). By providing a polarization controller, the SAI may be operable to measure polarization mode dispersion.
2. 3 Experiments
2.3.1 Single Mode Fiber
The dispersion properties of SMF28™ are well known and hence it was used to verify the theory of single arm interferometry. In this experiment a 39.5±0. 1 cm piece of the SMF28™ fiber was used in a SAI in order to characterize its dispersion.
Where λ0=1313 nm and So=0.086 ps/nm-km and D(λ) is measured in ps/nm-km.
The wavelength resolution of the measured dispersion curve, as determined by Eq. 22, is 2.4 nm. The measurable bandwidth according to Eq. 29 is 30 nm, which is the bandwidth actually used, as shown in
The simulated interference pattern is generated using Eq. 17 and the envelope of the interference pattern is generated using Eq. 18. In the simulation a test fiber (17) length of 0.395 m is assumed in order to match the experimental conditions. The path length of the air path (27) is determined via a calculation of the effective group index of the fiber was determined to be 1.472469 at the central wavelength, λo, via Eq. 46:
Where
κ(λ)=√{square root over (ncore(λ)2−neff(λ)2)}{square root over (ncore(λ)2−neff(λ)2)}
γ(λ)=√{square root over (neff(λ)2−ncladding(λ)2)}{square root over (neff(λ)2−ncladding(λ)2)} Eq. 47
Note that a is the core size of the fiber and J and K are Bessel functions of the first and second kind. The locations of equality in Eq. 46 determine the values of κ(λ) and γ(λ) as well as a mode of the fiber. The first of these modes is called the fundamental mode of the fiber. The values of ncore(λ) and ncladding(λ) are the index of bulk glass with the composition of the core and cladding respectively. The effective group index as a function of wavelength in SMF28™ fiber is determined using simulation.
2.3.2 Dispersion Compensating Fiber
As a second method of verification, dispersion was measured on a short piece of DCF, whose dispersion value is approximately one order of magnitude higher than that of SMF28™, and has an opposite sign. A 15.5±0.1 cm piece of DCF fiber was used, and the measurement results are given in
The standard deviation of the measured dispersion is calculated by taking the difference between the measured points and a linear fit and then determining the standard deviation of the difference. The standard deviation of the measured data is 0.99 ps/nm-km, which corresponds to a relative error of 0.9%. When multiplied by the length of the test fiber (17), this translates into a standard deviation of 0.00015 ps/nm.
Since DCF has negative dispersion values a procedure for determining the sign of the dispersion was developed. By examination of Eq. 22 repeated below for convenience
If the sign of the dispersion is negative then the location of the central wavelength will decrease as the path length of the air path (27) is increased. This is a quick method for determining the sign of the dispersion.
2.3.3 Twin-Hole Fiber
Twin-Hole Fiber (THF) has been used in fiber poling to facilitate parametric generation in fibers or making fiber-based electro-optic switching devices. In such nonlinear applications, dispersion of the fiber is an important parameter to be determined. The dispersion properties of THF, however, have never been reported. This is partly due to the lack of uniformity in the diameter of the THF along its length. The fiber has a 3-μm-diameter core and a numerical aperture that is higher than that of SMF28™. The cross section of a typical THF is shown in
The core is Ge-doped silica, and has an index similar to that of SMF28™. Therefore, the dispersion of THF is expected to be slightly lower than that of SMF28™. Without knowing the magnitude of the dispersion for THF, it is best to choose the largest length of THF available to increase the chance that the minimum bandwidth required for a measurement would fit in the available bandwidth of the tunable laser source. The largest length of THF available was 45±0.1 cm. This length of fiber is slightly longer than the length allowed by Eq. 37 but since the technique of wavelength windowing described in sections 4.5.1-4.5.3 was used the measurement of the envelope was still possible in this experiment.
The measurement results from the experiment on THF are given in
An important aspect of the previous three sections is the error associated with the measurement of each point in the dispersion parameter plots. The next section outlines the source and magnitude of the error associated with the measurement of the dispersion parameter.
In conclusion, the experimental results of Single Arm Interferometry confirm the accuracy of the present invention. They show that the dispersion parameter can be calculated from the envelope of the fringe pattern produced by the interference of 3 waves in a balanced SAI. The experiments on Single mode fiber (SMF28™) and Dispersion Compensating Fiber (DCF) were used to confirm the theory behind the technique and once the technique was confirmed it was used to measure the unknown dispersion parameter plot for THF. The length of twin-hole fiber used in the experiment was larger than allowed by Eq. 37 so the technique of wavelength windowing, described in sections 4.5.1-4.5.3, had to be used. This technique was shown theoretically and via simulation to extend the maximum length of fiber that can be characterized by this technique. Ultimately the largest length of fiber that can be characterized is limited by the largest air path (27) that can be produced in the experiment and the laser linewidth.
As mentioned, the resolution of balanced spectral interferometry, in particular, can be improved by replacing the combination broadband source and Optical Spectrum Analyzer with a tunable laser (31) and detector (23) system. Current tunable laser technology allows for a bandwidth of 130 nm and a 1 picometer resolution. This improves the range of fiber lengths that can be measured using this technique. Also of note is that the use of tunable lasers for dispersion measurement is becoming more widespread as they decrease in cost.
3: Benefits
3. 1 Significance to Research
The single arm interferometer is an alternative to the Michelson or the Mach-Zehnder configuration for interferometric measurements of the dispersion parameter. It is most useful for measurements of the dispersion parameter in short lengths of fiber. The technology may be used to eliminate the need for the arm balancing required by dual arm interferometers and by doing so allow for greater ease in the use of interferometric dispersion measurement techniques.
The new interferometer is significant for researchers since it can be studied and used alongside the earlier types of interferometers like the Michelson, the Mach-Zehnder and the Fabry Perot. This new interferometer provides researchers with another tool for studying dispersion in short length fibers and waveguides which will be useful in the development of specialty fibers. These specialty fibers require simple and accurate short length characterization since they are generally made in very small quantities and their geometry tends to vary as a function of position along the fiber.
3. 2 Significance to Industry
The new interferometer is significant to industry since it minimizes the need to compensate for unwanted reflections by eliminating the need for a coupler altogether. As a result this interferometer is a simpler (less expensive) interferometric dispersion measurement device capable of characterizing the dispersion of short length optical fiber. As a result it is a viable commercial competitor to the current Modulation Phase Shift (MPS) based devices currently on the market. The new interferometer, however, has an advantage over MPS based devices since it has the ability to measure short length fiber with high accuracy.
Also, since it can measure short lengths of fiber it has the ability for another type of measurement as well. Dispersion is a function of both material and dimensional (waveguide) properties of a fiber but if the dimensions, particularly the diameter of the fiber, vary then the dispersion will vary. If several small sections can be cut from various points on a long length fiber and the dispersion is measured in each of them then the variation in the dispersion can be plotted as a function of position in the fiber. This can then be directly related to the variation in the fiber diameter. The main point here is that a great deal of accuracy in measuring the fiber diameter can be achieved by measuring it indirectly via the dispersion and it would be an easy way for a fiber drawing company to perform quality control.
Greater commercial interest in this device will enable measurement of dispersion in smaller lengths of fiber since larger bandwidth tunable lasers will be developed. Also the advancement in the speed of the tunable laser and scanning process will make each measurement faster to obtain.
3.3 Module
One of the most interesting features of a single arm interferometer is the ease with which it can be built. This ease of construction lends itself very nicely to economical commercial assembly of a dispersion measurement device, such as an add-on module for a tunable laser system to allow it to make dispersion measurements. A conceptual design of such a module is illustrated in
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/CA2007/001824 | 10/15/2007 | WO | 00 | 7/5/2010 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2009/049393 | 4/23/2009 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
6137573 | Luke et al. | Oct 2000 | A |
6804008 | Morison et al. | Oct 2004 | B1 |
7227645 | Cyr | Jun 2007 | B2 |
Number | Date | Country | |
---|---|---|---|
20100296102 A1 | Nov 2010 | US |