Research Article | Open Access
Snezhana Georgieva Gocheva-Ilieva, Desislava Stoyanova Voynikova, Iliycho Petkov Iliev, "Modeling of Output Characteristics of a UV Cu+ Ne-CuBr Laser", Mathematical Problems in Engineering, vol. 2012, Article ID 420782, 21 pages, 2012. https://doi.org/10.1155/2012/420782
Modeling of Output Characteristics of a UV Cu+ Ne-CuBr Laser
This paper examines experiment data for a Ne-CuBr UV copper ion laser excited by longitudinal pulsed discharge emitting in multiline regime. The flexible multivariate adaptive regression splines (MARSs) method has been used to develop nonparametric regression models describing the laser output power and service life of the devices. The models have been constructed as explicit functions of 9 basic input laser characteristics. The obtained models account for local nonlinearities of the relationships within the various multivariate subregions. The built best MARS models account for over 98% of data. The models are used to estimate the investigated output laser characteristics of existing UV lasers. The capabilities for using the models in predicting existing and future experiments have been demonstrated. Specific analyses have been presented comparing the models with actual experiments. The obtained results are applicable for guiding and planning the engineering experiment. The modeling methodology can be applied for a wide range of similar lasers and laser devices.
During the process of engineering design, together with the purely experimental investigations, various methods for mathematical modeling are also applied. Standard mathematical models usually include differential and integral equations, optimization, or other types of mathematical problems through which the processes occurring in the laser medium are described and the behavior of the system is explored by simulating its various states [1–5]. The disadvantage of these methods is the difficulty in obtaining an explicit estimate of the final output characteristics such as laser output power, overall efficiency, or service life of the device. If experiment data from the process of developing laboratory and industrial prototypes is available, another approach is the application of statistical methods. These are useful in finding basic relationships within data for the specific type of device and the utilization of obtained results for guiding and planning of the engineering experiment. In the field of metal vapor lasers, traditional parametric statistical methods have been used in [6–9]. In [6, 9], MARS models are applied to different relationships between the data for CuBr lasers in the visible range. For ultraviolet lasers with copper bromide vapors, some particular MARS models for examining laser output power have been developed in [6, 10].
This paper conducts a statistical investigation of laser generation and service life of UV copper bromide vapor lasers. This new type of laser source is currently being developed at the Laboratory of Metal Vapor Lasers, Bulgarian Academy of Sciences [11–13]. The main experiment results obtained so far, used in this paper, have been published in [14–20].
The goal of this study is to determine the relationships between 9 input laser characteristics and two output characteristics: laser generation and laser service life by applying the MARS method based on available experiment data. The following problems are solved: (1) determining which input quantities (supplied power, geometric design of the tube, neon pressure, reservoirs temperature, etc.) influence significantly the change of the considered output characteristics; (2) defining adequate regression models describing explicitly the relationship between independent variables and dependent quantities; (3) application of the models for estimating known experiments; (4) application of the models for predicting experiments, in particular for guiding design and construction work when developing new UV copper bromide vapor lasers; (5) using the models to investigate the local influence of the main input laser quantities on output ones.
2. Subject of Study
We investigate experiment data for an ultraviolet copper ion excited copper bromide vapor laser (UV Cu+ Ne-CuBr laser). This ultraviolet laser is a promising innovative product with unique characteristics and stable operation. The first lasers of this type were constructed by the Laboratory of MVLs at the Institute of Solid State Physics of Bulgarian Academy of Sciences in 1999 [11, 14]. It is characterized by high quality of the laser beam and high output power. During the last decade, this type of laser has been the subject of intense experimental studies and its characteristics have been significantly improved (see [19, 20] and the cited therein detailed literature with experimental results).
Ultraviolet laser generation can be achieved using copper or gold ions. The typical technological problems of copper ion lasers—copper deposits on the inside walls of the tube and high discharge temperature—were solved thanks to the development of a new UV Ne-CuBr laser with a nanosecond pulsed longitudinal discharge. This laser emits in the deep ultraviolet spectrum on five spectral lines −248.6, 252.9, 259.7, 260.0, and 270.3 nm. Experimentally, a record average output power of 1.3 W has been obtained at all five lines and 0.85 W at the 248.6 nm line, but the laser emission was achieved only for a period of 20 hours, due to the technological problems [15, 16]. It has also been determined that adding small amounts of hydrogen (0.02–0.04 Torr) results in a twofold increase of laser generation . Further improvement by overcoming this problem was reported in .
This laser has large fields of application including medicine, microelectronics, photolithography, genetic engineering, and scientific research [23–25]. Due to the narrow emission range of just a few spectral lines and the high coherence of the beam, it is used for processing, which requires high resolution, such as recording information, fluorescence, high-precision drilling, cutting, cleaning, modification of newly developed materials, and so forth. Technological applications of the UV lasers created at the Laboratory of Metal Vapor Lasers are achieved. Laser-induced modification by the use of UV laser radiation with 248.6 nm wavelength has been performed in a conducting polymer. A significant growth in the electrocrystallization of copper on polymer layers was obtained. The refractivity of polymer layers in ultraviolet spectral regions was determined. Micronic holes of 10, 20, and 40 microns in size were drilled in glass, polymer, and ZnSe .
The goal of subsequent development of the examined UV lasers is to investigate the possibilities for improving its output characteristics: output power and service life.
3. Data Description
A total of 9 laser characteristics (independent variables, predictors) are studied, which are of physical significance and are considered independent. The variables we will examine are as follows: , mm: inside diameter of the laser tube; DR, mm: inside diameter of the internal rings in the ceramic tube or quartz tube insert; , cm: length of the active zone (distance between the electrodes); , kW: electric power supplied to the active volume, taking into account 50% losses of the total supplied power ; , Torr: neon gas pressure; , kHz: pulse rate frequency; , Torr: hydrogen gas pressure; , pF: equivalent capacity of the condensation battery; TR, °C: temperature of the copper bromide reservoir (equal to the temperature of the outer side of the quartz tube).
The dependent variables are laser output power (laser generation) , W, and laser service life , hours.
This study uses the data from experiments for 11 different laser devices, reported in [11–20]. The scatter plot of the relationships for each two of the eleven variables (9 independent and 2 dependent) is given in Figure 2.
The descriptive properties of the data indicate that applying the standard (conventional) parametric methods is not recommendable as the basic conditions for their validity are violated . More specifically, the two dependent variables are not normally distributed. There is also certain collinearity between independent variables. For that reason, in this paper, nonparametric statistical techniques are used, choosing in this case the MARS method.
4. MARS Method
The relatively new nonparametric MARS technique that performs well with complex data structures has been applied [21, 22, 26]. MARS allows building flexible models capable of describing local relationships in a multidimensional region of data, within which subregions are defined. For each of these multidimensional subregions, the most suitable type of relationship is chosen and a general model is built, which approximates as closely as possible the actual data in accordance with a specific criterion.
Let us denote the independent variables by and their dependent one by . We assume that are vectors with dimension . For our data , . The general form of a MARS model of the data for has the form where , , are the sought constant coefficients of the model, the number of functions included in the model, and basis functions of the selected type.
In the linear case, basis functions are presented in the form of one of the following two one-dimensional “mirror” functions: Here is a constant, called a knot of the basis function with values in the definition interval of the predictor variable . In the linear case, model (4.1) presented as a sum of basis functions of type (4.2) is plotted graphically in the form of a piecewise linear function.
When building nonlinear MARS models, basis functions can also include the products of two or more functions of the form (4.2) without repeating the index . The maximum number of multipliers is called order of interaction. The linear model has an order of interaction 1.
We have to note that unlike ordinary interpolation splines, in the MARS model (4.1), the knots of the basis functions, their number , the selection of functions, and the subintervals of determination are not known in advance. These are determined by different optimization conditions and prerequisites for estimating the proximity of the model to real data on .
When building the models, the researcher sets an initial maximum number of basis functions and the maximum order of interaction. The recommendation is . MARS procedures systematically select subintervals and subregions for each predictor, estimate the degree of significance to the model of a given basis function, and exclude those that do not contribute to the improvement of the model. Several different criteria are used to evaluate the model after each step. One criterion is the minimization of the sum of squares of errors: for a model with a number of current basis functions. Another basic criterion is the minimization of the coefficient of general cross-validation (GCV) [19, 20]. According to which “the best” model of models with basis functions is the one that minimizes the expression where , , , .
As a whole, the implementation of MARS method algorithms allows the automatic definition of those predictors that influence the examined dependent variable, as well as the degree of this influence, removing the statistically insignificant predictors. The final model usually requires that the residuals of the model are normally distributed.
An advantage of the method is that the resulting models are simpler in form and can easily be interpreted during their practical application.
Within this study only the best MARS models are presented. All nine variables are initially introduced as predictors. We use the notation for a model with initial maximum number of allowed basis functions (BFs) and up to th order of interactions, .
5. Building MARS Models of the Laser Generation of an UV Laser
First, we will present the best MARS models of laser generation .
The basic statistic figures of the constructed models are given in Table 2. We will describe in more detail two of these models: the best piecewise linear model with no interactions and the best model with first-order interactions.
5.1. Linear MARS Models of
Out of the linear models, we will present model (40;1). The figures of other linear models are given in the upper part of Table 2.
The generated best MARS model (40;1) includes the following 13 basis functions (BFs) with 8 predictors:
The obtained regression model for laser output power with these functions is
We calculate each of the functions in (5.1) consecutively. We have , , , and so forth. We substitute in model (5.2) obtaining the estimate The absolute error in case is about 60 mW, and the relative one is 7%. The overall relative error for the model is 8-9%, which is a satisfactory figure.
Model (5.1)-(5.2) is the best MARS model of a given type, which is selected so as to allow no overfitting of the model, as well as by using the measures of data fit SSE and GCV. The obtained basic statistics are given in Table 2. The model is significant at level 0.000.
The relative influence of individual predictors in model (5.2) is given in Table 3, column 2. It is apparent that the most significant variable is , whose influence is measured as 100%, and the influence of the others is calculated against it.
5.2. Nonlinear MARS Models
Of the second-order models we present the best model with up to 40 BF, which is denoted as model (40; 2). The obtained model includes the following 28 BFs:
The corresponding regression model for laser output power includes all of these functions with 8 independent variables. It has the following form:
The absolute error in case is about 44 mW, and the relative one is 5%.
A part of the partial contributions of pairs of predictors in model (5.5)-(5.6) to the value of are presented in Figures 3(a)–3(d). The biggest contribution of more than 6000 units is made by the interaction between and , which reaches the highest value in a large 2D region. Predictors and provide almost the same contribution. The other two interactions also have given effects.
5.3. Application of the Models for
For example, for the piecewise linear model (5.1)-(5.2), the graphs presented in Figure 3 show the local relationships between individual predictors and the dependent variable . Figure 3(a) indicates that the pressure of the buffer gas (neon) , which makes the biggest contribution to the model, should be taken within the interval (16, 21) Torr. For the inner diameter of the laser tube, the high values are achieved at mm (Figure 3(b)), which has been established experimentally in [16, 17]. The behavior of the pressure of hydrogen admixtures also fits quite well for all experiments, exhibiting maximum local influence in [0.02, 0.04] Torr (Figure 3(e)).
Analogically, for the nonlinear model (5.5)-(5.6), Figures 4(a)–4(d) show local relationships between the main pairs of significant input parameters in two-dimensional regions. The supplied electrical power makes the biggest contribution to the model and shows an increasing behavior by comparing Figures 4(a), 4(c) and 4(d). Figures 4(b) and 4(c), indicate that the neon pressure should be taken within the interval (15, 25) Torr. The hydrogen pressure gives an optimal contribution for [0.03, 0.04] Torr, according to Figures 4(d) and 4(b).
In addition to the mutual influence of predictors, it is also possible to plot slices in order to examine the local behavior of each of the two variables with one fixed at a value chosen by the researcher and the other changing within the whole interval of definition.
As an example, Figure 5(a) shows a slice from Figure 4(b) for the modeled behavior of the variable (neon gas pressure) for a fixed value of the hydrogen pressure Torr. In this case, there is a clearly identifiable maximum in the interval [11, 26] Torr. The slice in Figure 5(b) shows the distribution of the influence of the pair of predictors for a fixed Torr, depending on the increase of , as an element of Figure 4(c), which shows the maximum importance.
By carefully studying the resulting estimates of from a given model and the defined influences of the significant variables, as well as their interactions, it is possible to guide the experiment in order to improve the output characteristics.
In addition to the plots, the models can be used to examine the local influence of the main laser input quantities on output ones, as well as to predict new experiments. In order to demonstrate this, we will consider the local behavior of for the laser from  when varying the pressure of the neon buffer gas in the interval [15.5, 18.5] Torr. For we choose experiment data with fixed , , , , and . In Figure 6(a) continuous line indicates the experiment data, Predicted shows estimates from model (40;2).
Results for a future experiment prediction using model (40; 2) are given in Table 4. The following are fixed: mm, cm, Torr, Torr, kHz, °C, and pF. As expected, laser generation increases when input electric power is increased. In this table, the middle row #3 corresponds to an actual experiment with the measured value of mW, and the one predicted by the model is 1145 mW, which represents a relative error of under 5%.
6. MARS Models for Laser Service Life
This type of research has not been performed for metal vapor lasers so far. The models utilize as predictors all of the 9 input laser characteristics from Section 3 with response —laser service life, measured in hours. Again, our objective is to find a model with a sufficiently high coefficient of determination (over 95%) and the smallest possible mean-square error (under 5%) comparable with the experimental error.
The statistics of the obtained models and basic statistics are given in Table 5. The first best MARS model with 98% coefficient of determination is model (40;2).
6.1. Linear MARS Models of
Of the linear-type models, we will present model (30;1). The figures of other linear models are given in the upper part of Table 5.
The calculated best MARS model (30;1) includes the following 6 basis functions with 4 predictors:
The regression model for laser service life with these functions is
The relative influence of individual predictors in model (6.2) is given in Table 6, column 2. It is apparent that the most significant variables are , , , and (in descending order according to their influence). The behavior of the 4 main predictors in their intervals of determination and contribution in model (6.1)-(6.2) are shown in Figures 7(a)–7(d) in pure ordinal units.
6.2. Nonlinear MARS Models of
Of the second-order models we present best model with up to 30 BFs, which is denoted as model (30; 2). The obtained model includes the following 16 BFs: The corresponding regression model for laser service life includes 13 of these functions and 6 predictors. It has the following form:
6.3. Application of the Obtained MARS Models of
Model (6.3)-(6.4) analogically allows for estimation and prediction of the service life of the considered laser device depending on a given set of values of participating predictors. What is essential here is the relative influence of these predictors, which is to be taken into consideration by the design engineer when planning the experiment, in our case maintaining strictly the temperature of the reservoirs TR (equal to the temperature maintained at the outer wall of the laser tube (quartz or ceramic), and keeping the ratio between the pressure of the neon buffer gas and the bank of condenser within the established limits.
7. Comparison and Diagnosis of the Obtained MARS Models
Model comparison can be performed using the general statistical indices from Tables 2 and 5. It has to be noted that in the last column the predictors for the respective model are presented in descending order according to their relative importance.
The lowest values of and GCV are in linear models () both for and for . With the increase of the number of basis functions (BFs) the parameters of the models increase slightly. However, the investigation of the residuals of the linear models shows that they are not sufficiently adequate and that the relative error is comparatively high within 10–15%. This means that the relationship of the input independent variables with or is not linear.
For output power , Table 2 shows that second-order models from (40;2) exhibit a 98% coefficient of determination and that the participation of the predictors is stable. Since for (40;2) the GCV estimate is over 95%, this model is quite good. The next models (up to (60;2)) demonstrate practically the same fitting properties.
Of all models of , model (30;2) exhibits the best qualities.
Third-order models are comparable with the respective second-order models and demonstrate almost the same or slightly lower indices than second-order ones. Since these models are more complex in form, we conclude that the examined data demonstrate second degree local nonlinearities and it is best to describe these using second-order models.
In order to examine the diagnostics we have to note that all presented models are statistically significant at level . The coefficients of the models are significant with . The standardized error of the estimate is small. Model residuals are normally distributed.
Figure 9 shows the comparison of experimental data for laser output power with those calculated by model (40;2). It is observed that the model fits data quite well. For the coefficient of determination we have , which means that the model accounts for 98% of all data. The corresponding (see Table 2). The residuals of this model are normally distributed with .
This way, the two main criteria and statistical indices show very good predictive properties and goodness of fit of the constructed MARS model (40;2).
The fitting property of model (30;2) for against the experiment is given in Figure 10.
8. Physical Interpretation of the Results
The developed models correspond quite well with the experiment and reflect nonlinear local relationships in a multidimensional space of 9 variables.
Laser output power is most significantly influenced by input electric power, neon pressure, and hydrogen pressure. Of the geometric dimensions, output power is most dependent on the diameter of the rings DR, which for latest laser devices is equal to D—the diameter of the inner tube, with its optimal value between 7.1 and 7.3 mm.
The new results are those for the service life of the examined lasers, which is one of the most important issues for all types of devices. For the latest lasers  this period is quite acceptable, reaching 700–1000 hours. Furthermore, periodically a refreshment of laser tube is made, which involves shutting down the laser, cleaning, and changing the gas mixture. The obtained models unequivocally show that the main contribution to extending the service life of the laser is that of TR (inversely proportional—Figure 7(a)), C, and . Therefore, these characteristics need to be considered carefully when designing new lasers of the investigated type and their parameters should be constant or nearly constant in relation to those already established.
The basic results of the performed statistical modeling of laser output characteristics, laser generation and laser service life of a multiline ultraviolet copper ion vapor laser excited in a longitudinal pulsed discharge, have been obtained using nonparametric MARS models. The models demonstrate very good abilities and goodness of fit when predicting existing and future experiments.
It was determined that the best MARS models are nonlinear and contain second-order members. Of 9 input independent variables used as predictors, 8 influence laser output power with the most significant ones being the pressure of the applied electric power, the inner diameter of the laser tube (diameter of the rings), neon buffer gas pressure, and hydrogen pressure. Service life models are simpler and only 3 laser characteristics exhibit significant influence: temperature of the tube, neon gas pressure, and equivalent capacity of the condensation battery.
It is shown that the models can be used to estimate and predict local properties of laser generation and service life of the devices.
The techniques developed can also be further employed in the process of future industrial application of this type of lasers.
This paper is published in cooperation with Project of the Bulgarian Ministry of Education, Youth and Science, BG051PO001-3.3-05/0001 “Science and Business,” financed under Operational program “Human Resources Development” by the European Social Fund.
- M. J. Kushner and B. E. Warner, “Large-bore copper-vapor lasers: kinetics and scaling issues,” Journal of Applied Physics, vol. 54, no. 6, pp. 2970–2982, 1983.
- “Numerical modeling of low-temperature plasmas,” in Encyclopedia of Low-temperature Plasma, Series B, M. Ianus, Ed., vol. 7, Moscow, Russia, 2004.
- A. M. Boichenko, G. S. Evtushenko, and S. N. Torgaev, “Simulation of a CuBr laser,” Laser Physics, vol. 18, no. 12, pp. 1522–1525, 2008.
- B. L. Pan, G. Chen, B. N. Mao, and Z. X. Yao, “Kinetic process of UV Cu+ laser in Ne-CuBr longitudinal pulsed discharge,” Optics Express, vol. 14, no. 19, pp. 8644–8653, 2006.
- M. Bang-Ning, P. Bai-Liang, C. Li, W. Ya-Juan, and Y. Zhi-Xin, “Kinetic analysis of the factors limiting the output power of the Ne-CuBr UV laser,” Chinese Physics B, vol. 18, no. 4, pp. 1542–1546, 2009.
- S. G. Gocheva-Ilieva and I. P. Iliev, Statistical Models of Characteristics of Metal Vapor Lasers, Nova Science Publishers, New York, NY, USA, 2011.
- I. P. Iliev, S. G. Gocheva-Ilieva, D. N. Astadjov, N. P. Denev, and N. V. Sabotinov, “Statistical analysis of the CuBr laser efficiency improvement,” Optics and Laser Technology, vol. 40, no. 4, pp. 641–646, 2008.
- I. P. Iliev, S. G. Gocheva-Ilieva, D. N. Astadjov, N. P. Denev, and N. V. Sabotinov, “Statistical approach in planning experiments with a copper bromide vapor laser,” Quantum Electronics, vol. 38, no. 5, pp. 436–440, 2008.
- S. G. Gocheva-Ilieva and I. P. Iliev, “Parametric and nonparametric empirical regression models: case study of copper bromide laser generation,” Mathematical Problems in Engineering, vol. 2010, Article ID 697687, 15 pages, 2010.
- S. G. Gocheva-Ilieva and I. P. Iliev, “Modeling and prediction of laser generation in UV copper bromide laser via MARS,” in Proceedings of the 5th International Conference on Optics, Astrophysics and Astronomy (ICOAA '10), O. Martin, Ed., Advanced Research in Physics and Engineering Series, Mathematics and Computers in Science and Engineering, pp. 166–171, WSEAS Press, Cambridge, UK, 2010.
- N. Vuchkov, N. Sabotinov, and K. Temelkov, “Ultraviolet Copper laser,” Patent Reg. no. 103312, 06.04.1999, Published Applications for Inventions, BG 10/2000.
- N. K. Vuchkov, K. A. Temelkov, P. V. Zahariev, and N. V. Sabotinov, “Laser tube for ultraviolet copper laser,” Patent reg., no. 106759, 30.05.2002, Published Applications for Inventions BG 11/2003.
- N. K. Vuchkov, K. A. Temelkov, P. V. Zahariev, and P. V. N. V, “Laser tube for ultraviolet copper laser,” Bulgaria Patent No.64880, 03.08.2006.
- N. K. Vuchkov, K. A. Temelkov, and N. V. Sabotinov, “UV lasing on Cu+ in a Ne-CuBr pulsed longitudinal discharge,” IEEE Journal of Quantum Electronics, vol. 35, no. 12, pp. 1799–1804, 1999.
- N. K. Vuchkov, K. A. Temelkov, P. V. Zahariev, and N. V. Sabotinov, “Optimization of a UV Cu+ laser excited by pulse-longitudinal Ne-CuBr discharge,” IEEE Journal of Quantum Electronics, vol. 37, no. 4, pp. 511–517, 2001.
- N. K. Vuchkov, K. A. Temelkov, P. V. Zahariev, and N. V. Sabotinov, “Influence of the active zone diameter on the UV-ion Ne-CuBr laser performance,” IEEE Journal of Quantum Electronics, vol. 37, no. 12, pp. 1538–1546, 2001.
- N. K. Vuchkov, K. A. Temelkov, P. V. Zahariev, and N. V. Sabotinov, “Output parameters and a spectral study of UV Cu+ Ne-CuBr laser,” Optics and Laser Technology, vol. 36, no. 1, pp. 19–25, 2004.
- N. K. Vuchkov, K. A. Temelkov, P. V. Zahariev, and N. V. Sabotinov, “A new laser tube construction for the UV Cu+ Ne-CuBr laser,” in Laser and Laser Information Technologies, pp. 173–176, bgr, October 2003.
- N. K. Vuchkov, K. A. Temelkov, and N. V. Sabotinov, “Effect of hydrogen on the average output power of the UV Cu+ Ne-CuBr laser,” IEEE Journal of Quantum Electronics, vol. 41, no. 1, pp. 62–65, 2005.
- N. Vuchkov, “High discharge tube resource of the UV Cu+ Ne-CuBr laser and some applications,” in New Development in Lasers and Electric-Optics Research, pp. 41–74, Nova Science Publishers, New York, NY, USA, 2007.
- D. Steinberg, B. Bernstein, P. Colla, and K. Martin, MARS User Guide, Salford Systems, San Diego, Calif, USA, 2001.
- 2012, http://www.salfordsystems.com/mars.php.
- M. Ilieva, V. Tsakova, N. K. Vuchko, K. A. Temelkov, and N. V. Sabotinov, “UV copper ion laser treatment of poly(3,4-ethylenedioxythiophene),” Journal of Optoelectronics and Advanced Materials, vol. 9, no. 2, pp. 303–306, 2007.
- M. Ilieva, A. Stoyanova, V. Tsakova et al., “Effect of deep UV laser treatment on electroless silver precipitation in supported poly-3,4-ethylenedioxythiophene layers,” Journal of Optoelectronics and Advanced Materials, vol. 11, no. 10, pp. 1444–1447, 2009.
- J. Mouchovski, K. Temelkov, N. Vuchkov, and N. Sabotinov, “Calcium strontium fluoride crystals with different composition for UV-Laser application. Control of the growing rate and optical properties,” Comptes Rendus de L'Academie Bulgare des Sciences, vol. 62, no. 6, pp. 687–694, 2009.
- J. H. Friedman, “Multivariate adaptive regression splines,” The Annals of Statistics, vol. 19, no. 1, pp. 1–141, 1991.
Copyright © 2012 Snezhana Georgieva Gocheva-Ilieva et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.