Threshold Current Extraction Methods Compared

Scope

This article compares the four standard methods for extracting threshold current $I_{\text{th}}$ from a semiconductor laser diode LIV curve: the two-segment linear fit, the second-derivative peak, the first-derivative-of-log inflection point, and the fixed-power threshold. All four methods are applied to a single representative LIV dataset to demonstrate the numerical differences they produce. The conditions under which each method is preferred are discussed.

Why the choice matters

A semiconductor laser does not have a single, mathematically unambiguous threshold current. It has a transition region in which spontaneous emission gives way to stimulated emission, and any extraction method places the threshold somewhere within this region; different methods place it at different points. Thresholds quoted for similar devices are therefore comparable only when they were extracted the same way.

The implications:

• $I_{\text{th}}$ values from different methods are not interchangeable
• Comparisons across publications require matched methodology
• Extracted parameters that depend on $I_{\text{th}}$ (notably $T_0$, slope efficiency, wall-plug efficiency at fixed power above threshold) inherit the methodological choice
• Method selection should be reported alongside the value

The four methods reviewed here are the ones in common use. How far apart they fall depends on the width of the transition. On a clean curve from an ordinary edge-emitting laser, the three methods that locate the transition agree to within a small fraction of a milliampere, and the fixed-power threshold is offset from them by an amount set by its reference power; on a broad transition all four spread out. The worked example below shows both cases.

The transition region

Above threshold, the optical output power $P$ versus drive current $I$ is well-approximated as linear:

$$P(I)={\eta }_s(I-I_{\text{th}}),I\gg I_{\text{th}},$$

with slope efficiency ${\eta }_s$ and an extrapolated intercept $I_{\text{th}}$. Below threshold, $P$ is approximately linear in $I$ as well, with a much smaller slope set by spontaneous emission outcoupling.

In the transition region near $I_{\text{th}}$, $P(I)$ smoothly curves between the two linear regimes. The width of this transition is determined by the spontaneous emission factor $\beta$, the fraction of spontaneously emitted photons that couple into the lasing mode. In the rate-equation model given in that entry, the peak of $d^{2}P\mathrm{/}d{I}^2$ has a full width at half maximum of about $3\sqrt{\beta }{I}_{\text{th}}$. For a Fabry–Pérot semiconductor laser, with $\beta \sim {10}^{-5}$ to ${10}^{-4}$, that is 1 to 3% of the threshold current: 0.25 mA for $\beta ={10}^{-4}$ and a threshold of 8 mA. For high-$\beta$ devices such as microcavity and nanocavity lasers, the transition spreads over a range comparable to the threshold current itself.

The four extraction methods locate the threshold at different points within this transition.

Method 1: Two-segment linear fit

Two independent linear fits are performed: one to the sub-threshold region (gentle slope, near-zero output), and one to the above-threshold region (steep slope, linear output). The intersection of the two fit lines is taken as $I_{\text{th}}$.

Procedure:

1. Identify the lasing region, typically $I>1.2I_{\text{th}}^{\text{est}}$, where $I_{\text{th}}^{\text{est}}$ is a rough visual estimate. Fit a linear function $P_{\text{above}}(I)={\eta }_s(I-I_a)$ to this region.
2. Identify the sub-threshold region, typically $I<0.7I_{\text{th}}^{\text{est}}$. Fit a linear function $P_{\text{below}}(I)=aI+b$.
3. Solve for the intersection: $I_{\text{th}}=({\eta }_s{I}_a+b)\mathrm{/}({\eta }_s-a)$. In the limit $a\to 0$, this simplifies to $I_{\text{th}}=I_a+b\mathrm{/}{\eta }_s$.

For typical devices where the sub-threshold slope is much smaller than the above-threshold slope ($a\ll {\eta }_s$), the two-segment intersection is well-approximated as $I_{\text{th}}\approx I_a$, the linear extrapolation of the lasing region back to zero output power.

This method is the most common in practice and is the implicit assumption in many published laser characterization tools.

Advantages. Simple to compute, robust to measurement noise, intuitive interpretation, requires only modest LIV resolution near threshold.

Disadvantages. Sensitive to the choice of fitting ranges; substantial dependence on where the lasing region is taken to begin. For devices with curvature extending well above the kink (high-$\beta$ lasers, thermally compromised devices), the extracted $I_{\text{th}}$ depends strongly on the lower bound of the fitting range.

Method 2: Second-derivative peak

The optical output power is differentiated twice with respect to current. The threshold is taken at the current where $d^{2}P\mathrm{/}d{I}^2$ reaches its maximum value.

Mathematically this corresponds to the inflection point of the slope itself, the point of maximum curvature in the transition. For an ideal step-function transition between the two linear regimes, $d^{2}P\mathrm{/}d{I}^2$ would be a delta function at the kink; for a real device with finite transition width, it is a smooth peak.

Procedure:

1. Smooth the LIV data, typically with a Savitzky–Golay filter or a moving-average window.
2. Numerically differentiate twice with respect to $I$ (or fit a spline and take the analytical second derivative).
3. Identify the peak of $d^{2}P\mathrm{/}d{I}^2$ and locate its position $I_{\text{th}}$.

Advantages. No fitting range to choose: the method operates directly on the data. Mathematically well-defined and unambiguous. Particularly useful for devices with significant sub-threshold curvature where the two-segment method's fitting range becomes critical.

Disadvantages. Highly sensitive to measurement noise: taking two numerical derivatives amplifies noise by orders of magnitude. Requires current steps near threshold that are small compared with the width of the transition, which for an ordinary edge emitter is 1 to 3% of $I_{\text{th}}$; with coarser steps the peak can only land on a grid point unless it is refined by fitting a parabola through the highest point and its two neighbors. The smoothing filter parameters influence the extracted threshold value.

Method 3: First-derivative-of-log inflection

The natural logarithm of optical output power is differentiated once with respect to current. The threshold is taken at the inflection point of $\ln P$, where $d(\ln P)\mathrm{/}dI$ peaks, equivalently where the second derivative $d^2(\ln P)\mathrm{/}d{I}^2$ changes sign.

This method exploits the fact that $\ln P$ versus $I$ has a characteristic S-shape near threshold, with the inflection point identifying the steepest portion of the log-power transition. The method is particularly useful for low-power LED-like operation below threshold transitioning to lasing operation above.

Procedure:

1. Compute $\ln P_i$ for each measured point.
2. Smooth and numerically differentiate to obtain $d(\ln P)\mathrm{/}dI$.
3. Locate the maximum of $d(\ln P)\mathrm{/}dI$ within the transition region, which corresponds to the inflection point of $\ln P$. At the lowest currents, where $P$ rises from zero, $d(\ln P)\mathrm{/}dI$ is also large (it approaches $1\mathrm{/}I$ where $P\propto I$), so the search is confined to the region around the visible knee.

Advantages. Robust for measurements where the sub-threshold output is non-negligible. Less sensitive to noise than the second-derivative method (one fewer differentiation). Well-defined for high-$\beta$ devices where the transition is broad.

Disadvantages. Requires LIV data with positive output power throughout the measurement range. Measurements that drop to noise floor or detector dark in the sub-threshold region cannot be log-transformed. Sensitive to background offset; baseline subtraction is required before log transformation.

Method 4: Fixed-power threshold

The threshold is taken at the current that produces a fixed reference output power, typically 0.1 mW or 1 mW depending on device class.

Procedure:

1. Choose a reference power $P_{\text{ref}}$ appropriate to the device class.
2. Find the current $I_{\text{th}}=I(P=P_{\text{ref}})$, by interpolation if necessary.

Advantages. Trivial to compute. Reproducible: does not depend on fitting choices or differentiation parameters. Useful in production-test environments where consistency across thousands of devices is more important than absolute accuracy.

Disadvantages. Not a real threshold by any physical definition; the extracted $I_{\text{th}}$ depends on the choice of reference power. Cannot be compared across publications unless $P_{\text{ref}}$ is matched. Strongly device-dependent: a 0.1 mW reference is meaningful for a low-power telecom DFB but produces poorly-defined results for a 1 W pump diode.

Worked example

The example uses a computed curve rather than a measured one, so that the threshold is known exactly and each method can be checked against it. The curve comes from the rate-equation model in the spontaneous emission factor entry, for a laser with a threshold of 8.0 mA, a slope efficiency of 0.55 W/A and $\beta ={10}^{-4}$: the 25 °C values of the 1310 nm Fabry–Pérot laser in the T₀ extraction article. A background of spontaneous emission into the other modes, 0.2 µW per mA, is added; it follows the carrier density, so it stops rising at threshold. The model threshold, where the gain reaches the loss, is 8.0 mA. Drive current is swept from 0 to 30 mA in 0.5 mA steps.

Applying each method:

Two-segment linear fit. Linear fit on $I\in [12,26]$ mA: $P=0.550(I-7.995)$ in mW, with $R^2>0.9999$. Sub-threshold fit on $I\in [0,5.5]$ mA, below 0.7 $I_{\text{th}}$: $P=0.00036I-0.00012$ in mW. Intersection: $I_{\text{th}}=8.00$ mA.

Second-derivative peak. Differentiation with a 5-point Savitzky–Golay filter puts the peak of $d^{2}P\mathrm{/}d{I}^2$ at the 8.0 mA grid point. The transition is only 0.25 mA wide, narrower than the step, so the grid alone locates it to within ±0.25 mA; a parabola through the highest point and its two neighbors gives 8.00 mA.

First-derivative-of-log inflection. Within the transition region, $d(\ln P)\mathrm{/}dI$ peaks at the 8.0 mA grid point; the parabolic refinement gives 8.02 mA.

