Hindawi Publishing Corporation Mathematical Problems in Engineering Volume 2012, Article ID 420782,21 pages doi:10.1155/2012/420782

Research Article

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

Snezhana Georgieva Gocheva-Ilieva,1

Desislava Stoyanova Voynikova,1 and Iliycho Petkov Iliev2

1 Department of Applied Mathematics and Modeling, University of Plovdiv,

24 Tzar Assen Street, 4000 Plovdiv, Bulgaria

2 Department of Physics, Technical University of Sofia, Branch Plovdiv,

25 Tzanko Djusstabanov Street, 4000 Plovdiv, Bulgaria

Correspondence should be addressed to

Snezhana Georgieva Gocheva-Ilieva, snegocheva@yahoo.com Received 2 March 2012; Accepted 12 May 2012 Academic Editor: Hung Nguyen-Xuan

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.

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.

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 [15]. 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 [69]. 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

Table 1: Technical characteristics of a UV copper ion vapor laser [20].

Characteristic

Values

Radiation wavelength Operating mode Pulse rate frequency Pulse length Input electrical power Average output power Laser efficiency Total service life Active medium temperature Start time Buffer gas Neon Hydrogen admixture Active zone length Inside tube diameter

248.6, 252.9, 259.1, 260.0, 270.3 nm Periodic pulse, self-heating

Up to 1000 hours

1000 K 15-20 min 7-30 Torr 0-0.06 Torr 60-100 cm 4-26 mm

