# Improved Nonlinear Model Implementation for VCSEL Behavioral Modeling in Radio-Over-Fiber Links

Andrea Giovannini, Muhammad Usman Hadi, Lucas Iogna Prat, Najett Neji, Zerihun Gedeb Tegegne, Carlos Viana, Anne-Laure Billabert, Jean-Marc Laheurte, Jacopo Nanni, Diego Masotti, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Andrea Giovannini, Muhammad Usman Hadi, Lucas Iogna Prat, Najett Neji, Zerihun Gedeb Tegegne, et al. Improved Nonlinear Model Implementation for VCSEL Behavioral Modeling in Radio-Over-Fiber Links. Journal of Lightwave Technology, 2022, 40 (20), pp.6778-6784. 10.1109/JLT.2022.3195048. hal-03850780

HAL Id: hal-03850780

https://hal.science/hal-03850780

Submitted on 8 Dec 2022

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
# Improved Nonlinear Model Implementation for VCSEL Behavioral Modeling in Radio-Over-Fiber Links

Andrea Giovannini<sup>D</sup>, Muhammad Usman Hadi<sup>D</sup>, Lucas Iogna Prat, Najett Neji, Zerihun Gedeb Tegegne, Carlos Viana<sup>D</sup>, Anne-Laure Billabert, Jean-Marc Laheurte<sup>D</sup>, Jacopo Nanni<sup>D</sup>, Diego Masotti<sup>D</sup>, Giovanni Tartarini<sup>D</sup>, and Jean-Luc Polleux<sup>D</sup>

**Abstract**—This paper proposes an improved model for the nonlinear behavior of the components of Radio-over-Fiber (RoF) links focusing here on the laser source. In a previous work, it was shown with a flexible new opto-microwave design approach a first behavioral model of RoF links that exploited the hyperbolic tangent function (the so-called *tanh* model) to model the high-power nonlinearities. Since this was not sufficient to simulate adequately all the nonlinearities, in this work, for the first time in the authors' knowledge, the improved nonlinear model based on Cann's model is developed, on an electro optic transducer. Then the amplitude-to-amplitude conversion of the link is better taken into account with the introducing of nonlinear properties of the link at high input power in the model. The new transfer equation includes several parameters that can be adjusted to adapt the model to any of the components, especially the photo-emitters under study, which are the main cause of this type of nonlinear behavior of the system. In particular, this improved model is applied on the state-of-the-art 850 nm GaAs Vertical Cavity Surface Emitting Laser (VCSEL), simulating several Figures of Merit such as the Error-Vector Magnitude (EVM) of a Long Term Evolution (LTE) signal with a Quadrature Phase Shift Keying (QPSK) and the output amplitude-modulated power versus input amplitude-modulated power, also called AM/AM curve. This new proposed nonlinear model outperforms other existing behavioral models and its non-linear simulation better fits with the measurements. Finally, we

conclude that the performance improvement highly depends on the nonlinearity sharpness parameter, depending on the frequency and the bias current.

**Index Terms**—Behavioral model, improved cann's nonlinear model, microwave-photonics, radio-over-fiber (RoF), vertical cavity surface emitting laser (VCSEL).

## I. INTRODUCTION

A s a consequence of the increasing demand of high-speed and broadband wireless access, a trend of developing Radio-over-Fiber (RoF) technologies has appeared as a solution to reduce the infrastructure costs [1]–[5]. Indeed, at very high frequencies, the coverage distance is limited up to several meters due to the substantial signal attenuation in these frequency bands. Hence, radio-over-fiber (RoF) technologies are a good alternative to improve the coverage of Home Area Networks (HAN) [6] and can be privileged in a larger scale of applications as in stadium or tunnel.

To fulfill application requirements, complex modulation schemes are transmitted through an optical source into the fiber and re-emitted after optical detection into the final radio-cell. It is therefore important to preserve the quality of the radio signal through the RoF link and to properly design it in terms of frequency, noise and non-linearity behavior.

Behavioral modeling is in this regard an important approach that helps to provide a simple analysis of the component behavior, laser or photodetector for instance, into a system, or vice-versa to determine what are the minimum requirements for the device to operate into a given system. They opposed themselves to compact-circuit models or rate equations models [7] which requires a precise knowledge of the device physics. As part of behavioral models, we distinguished those models based on pure analytical equations such as Volterra series or memory polynomial such as those applied to lasers in [8] and [9] and those based on analytical equations to represent the system metrics of the device. We focus our interest on the latter behavioral models category that refers to the figures of merit of the device that provide a more flexible approach.

Many approaches to develop figures of merit taking into account the optical devices characteristics have been proposed in [1]–[22]. In [6], a complete set of Opto-Microwave (OM)