Fixed-power threshold at $P_{\text{ref}}=0.1$ mW. Linear interpolation between the 8.0 and 8.5 mA points gives 8.11 mA. The model curve itself reaches 0.1 mW at 8.14 mA; the difference is the error of interpolating linearly across the knee.

Summary:

On this sharp transition the three methods that locate the transition agree with the model threshold to within 0.02 mA. The fixed-power threshold lies 0.14 mA above it, because 0.1 mW is reached only after lasing has begun. The derivative methods depend on the grid: starting the same 0.5 mA sweep at a different offset moves their grid-point peaks anywhere from 7.75 to 8.25 mA, while the parabolic refinement stays within 0.03 mA of 8.00.

A broad transition. The same calculation with $\beta ={10}^{-2}$, as in a microcavity laser, and 0.05 mA steps spreads the methods out:

With $\beta ={10}^{-2}$ the emission into the lasing mode below threshold is a hundred times stronger, and the curve bends over a wide range of current. The lasing fit on 12–26 mA extrapolates to 7.78 mA, and the two-segment intersection comes back up to 7.98 mA only because the sub-threshold line is now steep enough to move it. The fixed-power threshold falls below the model threshold altogether, since 0.1 mW is reached at 71% of it. The model threshold is 8.0 mA in both cases.

For the sharp transition, the lasing-line intercept, the two-segment value and the second-derivative peak all coincide with the model threshold, and the slope efficiency extraction article obtains the same intercept from the same curve. Measured curves can have rounder knees than the model's, from effects it leaves out, such as current spreading, several longitudinal modes reaching threshold at slightly different currents, and self-heating; the rounder the knee, the further apart the methods fall.

Choosing among methods

The four methods locate the threshold at slightly different physical features of the transition. For a given application:

For new measurements, the two-segment linear fit is the most defensible default. It is the implicit assumption in much of the published literature on slope efficiency and characteristic temperature, and the extracted $I_{\text{th}}$ value enters cleanly into the corresponding formulas. The second-derivative method is the appropriate companion measurement when the visible kink location is the property of interest.

The fixed-power method should be reported as a fixed-power threshold and not labeled simply as "$I_{\text{th}}$". The value depends on $P_{\text{ref}}$, and the convention is publication-dependent.

Sources of error common to all methods

Independent of the extraction method, several measurement issues affect all four methods.

Self-heating in CW measurements. At currents approaching $I_{\text{th}}$, self-heating elevates the active region temperature, increasing $I_{\text{th}}$. For CW measurements, the threshold is therefore somewhat higher than the value at the heatsink temperature. Heating at higher currents also bends the lasing line, which moves a fitted intercept in the opposite direction, as the pulsed and CW comparison shows. Pulsed measurement at $\le 0.1\mathrm{\%}$ duty cycle largely removes both effects.

Optical detector noise floor. Sub-threshold output power is often below the detector noise floor, which puts $P=0$ in the sub-threshold region of the LIV data. The two-segment fit handles this correctly (the sub-threshold fit slope approaches zero). The first-derivative-of-log method requires the noise floor to be subtracted before log transformation.

Background or stray light reaching the detector. Stray light produces a power-independent baseline offset in the LIV. The two-segment fit absorbs it in the sub-threshold intercept, and the second-derivative peak is unaffected because differentiation removes a constant; the first-derivative-of-log and fixed-power methods need the baseline subtracted first.

Insufficient current resolution near threshold. Methods 2 and 3 need current steps near threshold that are a fraction of the width of the transition, which for an ordinary edge emitter is 1 to 3% of $I_{\text{th}}$; steps of 0.5% of $I_{\text{th}}$ or less resolve it. With coarser steps, such as the 0.5 mA (6%) steps of the worked example, the peak can only land on a grid point unless it is refined with a parabola.

Validation

For all methods, extracted $I_{\text{th}}$ should lie within the apparent transition region of the LIV, between where the sub-threshold and above-threshold linear approximations clearly hold. A threshold value outside the visible transition indicates extraction failure.

The three methods that locate the transition should agree to within about the width of the transition, a few tenths of a milliampere for an ordinary edge emitter. The fixed-power threshold is compared only with other fixed-power thresholds at the same $P_{\text{ref}}$. Larger disagreement indicates a broad transition, measurement noise, or an artifact in the dataset.

References

For the rate-equation description of the transition region and the spontaneous emission factor $\beta$, see L. A. Coldren, S. W. Corzine and M. L. Mašanović, Diode Lasers and Photonic Integrated Circuits, 2nd ed. (Wiley, 2012), chapter 2. For the high-$\beta$ regime, where the transition is broad and the threshold needs a more careful definition, see G. Björk, A. Karlsson and Y. Yamamoto, "Definition of a laser threshold," Physical Review A 50, 1675 (1994).

The T₀ extraction calculator fits the temperature dependence of thresholds extracted this way.