## 1.

## INTRODUCTION

Semiconductor lasers (SLs) are fundamental components of fiber based optical communication channels, due to their compact sizes, high reliability and low power operating conditions. Directly modulated (DM) SLs are the preferred choice in short range optical networks, e. g. Local Area Network (LAN) and Metropolitan Area Networks (MANs), due to the reduced system cost and complexity with respect to externally modulated laser sources [1]. Additionally, short and intense optical pulses with duration of tens of picoseconds can be generated with SLs when the driving conditions are properly set in order to work in gain switching [2] regime. Gain Switched (GS) SLs are practical pulse sources that find applications in a wide variety of fields, e. g. telemetry, optical ranging, 3D scanning, material processing and biomedical analysis.

A clear understanding of the basic dynamics processes in SLs is then relevant for the student audience of university photonics courses. In this work we present a web hosted educational tool for the simulation of DM SLs: the device model is based on a simple formulation of the well know Rate Equations (REs) of carrier and photon density and optical phase. The laser internal parameter and current driving conditions can be set by the user and the carrier density and laser output power are plotted as a function of time. Instantaneous frequency variations, typical of DM SLs, are numerically shifted to the audible frequency range and the resulting signal is sent to the computer loudspeakers, giving an intuitive description of the “chirp” phenomenon in SLs [3].

This work is organized as follows. In Section 2, the proposed simulation tool is described, giving details of the mathematical model in use, its numerical implementation and the user control interface. In Section 3, some examples of use are presented, in which the program is run in different conditions, changing the laser modulation parameters. Finally, in Section 4, the main conclusions of the proposed work are summarized.

## 2.

## SOFTWARE DESCRIPTION

## 2.1

### The Rate Equations Model

The light-matter interactions occurring in a SL can be modeled in a simple way by using the RE model [1, 4]. The RE model considers the interplay between carriers and photons in the laser cavity in terms of conservation equations. The number of electrons and holes is assumed to be the same, thus only one equation for carriers is considered in the model. The main assumption of the model is that carriers and photons are considered to be uniformly distributed in the longitudinal and transversal directions.

Assuming a single mode laser, the relationship between the carrier density, N(t), the photon density, S(t), and the optical phase, ϕ (t), can be expressed by the following equations:

where the descriptions of the symbols and their values are given in Tab. 1. Typical values for a 1550 nm Distributed Feed Back (DFB) laser [5] are used by default.

## Table 1.

Symbols and default values used in Eqs. 1-3.

1550 nm DFB Laser [5] | |||
---|---|---|---|

Symbol | Description | Value | Units |

q | Elementary charge | 1.602·10-19 | C |

Vact | Active volume | 2·10-11 | cm3 |

Γ | Confinement factor | 0.25 | - |

dG/dN | Differential gain | 1.7·10-16 | cm2 |

vG | Group velocity | 9.4·109 | cm/s |

N0 | Transparency carrier density | 1.6·1018 | cm-3 |

τN | Carrier lifetime | 1 | ns |

τP | Photon lifetime | 1 | ps |

ε | Nonlinear gain compression factor | 1.7·10-17 | cm3 |

β | Spontaneous emission coefficient | 1·10-4 | - |

α | Linewidth enhancement factor | 3.5 | - |

The RE model previously described has been implemented in Java® [6] programming language and the resulting applet is hosted by Universidad Politécnica de Madrid, in a web server at Tecnología Fotónica y Bioingeniería (www.tfo.upm.es/educational-tools). Equations 1-3 are solved with a Runge Kutta iterative algorithm for N(t), S(t) and ϕ(t), with temporal resolution of 1 ps for a vector of 10000 points (10 ns).

## 2.2

### The Graphic User Interface

The Graphic User Interface (GUI) is shown in Fig. 1. It is composed of four panels and three graphs. From the *Modulation Parameters* panel the user can choose the current waveform, the bias point, the peak-to-peak amplitude and the frequency of the modulating current. The waveform of the current can be selected by the user, either sinusoidal or rectangular.