Manuscript received 31 March 2022; revised 6 July 2022; accepted 22 July 2022. Date of publication 29 July 2022; date of current version 21 October 2022. This work was supported by the French FUI8 ORIGIN project and FUI MORF project and by the Université Franco Italienne - Università Italo Francese, through the sponsorship of the Joint Doctorate Project C3-2025 Efficient Wireless Coverage of In-Building Environments with Low Electromagnetic Impact within the Call Vinci 2019. (*Corresponding author: Andrea Giovannini*.)

Andrea Giovannini, Muhammad Usman Hadi, Jacopo Nanni, Diego Masotti, and Giovanni Tartarini are with the Dipartimento di Ingegneria dell'Energia Elettrica e dell'Informazione "Guglielmo Marconi", Università di Bologna, 40136 Bologna, Italy (e-mail: andrea.giovannini14@unibo.it; usmanhadi@ieee.org; jacopo.nanni3@unibo.it; diego.masotti@unibo.it; giovanni.tartarini@unibo.it).

Lucas Iogna Prat, Zerihun Gedeb Tegegne, Jean-Marc Laheurte, and Jean-Luc Polleux are with the Univ Gustave Eiffel, CNRS, ESYCOM, F-77454 Marne-la-Vallée, France (e-mail: lucas.iognaprat@edu.esiee.fr; zerihungedeb.tegegne@esiee.fr; jean-marc.laheurte@univ-eiffel.fr; jean-luc.polleux@esiee.fr).

Najett Neji is with the Université Paris Saclay, Univ Evry, 91025 Gif-sur-Yvette, France (e-mail: najett.neji@univ-evry.fr).

Carlos Viana is with the ICON Photonics, SAS, 77420 Champs-sur-Marne, France (e-mail: carlos.viana@icon-photonics.com).

Anne-Laure Billabert is with the Le Cnam, ESYCOM-CNRS, 75003 Paris, France (e-mail: anne-laure.billabert@lecnam.net).

Color versions of one or more figures in this article are available at https://doi.org/10.1109/JLT.2022.3195048.

Digital Object Identifier 10.1109/JLT.2022.3195048
quantities has been identified, which connects optical devices with microwave tools and concepts. These figures of merit, namely Opto-Microwave power, gain, noise figure, compression point (OM-P₁dB) and third-harmonic intermodulation product (OM-IP₃) have been defined, and behavioral models were proposed for the laser, the fiber and the photodetector. These models were implemented in a standard microwave systems simulator named ADS software from Keysight Technologies and efficiently represented the noise figure and the non-linearity of the link. However, some deviations to the measured non-linearity of the laser were observed [6].

In addition, several models have been developed in the last decades to describe the RoF links behavior [23], [24]. Due to the similarities between the large signal behavior of both lasers (e.g., VCSEL LIV curves [25], [26]) and power amplifiers (e.g., amplifier transfer function), some models have been developed, based on the existing literature about behavioral modeling of power amplifiers [27]–[35]. In [6], the large signal behavior is modeled by using the Hyperbolic Tangent (tanh) to fit the simulations with the measurements; an approach that was already proposed in the early 1990s. However, the main disadvantage of this method is that this function does not allow the user to adapt the model to any kind of device by changing the knee sharpness of the amplitude-to-amplitude (AM/AM) curve. In 1980 [36], Cann proposed an instantaneous nonlinearity model (in the context of radio-frequency power amplifiers) with adjustable knee sharpness; however, this model led to non-consistent results for two-tone intermodulation test, because of its nonanalyticity at zero. In 2012 [37], Cann proposed an improved nonlinearity model in order to overcome the pointed issues. Hence, in this paper, we adapt the Cann non-linear models of the RF power amplifiers [36] to photonic subsystems (such as lasers, but also photodetectors, and photonic link) as an improved solution for the behavioral model of RoF links.

The paper is organized as follows. In Section II, the RoF link system as well as the experimental setup for its characterization is described. In Section III, the behavioral modeling is developed. First a description of the hyperbolic tangent model is proposed. The Cann’s non-linearity behavior model is then proposed and adapted to RoF non-linear model. In Section IV, we present the complete opto-microwave model and we analyze its comparison to measurements including EVM for LTE signals transmission. Finally, in Section V, we summarize the main contributions and further work.

## II. ROF LINK DESCRIPTION AND EXPERIMENTAL SETUP

The concept of Radio-over-Fiber directly modulates the laser with an RF signal, transmits the modulated optical signal through the fiber, and detects its RF envelope through a photodetector. As described in [6], this scheme is suitable for home area networks (HAN) and enables the use of low-cost opto-electronic components such as Vertical Cavity Emitting Surface Lasers (VCSEL) together with multimode fiber, even up to millimeter-wave frequencies when using Intermediate-Frequency (IF) scheme. This paper however refers to a Direct-RoF scheme where the RF signal is within the bandwidth of the VCSEL and directly

<figure><img src="image_3.png" /><figcaption>Fig. 1. RoF link architecture.</figcaption></figure>

modulates the laser. This constitutes the core of both Direct and IF RoF systems.

