# Modeling of Output Characteristics of a UV Cu+ Ne-CuBr Laser

Authors: Snezhana Georgieva Gocheva-Ilieva; Desislava Stoyanova Voynikova; Iliycho Petkov Iliev
Journal: Mathematical Problems in Engineering
 (2012)

---

## Abstract

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.


---

## Body

## 1. Introduction

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.The examination has been conducted using MARS statistical software [21, 22].


## 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 [19]. Further improvement by overcoming this problem was reported in [20].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 [20].The schematic of the ultraviolet Ne-CuBr laser tube is shown in Figure1. Its main technical characteristics are given in Table 1.Table 1
Technical characteristics of a UV copper ion vapor laser [20].
Characteristic
Values
Radiation wavelength
248.6, 252.9, 259.1, 260.0, 270.3 nm
Operating mode
Periodic pulse, self-heating
Pulse rate frequency
15–25 kHz
Pulse length
20–50 ns
Input electrical power
1–2.5 kW
Average output power
0.2–1.3 W
Laser efficiency
<1%
Total service life
Up to 1000 hours
Active medium temperature
1000 K
Start time
15–20 min
Buffer gas Neon
7–30 Torr
Hydrogen admixture
0–0.06 Torr
Active zone length
60–100 cm
Inside tube diameter
4–26 mmFigure 1
Longitudinal principle scheme of the laser tube of a UV copper ion laser.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:D, mm: inside diameter of the laser tube; DR, mm: inside diameter of the internal rings in the ceramic tube or quartz tube insert; L, cm: length of the active zone (distance between the electrodes); PIN2, kW: electric power supplied to the active volume, taking into account 50% losses of the total supplied power PIN; PNE, Torr: neon gas pressure; PRF, kHz: pulse rate frequency; PH2, Torr: hydrogen gas pressure; C, 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)Pout, W, and laser service life LTime, hours.This study uses the data fromn=238 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.Figure 2
Scatter plot of all investigated data of UV Cu+ Ne-CuBr vapor lasers.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 [6]. 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 byX=(x1,x2,…,xp) and their dependent one by y=y(X). We assume that y,x1,x2,…,xp are vectors with dimension n. For our data n=238,   p=9. The general form of a MARS model y^[M] of the data for y has the form
(4.1)y^[M]=b0+∑j=1MbjBFj(X),
where b0, bj, j=1,2,…,M are the sought constant coefficients of the model, M the number of functions included in the model, and BFj(X) 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:(4.2)BFj(X)=max⁡(0,xk-ck)or   BFj(X)=max⁡(0,ck-xk).
Here ck is a constant, called a knot of the basis function with values in the definition interval of the predictor variable xk. 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 functionsBFj(X) can also include the products of two or more functions of the form (4.2) without repeating the index k. 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 M, 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 y(X).When building the models, the researcher sets an initial maximum numberM0 of basis functions and the maximum order of interaction. The recommendation is M0≥3p [21]. 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:
(4.3)SSE[Q]=∑i=1n(yi-y^[Q](Xi))2
for a model with a Q 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 M1 models with m basis functions is the one that minimizes the expression
(4.4)GCV[m]=∑i=1n(yi-y^[m](Xi))2n(1-C(m)/n)2,m=1,2,…,M1,
where C(m)=m+δ(m-1)/2, δ∈[2,3], Xi=(x1i,x2i,…,xpi), yi=y(x1i,x2i,…,xpi).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(M0,r) for a model with M0 initial maximum number of allowed basis functions (BFs) and up to rth order of interactions, r≥1.


## 5. Building MARS Models of the Laser Generation of an UV Laser