15-25 kHz 20-50 ns 1-2.5 kW 0.2-1.3 W

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 Figure 1. Its main technical characteristics are given in Table 1.

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); P[N2, kW: electric power supplied to the active volume, taking into account 50% losses of the total supplied power P[N; 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.

Zr oxide fibrous

insulation Window protectors

Quartz window /

Quartz tube / Ceramic tube insert

t-1-y——■ ' ' '

iQ [TTinnnq

Porous / Active zone CuBr Quartz reservoir t

copper Heater

electrode

Figure 1: Longitudinal principle scheme of the laser tube of a UV copper ion laser.

o o o o o

This study uses the data from n = 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.

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 by X = (x\,x2,...,xp) and their dependent one by y = y(X). We assume that y,x\,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

y[M]= bo + £ bj BF; (X), (4.1)

where bo, 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:

BFj (X) = max(0, xk - ck) or BFj (X) = max(0, ck - xk). (4.2)

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.

u OD 0 § OOŒOO OC B O OŒ5HD § 03DIJ30 tllMMWDKCOi) «■»»¡■h — CD O

Q o | OGD O O 0 fe 0 o O sjSffiD <JE) Ï5K) sW (gnnnD O ^ 0 O O

DR O | OOD O O 1 OC o b 8 ^^JED o ^SBEBD ((wmwm>) lip »inumra k O

i-J O 8 OOD O o 0 0 oo cÄüSo o 0 OSSD Sis0 © aMfcffh o àannBD o CD SssoB 0 O 8 (immm o OD 8o O O

1 i ¿98 O £ o ÎLoo t o o o a h 8 o Urn 0 1 <auu w <M § 1 £r o P? §

1 (a® ^iii K h h ¿à vL h |o 0

1 l 000 0 u 0 O 0 L OD i

S PH 1 000 J Sol II p 8 0 8 S> o £ L Iuuuid r l l

R h § OOD w | o 1 a®§n>| |c elle >° o ■«S № I© pjjgESBO TT £ lo e

h i Ami 1 i U II M

o gOQO o 0 o atskmo o£ h 8 o ma o q2böed 8° o o (fmmmmhflft o ännäö

C D DR L P[N2 Pne Ph2 Prf TR Pout Ltime

Figure 2: Scatter plot of all investigated data of UV Cu+ Ne-CuBr vapor lasers.

When building nonlinear MARS models, basis functions BFy (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 number M0 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:

SSE[q]= jr y - y[q](Xt))2 (4.3)

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

Table 2: Main statistics of the constructed best MARS models of laser output power Pout of UV copper ion lasers.

Model R2 MARSGCV R2 Sum of Squares BF in model Predictors

(30;1) 0.921 0.894 97.18 13 Pne, Pin2, L, Ph2, Prf, D, C, DR

(40;1) 0.921 0.892 97.18 13 PNE, PiN2, L, PH2, PRF, D, C, DR

(50;1) 0.922 0.890 96.83 14 Pne, Pin2, L, Ph2, Prf, D, C, DR

(30;2) 0.968 0.940 63.37 23 Pin2, Pne, Ph2, DR, L, Prf, C

(35;2) 0.977 0.951 55.80 27 Pin2, Pne, Ph2, DR, L, Prf, C

(40;2) 0.980 0.955 50.34 28 PiN2, PNE, PH2, DR, L, PRF, C, D

(45;2) 0.982 0.959 47.67 28 Pin2, Pne, Ph2, DR, L, Prf, C, D

(50;2) 0.984 0.957 45.50 31 Pin2, Pne, Ph2, DR, L, Prf, C, D

(30;3) 0.978 0.957 52.82 23 Pin2, DR, Pne, L, Ph2, TR, Prf, D, C

(35;3) 0.981 0.959 49.60 26 Pin2, DR, Pne, L, Ph2, TR, Prf, D, C

(40;3) 0.984 0.965 44.92 27 PiN2, DR, Pne, L, PH2, PRF, D, C, TR

(45;3) 0.984 0.964 44.92 27 PiN2, DR, Pne, L, Ph2, Prf, D, C, TR

(50;3) 0.985 0.965 44.42 27 ^N2, DR, Pne, L, Ph2, Prf, D, C, TR

which "the best" model of M1 models with m basis functions is the one that minimizes the expression

GCV[m] = Ell}yi -,y[m](X;))2, m = 1,2.....Ml, (4.4)

n(1 - C(m)/n)

where C(m) = m + S(m - 1)/2, 6 e [2,3], X, = (x^x^,... ,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 Mo 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 generation Pout.

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 Pout

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:

BF1 = max(0, DR - 5.7), BF3 = max(0,F\N2 - 525), BF4 = max(0,Pm - 0.03), BF5 = max(0,0.03 - Pm), BF8 = max(0,PRF - 19.5), BF9 = max(0,L - 90), BF11 = max(0, C - 372.184), (5.1)

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 power Pout with these functions is Pout = 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:

X128 = (D = 7.1 mm, DR = 7.1 mm, L = 86.5 cm, PN = 1500 W, PNE = 16.7 Torr,

Ph2 = 0.03 Torr, Prf = 19.5 kHz, C = 735 pF, 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

Pout^ = 840 mW.

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 TH2, whose influence is measured as 100%, and the influence of the others is calculated against it.

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.

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:

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/fkE - 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 power Pout includes all of these functions with 8 independent variables. It has the following form:

Pout = 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 (5.6)

+ 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

Pout^ = 856.4mW. (5.7)

The absolute error in case i = 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 PNEE 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.

6 8 10 12 14 16 18 20 22 24 26 D (mm)

3 4 5 6 7 8 9 10 11 12 500 600 700 800 900 1000 DR (mm) PnN2 (W)

"i i i i i T

0.02 0.04 0.06 0.08

Ph2 (Torr) (e)

"1 I I I I I I I I I I T 19 20 21 22 23 24 25 26

Prf (kHz) (f)

Figure 3: Continued.

(g) (h) Figure 3: Graphs of the main predictors in model (5.1)-(5.2) showing their relationship with the response

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); %

PIN2 97.3% 100%

Pne 100% 68.6%

Ph2 66.8% 59.7%

DR 19.3% 40.4%

L 76.7% 35.4%

Prf 55.2% 28.3%

C 30.9% 22.5%

D 53.0% 13.7%

TR — —

5.3. Application of the Models for Pout

In (5.4) and (5.7), it was demonstrated how the predicted values of the model Pout 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 Pout. 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)).

Table 4: Calculated values for predicted future experiments using MARS model (5.5)-(5.6).

Pin , W

Pout, mW-predicted value

#1 #2 #3

1500 1550 1600 1650 1700

1021 1083 1145 1208 1270

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

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 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 {PNra/PW} for a fixed PNE = 19.5 Torr, depending on the increase of P[N2, as an element of Figure 4(c), which shows the maximum importance.

By carefully studying the resulting estimates of Pout 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 Pout 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 = 1300 W we choose experiment data with fixed D = DR = 7.1 mm, L = 86.5 cm, PH2 = 0.03 Torr, PRf = 19.5 kHz, and C = 735 pF. 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: D = DR = 7.1mm, L = 86.5 cm, PNE = 16.7Torr, PH2 = 0.03 Torr, PRf = 25 kHz, TR = 490°C, and C = 372 pF. As expected, laser generation Pout increases when input electric power P\N 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%.

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 LTime—laser service life, measured in hours. Again, our objective is to find a model with a sufficiently high

10000 8000

% 6000 4000 2000

10000 8000 6000 4000 2000

4000 2000

4000 2000

^ 1000 ^ 30 800

4000 3000

4000 3000 2000 1000

600 600

0.05 ^fo

200 200

Figure 4: Graphs of contribution of the selected pairs of predictor variables in MARS model with (5.5)-(5.6) in ordinal units.

Slice PH2 = 0.0305882: pure ordinal

Pne (Torr) (a)

Slice Pne = 19.5294: pure ordinal

~1--1-1-r-1-r

0 10 20 30 40

-1-1-1-1-1-1-r

500 600 700 800 900 1000

PIN2 (W) (b)

Figure 5: View of slices from Figures 4(b) and 4(c) for the influence of fNE and pN2 on Pout in pure ordinal

15 15.5 16 16.5 17 17.5 18 18.5 19

íNe (Torr)

-♦— Experiment Prediction

Figure 6: Comparison of output power Pout with values predicted by model (40;2) at varying PNE for PIN = 1300 W.

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 LTime

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:

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 - P¡N2), BF13 = max(0, TR - 550).

The regression model for laser service life LTime with these functions is

LTime = 693.836 - 12.8726BF1 - 1.47903BF3 - 397.699BF5 + 18.566BF6 + 3.7592BF9 + 9.86482BF13.

Table 5: Main statistics of the constructed best MARS models of laser service life LTime of UV copper ion lasers.

Model R2 MARS GCV R2 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, Pm

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 Pn2 (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 LTime

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:

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,

(a) (b)

(c) (d)

Figure 7: Graphs of the main predictors in model (30;1) with (6.1)-(6.2) showing their relationship with the response LTime.

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 LT;me includes 13 of these functions and 6 predictors. It has the following form:

LTime = -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 Figures 8(a)-8(b).

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%

Pne 20.1% 39.2%

C 46.7% 30.8%

PRF — 2.5%

DR — 2%

P[N2 10.8% 1.6%

800 600 ; 400

% 200 O

O 500'

^ 1000 800

700 Crpp,

400 ¿j |

400 300 200

300 200 100 0

600 500 400

Figure 8: Graphs of the local mutual contribution of the main predictors to model (30;2).

6.3. Application of the Obtained MARS Models of LTime

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 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 R2 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 1015%. This means that the relationship of the input independent variables with Pout or Liime is not linear.

For output power Pout, 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 of LTime, 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 P = 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.

Figure 9 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 GCV R2 = 0.95 (see Table 2). The residuals of this model are normally distributed with N(0, 0.998).

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 LT¡me 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 [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

1200 -

1000 -

400 - O (GTfW

o/Ms®

0 R2 linear = 0.98

-200 -

0 250 500 750 1000 1250

Figure 9: Values of the experimental Pout against the predicted 5.5 5.6

5% confidential interval.

Figure 10: Values of the experimental LTime against the predicted Lnme by model (30;2) with (6.3)-(6.4) with 5% confidential interval.

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.

Acknowledgments

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.

References

[1] 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.

[2] "Numerical modeling of low-temperature plasmas," in Encyclopedia of Low-temperature Plasma, Series B, M. Ianus, Ed., vol. 7, Moscow, Russia, 2004.

[3] 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.

[4] 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.

[5] 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.

[6] S. G. Gocheva-Ilieva and I. P. Iliev, Statistical Models of Characteristics of Metal Vapor Lasers, Nova Science Publishers, New York, NY, USA, 2011.

[7] 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.

[8] 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.

[9] 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.

[10] 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.

[11] N. Vuchkov, N. Sabotinov, and K. Temelkov, "Ultraviolet Copper laser," Patent Reg. no. 103312, 06.04.1999, Published Applications for Inventions, BG 10/2000.

[12] 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.

[13] 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.

[14] 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.

[15] 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.

[16] 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.

[17] 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.

[18] 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.

[19] 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.

[20] 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.

[21] D. Steinberg, B. Bernstein, P. Colla, and K. Martin, MARS User Guide, Salford Systems, San Diego, Calif, USA, 2001.

[22] 2012, http://www.salfordsystems.com/mars.php.

[23] 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.

[24] 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.

[25] 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.

[26] J. H. Friedman, "Multivariate adaptive regression splines," The Annals of Statistics, vol. 19, no. 1, pp. 1-141, 1991.

Copyright of Mathematical Problems in Engineering is the property of Hindawi Publishing Corporation and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.