The typical architecture of the RoF link composed of a laser, an optical fiber and a photodetector is shown in Fig. 1.

For this study, the link is composed by a VCSEL laser, a single mode optical fiber connected to a variable optical attenuator and an InGaAs-based high speed photodiode (1414-50 New Focus). The laser is a 10 GHz GaAs VCSEL die, under probe, operating at 850 nm from Philips Technologies ULM Photonics GmbH. The purpose of the variable optical attenuator is to keep the fiber’s attenuation at a constant level, in order to easily compare the measurements at different biasing conditions. In particular, this total attenuation is set to be 5 dB.

The laser’s optical aperture diameter is 8 µm and the technology is based on mesa etching and subsequently laterally oxidized current apertures [32]–[35]. The die is biased through a 100 µm pitch RF GSG probe. The laser output beam is coupled to a lensed multimode fiber probe whose position is precisely fixed above the die with a nanopositioner to optimize the coupling efficiency. The photodetector is used as a reference photodetector for the extraction of the VCSEL performance. It is highly linear and shows a flat frequency response up to 25 GHz, with a responsivity of 0.22 A/W at 850 nm, under short-circuit conditions.

The optical fiber presents a defined angle of 6° with respect to the orthogonal direction of the VCSEL die in order to avoid parasitic optical cavities into the laser, due to backscattered light. The fiber output extremity is connected to the photoreceiver. The output signal is then connected to an Electrical Spectrum Analyzer (ESA).

The characteristics of the VCSEL laser and the New Focus photodiode used for the experimental measurements are summarized in Table I, where IP₁dB is the input 1 dB compression point, RIN is the Relative Intensity Noise of the laser, and fₐ₃dB is the cut-off frequency.

## III. BEHAVIORAL MODELLING OF ROF LINK

### A. Tanh Model

In our previous work [6], the input-output characteristics of both laser and photodetector were modelled by using the hyperbolic tangent as the transfer function to fit the simulations with the measurements. In that case, the definition of the equivalent
<table><caption>TABLE I
LD AND PD SPECIFICATIONS FOR OUR TEST BENCH</caption><thead><tr><th>Parameter</th><th>VCSEL LD</th><th>New Focus PD</th><th>Units</th></tr></thead><tbody><tr><td>IP<sub>1dB</sub></td><td>6</td><td>15</td><td>dBm</td></tr><tr><td>Gain</td><td>-27</td><td>-19</td><td>dB</td></tr><tr><td>Wavelength</td><td>850</td><td>850</td><td>nm</td></tr><tr><td>RIN</td><td>-147.5</td><td>-</td><td>dB/Hz</td></tr><tr><td>f<sub>.3dB</sub></td><td>12.3</td><td>25</td><td>GHz</td></tr></tbody></table>

optical current was given in (1):

$$I_{out} = I_{S} \cdot \tanh\left(\frac{G \cdot I_{in}}{I_{S}}\right) \quad (1)$$

Where $I_{out}$ is the component's output current, i.e., in our case the equivalent optical current at the laser's output $I_{opt,LD}$ as defined later in (2), $I_{S}$ is the limit value of the output signal which provides the saturation effect, G is the opto-microwave gain in the linear region, and $I_{in}$ is the injected input current to the component, i.e., the electrical injection current of the laser. In addition, $I_{opt,LD}$ is given by (2):

$$I_{opt,LD} = \alpha_{c} P_{opt,LD} \quad (2)$$

Where $\alpha_c$ is the conversion coefficient of the optical power to electrical current (considered equal to 1 A/W in our study). The opto-microwave (OM) gain of the laser, denoted $G_{LD}^{OM}$, is expressed as the ratio of the laser's equivalent OM output power to the electrical power at its input. When considering a 50 Ohm laser, the laser OM gain is expressed by (3):

$$G_{LD}^{OM} = \frac{P_{opt,LD}}{P_{in,LD}} = S_{LD}^{2} \quad (3)$$

Where $S_{LD}$ is the laser's slope efficiency. Similarly, the OM gain of the photodiode, denoted $G_{PD}^{OM}$, is expressed as the ratio of the output electrical power provided by the photodiode to the equivalent OM power delivered to its input. When considering a matched photodiode, it is given by:

$$G_{PD}^{OM} = \frac{P_{out,PD}}{P_{opt,PD}} = R_{PD}^{2} \quad (4)$$

Where $R_{PD}$ is the photodiode's responsivity. This model allows to control the output saturation level and the slope of the curve by changing the parameters $I_S$ and $G$ respectively. However, it is not possible to control the sharpness of the knee. For this reason, it was imperative to find a new model capable of simulating this behavior.

## B. First Cann's Nonlinear Model and Its Limit

In [36], Cann proposed an instantaneous nonlinear model for power amplifiers to analyze their non-linearity focusing on a variable named “knee sharpness”. This model was based on the frequency response of a RC high pass filter, that better implements the desired shape of the equivalent output optical current versus input current.