First, we will present the best MARS models of laser generationPout.The basic statistic figures of the constructed models are given in Table2. 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.Table 2
Main statistics of the constructed best MARS models of laser output powerPout  of UV copper ion lasers.
Model
R
2
MARS GCVR2
Sum of Squares
BF in model
Predictors
(30;1)
0.921
0.894
97.18
13
P
NE, PIN2, L, PH2, PRF, D, C, DR
(40;1)
0.921
0.892
97.18
13
P
NE, PIN2,
L,
PH2,
PRF,
D,
C, DR
(50;1)
0.922
0.890
96.83
14
P
N
E, PIN2, L, PH2, PRF, D, C, DR
(30;2)
0.968
0.940
63.37
23
P
IN
2, PNE, PH2, DR, L, PRF, C
(35;2)
0.977
0.951
55.80
27
P
IN
2, PNE, PH2, DR, L, PRF, C
(40;2)
0.980
0.955
50.34
28
P
IN2,
PNE,
PH2, DR,
L,
PRF,
C, D
(45;2)
0.982
0.959
47.67
28
P
IN
2, PNE, PH2, DR, L, PRF, C, D
(50;2)
0.984
0.957
45.50
31
P
IN
2, PNE, PH2, DR, L, PRF, C, D
(30;3)
0.978
0.957
52.82
23
P
IN
2, DR, PNE, L, PH2, TR, PRF, D, C
(35;3)
0.981
0.959
49.60
26
P
IN
2, DR, PNE, L, PH2, TR, PRF, D, C
(40;3)
0.984
0.965
44.92
27
P
IN2, DR,
PNE,
L, PH2,
PRF,
D,
C, TR
(45;3)
0.984
0.964
44.92
27
P
IN
2, DR, PNE, L, PH2, PRF, D, C, TR
(50;3)
0.985
0.965
44.42
27
P
IN
2, DR, PNE, L, PH2, PRF, D, C, TR


### 5.1. Linear MARS Models ofPout

Out of the linear models, we will present model (40;1). The figures of other linear models are given in the upper part of Table2.The generated best MARS model (40;1) includes the following 13 basis functions (BFs) with 8 predictors:(5.1)BF1=max⁡(0,DR-5.7),BF3=max⁡(0,PIN2-525),BF4=max⁡(0,PH2-0.03),BF5=max⁡(0,0.03-PH2),BF8=max⁡(0,PRF-19.5),BF9=max⁡(0,L-90),BF11=max⁡(0,C-372.184),BF12=max⁡(0,D-9.5),BF14=max⁡(0,L-86.5),BF16=max⁡(0,PNE-19.3),BF18=max⁡(0,PNE-18.6),BF20=max⁡(0,PNE-20),BF26=max⁡(0,D-7.1).The obtained regression model for laser output powerPout with these functions is
(5.2)P^out=403.83+49.71BF1+0.901119BF3-11325.2BF4-9840.28BF5+23.5177BF8-79.7689BF9+0.453381BF11+218.763BF12+42.4055BF14-1007.34BF16+514.971BF18+480.904BF20-215.734BF26.The model (5.1)-(5.2) can easily be used to estimate or predict laser generation. For example, let us take case i=128, with output laser power Pout128=900 mW [20] at the following values of the independent variables:
(5.3)X128=(D=7.1mm,DR=7.1mm,L=86.5cm,PIN=1500W,PNE=16.7Torr,PH2=0.03Torr,PRF=19.5kHz,C=735pF,Tr⁡=495°C).We calculate each of the functions in (5.1) consecutively. We have BF1=1.4, BF3=1500/2-525=225, BF4=BF5=BF8=BF9=0, and so forth. We substitute in model (5.2) obtaining the estimate
(5.4)P^out128=840mW.
The absolute error in case i=128 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 PH2, whose influence is measured as 100%, and the influence of the others is calculated against it.Table 3
Relative variable importance in the considered MARS model (40;1), described by (5.1)-(5.2), and model (40;2), described by (5.5)-(5.6), for laser output power  Pout  of UV copper ion laser.
Variable
In model (40;1), %
In model (40;2); %
P
IN
2
97.3%
100%
P
NE
100%
68.6%
P
H
2
66.8%
59.7%
DR
19.3%
40.4%
L
76.7%
35.4%
P
RF
55.2%
28.3%
C
30.9%
22.5%
D
53.0%
13.7%
TR
—
—The local behavior of the main predictors and their contribution to model (5.1)-(5.2) are shown in Figures 3(a)–3(h) in pure ordinal units.Graphs of the main predictors in model (5.1)-(5.2) showing their relationship with the response Pout.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)


