The present invention relates generally to a system and a method for creating and utilizing “multivariate paths” with multi-dimensional control to attain single mode operation in an electromagnetic radiation source for use in optical coherence tomography (OCT), optical frequency domain reflectometry (OFDR), spectroscopy, remote sensing, Light Detection and Ranging (LIDAR), and testing of optical components, etc.
Lasers are sometimes operated as wavelength-sweeping devices in remote sensing, LIDAR, and optical coherence tomography applications, as well as to test telecommunications components, among other applications. Discontinuities in the wavelength sweep or operation of the wavelength outside of a single mode can significantly affect the application that the laser is being used in. For example, the shape of a molecular gas absorption feature may be distorted by a discontinuity—a forward or backward jump—or operation of the wavelength outside of a single mode in the wavelength sweep of the laser. In another example, wavelength discontinuities or operation of the wavelength outside of a single mode can reduce the signal-to-noise of an OCT measurement of tissue, or a LIDAR (Laser RADAR) measurement using the Frequency Modulated Continuous Wave (FMCW) technique or more complex lidar/radar waveforms. Thus, it is usually desirable to eliminate from swept-wavelength lasers wavelength discontinuities, wavelength non-linearity, and operation of the laser outside of a single mode.
Attempts in the prior art to maintain single mode operation and to control a sweep profile for a laser are numerous, but unsatisfactory. Although it may be possible to carefully maintain single mode operation and control a sweep profile using a calibrated parameter set at an initial point in time (often referred to as pre-emphasis or predistortion), the passage of time or changes in, e.g., temperature or humidity will create discontinuities, non-linearities, and cause operation outside of a single mode. For example, mechanically-tuned external cavity lasers operate in near continuous single-mode using an external cavity mechanism coupled with a gain medium. Another example is the use of pre-emphasis to control a semiconductor laser.
In a typical single mode tunable laser, there are two key elements; a method for changing the wavelength (either internal to the cavity, or external), and another for altering the cavity length to optimize side mode suppression and to maintain single mode operation. In an analogous tunable laser, known as Littman-Metcalf configuration, there is a specific mechanical configuration which constrains the change in wavelength to happen coincident with a commensurate change in cavity length, thus maintaining single mode operation. In these mechanical systems, there is a mechanical construction which constrains the mechanical “path” which is traversed to one in which the wavelength increase is linear, and the path length difference is simultaneously changed in concert with the wavelength increase to maintain single mode operation with good SMSR. Further, in a laser that uses external feedback, the current through the gain medium is another parameter to optimize the SMSR or the wavelength tuning, making such lasers tunable via a multivariate parameter space.
Operation of the mechanical laser is maintained through accurate, tightly-toleranced components and precision alignment of the cavity, or using real-time elements such as piezoelectric transducers that adjust the cavity length in real-time. Other laser configurations use an intra-cavity element. Over time, however, the alignment of the laser degrades or the components wear, which may cause changes in the sweep profile versus time and operation outside of a single mode. As the ambient temperature, humidity, or pressure change, the alignment can degrade, which can also cause changes in the sweep profile versus time and operation outside of a single mode. Vibrations external or internal to the laser may also misalign the cavity, which again may cause changes in the sweep profile versus time and operation outside of a single mode.
Even in lasers with stable cavities, it is difficult to create wavelength sweeps (for example, using pre-emphasis) without wavelength discontinuities. Monolithically-constructed semiconductor lasers, or non-semiconductor monolithic lasers in general, are a class of single-mode laser for producing swept wavelengths. Monolithic semiconductor lasers include several sections or segments in the semiconductor, which serve, for example, as adjustable cavity mirrors, laser gain, coupled cavities, cavity phase and (optionally) external amplification. Examples are Vertical Cavity Surface Emitting Lasers (VCSELs), VCSELs with Micro-electromechanical systems (MEMS) tuning structures, Vernier-tuned Distributed Bragg Reflector (VT-DBR) lasers, Vernier-tuned ring lasers, Y-branch lasers, coupled cavity lasers, discrete mode lasers, injection-locked or externally-stabilized lasers, Super-Structure Grating Distributed Bragg Reflector (SSGDBR) lasers and similar devices. Because these lasers are typically monolithic with no moving parts (excepting the MEMs devices), their cavities and associated optics are extremely stable and can operate in single-longitudinal mode with narrow linewidth and long coherence length. Tunable monolithic lasers of this class require multiple control signals to tune the wavelength, presenting a challenge to creating wavelength sweeps without wavelength discontinuities.
There is a need for a system and method for determining and controlling an electromagnetic radiation source to generate rapid, highly repeatable, time-dependent, prescribed multi-wavelength patterns conformal to an optimized single mode operation, over the passage of time and without user intervention.
It is desirous to produce a multi-wavelength, time-sequenced pattern from a laser system that is conformal to an optimized single mode operation. The ratio between the main spectral mode amplitude and the most dominant side mode amplitude—i.e., the side mode suppression ratio (SMSR)—is used herein to characterize whether the laser is operating within a single mode. One such multi-wavelength pattern is known as a wavelength sweep. A wavelength sweep, as used in this disclosure, refers to a continuous (or step-wise continuous) movement over time from one wavelength to another, in a prescribed (and often preferably in a linear, monotonic) fashion. A wavelength sweep can represent the entire desired waveform, or it can represent one or more elements of a more complex waveform.
As would be understood by one of ordinary skill in the art, the optical frequency and optical wavelength are inversely proportional to each other and a wavelength sweep or movement in wavelength is also an optical frequency sweep or movement in optical frequency.
Another multi-wavelength pattern is referred to as wavelength hopping. Wavelength hopping, as used in this disclosure, is similar to wavelength stepping (as described below) but with much shorter periods of time for the continuous operation, or dwells, at any given wavelength. Hopping refers to a sharp transition (short time interval) between two distinct wavelengths. A sequential hopping pattern with small wavelength differences represents one method for obtaining a sweep, for instance. For example,
Yet another multi-wavelength pattern for a laser is known as wavelength stepping. Wavelength stepping, as used in this disclosure, refers to continuous operation at one wavelength, followed by continuous operation at another wavelength, possibly followed by continuous operation at a third wavelength, etc.
In
Yet another sweep pattern that can be realized with this method is wavelength modulation.
Finally, the wavelength sweep pattern can deviate substantially from a linear trend, and can instead be quite nonlinear in nature. An example of a non-linear pattern is shown in
A laser may also undergo a wavelength hop at a non-prescribed wavelength, and at a generally indeterminate time, to a generally indeterminate (but nearby) wavelength. This undesirable form of a wavelength hop is known in the industry as a ‘mode hop’. The interaction of multiple input parameters determines the output wavelength of a laser. Preventing such undesirable mode hops in the laser output thus requires understanding how to adjust the multiple parameters across the wavelength set, or wavelength range. A further complication to creating and maintaining a wavelength pattern is that time and temperature affect tuning in meaningful amounts. This means that for precise tuning accuracy, the Parameter Sweep Paths (PSPs) must also change over time. It is an aspect of this invention to adaptively measure and precisely control these PSPs as they change over time due to any factors, and without user intervention.
As described above, multiple parameters interact to control a laser output. Describing the relationship between the parameters is difficult due to the number of parameters, the interactions between the parameters, the dependence of the parameters on environmental factors, and the effect of the components' age on laser output.
The output wavelength is determined based on the parameters and the external factors mentioned above, through the expression:
where λ is the wavelength at a particular time tj and where ci are the m parameters, which interact to influence the desired output wavelength. The parameters ci might include, e.g., for a semiconductor monolithic tunable laser source (SMTLS), controllable factors such as front mirror current, back mirror current, phase current, gain, semiconductor optical amplifier current, and external factors such as feedback wavelength, temperature and humidity.
As described above, the degree to which a laser operates in a single longitudinal mode of the laser cavity is often measured by the SMSR of the laser. The SMSR is also governed by the multiple control parameters of the laser, as well as the external factors mentioned above, through the expression:
Attempts to control a laser by developing an overarching model that relates the control parameters to the laser output and holding the external factors constant (so called ‘pre-emphasis’) has not been successful in achieving the desired results, because of three factors. First, external factors cannot be held sufficiently constant. Second, because a laser is a non-linear system, the control parameters ci may have a tremendously complex interaction in determining the wavelength, SMSR, power or other performance metric of the laser. Third, the time-dependency necessary to control the laser to achieve specific outcomes is very complicated, and changes over time as the external factors (temperature, humidity, vibration) change and as the laser or its control systems age.
Other methods for obtaining longer-term control over a wavelength sweep, at least for a linear sweep, have been disclosed previously. Those methods involve an opto-electronic closed-loop approach (optical phased locked loops, OPLLs) and are not the subject of this invention. Examples of such approaches are described in U.S. Pat. Nos. 5,751,747; 6,034,976; 8,175,126; 7,995,630; and 8,913,636. Such methods produce good linearity, but suffer from speed limitations, require additional hardware, have complex closed loop operation, and are extremely difficult to adapt to non-linear types of waveforms such as frequency stepping and hopping. A limitation of the prior art methods is that when the relationship between laser control parameter and laser performance becomes sufficiently non-linear, the closed loop approach breaks down. Another limitation of prior art methods of closed loop control is that as the laser's response to the control parameter changes over time (due to aging of the laser, aging of the control electronics or mechanics, or due to changes in the environment of the laser) the closed-loop approach—which relies on a prior knowledge and stability of the control response relationship—breaks down.
In one embodiment, it is an aspect of this invention that control over the laser's PSPs does not require analog closed loop operation via a phase locked loop. It is another aspect of the present invention that control over the laser's PSPs may be not limited by feedback time of the control circuit, the measurement electronics, or the loop time. This makes certain embodiments of the present invention especially applicable to sweeps with wavelength versus time profiles where the waveform time is very short (less than about 0.5-10 microseconds) and the changes in wavelength during the waveform may happen very quickly (such as within nanoseconds), as may be encountered in the fields of use listed above.
Rather than trying to determine the functions F(ci,tj) or G(ci,tj), which can be very difficult or potentially intractable, it is only necessary to determine the limited set of potential solutions that satisfy the conditions of a desired set of wavelengths and minimum SMSR over a range of wavelengths. If the parameters are viewed as a multivariate space with each point in the space corresponding to a parameter combination causing a laser output having a corresponding wavelength and SMSR, there are only certain “multivariate paths” through the multivariate space that will maintain the single mode operation with optimal (e.g., maximal) SMSR. A multivariate path can be determined by identifying sub-paths that satisfy the conditions of a desired set of wavelengths and minimum SMSR over a subrange of the range of wavelengths. The sub-paths can be combined and traversed along in a prescribed manner (for example, in linear increasing wavelength or linear increasing optical frequency) as a multivariate path through the multivariate space defined by the parameters (=a wavelength sweep). Alternatively, the sub-paths can be traversed using a combination of wavelength sweeping (e.g., where the sweeping may be performed continuous or through a series of discrete steps), wavelength sub-sweeps (e.g., where the sweeping may be performed continuous or through a series of discrete steps), wavelength stepping, and wavelength hopping. Further, these combinations can involve sequential or adjacent wavelength transitions, or they may involve random selections of time-sequenced wavelengths within the range. Further, the instantaneous wavelengths between steps or hops can be either step wise (e.g., visually representing a staircase) or a continuous smooth monotonic transition from point to point.
The present invention provides a system and method for controlling a monolithic laser source to produce a single mode operation having an optimized side-mode suppression ratio over a set of wavelengths by adaptively determining a multivariate path for transitioning from one wavelength to another. In this way, and without user intervention, the invention may periodically determine and update an appropriate set of sub-paths (PSPs) for transitioning between wavelengths that satisfy a condition for a desired set of wavelengths and a minimum side mode suppression ratio rather than attempting to (e.g., a-priori) determine the relationship between the inputs, output wavelength, and side mode suppression ratio over all possible conditions.
According to one aspect, a method is presented for controlling an electromagnetic radiation source to produce a single mode operation having an optimized side-mode suppression ratio over a set of wavelengths within a prescribed temporal profile. The electromagnetic radiation source is configured to output electromagnetic radiation at a given wavelength based upon parameters. The method includes determining a set of parameter combinations that satisfy a condition for a desired set of wavelengths and a minimum side mode suppression ratio over the range of wavelengths. The set of parameter combinations define sub-paths for transitioning from one wavelength to another wavelength. Combinations of select sub-paths provide a multivariate path (PSP) for transitioning over the range of wavelengths. The method also includes adaptively controlling the laser, over time and other conditions, to emit electromagnetic radiation over the range of wavelengths by traversing the multivariate path in a desired manner. The desired manner is prescribed and time-sequenced, may include sweeps, hops, or steps, may be repeating in time, and the laser outputs the electromagnetic radiation in a single mode having the optimized side-mode suppression ratio over the range of wavelengths.
A number of features are described herein with respect to embodiments of the invention; it will be appreciated that features described with respect to a given embodiment also may be employed in connection with other embodiments.
The invention includes the features described herein, including the description, the annexed drawings, and, if appended, the claims, which set forth in detail certain illustrative embodiments. These embodiments are indicative, however, of but a few of the various ways in which the principles of the invention may be employed.
In the annexed drawings:
The present invention provides a method for adaptively controlling an electromagnetic radiation source to produce a single mode operation having an optimized side-mode suppression ratio (SMSR) over a set of wavelengths within a prescribed temporal profile. The electromagnetic radiation source is configured to output electromagnetic radiation at a given wavelength based upon parameters. The method may include periodically determining a set of parameter combinations without user intervention that satisfy a condition for a desired set of wavelengths and a minimum side mode suppression ratio over the range of wavelengths. The set of parameter combinations define sub-paths for transitioning from one wavelength to another wavelength. Combinations of select sub-paths provide a multivariate path for transitioning over the range of wavelengths. The method also includes controlling the laser to emit electromagnetic radiation over the range of wavelengths by traversing the multivariate path in a desired manner. The desired manner is prescribed and time-sequenced, may include sweeps, hops, or steps, may be repeating in time, and the laser outputs the electromagnetic radiation in a single mode having the optimized side-mode suppression ratio over the range of wavelengths.
An exemplary system 10 in accordance with aspects of the present disclosure is illustrated in
As will be understood by one of ordinary skill in the art, the controller 12 may have various implementations. For example, the controller (also referred to as circuitry) may include a processor or any other suitable device, such as a programmable circuit, integrated circuit, memory and I/O circuits, an application specific integrated circuit, microcontroller, complex programmable logic device, other programmable circuits, or the like. The controller may also include a non-transitory computer readable medium, such as random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), or any other suitable medium. Instructions for performing the method described below may be stored in the non-transitory computer readable medium and executed by the processor.
The controller 12 may, in additional to controlling the electromagnetic radiation source 14, function to determine a multivariate path for driving the electromagnetic radiation source 14 to output electromagnetic radiation 20 transitioning over a range of wavelengths. The controller controls the electromagnetic radiation source 14 to emit electromagnetic radiation 20 over the range of wavelengths by traversing the multivariate path. The controller 12 may also control one or more output devices to provide information related to the system and/or method to an associated user and/or cause information to be output to another device for additional processing or display, for example.
The system 10 may also include an optical system 16 and an electronic system 18. The optical system 16 may receive a portion of the outputted electromagnetic radiation 20 and generate signals or data in response to the wavelength, SMSR, optical power, or other properties of the electromagnetic radiation 20. The optical system 16 may include an interferometer. The optical system 16 may receive electromagnetic radiation output by the electromagnetic radiation source 14 even when the electromagnetic radiation is prevented from exiting the electromagnetic radiation source 14 by a shutter (e.g., a structure configured to prevent electromagnetic radiation from exiting the electromagnetic radiation source 14).
The optical system 16 may be coupled to the controller and the controller may analyze the data or signals for adaptively determining the wavelength and SMSR of the output electromagnetic radiation 20 in response to the parameters supplied to the electromagnetic radiation source 14. The electronic system 18 may measure the parameters at the electromagnetic radiation source 14. For example, the electronic system 18 may measure the certain electrical signals which are secondarily correlated to pertinent optical characteristics of the output of the electromagnetic radiation source 14. The electronic system 18 may be coupled to the controller 12 and the controller 12 may analyze the parameters as received by the electromagnetic radiation source 14 to determine the performance factors at the electromagnetic radiation source 14 and to detect defects in the current multivariate path being operated, and to periodically update the multivariate path.
An exemplary multivariate path 40 for an exemplary VT-DBR semiconductor laser is shown in
As will be understood by one of ordinary skill in the art, the multivariate path is not limited to two controllable parameters. Rather, the multivariate path may comprise a single controlled parameter (e.g., a current or voltage) and one or more uncontrolled parameters (e.g., humidity and/or temperature).
With reference to
As described above, the wavelength and SMSR of the electromagnetic radiation 20 output by the electromagnetic radiation source 14 for a given set of parameters may change with time due to changes in humidity, temperature, equipment misalignment, component damage, component misalignment, etc. For this reason, a multivariate space 42 may be generated if a multivariate space 42 has not been recently generated. The multivariate space 42 may be generated just prior to the determination of a multivariate path 40 and the generation of a sweep with the electromagnetic radiation source 14. A multivariate space 42 may also be generated if the wavelength and SMSR have not been measured for parameter combinations in the multivariate space 42.
A multivariate space 42 may be generated by supplying a parameter combination corresponding to each point in the multivariate space 42 and measuring the output wavelength and SMSR of the electromagnetic radiation source 14. A multivariate space 42 may also be generated by measuring (e.g., with the optical system 16) the output wavelength and SMSR of the electromagnetic radiation source 14 for a subset of the parameter combinations in the multivariate space 42. The measured wavelength and the measured side mode suppression ratio may then be incorporated into the multivariate space at the point in the multivariate space designated by the parameter combination.
The parameter combinations included in the multivariate space may be set by default or user defined. Also, when the electromagnetic radiation source is being controlled to emit electromagnetic radiation by traversing the multivariate path in the desired manner, the emitted electromagnetic radiation may be measured by the optical system 16 and used to update the multivariate space. That is, when a given parameter combination is already incorporated into the multivariate space, the multivariate space may be updated to include the measured wavelength and the measured side-mode suppression ratio.
The generation of the multivariate space may be performed according to different criteria. For example, the generation of the multivariate space may performed when the wavelength and the side-mode suppression ratio for a particular supplied parameter combination has not been measured for a given duration of time (e.g., one hour, 24 hours, or one month) or when a user requests generation of the multivariate space. Alternatively or additionally, the generation of the multivariate space may be performed when an ambient temperature at which the wavelength and the side-mode suppression ratio for a given parameter combinations was measured changes by at least a predetermined amount (e.g., five degrees Celsius or ten degrees Celsius). Alternatively or additionally, the multivariate space may be generated when a specified set of parameter combinations does not satisfy the condition for the desired set of wavelengths within a desired tolerance (e.g., within 1%, 3%, 5%, or 10%). As an example, the multivariate space may be generated when the measured side-mode suppression ratio and/or the time sequence of wavelengths emitted by the electromagnetic radiation source deviates by more than the desired tolerance.
Generating of the multivariate space may be performed by measuring with the optical system 16 the emitted electromagnetic radiation over the range of wavelengths when traversing the multivariate path in the desired manner. Alternatively or additionally, generating of the multivariate space may be performed during a time duration when the electromagnetic radiation source is not outputting electromagnetic radiation. For example, a shutter of the electromagnetic radiation source 14 may be closed such that electromagnetic radiation is not emitted externally by the electromagnetic radiation source 14. While shuttered, parameter combinations may be supplied to the electromagnetic radiation source 14 and the output of the electromagnetic radiation source 14 may be measured with the optical system 16. That is, even though the electromagnetic radiation is not being emitted to the external environment by the electromagnetic radiation source 14, the generated electromagnetic radiation may be received by the optical system 16. For example, the wavelength and the side-mode suppression ratio of the electromagnetic radiation generated by the electromagnetic radiation source 14 may be measured by the optical system 16.
As another example, the multivariate space 42 of
In process block 112, a set of parameter combinations are determined that satisfy (1) a condition for a desired set of wavelengths in the range of wavelengths and (2) a minimum SMSR over the range of wavelengths. The set of parameter combinations define sub-paths for transitioning from one wavelength to another in the wavelength range. For example, a parameter combination—i.e., a point—in the multivariate space 42 is determined that results in the output of electromagnetic radiation 20 having a wavelength A and a minimum SMSR over the range of wavelengths. A set of such parameter combinations can be determined that define a sub-path for transitioning from one wavelength to another wavelength in the range of wavelengths. In
In process block 114, sub-paths are selected from the available sub-paths. Sub-paths may be selected using any suitable criteria. For example, sub-paths may be selected based on the average SMSR of the sub-path, maximum SMSR of the path (e.g., only sub-paths with a maximum SMSR below a threshold are selected), or a combination thereof. Sub-paths may also be excluded from selection if the sub-path corresponds to a range of wavelengths smaller than a range threshold. In process block 116, the selected sub-paths are combined to provide a multivariate path for transitioning over the range of wavelengths. Although not easily represented on paper, a very small subset of all combinations of the parameters form distinct, discontinuous, paths through the multivariate space. Following these paths (e.g., combining the sub-paths into a multivariate path), it is possible to progress through the desired wavelengths in a desired fashion while satisfying the other conditions. The selected sub-paths may be combined by combining the parameter combinations in the select sub-paths. If the selected sub-paths overlap in wavelength range, the overlapping sub-paths may be shortened to remove the overlap. The sub-paths may be shortened by keeping only the overlapping parameter combinations with the highest SMSR.
Finally, in process block 118, the electromagnetic radiation source 14 is adaptively controlled by the controller 12, optical system 16, and electronic system 18 to emit electromagnetic radiation 20 over the range of wavelengths by traversing the multivariate path in a desired manner. The desired manner may for example be a linear sweep, and may therefore be linear increasing wavelength or optical frequency output in a single mode having the optimized side-mode suppression ratio over the range of wavelengths.
Following process block 118, the method may return to optional process block 110 or process block 112. This process flow provides a closed loop around the generation of the sweep and allows for the use of multiple parameters non-linear and even non-monotonic relationships between the control parameters and the optical frequency and SMSR of the emitted electromagnetic radiation. In this way, the method is not limited to a single control parameter combination with a clear and static functional relationship to the behavior of the electromagnetic radiation source 14 or the emitted electromagnetic radiation. By not limiting the method to a static control parameter combination, the method accounts for the effects of external temperature, strain, pressure, humidity, etc., and by the aging of the device, or of the control electronics that changes the relationship between the control parameter combinations and the properties of the emitted electromagnetic radiation.
Another example of a method for controlling an electromagnetic radiation source 14 to produce a single mode operation having an optimized SMSR over a set of wavelengths within a prescribed temporal profile is shown in
In process block 214, a set of parameter combinations are determined for each subrange of wavelengths that satisfy (1) a condition for a desired set of wavelengths in the subrange of wavelengths and (2) a minimum SMSR over the subrange of wavelengths. The set of parameter combinations define sub-paths for transitioning from one wavelength to another wavelength within the subrange of wavelengths. Sub-paths for transitioning within a subrange of wavelengths may be determined in the same manner as described previously.
In process block 216, additional untested parameter combinations may optionally be identified and added to the multivariate space 42. For example, a sub-path may be identified that extends into a region of multivariate space that does not contain data regarding output wavelength and SMSR. As another example, a range of parameter values may be identified as possibly connecting two unconnected sub-paths. If untested parameter combinations are identified, the previously untested parameter combination may be input to the electromagnetic radiation source 14, the output wavelength and SMSR measured, and this data added to the multivariate space 42.
In process block 218, sub-paths are selected from the sub-paths as described previously. In process block 220, the selected sub-paths are combined to provide a multivariate path for transitioning over the range of wavelengths. In process block 222, the electromagnetic radiation source 14 is controlled to emit electromagnetic radiation 20 over the range of wavelengths by traversing the multivariate path in a desired manner. The desired manner may be any form of wavelength sweep, as defined and described in additional detail below.
Following process block 222, the method may return to optional process block 210 or process block 212. As described above regarding method 100, this process flow provides a closed loop around the generation of the sweep and allows for the use of multiple parameters non-linear and even non-monotonic relationships between the control parameters and the optical frequency and SMSR of the emitted electromagnetic radiation.
In
In
Although the invention is shown and described with respect to illustrative embodiments, it is evident that equivalents and modifications will occur to those persons skilled in the art upon the reading and understanding hereof. The present invention includes all such equivalents and modifications and is limited only by the scope of the claims if appended hereto.
This application claims the benefit of Application No. 62/417,471 filed on Nov. 4, 2016, which is herein incorporated by reference in its entirety.
Number | Name | Date | Kind |
---|---|---|---|
9455549 | Minneman | Sep 2016 | B2 |
20120320368 | Jiao | Dec 2012 | A1 |
20140341244 | Minneman | Nov 2014 | A1 |
Number | Date | Country | |
---|---|---|---|
62417471 | Nov 2016 | US |