# The Linear Fitting Method for Model of PIN Array Receiver in Space Optical Communication

Authors: H. Qinggui
Journal: International Journal of Optics
 (2022)

---

## Abstract

We have proposed PIN array as receiver in space optical communication and adopted the double diode model, the widely accepted solar energy output theory, to study its characteristics of output current and voltage. In order to make the calculation simplified, the linear fitting method is put forward. Later, we adopt both the traditional method Newton-Raphson and the new linear fitting method to calculate the values of both maximum power and the corresponding voltage of 3 × 3 PIN array. The comparison result shows the calculated values with two methods are good consistent. It validates the feasibility of the new linear fitting method. In the next step, the experiment has been carried out. The experimental result validates the reasonability of adopting the double diode model to study PIN array. At the same time, it validates the feasibility of the new linear fitting method again.


---

## Body

## 1. Introduction

The combination of free space optical communication and ordinary optical fiber communication is considered to be a hopeful method to solve “the Last Mile Issue of Communication.” As to the combination, the coupling of space light to optical fiber is the key issue. In order to address this question, many scholars have conducted extensive research [1–7]. In 2013, Zhang et al. set up the space optical coupling alignment platform; they adopted both the position sensor and fast reflector to track the spatial light spot [1]. In 2014, Luo et al. used adaptive fiber couplers (AFC) to establish a closed loop control system to improve the coupling efficiency [2]. In 2014, Lei et al. adopted both the stepper motor and self-focusing lens to find the best coupling position in fiber array [3]. In reference [4], the authors present a new cladding design for high birefringence and low confinement loss photonic crystal fibers (PCFs) using a full-vector finite element method with anisotropic perfectly matched boundary layer. In reference [5], the authors propose a novel high-birefringence index-guiding photonic crystal fiber (PCF). This PCF is composed of a solid silica core and a cladding with two differently sized squeezed elliptical air-holes. In our previous works, we have proposed the conical fiber receiver to improve coupling efficiency and have studied its features from various perspectives [6–9]. Generally speaking, in these researches, the single optical fiber is used as the receiving device. The receiving area of the single optical fiber is too small; it brings difficulty for the coupling.In recent years, some scholars have begun to apply PIN photodiodes in FSO communication systems. In [10], the authors design and demonstrate an integrated receiver for a visible light communication (VLC) system based on LED and an array of silicon PIN diode detectors. In [11], based on PIN photodiode receivers, the authors employ both asymmetric density evolution technique and an optimization scheme combined search algorithms to optimize low-density parity-check (LDPC) codes for FSO systems. In [12], with the cascaded preequalization circuit, and a differential outputs PIN receiver, the authors proposed a cascaded amplitude equalizer for high-speed visible light communications (VLC) system. In our previous work, we have proposed 3×3 PIN diode array to receive the spatial light, each PIN photodiode could independently receive the space laser and convert it into electrical signals [13].The above literatures have studied PIN photodiode from various perspectives but failed to consider the effect of uneven illumination. In PIN diode array, due to the different light intensities, different PIN diode produces different amount of photocurrent. In order to track the maximum output power, we should build mathematical models to analyze the characteristics of output current and voltage. But such research is relatively lacking. For the question of how to track the maximum power, there are many interesting studies in the field of solar energy technology. They discovered the single photovoltaic cell can be regarded as a photoelectric cell. Its current and voltage characteristics can be described with double diode model. Then, the double diode equation in photovoltaic model is used to describe the current and voltage characteristics of PIN diode. Later, Newton-Raphson method is adopted to resolve the resolution. In this work, we attempt to adopt this method to study PIN array receiver in space optical communication. This is our innovation. Considering that the calculation of Newton-Raphson method is too complicated, we adopt the linear fitting method to do the calculation work. This is our innovation as well. At the same time, the experiment is carried out whose results prove the proposed linear fitting method is reasonable.


## 2. The Equivalent Circuit Model

In the photovoltaic cell, P-N junction absorbs light energy and turns it into electrical energy. In order to prevent current backflow, the bypass diode is required in the circuit. Based on engineering practice and theoretical analysis, the single photovoltaic cell can be regarded as a photoelectric cell. Its current and voltage characteristics can be described with double diode model [14–16]. For single PIN diode, it is similarly true [17, 18]. We can use the double diode equation in photovoltaic model to describe the current and voltage characteristics of PIN diode, as shown in Figure 1. Normally, the value of Rsh is large, its impact on the open-circuit voltage and short-circuit current is little; we can ignore its influence. Then, the output current of single PIN diode I is given by [19, 20].(1)I=Iph−Id,where Iph is the photocurrent of single PIN diode; and Id is the diode current; it is given by [19, 20](2)Id=IoexpqkATU+IRS−1,where Io is the reverse saturation current; it could be determined by the following equation [19, 20]:(3)Io=IrrTTr3expqEgAk1Tr−1T.Figure 1
The equivalent circuit of a single PIN diode.In the above two equations (2) and (3), q is the electronic charge (1.6 × 10−19C), k is the Boltzmann constant (it is 1.38 × 10−23 J/K); T is the temperature of single PIN diode, Tr is the temperature under standard test condition, Irr is the reverse saturation current on the temperature Tr, Rs is the equivalent series resistance (usually, it is hundred milli ohms or dozen milli ohms [21, 22]), A is the ideal factor of the single PIN diode, and Eg is the band gap energy.For the PIN diode array, as shown in Figure2, to assume the number of PIN diodes in series is Ns, that in parallel is Np, then the total number n is Ns×Np, we can gain the array’s V-I characteristic equation [21, 22]:(4)I=NpIph−NPIoexpqkATUNs+IRsNp−1.Figure 2
The equivalent circuit model of PIN array.