Equation (5) shows the Cann's first non-linear model, where y represents the equivalent output optical current, L is the output limit level, x is the input current and s is the knee sharpness variable.

$$y(x) = \frac{L}{\left(1 + \left(\frac{L}{x}\right)^s\right)^{1/s}} \quad (5)$$

This model is efficient for many reasons. First, it has a linear behavior for small input signals, and it tends to an asymptotic defined limit level for large input signals. In addition, it is an economical model to compute (closed form and analytic), capable of modeling various knee sharpness values by changing a single parameter. However, this model gives incorrect results for the third-order intermodulation product (IP<sub>3</sub>), since it has a nonphysical behavior for the large input signal regime. Because of the nonphysical behavior, the Cann model gives imprecise predictions for nonlinear products [38].

## C. Improved Cann's Nonlinear Model

In 2012, Cann developed an improved model where the limitations of its first nonlinear model were solved [37]. The solution for the equivalent optical current in this model is given by (6), where the variables in this equation can be related to our RoF system. Hence, y represents the output equivalent optical current of the laser, x is the input current of the laser, g is a gain parameter well approximated by the slope efficiency of the laser for sufficiently high k, L is the output limit level (saturation level) and k is the sharpness parameter.

$$y(x) = \frac{L}{k} \ln \left( \frac{1 + e^{k \left( \frac{g x}{L} + 1 \right)}}{1 + e^{k \left( \frac{g x}{L} - 1 \right)}} \right) - L \quad (6)$$

Cann demonstrated that this model suitably behaves for two-tone intermodulation products [37], and this outperforms his first model. While the other parameters of the model (i.e., gain and saturation level) can be easily measured, the sharpness parameter k is apparently determinable only with heuristics. A deeper analysis of the model allows to relate the parameter k to the input IP<sub>3</sub> (IIP<sub>3</sub>) and the input P<sub>1dB</sub> (IP<sub>1dB</sub>). In fact, as mentioned before, the gain parameter g in (6) is formerly not exactly equal to the slope efficiency of the laser, but is linked to it through (7):

$$g = \frac{1 + \cosh(k)}{\sinh(k)} S_{LD} \quad (7)$$

This relation is obtained by considering that the current linear gain (i.e., $\sqrt{G_{LD}^{OM}}$) must be equal to the first derivative of (6) computed in 0. The multiplying factor depends only on k and tends to 1 for increasing values of k. If k > 3 then g > 0.9 $S_{LD}$, so for k sufficiently high g = $S_{LD}$ is a good approximation. If we derivate up to the 3rd order and consider the IP<sub>3</sub> condition, then the saturation level L is linked to the IIP<sub>3</sub> through

$$L = kg\sqrt{\frac{3}{2} \cdot \frac{IIP_3}{1 + \cosh(k)}} \quad (8)$$

This last relation can be used to determine k and it requires to measure the saturation current of the laser. An alternative
<figure><img src="image_4.png" /><figcaption>Fig. 2. Laser (green box) internal model.</figcaption></figure>

is to consider the P<sub>1dB</sub> condition in order to obtain a relation dependent only on *k*. Assuming that in the P<sub>1dB</sub> biasing point the power associated to non-linear harmonics is negligible, then the following relation is found:

$$ \left[ \ln \left( \frac{1 + e^{\frac{4}{3} \frac{IP_{1dB}}{IIP_3} (1+\cosh(k)) + k}}{1 + e^{\frac{4}{3} \frac{IP_{1dB}}{IIP_3} (1+\cosh(k)) - k}} \right) - k \right] = \sqrt{\frac{4}{3(1 - \cosh(k))} \frac{IP_{1dB}}{IIP_3}} 10^{-1/20} \quad (9) $$

Equation (9) finally depends only on *k* and on the ratio between *IP*<sub>1dB</sub> and *IIP*<sub>3</sub>. Isolating *k* from (9) is rather challenging, but through numerical interpolation it is possible to find an appropriate polynomial approximating the curve. The Cann's model has been validated on commercial electrical amplifier and soft clipping limiters. Owing to this result, the choice has been done to exploit this Cann model to model an electro-optic transducer with a behavior approach. The steps and results are presented in the next paragraph.

### D. Electro-Optic Modelling Based on Cann's Model

The Cann's improved model is then introduced into the RoF model developed in [6] using the microwave circuit and system Keysight Technologies simulation software. The laser, the optical fiber and the photodetector were implemented through the use of Symbolically Defined Devices (SDD), which was introduced in its respective internal blocks. In each SDD cell, we define the transfer function of our model (only for the laser and photodetector since the optical fiber is considered as a linear system).

The laser model is implemented following the equivalent scheme in Fig. 2. A current generator with an internal 50 Ω resistance R<sub>G,in</sub> provides a current I<sub>G,in</sub> to the circuit. The noise sources of the laser are implemented into the laser model through a noise voltage source V<sub>N</sub> in μV/√Hz describing the RIN. The RIN function is implemented using the conventional frequency response given in [39] and using equation (10).

