System and method for receiver sensitivity improvement

11212011 · 2021-12-28

Assignee

Inventors

Cpc classification

International classification

Abstract

A system and method for ultrashort signal detection adds an optical weighting element upstream of a detector within a direct detection receiver. The optical weighting element is configured to generate an optical pulse that closely matches at least one ultrashort pulse within the input signal so that portions of the input signal that are nonoverlapping with the at least one ultrashort pulse are rejected.

Claims

1. An optical receiver for ultrafast signal detection, the receiver comprising: a direct detection receiver configured for receiving an input signal generated by a pulse source, the input signal comprising at least one ultrashort pulse within signal noise, the ultrashort pulse having an ultrashort pulse period within a time domain; and an optical weighting function implemented by a four-wave mixer disposed upstream of the direct detection receiver, the optical weighting function configured to impose optical weighting on different portions of the input signal, wherein the input signal is filtered by an optical bandpass filter and combined with a weighting pulse from the pulse source that has a pulse width matching the ultrashort pulse period, wherein one or more portions of the input signal within the ultrashort pulse period are multiplied by a non-zero function and other portions of the input signal that are outside of the ultrashort pulse period are multiplied by zero, and wherein the combined signal is guided into a highly non-linear fiber (HNLF) to generate an idler.

2. The optical receiver of claim 1, wherein the optical weighting function operates according to the relationship f ( t ) = { s ( t ) , t 0 - τ s 2 + nT < t < t 0 + τ s 2 + nT 0 , otherwise , where τ.sub.s is a weighting function duration, t.sub.0 is a pulse center, T is the ultrashort pulse period, and n is an integer.

3. An optical receiver for ultrashort signal detection, the receiver comprising: an optical weighting function implemented by one of an amplitude modulator and a four-wave mixer disposed upstream of an impulse response detector within a direct detection receiver, the optical weighting function configured to generate an optical pulse that closely matches an ultrashort pulse period of at least one ultrashort pulse within an input signal from a pulse source so that portions of the input signal within the pulse period are multiplied by a non-zero function while portions of the input signal that are nonoverlapping in time with the ultrashort pulse period are multiplied by zero.

4. The optical receiver of claim 3, wherein the optical weighting function operates according to the relationship f ( t ) = { s ( t ) , t 0 - τ s 2 + nT < t < t 0 + τ s 2 + nT 0 , otherwise , where τ.sub.s is a weighting function duration, t.sub.0 is a pulse center, T is the ultrashort pulse period, and n is an integer.

5. The optical receiver of claim 3, wherein the optical weighting function is implemented using a four-wave mixer, wherein the input signal is filtered by an optical bandpass filter and combined with a weighting pulse obtained from the pulse source, the weighting pulse having a pulse width that matches the ultrashort pulse period, wherein the combined signal is guided into a highly non-linear fiber (HNLF) to generate an idler.

6. A method for increasing sensitivity for detection of ultrafast optical signals, the method comprising: inserting one of an amplitude modulator and a four-wave mixer configured to implement an optical weighting function upstream of a detector within a direct detection receiver, wherein the optical weighting function is configured to generate an optical pulse that closely matches an ultrashort pulse period of at least one ultrashort pulse within an input signal from a pulse source so that portions of the input signal within the ultrashort pulse period are multiplied by a non-zero function while portions of the input signal that are nonoverlapping in time with the pulse period are multiplied by zero.

7. The method of claim 6, wherein the optical weighting function operates according to the relationship f ( t ) = { s ( t ) , t 0 - τ s 2 + nT < t < t 0 + τ s 2 + nT 0 , otherwise , where τ.sub.s is a weighting function duration, t.sub.0 is a pulse center, T is the pulse period, and n is an integer.

8. The method of claim 6, wherein the optical weighting function is implemented using a four-wave mixer, wherein the input signal is filtered by an optical bandpass filter and combined with a weighting pulse obtained from the pulse source, the weighting pulse having a pulse width that matches the ultrashort pulse period, wherein the combined signal is guided into a highly non-linear fiber (HNLF) to generate an idler.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1A is a diagrammatic view of the operating principle of an optical weighting receiver according to embodiment of the invention compared against a conventional direct detection receiver; FIG. 1B is enlarged view comparing simulated detector outputs of the OWR and DDR of FIG. 1A.

(2) FIGS. 2A-2F illustrate estimated error probability (P.sub.e) for simulated receiver parameters. FIGS. 2A and 2B illustrate a comparison of the estimated P.sub.e for a direct detection receiver (DDR) (FIG. 2A) and an OWR according to the present invention (FIG. 2B). FIGS. 2C and 2D show a simulated comparison of the estimated error probability for narrow-band DDR and narrow-band OWR, respectively. FIG. 2E illustrates the estimated P.sub.e for an OWR with a narrow-band detector and a constant optical filter bandwidth of 300 GHz where the gating pulse is varied between 0.8% and 50%. FIG. 2F shows estimated P.sub.e for an OWR with a narrow-band detector and a constant optical filter bandwidth of 300 GHz where the gating pulse width is comparable to the signal pulse width.

(3) FIGS. 3A and 3B illustrate exemplary experimental setups for a DDR and OWR, respectively.

(4) FIGS. 4A and 4B are plots of P.sub.e measurement as a function of optical and electrical filter bandwidth for a DDR and OWR, respectively.

(5) FIG. 5A is a plot of the receiver operating characteristics (ROC) for a DDR and OWR at different signal powers; FIG. 5B is an enlarged view of the area within the circle of FIG. 5B.

(6) FIGS. 6A and 6B illustrate simulation setups of an OWR and DDR, respectively.

DETAILED DESCRIPTION OF EMBODIMENTS OF THE INVENTION

(7) Previous theoretical and experimental work relates receiver performance to three parameters: the optical filter bandwidth, the electrical filter bandwidth, and the signal duty cycle. Matching bandwidth parameters to the input pulse improves receiver sensitivity for a given optical and an electrical filter transfer function. However, in many applications involving picosecond and femtosecond pulses this task remains unattainable. The inventive approach of decoupling the receiver performance from the electrical bandwidth in the described setting can be numerically tested by contrasting direct detection and the OWR scheme and is described herein.

(8) FIG. 1A illustrates the fundamental difficulties in detection of an ultrashort signal 102 with a conventional direct detection receiver (DDR) 152 and the improvement obtained with the optical weighting receiver (OWR) 122 according to the present invention. In FIG. 1A, input pulse 102 (in frequency domain (ƒ) is shown in the time domain (t) as a short pulse 112 (E.sub.signal) within background noise 114 (E.sub.noise), which is incident onto direct detection receiver 152 with an impulse response h.sub.e(t) 126. After integration 128, the detection current 134, 136 is represented as the integral of the sum of the signal and the noise power over the detector impulse response time h.sub.e(t). In the ideal (matched-detection) case, the electrical weighting function 126 matches the signal and only the noise that coincides with the signal reaches the integrator 128 and the noise is weighted exactly by the signal profile. In practice, however, the conventional electrical weighting function h.sub.e(t) (direct-detection case, in the time domain) is much longer than the ultrashort signal, mandating accumulation of excess noise much beyond the signal time extent, and loss of detection fidelity.

(9) For the optical weighting receiver 122, the short pulse 112 hidden in noise 114 is optically gated prior to detection 124 (ƒ(t)). Optical weighting 124 rejects optical noise beyond the signal temporal duration, thus preempting the mixing of the signal with the portion of the noise that does not overlap in time with signal to yield pulse 104 in the frequency domain, pulse 106 in the time domain. In practical terms, the new receiver circumvents the electrical bandwidth limitation necessarily imposed on photon-electron conversion by rejecting the noise beforehand. The remaining operations of squaring and electrical weighting (impulse response h.sub.e(t)) 126 and integration 128 produce detection current 130, 132, with vastly improved discrimination. FIG. 1B provides an enlarged image of the resulting detection currents for the DDR (134, 136) and OWR (130, 132) to clearly show the enhanced signal within the noise.

(10) The detection process is represented as the convolution of the optical power with the detection system impulse response. The result of the photon-to-electron conversion in the direct detection process is the current I.sub.dd(t), typically represented as a product of the responsivity parameter, R, and a convolution of the optical power, |E(t)|.sup.2 with the detection system impulse response h(t):

(11) I dd ( t ) = R - .Math. E s ( τ ) + E n ( τ ) .Math. 2 h ( τ - t ) d τ = R - .Math. E s ( τ ) .Math. 2 h ( τ - t ) d τ + R - ( .Math. E n ( τ ) .Math. 2 + 2 Re ( E s ( τ ) E n * ( τ ) ) ) h ( τ - t ) d τ ( 1 )
and |E(t)|=|E.sub.s(t)+E.sub.n(t)|, where E.sub.s.n are the signal and noise fields respectively. The OWR detection process can be mathematically represented in the following manner:
I.sub.OWR(t)=R∫.sub.−∞.sup.∞|E.sub.s(τ)h.sub.g(τ)+E.sub.n(τ)h.sub.g(τ)|.sup.2h(τ−t)dτ=R∫.sub.−∞.sup.∞|E.sub.s(τ)|.sup.2h.sub.g(τ)h(τ−t)dτ+R∫.sub.−∞.sup.∞(|E.sub.n(τ)|.sup.2+2Re(E.sub.s(τ)E.sub.n*(τ)))h.sub.g(τ)h(τ−t)  (2)
where ƒ(t) is the optical weighting function, given as

(12) f ( t ) = { s ( t ) , t 0 - τ s 2 + nT < t < t 0 + τ s 2 + nT 0 , otherwise , ( 3 )
where τ.sub.s is the weighting function duration, t.sub.0 is the pulse center, T is the pulse period, and n is an integer. In the ideal (matched-detection) case, the electrical weighting function matches the signal and only the noise that coincides with the signal reaches the integrator and the noise is weighted exactly by the signal profile. In practice, however, the conventional electrical weighting function in the direct detection case, in the time domain, is much longer than the ultrashort signal, resulting in the accumulation of excess noise beyond the signal time duration, severely degrading the detected signal fidelity.

(13) The inventive approach deals with the noise addition problem by rejecting optical noise outside of the signal time span, before detection. In an exemplary embodiment, this method is implemented by adding optical weighting to reject optical noise beyond the signal temporal duration to avoid mixing of the signal with the portion of the noise that does not overlap in time with the signal. In practical terms, the OWR circumvents the electrical bandwidth limitation that is necessarily imposed on photon-electron conversion by rejecting the noise beforehand.

(14) The benefits of the OWR technique were first modeled and compared to direct detection in detection of an ultrashort signal represented by a Gaussian 1.5 ps long noisy signal repeating over a 100 ps period. Detected signal fidelity was characterized by calculating the probability of error (P.sub.e), using ta Bit Error Ratio (BER) procedure derived from well-established telecommunication characterization methods. Specifically, the widely-used system characterization approach described by Winzer and Kalmar (“Sensitivity enhancement of optical receivers by impulsive coding”, Journal of Lightwave Technology 17(2), 171-177 (1999)) was adopted to test the probability of error (P.sub.e) driven by two physical parameters: optical filter bandwidth and electrical filter bandwidth. The evaluation was performed using commercial simulation software from VPIphotonics Inc. (Norwood, Mass.). Numerical simulations were performed using the setups shown in FIGS. 6A and 6B, representing the OWR and DDR, respectively. The OWR receiver performance was tested by detecting the 1.5 ps signal 602 with 10 GHz repetition rate. Signal pulses with 0.14 mW peak power were amplitude modulated by amplitude modulator (AM) 604 with a pseudo-random bit sequence 606 and combined with white noise from optical noise source (ONS) 608 with power spectral density of 1.92×10.sup.−17 W/Hz. The resulting optical-signal-to-noise-ratio of 0.3 is estimated at operation 609 by measuring the average signal and the noise power independently. The combined fields are next passed through a Gaussian-shaped optical-band-pass-filter (OBPF) 610 and either detected directly at detector 614 or optically gated at amplitude modulator 612 with 1.2 ps pulses before detection. Pulse detection was performed by the square-law detector and both a dark current of 0.1 nA and shot noise were added from electrical noise source (ENS) 616 in the detection process. The total current is passed through a fourth-order Bessel-shaped low-pass filter (LPF) 618 with bandwidth B.sub.e and guided to the P.sub.e test engine 620. The DDR shown in FIG. 6B has a similar architecture to that of the OWR of FIG. 6A, with the main difference being the absence of the amplitude modulator 612 that introduces optical weighting before the detector 614. Both schemes were tested under equivalent conditions.

(15) The inventive approach to achieving optical weighting of the input signal is not limited to use of the amplitude modulator (AM) arrangement as in the example shown in FIG. 6B. In fact, the experimental set-up described with reference to FIG. 3B, described in more detail below, utilizes a four-wave mixing (FWM) process within a highly non-linear fiber (HNLF) to generate an idler capable of achieving the appropriate optical weighting. The FWM approach may actually be preferable to the AM method for speed purposes. Selection of other high speed optical processes for generating a matched weighting pulse will be apparent to those of skill in the art.

(16) Results of the simulations are shown in FIGS. 2A-2F. FIG. 2A represents the ideal, wide-band DDR with a minimum P.sub.e of 0.00012. For purposes of this description, this value is considered as optimal. The P.sub.e associated with the OWR scheme (FIG. 2B) also reaches the optimal P.sub.e level. However, in a more realistic setting, employing a narrow-band detector, the DDR performance deteriorates notably. This scenario becomes apparent in FIGS. 2C-2D, where the minimum P.sub.e achieved with the OWR method remains an order of magnitude lower than that of the DDR. Further, the OWR performance remains electrical-filter-bandwidth-independent, allowing the scheme to attain the optimal P.sub.e magnitude with the narrow-band receiver. Note that in the absence of optical weighting, no amount of electrical filtering, post square-law detection, can bring the performance below the level of 10.sup.−3.

(17) To evaluate the OWR's performance dependence on the gating pulse width, the electrical filter bandwidth was varied in the 10 GHz to 40 GHz range and the gating pulse width 0.8 to 32 times that of the input pulse. As expected, the OWR performance depends strongly on the gating pulse width, demonstrated by significant P.sub.e variations depicted in FIGS. 2E-2F. (Note that FIG. 2F is an expansion along the y-axis of the area within the black rectangle at the bottom of FIG. 2E.) Plotted P.sub.e maps point to an optimal optical weighting design—matching the weighting pulse to the signal directly optimizes the receiver performance.

EXAMPLE 1

Experimental Verification

(18) The experimental setup for the verification of the benefits of the inventive OWR approach, emulating the standard telecommunication pre-amplified receiver architecture, is diagrammatically illustrated in FIGS. 3A and 3B for a DDR and OWR, respectively.

(19) Probability of error (P.sub.e) measurement was performed using signal pulses 304 of 1.5 ps temporal width and 10 GHz repetition rate, centered at 1557.4, having sinusoidal amplitude and phase modulation, produced by pulse source 302 and compressed in a dispersive fiber 306. The pulse source 302 was divided to create the signal and the weighting pulse source 302′ (seen in FIG. 3B). The signal was amplitude modulated at amplitude modulator 308 with a 2.sup.7-1-long pseudo-random bit sequence 310 and then frequency converted to 1544 nm in the highly nonlinear fiber (HNLF) by means of the four-wave mixing process. The HNLF was characterized by the following global parameters: 182 m fiber length, 1555 nm average zero dispersion wavelength, 0.029 ps/(nm.sup.2km) dispersion slope and 28 W.sup.−1km.sup.−1 nonlinear parameter. The signal was next attenuated at attenuator 312 to 24.2 nW and then amplified in an Erbium-doped fiber amplifier (EDFA) 314. The signal average power is monitored at the EDFA input. The value reported in this measurement corresponds to the power before the amplifier front loss—this loss is 0.5 dB (0.5 dB insertion loss corresponds to the loss of two components before the Erbium fiber: the wavelength division multiplexer for combining the pump and signal, and the isolator. The splice loss is not taken into account. The bandwidth of optical bandpass filter (OBPF) 316 at the amplifier output is varied, having the following values: 0.6 nm, 1 nm, 2 nm, 5 nm and 7 nm.

(20) In case of the optical weighting receiver of FIG. 3B, the amplified signal is filtered at filter 316 after which a weighting pulse is generated in an optical fiber using a four-wave mixer (FWM). In this approach, the filtered signal is combined with the weighting pulse from source 302′ at combiner 330, which acts as the pump wave 350, then guided into a 20 m long HNLF 320. The generated idler, centered at 1571 nm, is then filtered using a 20 nm OBP filter 321 to remove the weighting pump and the signal. The idler 332 is then detected using a 40 GHz photodiode 318, followed by a 33 GHz electrical amplifier 322, a low-pass filter 324 and an acquisition system 326, which in the test set-up was a real-time oscilloscope. The electrical filter bandwidth is varied at the oscilloscope, having the following set of values: 11 GHz, 20 GHz, 25 GHz, 30 GHz, and 33 GHz. The sampling rate of the real-time oscilloscope was set to 100 GS/s and 10.sup.5 points are collected per acquisition, corresponding to 10.sup.4 bits. The reported P.sub.e value is averaged over six measurements.

(21) In the direct detection case, the amplified signal is detected directly by the same detection system as described above. In both the direct detection receiver scheme and the OWR scheme, a single polarization is detected at the amplifier output, by adding a polarization beam splitter after the amplifier. Both measurements are performed under the same conditions and the P.sub.e averaged over six measurements. The P.sub.e map is represented in FIGS. 4A and 4B for the DDR and OWR architectures, respectively.

(22) Measurements represented in FIGS. 4A and 4B confirm the significant improvement provided by the OWR approach. As expected, the OWR scheme outperforms the standard direct detection receiver, demonstrated by the seventeen-fold improvement in the P.sub.e measurement. Specifically, a minimum P.sub.e of 0.0032 is measured with the DDR scheme, while an order of magnitude lower minimum P.sub.e of 0.00018 is assessed with the OWR scheme under equivalent conditions. Still, some electrical filter bandwidth dependence is observed in the OWR measurement, a trend contributed to the sub-optimal gating pulse width and the finite gating pulse extinction ratio.

(23) Finally, the two detection schemes are contrasted by measuring the receiver operating characteristics (ROC) curve at different average signal power levels <P.sub.s>. FIGS. 5A and 5B represent the dark count probability, measured as a function of the detection efficiency by varying the detection threshold. Measurements were performed for a different average signal power at the amplifier input, ranging from 1.54 nW to 23.99 nW (corresponding to 1.2 to 18.6 average number of photons per bit). In the ROC plots, the DDR results are represented by the “°” marker while the OWR results are represented by “*” marker. The curve pairs for the DDR and OWR are numbered from right to left as indicated in Table 1 below:

(24) TABLE-US-00001 Curve Pair No. Signal Power <P.sub.S> 1 23.99 nW 2 12.19 nW 3  6.15 nW 4  1.54 nW
FIG. 5B is an enlargement of the area in the circle within FIG. 5A, showing the detailed results for the 1.54 nW signal power (curve pair no. 4).

(25) The dark count probability is measured as the probability of sending the logical “zero” state and detecting the “one” state, P(1|0); the detection efficiency is measured as the probability of sending the logical “one” state and detecting the logical “one” state, P(1|1). In this measurement, the sampling rate was set to 100 GS/s, as in the case of the P.sub.e measurement. The number of bits per acquisition was 10.sup.4 and was chosen for computational expediency. The ROC curves are assembled by taking multiple single-shot measurements at distinct points in time to increase the total number of points. The total number of points was greater than 10.sup.5.

(26) The average photon number per bit is estimated from the average power, measured at the amplifier input. The average photon number per bit <N.sub.bit> of the classically attenuated signal is calculated with the following equation:

(27) .Math. N bit .Math. = .Math. P .Math. T hf
where <P> is the average signal power, T is signal period, h is Planck's constant and ƒ is signal frequency. The average photon number per pulse is twice the average photon number per bit, where the factor of two originates from the PRBS modulation, having the same probability of 0 and 1.

(28) These measurements again clarify superiority of the inventive OWR approach, demonstrated by more than 5-fold detection efficiency increase for a fixed dark count probability of 0.001 and average signal power at the amplifier input of <P.sub.s>=1.54 nW. The OWR scheme shows unprecedented performance in direct detection of ultrashort signal obscured by noise. These results represent an important milestone in optimizing ultrashort signal detection.

(29) The OWR approach described herein has been shown to outperform a typical narrow-band receiver in ultrashort pulse detection. Specifically, experimental results demonstrate the viability of decoupling the receiver performance from its electrical bandwidth in the OWR architecture, allowing the scheme to reach the optimal P.sub.e value with narrow-band receiver. Experimental measurements are consistent with theoretical predictions; seventeen-fold improvement in the P.sub.e is achieved with the OWR scheme, compared to a direct detection scheme. Reconcilable results are observed in the ROC measurement, where more than 5-fold increase in the detection efficiency is evaluated for the average input signal power of 1.54 nW at the amplifier input.

(30) The inventive OWR approach shows immense potential in achieving optimal ultrafast pulse detection, with a plethora of applications across disciplines. Specifically, besides typical telecommunication receiver architecture demonstrated in this work, many applications in medicine and biology could largely benefit from the OWR receiver. Recent breakthroughs in histopathological tissue imaging provide one example of such an application that is limited by its reliance on existing direct detection receiver architecture for ultrafast pulse detection. Optical coherence tomography-based imaging is another example of an application in which a coherent detection scheme is implemented for ultrashort pulse detection. While coherent detection technique is inherently more sensitive than the direct detection method, it is also a subject to strict arm balancing requirements that are often difficult to meet. The OWR scheme, thus, has potential to improve the detection sensitivity of the misbalanced coherent receiver. Such improvement could significantly increase the acquisition rate and possibly the depth range in optical coherence tomography.

(31) In summary, receiver sensitivity improvement is achieved by employing optical weighting prior to detection of the picosecond-long signal obscured by noise. Findings indicate the possibility of decoupling receiver performance from its electrical bandwidth, ultimately optimizing the probability of error. Such sensitivity improvement is beyond the reach of the standard, narrow-band, direct detection receiver technologies. Although the preceding description involves sensitive detection of picosecond-long pulses, the technique applies to pulses of arbitrary duration. Indeed, the technological advancement represented by the OWR approach provides a path to a new class of opto-electronic detector technology with a potential to bring profoundly transformative effects across many disciplines.

(32) As will be apparent to persons of skill in the art, modifications and variations to the specific detector types and architectures described in the examples may be made without deviating from the inventive approach of introducing an optical weighting element and operation prior to detection of ultrafast signals. Accordingly, the scope of the invention is not intended to be limited to the specific embodiments described herein, but is intended to encompass all optical detector arrangements that employ optical weighting to preempt mixing of noise within the input signal with the portion of the signal corresponding to the ultrashort pulse.

REFERENCES

(33) [1] Khilo, A., et al., “Photonic ADC: overcoming the bottleneck of electronic jitter”, Optics Express 20(4), 4454-4469 (2012). [2] Koos, C., et al., “All-optical high-speed signal processing with silicon-organic hybrid slot waveguides”, Nature Photonics 3, 216-219 (2009). [3] Kuo, B. P.-P., et al., “Transmission of 640-Gb/s RZ-OOK channel over 100-km SSMF by wavelength-transparent conjugation”. Journal of Lightwave Technology 29(4), 516-523 (2011). [4] Sibbett W., et al., “The development and application of femtosecond laser systems,” Optics Express 20(7), 6989-7001 (2012). [5] Konig, K, et al., Clinical two-photon microendoscopy,” Microscopy Research and Technique 70, 398-402 (2007). [6] Barretto, R. P. J., et al., “Time-lapse imaging of disease progression in deep brain using fluorescence microendoscopy,” Nature Medicine 17, 223-228 (2011). [7] Mejan, G. et al., “Remote detection and identification of biological aerosols using femtosecond terawatt lidar system,” Applied Physics B 78, 535-537 (2004). [8] Udem, T. H., et al., “Optical frequency metrology,” Nature 416, 233-237 (2002). [9] McDonoigh, R. N.; and Whalen, A. D. Detection of signals in noise (Academic Press, San Diego, 1995). [10] Kay, S. Fundamentals of statistical signal processing, Volume II: Detection theory (Prentice Hall, Upper Saddle River, 1998). [11] Winzer, P. J.; and Kalmar, A., “Sensitivity enhancement of optical receivers by impulsive coding,” Journal of Lightwave Technology 17(2), 171-177 (1999). [12] Chen, Y., et al., “Reducing ASE effect in coherent detection by employing double-pass fiber preamplifier and time-domain filter,” IEEE Journal of Quantum Electronics 45(10), 1289-1296 (2009). [13] Sanjoh, H., et al., “Time-domain filtering to reduce signal-spontaneous beat noise for a pulse-type multiwavelength light source,” IEEE Photonics Technology Letters 15(5), 757-759 (2003). [14] Derickson, D.; and Muller. M. Digital Communications Test and Measurement (Prentice Hall, Boston, 2008). [15] Slavik, R., et al., “All-optical phase and amplitude regenerator for next-generation telecommunications systems.” Nature Photonics 4, 690-695 (2010). [16] Otsuji, T., et al., “10-80-Gb/s highly extinctive electrooptic pulse pattern generation,” IEEE Journal of Selected Topics in Quantum Electronics 2(3), 643-649 (1996). [17] Andrekson, P. A.; and Westlund, M. Nonlinear optical fiber based high resolution all-optical waveform sampling. Laser & Photonics Reviews 1(3), 231-248 (2007). [18] Agrawal, G. P. Nonlinear Fiber Optics (Academic Press, San Diego, 2007). [19] Osche, G. R. Optical detection theory for laser applications (Wiley, Hoboken, 2002). [20] Tu, H, et al., “Stain-free histopathology by programmable supercontinuum pulses,” Nature Photonics 10, 534-540 (2016). [21] Huang, D., et al., Optical coherence tomography. Science 254, 1178-1181 (1991).