### 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:(5.5)BF1=max⁡(0,DR-5.7),BF2=max⁡(0,5.7-DR),BF3=max⁡(0,PIN2-525),BF4=max⁡(0,PH2-0.03),BF5=max⁡(0,0.03-PH2),BF8=max⁡(0,PNE-19.3)BF3,BF12=max⁡(0,PNE-20)BF3,BF14=max⁡(0,PRF-20)BF3,BF15=max⁡(0,20-PRF)BF3,BF16=max⁡(0,L-90)BF3,BF17=max⁡(0,90-L)BF3,BF18=max⁡(0,PIN2-650)BF5,BF19=max⁡(0,650-PIN2)BF5,BF21=max⁡(0,21-PRF)BF5,BF22=max⁡(0,PNE-13.5)BF1,BF23=max⁡(0,13.5-PNE)BF1,BF24=max⁡(0,C-372.184)BF1,BF25=max⁡(0,L-90)BF1,BF27=max⁡(0,PIN2-670)BF1,BF28=max⁡(0,670-PIN2)BF1,BF29=max⁡(0,PNE-18.3)BF3,BF31=max⁡(0,PNE-19.75)BF3,BF33=max⁡(0,PNE-17.5),BF34=max⁡(0,17.5-PNE),BF35=max⁡(0,PNE-5.5)BF2,BF36=max⁡(0,D-4)BF2,BF37=max⁡(0,PNE-21.5)BF5,BF39=max⁡(0,PNE-20).The corresponding regression model for laser output powerPout includes all of these functions with 8 independent variables. It has the following form:
(5.6)P^out=687.002-5793.51BF2+1.48468BF3-10671BF4-15811.5BF5-3.04711BF8+0.957235BF10+4.28622BF12+0.164663BF14-1.38904BF15-0.110126BF16-0.0553432BF17-59.017BF18+245.219BF19+3593.7BF21-26.8772BF22+7.20284BF23+0.348349BF24+22.0488BF25+0.266576BF27-0.814258BF28+1.31689BF29-3.51103BF31-52.8455BF33-63.9875BF34+368.813BF35+537.251BF36-12825.2BF37-180.136BF39.Likewise, model (5.5)-(5.6) can be used to estimate and predict output power. For example, in the same case i=128, with the data (5.3), we obtain the estimate
(5.7)P^out128=856.4mW.The absolute error in casei=128 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 Pout are presented in Figures 3(a)–3(d). The biggest contribution of more than 6000 units is made by the interaction between PNE and PH2, which reaches the highest value in a large 2D region. Predictors C and DR provide almost the same contribution. The other two interactions also have given effects.


### 5.3. Application of the Models forPout

In (5.4) and (5.7), it was demonstrated how the predicted values of the model P^out are found. The properties of approximation of the first- and second-order models can be easily examined using graphs.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 P^out. Figure 3(a) indicates that the pressure of the buffer gas (neon) PNE, 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 4<D<10 mm (Figure 3(b)), which has been established experimentally in [16, 17]. The behavior of the pressure of hydrogen admixtures PH2 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 PIN 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 PNE 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).Graphs of contribution of the selected pairs of predictor variables in MARS model with (5.5)-(5.6) in ordinal units.
(a)
(b)
(c)
(d)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, Figure5(a) shows a slice from Figure 4(b) for the modeled behavior of the variable PNE (neon gas pressure) for a fixed value of the hydrogen pressure PH2=0.03 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 {PNE,PIN2} for a fixed PNE=19.5 Torr, depending on the increase of PIN2, as an element of Figure 4(c), which shows the maximum importance.View of slices from Figures4(b) and 4(c) for the influence of PNE and PIN2 on Pout in pure ordinal units.
(a)
(b)By carefully studying the resulting estimates ofP^out 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 ofPout for the laser from [18] when varying the pressure of the neon buffer gas PNE in the interval [15.5, 18.5] Torr. For PIN=1300W we choose experiment data with fixed D=DR=7.1mm, L=86.5cm, PH2=0.03Torr, PRF=19.5kHz, and C=735pF. In Figure 6(a) continuous line indicates the experiment data, Predicted shows estimates from model (40;2).Figure 6
Comparison of output powerPout with values predicted by model (40;2) at varying PNE for PIN = 1300 W.Results for a future experiment prediction using model (40; 2) are given in Table4. The following are fixed: D=DR=7.1 mm, L=86.5 cm, PNE=16.7 Torr, PH2=0.03 Torr, PRF=25 kHz, TR=490°C, and C=372 pF. As expected, laser generation Pout increases when input electric power PIN is increased. In this table, the middle row #3 corresponds to an actual experiment with the measured value of Pout=1200 mW, and the one predicted by the model is 1145 mW, which represents a relative error of under 5%.Table 4
Calculated values for predicted future experiments using MARS model (5.5)-(5.6).
P
IN
, W
P
out, mW-predicted value
#1
1500
1021
#2
1550
1083
#3
1600
1145
#4
1650
1208
#5
1700
1270