$$ RIN(f) = A \cdot \frac{4\pi^2 \cdot f^2 + B}{16\pi^4(f_R^2 - f^2)^2 + 4\pi^2 \cdot \gamma^2 \cdot f^2} \quad (10) $$

Where A and B are device-dependent constants, *f* is the modulation frequency, *f*<sub>R</sub> is the relaxation oscillation frequency and *γ* the damping factor of the intrinsic laser.

I<sub>in</sub> is the laser injection current while R<sub>in</sub> is the input resistance of the laser, which can be conveniently mismatched with respect to the generator if a reflection needs to be taken into account.

However, in this work we assume a matched input port. On the output side, a generator with an internal 50 Ω resistance R<sub>G,out</sub> provides a current I<sub>G,out</sub> based on (6) and implementing also a frequency dependent equation and non-linear analytical equation according to [23]. The output of the device I<sub>out</sub> represents the output optical power, P<sub>opt,LD</sub>, as an equivalent electrical current noted I<sub>opt,LD</sub> feeding a matched 50 Ω load R<sub>L</sub> (in this case it corresponds to the input impedance of the optical fiber OM block). The laser overall modulation frequency response H(*f*) is implemented through a common third order system, given by equation (11):

$$ H(f) = \frac{1}{1 + j \cdot \frac{f}{f_P}} \cdot \frac{f_R^2}{f_R^2 - f^2 + j \cdot f \cdot \frac{\gamma}{2\pi}} \quad (11) $$

Where *f*<sub>P</sub> is the electrical parasitic cutoff frequency due to the access lines. The large signal behavior is then implemented into the model by introducing the modified Cann model equation (6), where *y* is the output current I<sub>opt,LD</sub> of the laser, *x* is the injected input current into the laser active part, and *g* is the small-signal responsivity divided by 1A/W to suit the conversion P<sub>opt,LD</sub> to I<sub>opt,LD</sub>. The parameters *L*, *k* and *g* are derivated considering the P<sub>1dB</sub> and IP<sub>3</sub> as found in (7), (8) and (9). The optical fiber and photodetector model's definitions are described in [6]. The optical fiber model includes the coupling losses between the laser and the fiber. Its attenuation value is slightly adjusted to fit with the measurements, as misalignments may vary. Considering the short length of the fiber, the frequency dependent term of the attenuation due to modal dispersion has not been considered. This device is considered linear and noiseless. The P<sub>1dB</sub> value of the photodiode model is set to a high value, because the nonlinearities of the device are negligible in our case study. To be noted that the model formerly transduces an electrical current into an optical power, represented by the equivalent OM current.

## IV. STUDY OF THE ROF LINK: SIMULATIONS WITH THE IMPROVED CANN'S MODEL COMPARED TO ROF MEASUREMENTS

### A. System Simulation: Gain, AM/AM and EVM Curves

By using the previously described model, we investigate several figures of merit (OM gain, AM/AM curve, and Error Vector Magnitude (EVM) performance) and compare simulation results with measurements. The data of the Gain and the AM/AM curves is relative to the input electrical power of the laser to the output electrical power at the photodiode, with a maximum RF input power of 15 dBm. For this, we perform a measurement of a RoF link that is constituted by the VCSEL laser, multimode optical fiber (MMF) and New Focus photodiode directly modulated at a carrier frequency of 2 GHz. By knowing the responsivity of the photodiode and the attenuation of the fiber plus the optical attenuator, it is possible to extract the laser's output OM power, useful for the behavioral model. Regarding the EVM measurements, a LTE signal is transmitted through the RoF link with a Quadrature Phase Shift Keying (QPSK) modulation at a center frequency of 2 GHz and a bandwidth of 20 MHz. A Vector Signal Generator (VSG) Agilent model N5182B MXG X-Series
<figure><img src="image_5.png" /><figcaption>Fig. 3. Measured and simulated opto-microwave gain (at 2 GHz using the VCSEL – 8 mA biasing).</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>Fig. 4. Measured and simulated AM/AM curves (at 2GHz using the VCSEL – 8mA DC biasing). The simulation considers first and improved Cann’s model as well as tanh model.</figcaption></figure>

is exploited to generate the signal emulated through a domestic MATLAB code. The VSG-generated signal directly modulates the VCSEL laser. The photo-detected signal is amplified and captured with a Vector Signal Analyser from Rhode&Schwartz. The captured signal is then processed offline in MATLAB and demodulated, followed by EVM estimation.

Fig. 3 presents the comparison of the simulated OM gain (solid lines) with the measured OM gain of the proposed RoF link at 2 GHz where the VCSEL is biased at 8 mA. It can be observed that the simulated curves are fitting with the measurements data.

