Photonica
← Articles
Procedure · Laser characterization

Slope Efficiency and Differential Quantum Efficiency from LIV Measurements

Procedure for extracting slope efficiency and external differential quantum efficiency from laser diode LIV data, with discussion of measurement geometry, facet-collection considerations, and common extraction errors.

Published May 8, 20268 min read

Scope

This article describes the procedure for extracting external slope efficiency ηs\eta_s and external differential quantum efficiency ηd\eta_d of a semiconductor laser diode from light–current–voltage (LIV) measurements. Coverage is restricted to single-mode and multi-mode edge-emitting devices in continuous-wave or low-duty-cycle pulsed operation. Internal quantum efficiency separation by inverse-length method is treated separately.

Definitions

Slope efficiency ηs\eta_s is defined as the local derivative of optical output power PoutP_\text{out} with respect to drive current II in the lasing region of the LIV curve:

ηs  =  dPoutdII>Ith,\eta_s \;=\; \left. \frac{d P_\text{out}}{d I} \right|_{I > I_\text{th}},

with units of W/A or mW/mA. It is determined by linear regression on the LIV curve above threshold.

External differential quantum efficiency ηd\eta_d is the dimensionless ratio of photons emitted from the device facet per electron injected above threshold:

ηd  =  (PoutPout(Ith))/hν(IIth)/q  =  qhνηs  =  qλhcηs,\eta_d \;=\; \frac{(P_\text{out} - P_\text{out}(I_\text{th}))/h\nu}{(I - I_\text{th})/q} \;=\; \frac{q}{h\nu} \, \eta_s \;=\; \frac{q \lambda}{h c} \, \eta_s,

where qq is the electron charge, hνh\nu is the photon energy, λ\lambda is the emission wavelength, hh is Planck's constant, and cc is the speed of light in vacuum.

For practical conversion at common telecom and laser wavelengths, the factor qλ/(hc)=λ/(1.2398 μmV)q\lambda/(hc) = \lambda/(1.2398~\mu\text{m} \cdot \text{V}) gives:

Wavelengthqλ/(hc)q\lambda/(hc)ηd\eta_d for ηs=0.5\eta_s = 0.5 W/A
850 nm0.686 A/W34.3%
1310 nm1.057 A/W52.8%
1550 nm1.250 A/W62.5%
1064 nm0.858 A/W42.9%

Reported values of ηd\eta_d above unity are physical for devices where both facets are collected and the round-trip differential efficiency exceeds the front-facet collection. Single-facet values are bounded by unity.

Distinction from other efficiencies

Slope efficiency must not be confused with:

  • Wall-plug efficiency ηWPE=Pout/(IV)\eta_\text{WPE} = P_\text{out} / (I V), which is the total electrical-to-optical conversion at a specific operating point. ηWPE\eta_\text{WPE} depends on both the slope and the absolute operating current, voltage, and threshold; it is a different parameter and cannot be derived from ηs\eta_s alone.
  • Internal quantum efficiency ηi\eta_i, the fraction of injected carriers that produce stimulated photons inside the cavity. ηi\eta_i is extracted from inverse-length measurements (1/ηd1/\eta_d vs. LL) over a set of devices with varying cavity lengths and is not obtainable from a single-device LIV.
  • Power conversion efficiency at threshold or below, which is not defined in the strict sense — slope efficiency describes only the lasing region.

Required measurements

The measurement configuration is the same as for any LIV characterization:

FunctionComponentSpecification
Current sourceSMU or pulsed current driver0.1%\leq 0.1\% resolution; rise time <100<100 ns for pulsed
Voltage measurement4-wire Kelvin sense100\leq 100 μV resolution
Optical detectionCalibrated power meter or integrating sphereLinear over the lasing power range; spectral calibration at the test wavelength
Temperature controlTEC + thermistor±0.1\pm 0.1 K stability

Two geometric configurations are common and produce different ηs\eta_s values that must not be compared without correction.