## 5.1. Linear MARS Models ofPout

Out of the linear models, we will present model (40;1). The figures of other linear models are given in the upper part of Table2.The generated best MARS model (40;1) includes the following 13 basis functions (BFs) with 8 predictors:(5.1)BF1=max⁡(0,DR-5.7),BF3=max⁡(0,PIN2-525),BF4=max⁡(0,PH2-0.03),BF5=max⁡(0,0.03-PH2),BF8=max⁡(0,PRF-19.5),BF9=max⁡(0,L-90),BF11=max⁡(0,C-372.184),BF12=max⁡(0,D-9.5),BF14=max⁡(0,L-86.5),BF16=max⁡(0,PNE-19.3),BF18=max⁡(0,PNE-18.6),BF20=max⁡(0,PNE-20),BF26=max⁡(0,D-7.1).The obtained regression model for laser output powerPout with these functions is
(5.2)P^out=403.83+49.71BF1+0.901119BF3-11325.2BF4-9840.28BF5+23.5177BF8-79.7689BF9+0.453381BF11+218.763BF12+42.4055BF14-1007.34BF16+514.971BF18+480.904BF20-215.734BF26.The model (5.1)-(5.2) can easily be used to estimate or predict laser generation. For example, let us take case i=128, with output laser power Pout128=900 mW [20] at the following values of the independent variables:
(5.3)X128=(D=7.1mm,DR=7.1mm,L=86.5cm,PIN=1500W,PNE=16.7Torr,PH2=0.03Torr,PRF=19.5kHz,C=735pF,Tr⁡=495°C).We calculate each of the functions in (5.1) consecutively. We have BF1=1.4, BF3=1500/2-525=225, BF4=BF5=BF8=BF9=0, and so forth. We substitute in model (5.2) obtaining the estimate
(5.4)P^out128=840mW.
The absolute error in case i=128 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 PH2, whose influence is measured as 100%, and the influence of the others is calculated against it.Table 3
Relative variable importance in the considered MARS model (40;1), described by (5.1)-(5.2), and model (40;2), described by (5.5)-(5.6), for laser output power  Pout  of UV copper ion laser.
Variable
In model (40;1), %
In model (40;2); %
P
IN
2
97.3%
100%
P
NE
100%
68.6%
P
H
2
66.8%
59.7%
DR
19.3%
40.4%
L
76.7%
35.4%
P
RF
55.2%
28.3%
C
30.9%
22.5%
D
53.0%
13.7%
TR
—
—The local behavior of the main predictors and their contribution to model (5.1)-(5.2) are shown in Figures 3(a)–3(h) in pure ordinal units.Graphs of the main predictors in model (5.1)-(5.2) showing their relationship with the response Pout.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)