Based on the nonlinear model described above, we also perform the simulations at 2 GHz based on first and improved Cann’s nonlinear model as well as tanh model as presented in Fig. 4. It can be noted that the new models are properly fitting with the measurement and behave much better than the tanh model. Indeed, it is impossible to change the shape of the simulated tanh model. Nevertheless, using the improved Cann’s model, we can control the shape only by changing the value of the parameter *k* (see equation (6)).

Now we are interested in the Error Vector Magnitude (EVM) which is an important indicator of the RoF link performance. The measured EVM is compared with respect to the three models

<figure><img src="image_7.png" /><figcaption>Fig. 5. Measured and simulated EVM curves for the first test bench (at 2 GHz using the VCSEL – 8 mA biasing). (a) Measured EVM curve comparing with each simulated EVM curve. (b) Measured and simulated EVM curves altogether.</figcaption></figure>

(tanh, Cann’s and improved Cann’s model) as shown in Fig. 5(a). The Back-to-Back (B2B) EVM experimental results are also indicated in both figures. In Fig. 5(b) we can notice that the curve corresponding to the improved Cann model fits better to the measurement data compared to the tanh model especially at high input RF power. Similarly, the improved Cann model fits closer to the experimental data as compared to Cann model itself.

The parameters of the tanh and Cann’s models are determined either by measurements or by numerical estimation following a Least Absolute Error (LAE) method, which aims at minimizing the deviation between the model and the AM/AM measurements. For the estimation of *k* in the Improved Cann’s model (6), we numerically solved (9) and found the values of the ratio *IP*₁dB/*IIP*₃ for *k* ranging between 0.1 and 30, then fit the points with a 4th order polynomial. By measuring *IP*₁dB and *IIP*₃ we then extract *k* and determine *L* and *g* with (7) and (8).

## B. Sharpness Parameter Effect on the Model

Finally, we analyze the influence of the sharpness parameter on the link performance when the improved Cann’s nonlinear model is used. To fully test the model, we perform AM/AM measurement for full frequency band and VCSEL bias current sweep, respectively from 1 to 6 GHz and from 4 to 10 mA. Then the Sharpness parameter is extracted through a mathematical fit of the experimental output power vs input power curves (AM/AM) over the range of the sweep parameters (frequency and source current). The result for this Least Absolute Error (LAE) extraction provides the extracted *k* under the improved Cann’s model as shown in Fig. 6.
<figure><img src="image_8.png" /><figcaption>Fig. 6. Sharpness parameter as a function of the bias Current and the Frequency for improved Cann’s model.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>Fig. 7. The IP<sub>1dB</sub> compression point (in dBm) of the VCSEL as a function of biasing current function of frequency.</figcaption></figure>

At low VCSEL bias current, the sharpness parameter *k* values are lower. In this case, the operation point is very close to the threshold of the VCSEL, and this might be the cause of the nonlinearities. For higher biasing, the *k* parameter increases with respect to the frequency, which means that the curve becomes sharper. Hence, the transition between the linear and the saturation regions is faster.

The variation of the sharpness parameter can be explained by the change in the intrinsic dynamic behavior of the VCSEL as a function of the current bias. Fig. 7 depicts the experimental result of P<sub>1dB</sub> compression point of the VCSEL as a function of bias current and frequency. As the bias current of the VCSEL increases the relaxation oscillation frequency is also increased. This influences the nonlinearities of the VCSEL such as the P<sub>1dB</sub> compression point also.

## V. CONCLUSION

This paper has proposed a behavioral model for VCSEL lasers employed into RoF systems with a focus on the modeling of their non-linearity. It exploits a Cann and Improved-Cann model, which were known previously in the field of RF power amplifiers and transposed in this paper to the opto-microwave domain. The resulting model includes the frequency response, RIN modeling and non-linear behavior.

Comparisons between the numerical simulations and measurements of the AM/AM curves over the range 1-10GHz and EVM data in the case of an LTE signal at 2GHz were provided

for a VCSEL RoF link. The non-linearity sharpness parameter *k* is key to represent the non-linearity of the device. Its extraction over the full 1-6GHz has been provided showing that the parameter is highly dependent on the frequency and biasing current.

Behavioral modeling is in this regard an important approach that helps to provide a simple analysis of the component behavior into a system, or vice-versa to determine what are the minimum requirements in terms of compression point, non-linearity sharpness, bandwidth or gain for the device to operate into a given system.

While previous *tanh* model did not evaluate properly the region of EVM non-linearity, the proposed model provides a convenient tool that helps to predict the maximum *k* value for a device to operate with a given maximum EVM for example. It helps the selection of a given VCSEL of the prediction of the performance of a given laser into an overall system.

To improve the Cann’s nonlinear model, and as further work, we are focusing on providing analytical expressions of some OM figures of merit. We are more particularly interested on the 1dB compression point and the third order intermodulation product. These parameters would depend on the saturation value and the sharpness variable.

## REFERENCES

[1] A. Stohr et al., "60 GHz radio-over-fiber technologies for broadband wireless services," *J. Opt. Commun. Netw.*, vol. 8, no. 5, pp. 471-487, May 2009.

[2] B. Charbonnier et al., "Photonics for broadband radio communications at 60 Hz in access and home networks," in *Proc. Int. Topical Meeting Microw. Photon.*, 2008, pp. 5-8.

[3] M. Huchard, M. Weiss, A. Pizzinat, S. Meyer, P. Guignard, and B. Charbonnier, "Ultrabroadband wireless home network based on 60-GHz WPAN cells interconnected via RoF," *J. Lightw. Technol.*, vol. 26, no. 15, pp. 2364-2372, Aug. 2008.

[4] A. Ng'oma and M. Sauer, "Radio-over-fiber technologies for high data rates wireless applications," in *Proc. IEEE Sarnoff Symp.*, 2009, pp. 1-6.

[5] A. Ng'oma, G. -J. Rijckenberg, and A. M. J. Koonen, "Building extended-reach radio-over-fiber links by exploiting optical frequency multiplication's dispersion tolerance," in *Proc. Int. Microw. Symp.*, 2007, pp. 2179-2200.

[6] C. Viana, Z. G. Tegegne, J. Polleux, and C. Algani, "Flexible new opto-microwave design approach for radio-over-fiber applications: A case study of low cost 60 GHz VCSEL-based IF-RoF link," *IEEE Trans. Microw. Theory Techn.*, vol. 66, no. 9, pp. 4293-4305, Sep. 2018.

[7] P. V. Mena, S.-M. Kang, and T. A. DeTemple, "Rate-equation-based laser models with a single solution regime," *J. Lightw. Technol.*, vol. 15, no. 4, pp. 717-730, Apr. 1997, doi: 10.1109/50.566695.

[8] M. U. Hadi, J. Nanni, O. Venard, G. Baudoin, J. L. Polleux, and G. Tartarini, "Practically feasible closed-loop digital predistortion for VCSEL-MMF-based radio-over-fiber links," *Radioengineering J.*, Czech Slovak Tech. Universities, vol. 29, no. 1, pp. 37-43, 2020.

[9] M. U. Hadi et al., "Assessment of digital predistortion methods for DFB-SSMF-radio-over-fiber links linearization," *Microw. Opt. Technol. Lett.*, vol. 62, no. 2, pp. 540-546, Jul. 2020.

[10] C. H. Cox, III, G. E. Betts, and L. M. Johnson, "An analytic and experimental comparison of direct and external modulation in analog fiber-optic links," *MMT*, vol. 38, no. 5, pp. 501-509, 1990.

[11] C. H. Cox, "Gain and noise figure in analogue fiber-optic links," *IEE Proc. J. Optoelectron.*, vol. 139, no. 4, pp. 238-242, 1992.

[12] C. H. Cox, E. I. Ackerman, and G. E. Betts, "Relationship between gain and noise figure of an optical analog link," in *Proc. IEEE MTT-S Int. Microw. Symp. Dig.*, 1996, vol. 3, pp. 1551-1554.

[13] C. H. Cox, III, E. I. Ackerman, G. E. Betts, and J. L. Prince, "Techniques and performance of Intensity-Modulation Direct-Detection analog optical links," *IEEE Trans. Microw. Theory Techn.*, vol. 45, no. 8, pp. 1375-1383, Aug. 1997.
[14] E. I. Ackerman and C. H. Cox, III, "State of the art in analog fiber-optic link technology," in *Proc. URSI Int. Symp. Signals, Syst. Electron.*, 1998, pp. 372-377.

[15] E. I. Ackerman, C. Cox, G. Betts, H. Roussell, F. O'Donnell, and K. Ray, "Input impedance conditions for minimizing the noise figure of an analog optical link," *IEEE Trans. Microw. Theory Techn.*, vol. 46, no. 12, pp. 2025-2031, Dec. 1998.

[16] E. I. Ackerman, "The effective gains and noise figures of individual components in an analog photonic link," in *Proc. Int. Topical Meeting Microw. Photon.*, 2003, pp. 369-372.

[17] C. H. Cox, E. I. Ackerman, G. E. Betts, and J. L. Prince, "Limits on the performance of RF-over-fiber links and their impact on device design," *IEEE Trans. Microw. Theory Techn.*, vol. 54, no. 2, pp. 906-920, Feb. 2006.

[18] E. I. Ackerman et al., "RF-over-fiber links with very low noise figure," *J. Lightw. Technol.*, vol. 26, no. 15, pp. 2441-2448, Aug. 2008.

[19] H. A. Haus, "Noise figure definition valid from RF to optical frequencies," *IEEE J. Sel. Topics Quantum Electron.*, vol. 6, no. 2, pp. 240-247, Mar./Apr. 2000.

[20] D. M. Baney, P. Gallion, and R. S. Tucker, "Theory and measurement techniques for the noise figure of optical amplifiers," *Opt. Fibre Technol.*, vol. 6, no. 2, pp. 122-154, 2000.

[21] R. S. Tucker and D. M. Baney, "Optical noise figure: Theory and measurements," in *Proc. Opt. Fiber Commun. Conf.*, 2001, vol. 3, Paper W11.

[22] I. Frigyes, B. G. Habermajaer, A. Molnar, A. Seeds, and F. Som, "Noise and loss characteristics of microwave direct modulated optical links," in *Proc. 27th Eur. Microw. Conf.*, 1997, vol. 1, pp. 309-317.

[23] G. Plantier, C. Bes, and T. Bosch, "Behavioral model of a self-mixing laser diode sensor," *IEEE J. Quantum Electron.*, vol. 41, no. 9, pp. 1157-1167, Sep. 2005.

[24] N. J. Gomes, P. P. Monteiro, and A. Gameiro, "Compensation of impairments in the radio over fiber infrastructure," in *Next Generation Wireless Communications Using Radio over Fiber*. Wiley, 2012, pp. 211-245, doi: 10.1002/9781118306017.ch10.

[25] K.-L. Chi et al., "Single-mode 850 nm VCSELs for 54 Gbit/sec on-off keying transmission over 1km multi-mode fiber," *IEEE Photon. Technol. Lett.*, vol. 28, no. 12, pp. 1367-1370, Jun. 2016.

[26] J. Estaran et al., "Quad 14 Gbps L-band VCSEL-based system for WDM migration of 4-lanes 56 Gbps optical data links," *Opt. Exp.*, vol. 20, no. 27, pp. 28524-28531, Dec. 2012.

[27] C. Crespo-Cadenas, J. Reina-Tosina, and M. J. Madero-Ayora, "Volterra behavioral model for wideband RF amplifiers," *IEEE Trans. Microw. Theory Techn.*, vol. 55, no. 3, pp. 449-457, Mar. 2007.

[28] A. Zhu and T. J. Brazil, "Behavioral modeling of RF power amplifiers based on pruned Volterra series," *IEEE Microw. Wireless Compon. Lett.*, vol. 14, no. 12, pp. 563-565, Dec. 2004.

[29] D. Mirri, G. Luculano, F. Filicori, G. Pasini, G. Vannini, and G. P. Gabriella, "A modified volterra series approach for nonlinear dynamic systems modeling," *IEEE Trans. Circuits Syst. I: Fundam. Theory Appl.*, vol. 49, no. 8, pp. 1118-1128, Aug. 2002.

[30] A. Zhu, J. C. Pedro, and T. J. Brazil, "Dynamic deviation reduction-based volterra behavioural modeling of RF power amplifiers," *IEEE Trans. Microw. Theory Techn.*, vol. 54, no. 12, pp. 4323-4332, Dec. 2006.

[31] M. Isaksson and D. Ronnow, "A Kautz-Volterra behavioral model for RF power amplifiers," in *Proc. IEEE MTT-S Int. Microw. Symp. Dig.*, 2006, pp. 485-488.

[32] A. Zhu and T. J. Brazil, "RF power amplifier behavioral modeling using volterra expansion with laguerre functions," in *Proc. IEEE MTT-S Int. Microw. Symp. Dig.*, 2005, pp. 963-966.

[33] H. Ku, M. D. McKinley, and J. S. Kenney, "Extraction of accurate behavioral models for power amplifiers with memory effects using two-tone measurements," in *Proc. IEEE MTT-S Int. Microw. Symp. Dig.*, 2002, vol. 1, pp. 139-142.

[34] P. M. Cabral, J. Pedro, and N. B. Carvalho, "Modeling nonlinear memory effects on the AM/AM, AM/PM and two-tone IMD in microwave PA circuits," *Int. J. RF Microw. Comput.-Aided Eng.*, vol. 16, no. 1, pp. 13-23, 2006.

[35] H. Wang and P. Liu, "Behavioral models for power amplifiers using a difference-frequency dual-signal injection method," in *Proc. PIERS*, Mar. 2010, pp. 950-953.

[36] A. J. Cann, "Nonlinearity model with variable knee sharpness," *IEEE Trans. Aerosp. Electron. Syst.*, vol. AES-16, no. 6, pp. 874-877, Nov. 1980.

[37] A. J. Cann, "Improved nonlinearity model with variable knee sharpness," *IEEE Trans. Aerosp. Electron. Syst.*, vol. 48, no. 4, pp. 3637-3646, Oct. 2012.

[38] S. L. Loyka, "On the use of cann's model for nonlinear behavioral-level simulation," *IEEE Trans. Veh. Technol.*, vol. 49, no. 5, pp. 1982-1985, Sep. 2000.

[39] C. J. O'Brien, M. L. Majewski, and A. D. Rakic, "A critical comparison of high-speed VCSEL characterization techniques," *J. Lightw. Technol.*, vol. 25, no. 2, pp. 597-605, Feb. 2007.