## 3. The Linear Fitting Method for PIN Diode Model

Formulas (1) and (4) are general mathematical models. However, they belong to exponential transcendental equations. Usually, we should adopt Newton-Raphson method to calculate it. But, this method needs large computation costs [23, 24]. Here, we put forward the new linear fitting method to simplify the calculation process.With the MATLAB software, we can adopt the traditional method to figure out the relationship between the light intensitySi and the voltage under maximum output power condition (we call it as maximum power voltage Vmpp). To assume the PIN array is 3×3, the ideal factor A is 1.2, and the band gap energy Eg is 1.12, the temperature unchanged (T is 25oC), then we can gain the relationship between Vmpp and Si, as shown in the actual curve in Figure 3.Figure 3
The actual curve and fitting curve betweenVmpp and Si.As shown in the actual curve in Figure3, we can use the following linear equation to describe the relation between Vmpp and Si [25, 26].(5)Vmpp′=k1lnSiSir+b1,where Vmpp′ is the fitting maximum power voltage, and Sir is the standard light intensity (100 mW/cm2). This equation’s meaning is that we can use this fitting equation to describe the relation between Vmpp and Si, as the fitting curve shown in Figure 3.To assume light intensity maintained unchanged, for example, it is 20 mW/cm2. We use the traditional Newton-Raphson method to resolve the mathematical model formula from (1) to (4). Then, we can figure out the relationship between maximum power voltage Vmpp and temperature T, as shown in Figure 4. From those data points in the figure, we use equation (6) to describe the relation between the fitting maximum power voltage Vmpp′ and the maximum power voltage Vmpp,(6)Vmpp′=Vmpp+kT25−T,where kT is the temperature coefficient.Figure 4
The relations betweenVmpp and T under different Si.In the next step, whenSi are 40, 60, 80, and 100 mW/cm2, we can repeat the above calculation process to gain the different fitting lines, as shown in Figure 4. We can see that, when Si is different, the influence of temperature on the relations between Vmpp′ and Vmpp is different. Then, we use equation (7) to describe the relation between the temperature coefficient kT and the light intensity Si.(7)kT=k2Si+b2.Similarly, whenSi are 20, 40, 60, 80, and 100 mW/cm2, we can gain the relation curves between the maximum output power Pmpp and Vmpp, as shown in Figure 5. From the figure, we can also discover that when Si is maintained unchanged, we can use the following linear equation to describe the relation between Pmpp and Vmpp.(8)Pmpp=kvVmpp+bv.Figure 5
The relations betweenVmpp and Pmpp under different Si.In the above equation, the coefficientskv and bv are related with Si; we can look them as the functions of the light intensity Si.(9)kv=k3Si+b3,bv=k4Si+b4.From the above process, it can be seen that we can use the fitting linear equations to describe the relation betweenPmpp, Vmpp, and Si. The fitting linear equations can make the solving process simplified.


## 4. The Calculation Comparison between the Linear Fitting Method and Traditional Newton-Raphson Method