The laser parameters can be introduced by the user in the *Laser Parameters* panel, placed on the right side of the GUI. Values shown in Tab. 1 are considered by default. Carrier density threshold (N_{TH}) and threshold current (I_{TH}) are calculated as N_{TH} = N_{0} + (τ_{P}·dG/dN·v_{G}Γ)^{-1} and I_{TH} = N_{TH}·q·V_{act}·τ_{N}^{-1}, by solving Eqs. 1-2 in steady state regime at threshold.

By pressing the *Run* button in the *Modulation Parameter* panel, the program solves Eqs. 1-3, and returns the values of carrier density, photon density and optical phase at each temporal instant. The instantaneous frequency, *v*(t), is calculated from the temporal derivative of the phase, as *v*(t) = (2π)^{-1}d/dt(ϕ(t)). Three graphs on the left side of the GUI show the results of the simulation. From up to down, the temporal profiles of injected current, carrier density and output power are shown. Current and carrier density are plotted as a function of I_{TH} and N_{TH}, respectively.

From the laser intensity and chirp, the software builds an amplitude and frequency modulated signal of the form A(t)·cos(2π(f_{0} + f(t)), where A(t) is proportional to S(t), f(t) is the frequency *v*(t) shifted into the audible frequency range. This conversion is obtained with a constant factor of 10^{9}, changing nanoseconds into seconds and GHz into Hz. The Sound control is placed at the bottom of the right side of the GUI and allows switching on and off the reproduced sound. The Save Data panel allows the user to export the simulation results to a text file, where the time, current, number of carriers, output power and instantaneous frequency vectors are stored in a tab separated format.

## 3.

## EXAMPLES

In this Section, we present some simulation results that can be useful for understanding the dynamics of a DM SL. We show the dependence of the laser output, in terms of intensity and instantaneous frequency, on its internal parameters and on the modulation conditions, i. e. current waveform, current amplitude and bias point. In the proposed examples, attention is given both to the graphical representation of the results and on the impact of the audible output on the user understanding of the simulated phenomenon.

These examples are directed to a student audience with previous basic notions on the theory of laser dynamics. Here we explain how to use the proposed simulation tool and its potential for a better learning experience, giving a brief theoretical description of the simulated phenomena. Results are exported with the Save Data button and presented in figures.

## 3.1

### Laser Switch On

We first consider the simple case in which the SL is switched on with a rectangular current function. In order to have this, in the *Modulation Parameters* panel, the rectangular current waveform is selected, I_{BIAS} is set to 1.5 I_{TH}, the amplitude is set to 3 I_{TH} and the frequency is set to 0.25 GHz (4 ns period with 50 % duty cycle).

This corresponds to inject the laser with a current step function with low level equal to 0 and a high level current of 3 I_{TH}. The laser parameters are set to the default values.

After the Run button is pressed, the program takes some seconds to complete the simulation and the results are finally shown in the GUI graphs. Current, carrier density, output power and instantaneous frequency are exported with the *Save Data* button and results are plotted in Fig. 2 (a) and (b).

When the rectangular current is applied, the carrier density grows as 1- exp(-t/τ_{N}) from the low level steady state to N_{TH}, as expected from theory [1, 4]. N(t) undergoes damped oscillations around the threshold density value before reaching the steady state value clamped at N_{TH}. This corresponds to the relaxation oscillations of the output power, due to carrier induced gain modulation, and of the instantaneous frequency, due to the carrier induced index modulation.

The user can appreciate in an intuitive manner, the coupled amplitude/frequency modulation, by pressing the Sound in the GUI. As expected, when the laser is off, for I(t) = 0, no sound is produced, then a chirped sound is emitted, when power and instantaneous frequency oscillate, and finally a single frequency tone is produced, corresponding to the optical carrier of the laser shifted into the audible band.

## 3.2

### Direct Modulation

In this example, we consider the direct modulation of the laser with a periodic rectangular current waveform, with the low and high current level both above the threshold current value I_{TH}. In the Modulation Parameters panel, the rect current waveform is selected, BIAS is set to 4 I_{TH}, the amplitude is set to 2 I_{TH} and the frequency field is set to 0.5 GHz. After pressing the Run button, the results are plotted in the GUI graphs. In this case, the well known phenomenon of transient and adiabatic chirp in DM SLs is clearly appreciated after listening to the audible output. The transient chirp is usually referred to the instantaneous frequency variation associated to the transition between the two levels of the injected current, i. e. the damped oscillations of *v*(t) which follows from low to high and from high to low current level transitions. The adiabatic chirp is the frequency difference between the frequency produced during the high level and the frequency produced during the low level of injected current. The adiabatic chirp is directly related to the gain compression factor ε. As shown in Fig. 3 (a) and (b), the carrier density is not clamped to the threshold value NTH as in the ideal case (ε = 0).

In the proposed simulation, the audio output is composite of a single higher tone during the high level state, a lower single tone during the low level state (corresponding to the adiabatic chirp) and a chirped sound when the current switches from low to high level and viceversa (corresponding to the transient chirp). This is also appreciated from the instantaneous frequency temporal profile, where the frequency difference between low and high level of the injected current due to the adiabatic chirp, is clearly shown in Fig. 4 (a) and (b). If the gain compression factor ε is set to zero, the adiabatic chirp signature disappears from the audio output and the graphs, as expected from theory.

## 4.

## CONCLUSIONS

We have presented novel educational software for the simulation of the dynamics of directly modulated semiconductor lasers. The proposed tool is directed to undergraduate students of photonics and optical communications courses.

The software is based on the well known RE model for SLs, and it is numerically implemented on a Java® applet. The GUI allows the student to choose the laser internal parameters and the modulating conditions. The simulation results are plotted and the amplitude/frequency modulation typical of current injected SLs is shifted to the audible frequency range for an intuitive understanding of the phenomenon. Several examples are given, i. e. the laser turn on and laser direct modulation, and results are discussed for different modulation parameters.

The proposed software allows the user to explore and test his knowledge on SL dynamics with an interactive tool by changing the operating condition of the simulated device and see and listen to the consequent results in the laser output.

## ACKNOWLEDGEMENTS

This work has been supported by Universidad Politécnica de Madrid, through the project “Desarrollo de técnicas de aprendizaje y metodología activas para la impartición de la Fotonica”. We gratefully acknowledge the support of the Ministerio de Economia y Competitividad (Spain) through project TEC2012-38864-C03-02.

## REFERENCES

G. P. Agrawal, [Fiber-Optic Communication Systems], John Wiley & Sons, Inc., 2002 https://doi.org/10.1002/0471221147Google Scholar

K. Y. Lau, “Gain switching of semiconductor injection lasers” Appl. Phys. Lett. 52 (4), 257–259 (1988) https://doi.org/10.1063/1.99486Google Scholar

A. Consoli, B. Bonilla, P. R. Horche and I. Esquivias “Sound of lasers (SOL) - An audiovisual approach to semiconductor laser dynamics”, Proc. of the 15th International Conference on Interactive Collaborative Learning (ICL), 26 - 28 September 2012, Villach, Austria. https://doi.org/10.1109/ICL.2012.6402037Google Scholar

L. A. Coldren and S. W. Corzine, [Diode Lasers and Photonic Integrated Circuits], John Wiley & Sons, Inc., 1995.Google Scholar

John C. Cartledge and R. C. Srinivasan, “Extraction of DFB Laser Rate Equation Parameters for System Simulation Purposes”, J. Lightwave Technol., vol. 15, iss. 5, pp. 852–860, 1997. https://doi.org/10.1109/50.580827Google Scholar