Single-facet butt collection. A photodiode or fiber is positioned close to one facet of the laser. Captures only one direction of emission and is sensitive to far-field divergence. For Fabry–Pérot devices with symmetric cleaved facets, the front-facet measurement captures approximately half of the total emitted power. For HR/AR-coated devices, the front-facet measurement captures the asymmetric majority of the emitted power.

Integrating sphere collection. The laser is placed at the entrance of an integrating sphere. Captures all emitted power independent of beam divergence. Required for devices with elliptical or non-Gaussian emission, for accurate ηd\eta_d extraction, and for any device with HR/AR coatings where front-facet collection is geometry-dependent.

The collection geometry must be reported alongside the ηs\eta_s value. Direct comparison of single-facet ηs\eta_s across devices is valid only when the facet geometry is identical.

Extraction procedure

1. Acquire LIV with adequate dynamic range

The LIV sweep must include current values extending well above IthI_\text{th} — typically I2IthI \geq 2 I_\text{th} to 3Ith3 I_\text{th}, and ideally up to the rated operating current — so that the lasing region is well-sampled. Current step size in the lasing region should be small enough to give at least 30 data points in the linear fitting range; for a typical telecom DFB with Ith=8I_\text{th} = 8 mA and operation to 80 mA, a step of 0.5–1 mA is sufficient.

For high-power devices, pulsed measurement at 0.1%\leq 0.1\% duty cycle eliminates self-heating-induced curvature in the LIV that would otherwise depress the apparent slope.

2. Identify the lasing region

The lasing region is the range where d2P/dI20d^2 P / d I^2 \approx 0. Three behaviors define its endpoints:

Lower bound. Immediately above threshold the LIV curvature is nonzero due to the gradual transition from spontaneous to stimulated emission. The lower bound of the linear region is conservatively taken at I1.2IthI \geq 1.2 \, I_\text{th}.

Upper bound. At high currents, three mechanisms reduce slope:

  • Thermal rollover. Self-heating reduces gain. Visible as smooth concave-down curvature.
  • Carrier leakage. Carriers escape over the heterojunction barriers. Material- and structure-dependent.
  • Spatial hole burning and kinks. Mode redistribution produces step changes in slope.

The upper bound is taken at the current where the LIV begins to deviate from linear by more than a stated tolerance (commonly 1% of expected linear output, or where the second derivative magnitude exceeds 1% of the slope).

3. Fit

Perform a least-squares linear fit on (I,Pout)(I, P_\text{out}) over the identified lasing region:

Pout(I)  =  ηs(IIth),P_\text{out}(I) \;=\; \eta_s \, (I - I_\text{th}^*),

where IthI_\text{th}^* is the linear-extrapolated threshold from this fit (not necessarily the same as the threshold extracted from the kink or second-derivative method). The slope of the fit is ηs\eta_s.

Report:

  • ηs\eta_s in W/A or mW/mA
  • The fitting range [Imin,Imax][I_\text{min}, I_\text{max}]
  • The number of data points NN
  • The fit residual standard deviation σr\sigma_r
  • R2R^2 of the fit

A high-quality LIV in the lasing region should yield R2>0.999R^2 > 0.999 and residuals randomly distributed with no trend.

4. Convert to differential quantum efficiency

For a device with center wavelength λ\lambda:

ηd  =  ηsqλhc.\eta_d \;=\; \eta_s \cdot \frac{q \lambda}{h c}.

If the emission spectrum changes significantly over the fitting range (e.g., mode-hopping), use the peak wavelength averaged over the same current range, or extract ηd\eta_d at fixed wavelength using a wavelength-resolved measurement.

Worked example

Using the same 1310 nm Fabry–Pérot InP device from the T₀ extraction article, the LIV at 25 °C gives:

Current (mA)Power (mW)
8.70.0 (threshold)
12.02.45
14.03.50
16.04.55
18.05.60
20.06.65

Linear least-squares fit over I[12,20]I \in [12, 20] mA yields:

ηs  =  0.525 W/A,Ith  =  7.33 mA,R2=0.9999.\eta_s \;=\; 0.525 \text{ W/A}, \qquad I_\text{th}^* \;=\; 7.33 \text{ mA}, \qquad R^2 = 0.9999.