## 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:(5.5)BF1=max⁡(0,DR-5.7),BF2=max⁡(0,5.7-DR),BF3=max⁡(0,PIN2-525),BF4=max⁡(0,PH2-0.03),BF5=max⁡(0,0.03-PH2),BF8=max⁡(0,PNE-19.3)BF3,BF12=max⁡(0,PNE-20)BF3,BF14=max⁡(0,PRF-20)BF3,BF15=max⁡(0,20-PRF)BF3,BF16=max⁡(0,L-90)BF3,BF17=max⁡(0,90-L)BF3,BF18=max⁡(0,PIN2-650)BF5,BF19=max⁡(0,650-PIN2)BF5,BF21=max⁡(0,21-PRF)BF5,BF22=max⁡(0,PNE-13.5)BF1,BF23=max⁡(0,13.5-PNE)BF1,BF24=max⁡(0,C-372.184)BF1,BF25=max⁡(0,L-90)BF1,BF27=max⁡(0,PIN2-670)BF1,BF28=max⁡(0,670-PIN2)BF1,BF29=max⁡(0,PNE-18.3)BF3,BF31=max⁡(0,PNE-19.75)BF3,BF33=max⁡(0,PNE-17.5),BF34=max⁡(0,17.5-PNE),BF35=max⁡(0,PNE-5.5)BF2,BF36=max⁡(0,D-4)BF2,BF37=max⁡(0,PNE-21.5)BF5,BF39=max⁡(0,PNE-20).The corresponding regression model for laser output powerPout includes all of these functions with 8 independent variables. It has the following form:
(5.6)P^out=687.002-5793.51BF2+1.48468BF3-10671BF4-15811.5BF5-3.04711BF8+0.957235BF10+4.28622BF12+0.164663BF14-1.38904BF15-0.110126BF16-0.0553432BF17-59.017BF18+245.219BF19+3593.7BF21-26.8772BF22+7.20284BF23+0.348349BF24+22.0488BF25+0.266576BF27-0.814258BF28+1.31689BF29-3.51103BF31-52.8455BF33-63.9875BF34+368.813BF35+537.251BF36-12825.2BF37-180.136BF39.Likewise, model (5.5)-(5.6) can be used to estimate and predict output power. For example, in the same case i=128, with the data (5.3), we obtain the estimate
(5.7)P^out128=856.4mW.The absolute error in casei=128 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 Pout are presented in Figures 3(a)–3(d). The biggest contribution of more than 6000 units is made by the interaction between PNE and PH2, which reaches the highest value in a large 2D region. Predictors C and DR provide almost the same contribution. The other two interactions also have given effects.


## 5.3. Application of the Models forPout

In (5.4) and (5.7), it was demonstrated how the predicted values of the model P^out are found. The properties of approximation of the first- and second-order models can be easily examined using graphs.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 P^out. Figure 3(a) indicates that the pressure of the buffer gas (neon) PNE, 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 4<D<10 mm (Figure 3(b)), which has been established experimentally in [16, 17]. The behavior of the pressure of hydrogen admixtures PH2 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 PIN 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 PNE 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).Graphs of contribution of the selected pairs of predictor variables in MARS model with (5.5)-(5.6) in ordinal units.
(a)
(b)
(c)
(d)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, Figure5(a) shows a slice from Figure 4(b) for the modeled behavior of the variable PNE (neon gas pressure) for a fixed value of the hydrogen pressure PH2=0.03 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 {PNE,PIN2} for a fixed PNE=19.5 Torr, depending on the increase of PIN2, as an element of Figure 4(c), which shows the maximum importance.View of slices from Figures4(b) and 4(c) for the influence of PNE and PIN2 on Pout in pure ordinal units.
(a)
(b)By carefully studying the resulting estimates ofP^out 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 ofPout for the laser from [18] when varying the pressure of the neon buffer gas PNE in the interval [15.5, 18.5] Torr. For PIN=1300W we choose experiment data with fixed D=DR=7.1mm, L=86.5cm, PH2=0.03Torr, PRF=19.5kHz, and C=735pF. In Figure 6(a) continuous line indicates the experiment data, Predicted shows estimates from model (40;2).Figure 6
Comparison of output powerPout with values predicted by model (40;2) at varying PNE for PIN = 1300 W.Results for a future experiment prediction using model (40; 2) are given in Table4. The following are fixed: D=DR=7.1 mm, L=86.5 cm, PNE=16.7 Torr, PH2=0.03 Torr, PRF=25 kHz, TR=490°C, and C=372 pF. As expected, laser generation Pout increases when input electric power PIN is increased. In this table, the middle row #3 corresponds to an actual experiment with the measured value of Pout=1200 mW, and the one predicted by the model is 1145 mW, which represents a relative error of under 5%.Table 4
Calculated values for predicted future experiments using MARS model (5.5)-(5.6).
P
IN
, W
P
out, mW-predicted value
#1
1500
1021
#2
1550
1083
#3
1600
1145
#4
1650
1208
#5
1700
1270