For the coefficientsk1 and b1 in equation (5), we adopt the least-square method to compute their values [14–16]. The specific method is like the following steps. We select a lot of points on the actual curve in Figure 2, use Sii to represent their abscissa values, use Vmppi to represent their ordinate values, and use n to represent the total number of discrete points. Then, the coefficients k1 and b1 are given by(10)k1=∑Sii∑Vmppi−n∑SirVmppi∑Sii2−n∑Sii2,b1=∑SiiVmppi∑Sii−∑Vmppi∑Sii2∑Sii2−n∑Sii2.To assume the PIN array is3×3, the ideal factor A is 1.2, and the band gap energy Eg is 1.12, with the above least-square method, we can work out the values of the coefficients k1 and b1, which are k1=1.082 and b1=17.678, and in this way, we can gain the values of other coefficients, like the following:(11)k2=0.0,b2=0.076,k3=0.039,b3=−0.541,k4=−0.044,b4=8.540.In order to verify the effectiveness of this method, we make a comparison. Firstly, we adopt traditional Newton-Raphson method to calculate the values of both maximum power and the corresponding voltage, looking them asPmpp1 and Vmpp1. Then, we adopt new method to calculate the values of both maximum power and the corresponding voltage, looking them as Pmpp2 and Vmpp2. Later, the relative error δ is defined as the following equations.ForVmpp,(12)δ=Vmpp1−Vmpp2Vmpp1×100%.ForPmpp,(13)δ=Pmpp1−Pmpp2Pmpp1×100%.Finally, the result is shown in Table1.Table 1
Comparison between Newton-Raphson and proposed methods.
Temperature/°CLight intensity/(mW/cm2)Vmpp1/VVmpp2/VErrorδ/%Pmpp1/WPmpp2/W2510017.7517.680.3961.8162.511.13358016.3216.300.1245.7646.050.6506014.8714.611.732.3732,650.86604013.2513.011.8120.1020.311.04702011.8211.790.2510.0310.464.28From Table1, it can be seen that both Vmpp1 and Pmpp1 are coincident with Vmpp2 and Pmpp2, respectively. It validates the feasibility of the new linear fitting method. Specifically, for the relative error δ, it changes from 0.25% to 4.28%; this value is very small. In other words, the matching degree is high.In the view of the computational complexity, the traditional Newton-Raphson method needs to solve the transcendental equation. It includes the exponent arithmetic operations and needs multiple iterations. In a word, it needs much calculation cost. However, the linear fitting method does not need so much calculation. This method greatly simplifies the calculation process.


## 5. Experimental Test

In order to check the reasonability of the new linear method again, the author had a discussion with the technicians of Sichuan Coreway Optech Co., Ltd. At the same time, the test experiment is done.In our previous work, as shown in Figure6, we have produced 3×3 PIN array to do the experiment of space optical communication [13].Figure 6
The appearance of the encapsulated PIN array.At this time, we also use this PIN array to do the test experiment. In this experiment, we use TES-1339 illuminometer to test the light intensity. TES-1339 illuminometer is produced in Taiwan; it is widely used in various industries and fields [4–7]. The experiment starts at 6:00 a.m. and ends at 2:00 p.m. At the same time, we adopt the new linear fitting method to calculate the theoretical data. In order to compare the experiment data with the theoretical data, we define the error rate in the following equation.Error rate = (| measured data−theoretical data |/measured data) × 100%.Finally, the measurement results and the error rate are shown in Table2.Table 2
Comparison between experimental and proposed theoretical data.
TimeTemperature (°C)Light intensity (W/cm2)Power (W)Error rate (%)Measured dataSimulation data6:001474.34.79.37:0017168.38.97.28:0018.2199.29.53.29:0020.1189.09.22.210:0022201110.54.511:00232111.311.51.712:00252512.312.73.213:00272211.611.92.614:00282311.912.22.5From Table2, it can be seen that when the light intensity is low, the error rate is large. On the contrary, when the light intensity is high, the error rate is small. For example, when the light intensities are 7 W/cm2 and 16 W/cm2, the error rates are 9.3% and 7.2%, respectively. When the light intensities are 22 W/cm2 and 23 W/cm2, the error rates are 2.6% and 2.5%. This is because the losses of the instrument maintain unchanged. When the light intensity is lower, the fixed losses account for higher proportion. Then, the error is bigger. On the contrary, when the light intensity is higher, the fixed losses account for lower proportion. Then, the error is smaller.From Table2, it can be seen that the temperature also influences the error rate. When the temperature is low, the error rate is high. On the contrary, when the temperature is high, the error rate is low. For example, when the temperature is 14°C, the error rate is 9.3%. However, when the temperature is 28°C, the error rate is 2.5%. Maybe it is because of slightly underestimating the contribution of temperature in the double diode model. The new linear fitting method is deduced from the mathematical model. Thus, it produces the errors. In addition, the experimental data will be affected by the environmental parameters of last time. But this factor is not considered in the mathematical model. On the whole, the average error is about 4%; it shows the experimental result is in relatively good agreement with the numerical calculation result. The experiment proves the reasonability of the new linear fitting method again.


## 6. The Conclusion

The double diode model based on semiconductor diode theories have been widely used to extract the maximum power point of the solar PV array. In this paper, we use this model to study PIN array in space optical communication. Because the mathematical model contains transcendental equation, it makes the solving process complicated. In order to make the calculation simple, we put forward the new linear fitting method to do the calculation process. The comparison of calculation results has been conducted between the new and traditional methods. It shows the relative errorδ changes from 0.25% to 4.28% in Table 2; this value is very small. In other words, matching degree is high. It validates the feasibility of the new linear fitting method. In addition, the test experiment has been carried out. The experimental data are compared with the theoretical data. By the comparison, we find it is reasonable to adopt the double diode model to study PIN array. At the same time, the experiment validates the feasibility of the new linear fitting method once more. This work shows we can use the linear fitting method to trace the maximum power of PIN array in FSO communication systems.



---

*Source: 2712863-2022-03-22.xml*