The extrapolated threshold Ith=7.33I_\text{th}^* = 7.33 mA is slightly below the two-segment kink threshold of Ith=8.7I_\text{th} = 8.7 mA, reflecting the gradual curvature in the near-threshold region. Converting to differential quantum efficiency:

ηd  =  0.525 W/A×1.057 A/W  =  0.555  =  55.5%.\eta_d \;=\; 0.525 \text{ W/A} \times 1.057 \text{ A/W} \;=\; 0.555 \;=\; 55.5\%.

This is the single-facet external differential quantum efficiency. For a symmetric uncoated Fabry–Pérot device, the total ηd\eta_d collected from both facets would be approximately twice this, or 100%\sim 100\% — at the upper limit of the physical range and consistent with a low-loss device.

Sources of extraction error

The following dominate the spread of ηs\eta_s values reported in the literature for nominally identical devices.

Inconsistent fitting range. ηs\eta_s extracted on I[1.05Ith,1.5Ith]I \in [1.05 I_\text{th}, 1.5 I_\text{th}] differs from ηs\eta_s extracted on I[1.5Ith,3Ith]I \in [1.5 I_\text{th}, 3 I_\text{th}] by 5–15% for typical devices, due to the gradual transition region near threshold and to thermal rollover at high II. The fitting range must be reported with the value.

Self-heating in CW measurements. For high-power devices in CW mode, junction temperature rise during the sweep reduces gain at the end of the sweep relative to the start. The apparent ηs\eta_s is biased low. Magnitude depends on device thermal impedance, mounting, and sweep direction; for typical high-power devices, the bias is 5–20% relative.

Photodetector saturation or non-linearity. Power meters and photodiodes have linear ranges. Exceeding the linear range produces apparent slope reduction not related to the device. Verify the detector linear range covers the LIV power range, including margin for the high-current end.

Spectral mismatch in detector calibration. Power meter responsivity varies with wavelength. A meter calibrated at 1550 nm used at 1310 nm without recalibration can introduce 5–20% systematic error. Detectors must be calibrated at the actual measurement wavelength.

Far-field truncation. Single-facet butt-coupled photodiodes capture only the central portion of the laser far-field. For lasers with broad far fields or non-Gaussian emission, the captured fraction varies with operating conditions, biasing ηs\eta_s low and producing artificial current-dependence in the slope. Use an integrating sphere or a beam-expanded lens-coupled detector to eliminate.

Cleave or facet quality variation. ηd\eta_d depends strongly on facet reflectivity. Devices with poor or uneven cleaves show reduced and variable ηd\eta_d that does not reflect bulk device properties. Visual inspection of facets under high magnification before measurement identifies most facet issues.

Wavelength drift across fit range. For devices that mode-hop or have appreciable wavelength shift over the fitting range, the conversion factor qλ/(hc)q\lambda/(hc) changes slightly. For typical edge-emitters with α=0.30.5\alpha = 0.3{-}0.5 nm/^\circ C and modest self-heating, the wavelength-dependence error is 1%\sim 1\% and usually neglected.

Validation

Two checks confirm extraction quality.

The fit residuals should be randomly distributed about zero with no systematic trend versus current. Visible curvature in the residual plot indicates either an inappropriate fitting range or a real device non-linearity.

The extracted ηd\eta_d should fall within the literature range for the material system, wavelength, and device structure. Typical bounds:

Systemηd\eta_d range (single facet, uncoated FP)
InGaAsP/InP 1550 nm20–35%
InGaAsP/InP 1310 nm25–45%
InGaAs/GaAs 980 nm pump35–60%
AlGaAs/GaAs 808 nm30–50%
GaN/InGaN 405 nm20–40%

Values outside these ranges by more than 30% indicate either a high-performance specialty device or an extraction error.

References

For the textbook treatment of slope efficiency, differential quantum efficiency, and the inverse-length method for separating internal efficiency from mirror loss, see Coldren, Corzine, and Mašanović (2012), chapter 2. For the influence of self-heating on extracted slope and methods for thermal correction, see Piprek (2003), chapter 7. For ISO measurement methodology including detector calibration and far-field considerations, see IEC 61290-1.