## 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 Section3 with response LTime—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 Table5. The first best MARS model with 98% coefficient of determination is model (40;2).Table 5
Main statistics of the constructed best MARS models of laser service lifeLTime of UV copper ion lasers.
Model
R
2
MARS GCVR2
Std. Err. of the estimate
BF  in model
Predictors
(30;1)
0.991
0.989
22.53
6
TR,
C,
PNE,
PIN2
(40;1)
0.991
0.989
22.55
6
TR,C, PNE, PIN2, L
(50;1)
0.991
0.989
22.53
6
TR,C, PNE, PIN2
(30;2)
0.997
0.995
13.90
13
TR,
PNE,
C,
PRF, DR,
PIN2
(35;2)
0.997
0.995
13.90
13
TR,PNE, C, PRF, DR, PIN2,
(40;2)
0.997
0.995
13.90
13
TR,PNE, C, PRF, DR, PIN2,
(30;3)
1.000
0.999
0.14
21
TR,C, PNE, DR, PH2
(35;3)
1.000
0.999
0.33
24
TR,C, PNE, DR, PH2
(40;3)
1.000
0.999
0.33
24
TR,C, PNE, DR, PH2


### 6.1. Linear MARS Models ofLTime

Of the linear-type models, we will present model (30;1). The figures of other linear models are given in the upper part of Table5.The calculated best MARS model (30;1) includes the following 6 basis functions with 4 predictors:(6.1)BF1=max⁡(0,TR-500),BF3=max⁡(0,C-372.184),BF5=max⁡(0,8.4-PNE),BF6=max⁡(0,TR-560),BF9=max⁡⁡(0,650-PIN2),BF13=max⁡(0,TR-550).The regression model for laser service lifeLTime with these functions is
(6.2)L^Time=693.836-12.8726BF1-1.47903BF3-397.699BF5+18.566BF6+3.7592BF9+9.86482BF13.Model (6.1)-(6.2) can be used to estimate and predict the service life of the examined lasers as shown by the models of Pout. The model is significant at level 0.000.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 TR, C, PNE, and PIN2 (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.Table 6
Relative variable importance in the considered MARS models (30;1), described by (6.1)-(6.2), and (30;2) with (6.3)-(6.4), for laser service life  LTime of UV copper ion laser.
Variable
Model (30;1)
Model (30;2)
TR
100%
100%
P
NE
20.1%
39.2%
C
46.7%
30.8%
P
RF
—
2.5%
DR
—
2%
P
IN
2
10.8%
1.6%Graphs of the main predictors in model (30;1) with (6.1)-(6.2) showing their relationship with the response LTime.
(a)
(b)
(c)
(d)


### 6.2. Nonlinear MARS Models ofLTime

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:(6.3)BF1=max⁡(0,TR-500),BF3=max⁡(0,C-372.184),BF4=max⁡(0,PNE-8.4),BF5=max⁡(0,8.4-PNE),BF6=max⁡(0,TR-550)BF4,BF7=max⁡(0,550-TR)BF4,BF8=max⁡(0,PNE-5.5)BF1,BF10=max⁡(0,560-TR),BF12=max⁡(0,7.3-DR)BF10,BF13=max⁡(0,PNE-16.7)BF10,BF15=max⁡(0,PNE-16.5)BF10,BF17=max⁡(0,PIN2-525)BF3,BF18=max⁡(0,TR-490)BF5,BF21=max⁡(0,PIN2-700)BF10,BF23=max⁡(0,PIN2-750)BF10,BF25=max⁡(0,PRF-19.5)BF1.
The corresponding regression model for laser service life LTime includes 13 of these functions and 6 predictors. It has the following form:
(6.4)L^Time=-103.913-2.59996BF3-22.9161BF4-2.00281BF6+2.35338BF7+0.71495BF8-4.50066BF12+8.48605BF13-9.19289BF15+0.0105753BF17+0.456027BF18-0.112225BF21+0.104861BF23+1.12861BF25.The absolute error of model (6.3)-(6.4) is about 0.006%, and the relative one is 4%.The importance of the two main predictors is shown in Figures8(a)–8(b).Graphs of the local mutual contribution of the main predictors to model (30;2).
(a)
(b)


### 6.3. Application of the Obtained MARS Models ofLTime

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 PNE and the bank of condenser within the established limits.


## 6.1. Linear MARS Models ofLTime

Of the linear-type models, we will present model (30;1). The figures of other linear models are given in the upper part of Table5.The calculated best MARS model (30;1) includes the following 6 basis functions with 4 predictors:(6.1)BF1=max⁡(0,TR-500),BF3=max⁡(0,C-372.184),BF5=max⁡(0,8.4-PNE),BF6=max⁡(0,TR-560),BF9=max⁡⁡(0,650-PIN2),BF13=max⁡(0,TR-550).The regression model for laser service lifeLTime with these functions is
(6.2)L^Time=693.836-12.8726BF1-1.47903BF3-397.699BF5+18.566BF6+3.7592BF9+9.86482BF13.Model (6.1)-(6.2) can be used to estimate and predict the service life of the examined lasers as shown by the models of Pout. The model is significant at level 0.000.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 TR, C, PNE, and PIN2 (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.Table 6
Relative variable importance in the considered MARS models (30;1), described by (6.1)-(6.2), and (30;2) with (6.3)-(6.4), for laser service life  LTime of UV copper ion laser.
Variable
Model (30;1)
Model (30;2)
TR
100%
100%
P
NE
20.1%
39.2%
C
46.7%
30.8%
P
RF
—
2.5%
DR
—
2%
P
IN
2
10.8%
1.6%Graphs of the main predictors in model (30;1) with (6.1)-(6.2) showing their relationship with the response LTime.
(a)
(b)
(c)
(d)


## 6.2. Nonlinear MARS Models ofLTime

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:(6.3)BF1=max⁡(0,TR-500),BF3=max⁡(0,C-372.184),BF4=max⁡(0,PNE-8.4),BF5=max⁡(0,8.4-PNE),BF6=max⁡(0,TR-550)BF4,BF7=max⁡(0,550-TR)BF4,BF8=max⁡(0,PNE-5.5)BF1,BF10=max⁡(0,560-TR),BF12=max⁡(0,7.3-DR)BF10,BF13=max⁡(0,PNE-16.7)BF10,BF15=max⁡(0,PNE-16.5)BF10,BF17=max⁡(0,PIN2-525)BF3,BF18=max⁡(0,TR-490)BF5,BF21=max⁡(0,PIN2-700)BF10,BF23=max⁡(0,PIN2-750)BF10,BF25=max⁡(0,PRF-19.5)BF1.
The corresponding regression model for laser service life LTime includes 13 of these functions and 6 predictors. It has the following form:
(6.4)L^Time=-103.913-2.59996BF3-22.9161BF4-2.00281BF6+2.35338BF7+0.71495BF8-4.50066BF12+8.48605BF13-9.19289BF15+0.0105753BF17+0.456027BF18-0.112225BF21+0.104861BF23+1.12861BF25.The absolute error of model (6.3)-(6.4) is about 0.006%, and the relative one is 4%.The importance of the two main predictors is shown in Figures8(a)–8(b).Graphs of the local mutual contribution of the main predictors to model (30;2).
(a)
(b)


## 6.3. Application of the Obtained MARS Models ofLTime

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 PNE 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 Tables2 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 ofR2 and GCV R2 are in linear models (r=1) both for Pout and for LTime. 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 Pout or LTime is not linear.For output powerPout, 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 R2 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 ofLTime, 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 levelP=0.00000. The coefficients of the models are significant with P≤0.005. The standardized error of the estimate is small. Model residuals are normally distributed.Figure9 shows the comparison of experimental data for laser output power Pout with those calculated by model (40;2). It is observed that the model fits data quite well. For the coefficient of determination we have R2=0.98, which means that the model accounts for 98% of all data. The corresponding GCVR2=0.95 (see Table 2). The residuals of this model are normally distributed with N(0,0.998).Figure 9
Values of the experimentalPout against the predicted P^out by model (40;2) with (5.5)-(5.6) with 5% confidential interval.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) forLTime against the experiment is given in Figure 10.Figure 10
Values of the experimentalLTime against the predicted L^Time by model (30;2) with (6.3)-(6.4) with 5% confidential interval.


## 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 toD—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 [20] 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 PNE. 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.


## 9. Conclusion

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.



---

*Source: 420782-2012-07-18.xml*