USF Libraries
USF Digital Collections

Computer simulation of stand-off LIBS and Raman LIDAR for remote sensing of distant compounds

MISSING IMAGE

Material Information

Title:
Computer simulation of stand-off LIBS and Raman LIDAR for remote sensing of distant compounds
Physical Description:
Book
Language:
English
Creator:
Pliutau, Dzianis
Publisher:
University of South Florida
Place of Publication:
Tampa, Fla.
Publication Date:

Subjects

Subjects / Keywords:
Modified lidar equation
Atmospheric attenuation
Laser-induced breakdown spectroscopy
Raman spectroscopy
Dissertations, Academic -- Physics -- Masters -- USF   ( lcsh )
Genre:
bibliography   ( marcgt )
theses   ( marcgt )
non-fiction   ( marcgt )

Notes

Abstract:
Long range stand-off Raman and Laser-Induced Breakdown Spectroscopy (LIBS) lidar signal simulations have been carried out using a modified UV-visible atmospheric transmission program and a modified lidar equation. The Hitran-PC atmospheric transmission program which normally operates over the wavelength range of 400 nm to the far-IR was modified to provide UV atmospheric attenuation (200 nm -- 400 nm) using the optical cross section data contained in the HITRAN database. The two-way lidar equation was modified in order to simulate the one-way propagation of the Raman and LIBS spectral, and thus provide calculations of the expected Raman or LIBS signal as a function of range. Estimation of the LIBS and Raman spectral intensity was then calculated for several remote sensing cases. In particular, the atmospheric attenuation spectra generated with the modified Hitran-PC program were combined with the calculated LIBS and Raman lidar emission spectra at the remote excitation site using a modified lidar equation to determine for the first time to our knowledge the power and S/N ratio versus range of the LIBS and Raman Lidar complete spectrum as a function of wavelength in the UV -- IR region. Previous simulations had only made S/N versus range calculations at a single wavelength or for the total integrated emission. These results are important as they can be used for future design of stand-off LIBS and Raman lidar systems, and for comparisons with experimental measurements. In particular, we are planning to use our simulations for comparison of 266 nm excited LIBS and Raman lidar measurements of energetic compounds at ranges of a few tens of meters.
Thesis:
Thesis (M.S.)--University of South Florida, 2007.
Bibliography:
Includes bibliographical references.
System Details:
System requirements: World Wide Web browser and PDF reader.
System Details:
Mode of access: World Wide Web.
Statement of Responsibility:
Dzianis (Denis) Pliutau.
General Note:
Title from PDF of title page.
General Note:
Document formatted into pages; contains 76 pages.

Record Information

Source Institution:
University of South Florida Library
Holding Location:
University of South Florida
Rights Management:
All applicable rights reserved by the source institution and holding location.
Resource Identifier:
aleph - 001959234
oclc - 248647454
usfldc doi - E14-SFE0002208
usfldc handle - e14.2208
System ID:
SFS0026526:00001


This item is only available as the following downloads:


Full Text
xml version 1.0 encoding UTF-8 standalone no
record xmlns http:www.loc.govMARC21slim xmlns:xsi http:www.w3.org2001XMLSchema-instance xsi:schemaLocation http:www.loc.govstandardsmarcxmlschemaMARC21slim.xsd
leader nam Ka
controlfield tag 001 001959234
003 fts
005 20080912141434.0
006 m||||e|||d||||||||
007 cr mnu|||uuuuu
008 080912s2007 flu sbm 000 0 eng d
datafield ind1 8 ind2 024
subfield code a E14-SFE0002208
035
(OCoLC)248647454
040
FHM
c FHM
049
FHMM
090
QC21.2 (ONLINE)
1 100
Pliutau, Dzianis.
0 245
Computer simulation of stand-off LIBS and Raman LIDAR for remote sensing of distant compounds
h [electronic resource] /
Dzianis (Denis) Pliutau.
260
[Tampa, Fla.] :
b University of South Florida,
2007.
3 520
Long range stand-off Raman and Laser-Induced Breakdown Spectroscopy (LIBS) lidar signal simulations have been carried out using a modified UV-visible atmospheric transmission program and a modified lidar equation. The Hitran-PC atmospheric transmission program which normally operates over the wavelength range of 400 nm to the far-IR was modified to provide UV atmospheric attenuation (200 nm -- 400 nm) using the optical cross section data contained in the HITRAN database. The two-way lidar equation was modified in order to simulate the one-way propagation of the Raman and LIBS spectral, and thus provide calculations of the expected Raman or LIBS signal as a function of range. Estimation of the LIBS and Raman spectral intensity was then calculated for several remote sensing cases. In particular, the atmospheric attenuation spectra generated with the modified Hitran-PC program were combined with the calculated LIBS and Raman lidar emission spectra at the remote excitation site using a modified lidar equation to determine for the first time to our knowledge the power and S/N ratio versus range of the LIBS and Raman Lidar complete spectrum as a function of wavelength in the UV -- IR region. Previous simulations had only made S/N versus range calculations at a single wavelength or for the total integrated emission. These results are important as they can be used for future design of stand-off LIBS and Raman lidar systems, and for comparisons with experimental measurements. In particular, we are planning to use our simulations for comparison of 266 nm excited LIBS and Raman lidar measurements of energetic compounds at ranges of a few tens of meters.
502
Thesis (M.S.)--University of South Florida, 2007.
504
Includes bibliographical references.
516
Text (Electronic thesis) in PDF format.
538
System requirements: World Wide Web browser and PDF reader.
Mode of access: World Wide Web.
500
Title from PDF of title page.
Document formatted into pages; contains 76 pages.
590
Adviser: Dennis K. Killinger, Ph.D.
653
Modified lidar equation.
Atmospheric attenuation.
Laser-induced breakdown spectroscopy.
Raman spectroscopy.
690
Dissertations, Academic
z USF
x Physics
Masters.
773
t USF Electronic Theses and Dissertations.
4 856
u http://digital.lib.usf.edu/?e14.2208



PAGE 1

Computer Simulation of Sta nd-Off LIBS and Raman Lidar for Remote Sensing of Distant Compounds by Dzianis (Denis) Pliutau A thesis submitted in partial fulfillment of the requirements for the degree of Maser of Science Department of Physics College of Arts and Sciences University of South Florida Major Professor: Dennis K. Killinger, Ph.D. Nicholas Djeu, Ph.D. Myung K. Kim, Ph.D. Date of Approval: November 14, 2007 Keywords: modified lidar equation, atmosphe ric attenuation, lase r-induced breakdown spectroscopy, Raman spectroscopy Copyright 2007, Dzianis (Denis) Pliutau

PAGE 2

DEDICATION Dedicated to everyone who will ever open this thesis.

PAGE 3

ACKNOWLEDGEMENTS I would like to thank Dr. Denni s Killinger for being a respons ible advisor as well as for the opportunity to work on HITRAN-PC. I also thank Dr. Nicholas Djeu and Dr. Myung K. Kim for serving in my committee.

PAGE 4

i TABLE OF CONTENTS LIST OF TABLES iii LIST OF FIGURES iv ABSTRACT vii CHAPTER 1. INTRODUCTION 1 CHAPTER 2. INTRODUCTION TO LASER-INDUCED BREAKDOWN SPECTROSCOPY 5 2.1 Point LIBS 5 2.2 Laser sources used in LIBS 7 2.3 Remote LIBS 9 2.3.1 Remote LIBS using picoand nanosecond pulses 12 2.3.2 Conventional remote LIBS using femtosecond pulses 12 2.3.3 Filament-induced breakdown spectroscopy 14 CHAPTER 3. INTRODUCTION TO RAMAN SPECTROSCOPY 22 3.1 Point Raman Spectroscopy 22 3.2 Remote Raman Spectroscopy 24 3.3 Joint remote LIBS and Raman lidar measurements 30 CHAPTER 4. NEW ATMOSPHERIC TRANSMISSION COMPUTER CODES FOR UV-VISIBLE-IR WAVELENGTHS USEFUL FOR LIBS AND RAMAN LIDAR 32 4.1 Current atmospheric transmission simulations 32

PAGE 5

ii 4.2 Modification of Hitran-PC softwa re suitable for remote LIBS and Raman signal attenuation modeling 34 CHAPTER 5. MODIFICATION OF LI DAR EQUATION SUITABLE FOR LIBS AND RAMAN LIDAR 40 5.1 Current lidar equations 40 5.2 Modified LIDAR equation suitable for LIBS and Raman LIDAR 43 CHAPTER 6. ESTIMATED LIBS AND RAMAN LIDAR EMISSION 45 6.1 Estimated LIBS power spect rum at the excita tion site 45 6.2 Estimated Raman power spect rum at the excitation site 48 CHAPTER 7. LIBS AND RAM AN LIDAR SIMULATIONS 54 7.1 LIBS Lidar simulations 54 7.2 Raman Lidar Simulations 56 7.3 Comparison of Initial Simulations with Selected LIBS and Raman LIDAR studies 59 CHAPTER 8. SUMMARY AND FU TURE RESEARCH 62 LIST OF REFERENCES 63 APPENDIX A. MODIFICATION OF HITRAN-PC PROGRAM TO ADD UV CROSS SECTION ABSORBTION 68 APPENDIX B. PROCEDURE TO RESCALE THE DIGITIZED SPECTRUM TO THE RESOLU TION OF THE ATMOSPHERIC ATTENUATION SPECTRUM 72 APPENDIX C. MATLAB PROGRAM TO LINK THE ATMOSPHERIC ATTENUATION DATA, POWER SPECTRUM DATA AND PLOT 3D S/N VS. RANGE VS. WAVELENGTH FOR LIBS/RAMAN LIDAR 75

PAGE 6

iii LIST OF TABLES Table 2.1 Characteristic Raman shifts for selected gaseous, liquid and solid molecules. 25

PAGE 7

iv LIST OF FIGURES Figure 2.1 Diagram and a photo of a typical point LIBS set-up. (Adapted from D. Cremers, L. Radziemski, Handbook of Laser-Induced breakdown spectroscopy John Wiley & Sons, Ltd, 2006) .1 8 Figure 2.2 Schematic of LIBS lidar system. 10 Figure 2.3 LIBS spectra from stainl ess steel sample (liquid and solid at 1420C) obtained at a range of 7.5 m. Laser used: 220mJ per pulse, 1064nm, 5ns (From Palanco7). 11 Figure 2.4 R-LIBS spectra for a copper sample at 25m: (a) using 75 fs pulses, (b) using 200 ps pulses. Wavelengt hs of observed copper lines are indicated in (a), positions of observed atomic oxygen and nitrogen lines are indicated in (b) (From Rohwetter12). 13 Figure 2.5 R-FIBS spectrum taken for aluminum sample located 50m away. Laser used: 108mJ, 10fs, 800nm (From Liu3). 15 Figure 2.6 R-FIBS spectra of copper (a) and steel (b), measured at 90 m distance. Laser: 250mJ, 80fs, 800 nm (From Stelmaszczyk14). 16 Figure 2.7 Distance dependence of the range-corrected R-FIBS signal from the 521.8 nm line of copper (From Stemaszczyk14). 18 Figure 2.8 Extrapolation of the single-shot detecti on limit for the integral intensity of the detected LIBS signal (From Xu16). 19 Figure 2.9 Signal intensity extrapol ation of Al I: 396.15 nm line as a function of distance (From Liu3). 20 Figure 3.1 Schematic of a point Raman spectroscopy system. 23 Figure 3.2 Schematic of a typical remote Raman system. (Adapted from Chen5) 26 Figure 3.3 Portions of remote Raman spectra of ethyl be nzene and toluene in the spectral range 180 – 1800 cm-1 at 10m distance excited with a

PAGE 8

v 532nm, 35mJ, 8ns laser (From Sharma4). 28 Figure 3.4 A portion of remote Raman spectra naphthalene sample at 122 m distance measured with the direct ly coupled coaxial pulsed Raman system. Laser 532 nm, 25 mJ per pulse at the sample, 20 Hz, slit 100 m (From Sharma20). 29 Figure 3.5 (a) LIBS and (b) remote Raman spectra of a chalk sample, Raman spectrum of a calcite spectrum for reference. Laser: 532 nm, 35 mJ, 8ns (From Weins24). 31 Figure 4.1 File structure of the HITRAN database showing the large line-by-line molecular database and the smaller UV cross section database (Adapted from Rothman25). 33 Figure 4.2 Transmission spectrum of the atmosphere generated by HITRAN-PC program: US. Standard atmosphere, 2000 m path. 35 Figure 4.3 Composite transmission spect rum of the atmosphere in the UV region calculated from the HITRAN-2004 database UV cross-sections data for typical urban concentrations of BrO, N2O, NO2, NO3, O3 using the modified HITRAN-PC program (a – without Rayleigh scattering, b – with Rayleigh scattering). Path: 2000 m. Partial pressures: BrO (1-2ppt), N2O (0.3 ppm), NO2 (0.8 ppm), NO3 (0.35 ppb), O3 (0.5 ppm) 37 Figure 4.4 Combined spectra of Fig 4.2 and 4.3 for a 250 – 850 nm spectral range (2 km path). 38 Figure 4.5 Combined UV cross se ction and line-by-line spectra obtained with the Hitran-PC program for a 250 nm – 2.3 m region (2 km path). 39 Figure 5.1 Schematic of a typical lidar or differential-absorption lidar system (From Killinger28). 41 Figure 5.2 Pictorial schematic of lidar equation geometry showing backscatter of the lidar signa l and collection by a telescope. (From Killinger28). 42 Figure 6.1 Estimated power spectrum of a LIBS signal at the excitation site, Ps( ), for a peak laser pulse power of 1 kW and 1% LIBS efficiency calculated from the LIBS spectrum shown in Fig 2.5. 47 Figure 6.2 Estimated power spectrum of a Raman signal at the excitation site,

PAGE 9

vi Ps( ), for a peak laser power of 25 MW and 10-10 Raman efficiency, and laser excitation wavelength of 266 nm. Raman spectrum was calculated from Raman spectrum shown in Fig 3.4 (Naphthalene), obtained at 532 nm, no 1/ 4 correction was used. 49 Figure 6.3 Estimated Raman emission lin es for different laser excitation wavelengths and several different gas molecules located at the distant hard target (with -4 dependence taken into account). 51 Figure 6.4 Estimated Raman emission for the spectrum in Fig 6.2 (remote hard target of naphthalene ; 266 nm excitation) for different laser excitation at 266 nm, 532 nm, 1064 nm, and 1.54 m (without -4 dependence). 52 Figure 6.5 Composite plot of Raman emission for four different laser excitation wavelengths (Fig 6.4) and transmission spectrum of the atmosphere for a 2 km path (Fig. 4.5). 53 Figure 7.1 1 Calculated power spectr um of LIBS as a function of range (1 – 25 m to the target) and wave length for a peak laser pulse power of 1 KW, 1% LIBS efficiency, A = 0.01 m2 telescope area, K = 0.01 lidar efficiency. (Based upon Fig. 2.5 and power spectrum in Fig 6.1). 55 Figure 7.2 Calculated S/N spectrum of a LIBS signal as a function of range (10 – 15 km to the targ et) and wavelength for a peak laser pulse power of 1 KW and 1% LIBS efficiency, A = 0.01 m2 telescope area, K = 0.01 lidar e fficiency. (Based upon Fig. 2.5 and power spectrum in Fig 6.1). 57 Figure 7.3 Calculated power spectru m of a Raman signal as a function of range (1 – 25 m to the targ et) and wavelength for a peak laser pulse power of 25 MW and 10-10 Raman efficiency, telescope area A = 0.01 m2, lidar efficiency K = 0.01. (Based upon spectrum in Fig. 3.4 and power spectrum in 6.2 using 266 nm excitation). 58 Figure 7.4 Calculated S/N spectrum of a Raman signal as a function of range (100 – 300 m to the target ) and wavelength for a peak laser pulse power of 25 MW and 10-10 Raman efficiency, A = 0.01 m2 telescope area, K = 0.01 lidar efficiency. 60

PAGE 10

vii COMPUTER SIMULATION OF STANDOFF LIBS AND RAMAN LIDAR FOR REMOTE SENSING OF DISTANT COMPOUNDS Denis Pliutau ABSTRACT Long range stand-off Raman and Lase r-Induced Breakdown Spectroscopy (LIBS) lidar signal simulations have been carried out using a modified UV-visible atmospheric transmission program and a modified lidar equation. The Hitran-PC atmospheric transmission program which normally operate s over the wavelength range of 400 nm to the far-IR was modified to provide UV atmo spheric attenuation (200 nm – 400 nm) using the optical cross section data contained in the HITRAN database. The two-way lidar equation was modified in order to simulate the one-way propagation of the Raman and LIBS spectral, and thus provide calculations of the expected Raman or LIBS signal as a function of range. Estimation of the LIBS and Raman spectral intensity was then calculated for several remote sensing cases. In partic ular, the atmospheric attenuation spectra generated with the modified H itran-PC program were combined with the calculated LIBS and Raman lidar emission spectra at the remo te excitation site using a modified lidar equation to determine for the first time to our knowledge the power and S/N ratio versus

PAGE 11

viii range of the LIBS and Raman Lidar complete spectrum as a function of wavelength in the UV – IR region. Previous simulations ha d only made S/N versus range calculations at a single wavelength or for the total integrated emission. These results are important as they can be used for future design of stand-off LIBS and Raman lidar systems, and for compar isons with experimental measurements. In particular, we are planning to use our si mulations for comparison of 266 nm excited LIBS and Raman lidar measurements of energe tic compounds at ranges of a few tens of meters.

PAGE 12

1 CHAPTER 1. INTRODUCTION Laser-Induced-Breakdown Spectroscopy (LIBS)1-3 and Raman spectroscopy4-6 are optical spectroscopic techni ques that have been used for decades for the spectral identification and detection of a wide range of chemical species. LIBS uses a focused laser beam to produce a high-opt ical field induced dielectric breakdown at a target site, and the detection of the atomic or ionic emission from the excited species within the plasma. Raman spectroscopy often uses a hi gh power laser to exci te the weak Raman excitation emission lines from material, where the Raman emission is shifted in wavelength by a unique value that can be used to identify the material. Both of these techniques have been used recently for the remote or stand-off de tection of different materials.3-24 For example, stand-off LIBS and Ra man lidar systems are being considered for the remote detection of energetic materials. Because of these applications, it is impor tant to be able to compute the system Signal-to-Noise (S/N) value as a function of detection range for these kinds of LIBS and Raman lidar systems. Unfortunately, in the pa st most of these simulations have involved detailed studies of the LIBS or Raman spect ra at the emission point, but very little analysis of the propagation of the spectral si gnature through the atmosphere or the lidar detection of the signals has been done. In some past simulations, the LIBS or Raman lidar signal was calculated only at a single emission wavelength and often of either 1064 nm or 535 nm excitation.3, 16 Since the LIBS and Raman spectra can be complex

PAGE 13

2 consisting of tens of different emission lines it is important to c onsider all of these spectral lines in any remote sensing system sin ce often the ratio of su ch lines is used to determine the concentration of the remotely sensed material. Toward this end, we have studied the pr opagation and detection of a remotely emitted LIBS or Raman signal using modifi ed atmospheric transmission codes and databases, and use of a modified lidar e quation suitable for the LIBS and Raman standoff case. Because recent LIBS work has star ted to use deep-UV excitation, we have also expanded these studies to be useful in the UVvisible-mid-IR spectral range. Previous work was more directed toward the visibl e and near-IR. Our studies resulted in combining two different atmospheric tran smission simulations using the HITRAN line spectra database and the HITR AN UV cross section database.25 The LIBS or Raman lidar system was modeled using a remotely pr oduced LIBS plasma spark or a remotely generated Raman emission target area, and then back propagating the LIBS or Raman spectra toward the remote lidar telescope a nd detector system. Atmospheric absorption and the lidar 1/R2 decrease in the signal was taken into account. The resultant signal was compared to the Noise-Equivalent-Power (NEP) of the lidar detector, and the S/N as a function of range and as a function of wavelength was calculated. Several examples of Raman and LIBS emi ssion spectra were used to simulate the lidar system. For example, the LIBS spectra from an aluminum target had lines from 250 nm to 700 nm, and after propagating throug h the atmosphere showed a S/N value of about 10 out to ranges of near 10000 m. Th e Raman lidar analysis studied was more complex in that different excitation lase r wavelengths were studied including 1500nm, 1064 nm, 535 nm, and 266 nm. The Raman em ission for each of th ese excitations was

PAGE 14

3 calculated, and the resultant Raman em ission spectrum propagated through the atmosphere and lidar detection system. Our results show that different excitation wavelengths will produce different Raman lidar S/N values. These simulation models and computer codes will be useful in future application of LIBS and Raman lidar techniques to a wi de variety of different remote sensing applications. The results of our study are presented in the following chapters. Chapter 2 provides background information on point LI BS spectroscopy, and summarizes the current state and practical ranges of remote LI BS technique. Chapter 3 presents an overview of the use of point Raman spectrosc opy for chemical detection, as well as the current state of remote Raman measurement te chniques with particular attention paid to the detection of liquid and so lid samples. In Chapter 4 an overview of current atmospheric transmission simulation software is presented. In part icular, modifications of the Hitran-PC software package which us es the HITRAN line spectra database are described in order to simulate LIBS and Ra man signal attenuation by the atmosphere. Chapter 5 describes the standard LIDAR equa tion and covers the modifications to the LIDAR equation to model the transmission and detection of the remote LIBS and Raman LIDAR signals. In Chapter 6 estimates are made of the LIBS and Raman power spectra at the excitation site for eventual use with the modified LIDAR equation. In Chapter 7 the modified LIDAR equation is used in combination with the modified HITRAN generated atmospheric attenuation coefficients and a typical detector NEP in order to obtain the estimated power of LIBS and Rama n spectra at different wavelengths as a

PAGE 15

4 function of range. Several examples are presented based upon published LIBS emission spectra and realistic parameters. Finally, Ch apter 8 describes future work in this area including conducting remote se nsing LIBS and Raman lidar experiments and comparing these results to our simulations presented in this thesis.

PAGE 16

5 CHAPTER 2. INTRODUCTION TO LASER-INDUCED BREAKDOWN SPECTROSCOPY In this chapter an overview of point a nd remote LIBS techniques are presented. The main components of typical point and remo te LIBS setups are shown, and the current state of remote LIBS and vari ations of the method are summ arized. In addition, current studies to predict remote LIBS detection ranges are presented. 2.1 Point LIBS Atomic emission spectroscopy (AES) is a method of quantitative elemental analysis where the sample is evaporated, atomized and the spectrum of the excited atomized elements is recorded to determine the chemical elements present in the sample. Several excitation sources are typically used in atomic emission spectroscopy to form the plasma and influence the description of the techniques; these include InductivelyCoupled Plasma (ICP), Microwave-Induced Pl asma (MIP), and spark and laser-induced excitation.1 Laser-induced breakdown spectroscopy (LIB S), also called lase r-induced plasma spectroscopy (LIPS), is a met hod of atomic emission spectr oscopy (AES) utilizing laser radiation for vaporization, atomization and ex citation of a sample. Out of the four typical excitation sources mentioned above laser-induced breakdown is the only method that can be used for remote sensing at a distance.

PAGE 17

6 There are basically two ways to use LIBS for analysis One is point sensing, often done in the lab and at close ranges. The other is remote LIBS, conducted at ranges of 1 to 1000 m. Both of these are covered in the following. Points LIBS technique utilizes a power ful laser pulse which is focused at the sample surface in order to evaporate and atom ize atoms in the area where the laser pulse is focused. The spectrum produced by the vapor ized excited atoms and ions is gathered by the collection optics and sent to a spectrometer and a computer for further processing. Quantitative LIBS analysis is possible because the spectral wavelengths of the characteristic atomic and ionic emission lines for the individual elements are fixed and known. The intensities of these lines are us ed in order to determine the quantity of a particular element in the sample. However, quantitative analysis is complicated because the intensity of a particular line depends both on the quantity of the element being analyzed as well as on the general composition of the sa mple. This phenomenon is referred to as the chemical matrix effect. 1 Therefore accurate quantitative measurements are only possible if the standards used for the calibration of a LIBS system have a chemical composition close to the one being analyzed. This requirement is easily achievable for routine LIBS analysis at i ndustrial factories, fo r example for quality control of the manufacturing of steel. Unfortunately remote LIBS technique sometimes involves limited knowledge about the sample be ing analyzed, theref ore the quantitative accuracy of remote LIBS may be more li mited as compared to the routine LIBS applications in manufacturing environments. A typical point LIBS setup and a diagram are shown in Fig 2.1. The main components of most LIBS systems are a pulse d laser used to form the plasma, optical

PAGE 18

7 focusing system that focuses and directs the laser pulse to the samp le, light collecting system (lens, fiber optics) that collects the radiation from the microplasma and transports it to the detection system, delay generator to control the time between sending the laser pulse and the LIBS signal det ection, a spectrometer to determine intensities of the detected radiation at various wavelengths, a nd a computer and electronics to process the spectrum.1 2.2 Laser sources used in LIBS Typical laser pulse energies used in LIBS are 10 to 500 mJ with typical pulse widths of 5 to 20ns. Powerful laser pulses on the order of at leas t 5 MW are needed to form the plasma.1,2 The lasers most commonly used for LIBS are a Nd:YAG laser (1064 nm, 532 nm, 355 nm, 266 nm), XeCl (308 nm), KrF (248 nm) and ArF (194 nm) eximer lasers, and CO2 laser (10.6 m).1 Nd:YAG lasers generating nanosecond pulses are preferred for LIBS applications because these lasers are t echnologically well developed, compact, reliable, and provide high power pulses required for effective LI BS implementation. Injection seeded Nd:YAG lasers have also been inves tigated for the application to LIBS.1 The use of injection seeding reduces the linewidth of th e laser pulse, improves st ability and temporal profile. Even though nanosecond lasers are most commonly used in LIBS, picosecond and femtosecond lasers are also being evaluated. In particular, the use of femtosecond pulses

PAGE 19

8 Figure 2.1 Diagram and a photo of a typical point LIBS set-up (Adapted from D. Cremers, L. Radziemski, Handbook of LaserInduced breakdown spectroscopy John Wiley & Sons, Ltd, 2006).1

PAGE 20

9 provides a number of definite advantages to be outlined further in this thesis. The application of picoand femtosecond lasers to LIBS is currently limited due to the complexity of the equipment necessary to generate such short pulses. It is important to point out that some la ser wavelengths couple more readily into a specific material being analyzed compared to other laser wavelengths. Therefore, the optimum laser used for LIBS depends on the application and the la ser wavelength. 2.3 Remote LIBS While the point LIBS technique is comparatively well established and investigated, the remote LIBS technique at te ns of meters has only become a source of active research in the last 5 – 10 years. Recent publications repor t successful remote LIBS measurements at distances of up to 180 meters with potential distances of several kilometers, but these measurements were more qualitative than quantitative.3 A typical remote LIBS setup is given in Fig. 2.2. As can be seen, the components of a typical remote LIBS setup ar e the same as those used in point LIBS, except a telescope is used to transmit the focu sed laser beam and a telescope is used to collect the LIBS signal. Stand-off LIBS technique has been demons trated to be effective for the analysis of liquid steel,7 biological materials,8 cultural heritage,9 and remote atmospheric analysis.10 For example, typical remote LIBS spectra of liquid and solid stainless steel at 1420 C and 7.5 m range are shown in Fig 2.3.7 As can be seen from the figure, the LIBS emission from steel contains lines from Fe, Mn, Ti, Ni, Ca and Cr lines which are

PAGE 21

10 Figure 2.2 Schematic of LIBS lidar system.

PAGE 22

11 Figure 2.3 LIBS spectra from stainless st eel sample (liquid a nd solid at 1420C) obtained at a range of 7.5 m. Laser us ed: 220mJ per pulse, 1064nm, 5ns (From Palanco7; reprinted with permission).

PAGE 23

12 located in the UV region and their intensity de pends on the physical state of the material being tested 2.3.1 Remote LIBS using picoand nanosecond pulses LIBS measurements are usually carried out using nanosecond laser pulses at 1064 nm with a 5-20ns duration and energies on th e order of 10-500 mJ. The pulse power is usually on the order of 0.3 to 50MW with th e corresponding focused power densities of 3.8 to 6.4 GW/cm2 (for a 0.1 mm spot size diameter).1 Nanosecond lasers are preferred for remote LIBS applications because they are technologically well developed. Comp act and filed-deployable nanosecond laser systems are currently available for re mote LIBS with conventional focusing.1 Using nanosecond lasers ranges of up to 120 meters have been reported with remote LIBS.11 Picosecond (10-12 s) laser pulses have also been tested and compared to nanoand femtosecond regimes for the use in remote LIBS applications.11 2.3.2 Conventional remote LI BS using femtosecond pulses Femtosecond lasers with 10-15 s pulse durations have been investigated for use with remote LIBS application. Unfocused powers from these lasers may reach tens of terawatts. A comparative study of nano-, picoand femtosecond laser pulses for use with remote LIBS has been reported previously.12, 13 For example, a comparison spectra for a copper sample recorded at 25m range using 200 picosecond and 75 femtosecond pulses are presented in Fig. 2.4.12, 13 As can be seen from the figure, while

PAGE 24

13 Figure 2.4 R-LIBS spectra for a co pper sample at 25m: (a) using 75 fs pulses, (b) using 200 ps pulses. Wave lengths of observed copper lines are indicated in (a), positions of obs erved atomic oxygen and nitrogen lines are indicated in (b) (From Rohwetter12; reprinted with permission).

PAGE 25

14 the spectra are identical in the shorter wavelength range, the femtosecond LIBS signal shows an absence of oxygen and nitrogen em ission lines in the recorded spectrum.13 2.3.3 Filament-induced breakdown spectroscopy Because of high optical powers generated by femtosecond pulse lasers, a variant of LIBS is possible based on atmospheric fi lamentation produced by these pulses. The filamentation is the result of the Kerr effect as the inte nse laser pulse passes through the air resulting in the change of the medium refractive index. The induced refractive index change forms a lens acting to focus the light pulse. This approach is usually referred to as remote filament-induced breakdown spectroscopy (R-FIBS).14 The R-FIBS technique has been successfully applied fo r the analysis of a number of materials. For example, a typical R-FIBS spectrum of an aluminum samp le recorded at a 50m distance for a laser (Parameters: 108mJ, 10fs, 800 nm) is presented in Fig. 2.5.3 As can be seen, Al and Mg lines are clearly presented in the spectrum in the UV-VIS region. Also, spectra of Cu and Fe recorded at a 90m distance for a la ser (Parameters: 250mJ, 80fs, 800 nm) are presented in Fig. 2.6.14 As can be seen from the fi gure, Fe and Cu lines are clearly identified in the R-FIBS spectra. Femtosecond laser systems generating powers on the order of 1014 W needed for the filamentation are currently laboratory instruments due to their complexity and controlled operation requirements. An example of such a system is a “ Teramobile ” system.15 It should be added that filamenta tion requires careful control of the laser power and characteristics, and may be difficult to obtain consistent r eadings of the LIBS signal.

PAGE 26

15 Figure 2.5 R-FIBS spectrum taken for aluminum sample located 50m away. Laser used: 108mJ, 10fs, 800nm (From Liu3; reprinted with permission).

PAGE 27

16 Figure 2.6 R-FIBS spectra of copper (a) and steel (b), measured at 90 m distance. Laser: 250mJ, 80fs, 800 nm (From Stelmaszczyk14; reprinted with permission).

PAGE 28

17 It has been also suggested that the R-FI BS detected signal is not dependent on the sample distance, in contrast to the remote-LIBS excited by a focused beam.14 An experimental measurement as a function of range is shown in Fig. 2.7 along with a 1/R2 dependence. As can be seen, the uncertain ty of the signal is qu ite large and may fall within the 1/R2 dependence reported. As such, th e reported range independence of the detected signal in the filament-induced LIBS may be questionable and may require more experimental study.14 The introduction of the R-FIBS techni que suggests the possi bility of extending the range for remote LIBS detection to several kilometers.3, 16 Xu et. al.16 carried out calculations using the integrated intensity of the recorded signal and the noise of the detector to predict a possible detection range of 1 km. The result of their calculations is presented in Fig. 2.8. As can be seen, th e maximum detectable range for integrated intensity is estimated to be 1 km. Along these lines, Liu et. al.3 used a single Al 396.15 nm LIBS emission line signal and recorded the b ackground spectrum in order to determine the possible maximum detection range. The result of thei r calculation is presen ted in Fig. 2.9. As can be seen, the calculation suggested a possible detecti on range of 1.9 km. Both Xu et. al.3 and Liu et. al.16 took into account the 1/R2 attenuation of the emitted signal before reaching the detector. However, no atmospheric attenuation was taken into account as well as no analysis wa s carried out for different lines in the LIBS spectrum.

PAGE 29

18 Figure 2.7 Distance dependence of the range-corrected R-FIBS signal from the 521.8 nm line of copper (From Stemaszczyk14; reprinted with permission).

PAGE 30

19 Figure 2.8 Extrapolation of the single-shot detecti on limit for the integral intensity of the detected LIBS signal (From Xu16; reprinted with permission).

PAGE 31

20 Figure 2.9 Signal intensity extrapol ation of Al I: 396.15 nm line as a function of distance (From Liu3; reprinted with permission).

PAGE 32

21 With the possibility of extending the ra nge of stand-off LIBS measurements, the question of attenuation of the LIBS atomic emission spectrum by the atmosphere prior to reaching the detector becomes important. Unfortunately no lidar range calculations or simulations of remote-LIBS signal over the entire signal spectral range that ta ke into account atmospheric absorption have been carried out so far. A calculation of this kind is presented in this thesis.

PAGE 33

22 CHAPTER 3. INTRODUCTION TO RAMAN SPECTROSCOPY An overview of point and remote Ra man spectroscopy is presented in this chapter, along with typical e xperimental setups for both tech niques. The current state of remote Raman measurements is discussed 3.1 Point Raman Spectroscopy Raman spectroscopy has for many years been utilized as a complimentary tool to IR absorption spectroscopy for the investiga tion of vibrational and rotational spectra of molecules.6 The advantage of Raman spectrosc opy as compared to the regular IR absorption spectroscopy is the ability to inves tigate molecules without an initial dipole moment even though the intensity of Raman emission is comparatively week. Raman spectroscopy has become a useful analytical tool ever since th e invention of the laser. Similarly to LIBS, Raman spectroscopy may be subdivided into point Raman spectroscopy which is typically carried out in a laboratory environment at a close range and remote Raman measurements at a range of a meter to several kilometers. Both techniques are outlined in the following. A typical point Raman spectroscop y setup is presented in Fig. 3.1.6 As can be seen, a typical setup usually includes a CW or pulsed laser to excite the sample molecules, collection optics to collect the Ra man scattering signal, a detector (PMT with a photon counter or a modulator with a lock-in amplifier), an d a computer to process the spectrum and control the setup. It is importa nt to point out that a double monochromator

PAGE 34

23 Figure 3.1 Schematic of a poi nt Raman spectroscopy system.

PAGE 35

24 is usually used in order to eliminate the s cattered radiation of th e laser used for the excitation. 3.2 Remote Raman spectroscopy Remote Raman-scattering lidar technique was developed in the late 60s and originally demonstrated on wa ter vapor in the atmosphere.17 Ever since the introduction of this technique it has been successfully applied to the detect ion of several other atmospheric gases. Only major atmosphe ric gases may be measured with the Raman technique because the applicability of standoff Raman scattering technique is limited by the weak Raman scattering cross section. T hus, effective measurements of atmospheric gases are only possible for major atmospheric constituents, such as H2O, N2, O2, and CO2.18 Because the Raman cross section depends on the wavelength as -4, shorter wavelengths are preferred for Raman lidar systems. Since each molecule has a different Raman shift, one laser wavelength may be used for a large variety of species. Table 3.1 contains information about the characteristic Raman shifts of sel ected gaseous, solid and liquid species.19, 20 A typical remote Raman spectrosc opy setup is presented in Fig. 3.2.5 As can be seen, the difference between the point and re mote Raman systems is the introduction of a telescope to collect the weak Raman scattering signal. Stand-off Raman measurements may be car ried out at night or during the day, but there are differences in the components us ed to construct a remote Raman system depending upon the level of solar light that has to be reject ed in the optic al detection system. Nighttime measurements are best made using lasers operating around 300-350

PAGE 36

25 Molecule name Type of substance Raman shift, cm-1 N2 gas 2331 H2O gas 3460 O2 gas 1566 H2 gas 4155 O2 liquid 1552 N2 liquid 2326.5 Benzene liquid 992 CS2 liquid 655.6 Toluene liquid 1003 SiO2 solid 467 LiNbO3 solid / crystal 256 258 637 643 LiTaO3 solid / crystal 201 215 Table 3.1 Characteristic Raman shifts for selected gaseous, liquid and solid molecules.19, 20

PAGE 37

26 Figure 3.2 Schematic of a t ypical remote Raman system (Adapted from Chen5).

PAGE 38

27 nm since no sunlight is present. In the daytime, solar background radiation interferes with the measurement of Raman emission lin es for wavelengths longer than about 265 nm. Therefore daytime measurements are best made with lasers in the 260-265 nm spectral range using a solar-blind lidar.21 However, for laser wavelengths of about 265 nm or shorter, ozone and oxygen attenuation increases.22 Therefore, since both the ozone and oxygen absorb in the UV region and their concentrations vary, it is important to measure ozone and oxygen concentrations if quantitative remote Raman measurements are to be performed. Stand-off Raman techniques have now been applied for the dete ction of a variety of liquid and solid samples with current reported detection ranges of up to 533m.23 For example, a stand-off Raman spectrum of ethyl benzene and toluene recorded at a distance of 10 m using a 532 nm laser (parameters: 35mJ, 8ns) is presented in Fig 3.3.20 As can be seen, several distinct Raman scattering lin es are present in the spectrum. Another example is shown in Fig. 3.4 for a remote Ra man spectrum of naphthalene recorded at a distance of 122m using a 532 nm laser (parameters: 25mJ / pulse).20 As can be seen, the naphthalene spectrum is shown for a 1s and 10s signal integration time. A detailed review article on the current state of the remote Raman technique has recently been published.20 Measurements of cyclohexane and acetone trile with 266 nm excitation have also been carried out with possible pred icted detection rages of over 1000 m.23

PAGE 39

28 Figure 3.3 Portions of remote Raman spectra of ethyl be nzene and toluene in the spectral range 180 – 1800 cm-1 at 10m distance excited with a 532nm, 35mJ, 8ns laser (From Sharma4; reprinted with permission).

PAGE 40

29 Figure 3.4 A portion of remote Raman spectra naphthalene sample at 122 m distance measured with the directly coupled coaxial pu lsed Raman system. Laser 532 nm, 25 mJ per pulse at the sample, 20 Hz, slit 100 m (From Sharma20; reprinted with permission).

PAGE 41

30 3.3 Joint remote LIBS and Raman lidar measurements Simultaneous LIBS and Raman measurements have also been attempted recently.24 For example, the joint stand-off LI BS and Raman spectra of chalk has been measured and these results are presented in Fig. 3.5. While the spectra for the two techniques are different, the joint LIBS a nd Raman techniques are complimentary to each other. For example, there are several gr oups now conducting research on the increased detection probability of selected material s using both remote LIBS and Raman lidar.

PAGE 42

31 Figure 3.5 (a) LIBS and (b) remote Ra man spectra of a chalk sample, Raman spectrum of a calcite spectrum for re ference. Laser: 532 nm, 35 mJ, 8ns. ( From Weins24; re p rinted with p ermission ) (a) (b)

PAGE 43

32 CHAPTER 4. NEW ATMOSPHERIC TRANSMISSION COMPUTER CODES FOR UV-VISIBLE-IR W AVELENGTHS USEFUL FOR LIBS AND RAMAN LIDAR In this chapter an overview of curre nt atmospheric transmission models and simulation computer codes are given. A modification to the Hitran-PC atmospheric line spectra program is described which make s it possible to generate atmospheric transmission spectra suitable for remo te UV LIBS and Raman simulations. 4.1 Current atmospheric transmission simulations HITRAN is the largest and most used database containing spectral absorption information about major atmospheric consti tuents and was developed by the Air Force Research Labs over the past 30 years.25 It is extensively used by various atmospheric transmission simulation codes to quantitatively predict the transmission of the atmosphere as well as individual gases at different pressures and temperatures. The structure of the overall Hitran-2004 da tabase is presented in Fig. 4.1, but it consists mainly of the large HITRAN lin e-by-line database which contains over 2,000,000 individual line intensitie s and parameters for most of the gases in the atmosphere for the wavelength range of 400 nm to the microwave region. In addition, a smaller optical absorption cross section f ile contains parameter values for a few molecules in the UV (200 nm to 400 nm).

PAGE 44

33 Figure 4.1 File structur e of the HITRAN database s howing the large line-by-line molecular database and the smaller UV cross section database(Adapted from Rothman25).

PAGE 45

34 There are a number of atmospheric transm ission codes that use the HITRAN database, including the Air Force developed FA SCODE (PCLnWin) program and the USF developed HITRAN-PC program developed for use on a PC. Both FASCODE and HITRAN-PC use the line-by-line Hitran data base in order to generate atmospheric transmission spectra.25 These programs can be used to calcul ate the high resolution transmission spectrum of the atmosphere. For example, Fig. 4.2 shows the calculated transmission spectrum of the atmosphere for a 2 km horiz ontal path using the HITRAN-PC program. As can be seen from the transmission spectrum, the atmospheric attenuation is considerable for some wavelength regions and thus must be taken in to account if long range stand-off LIBS or Raman emission wave lengths coincide with these absorption regions. 4.2 Modification of Hitran-PC softwa re suitable for remote LIBS and Raman signal attenuation modeling The HITRAN-PC program was modified to also add the attenuation due to the UV cross section database to th at of the line-by-line HITRAN database. The number of UV cross section molecules in the HITRAN da tabase for typical atmospheric gases is currently limited to N2O, NO2, NO3, O3 and BrO. Atmospheric attenuation coefficients in the 0.25 – 0.4 m were calculated using the modified Hitran-PC software and the cros s section data. For example, a 2000 m path composite UV atmospheric transmission sp ectrum for available UV molecules and

PAGE 46

35 Figure 4.2 Transmission spectrum of the atmosphere generate d by HITRAN-PC program: US. Sta ndard atmosphere, 2000 m path.

PAGE 47

36 typical urban gas concentra tions is presented in Fig 4.3 with both Rayleigh scattering attenuation component and without it. As can be seen, the primary contribution to absorption is due to the ozone in the shor t wavelength region. Figure 4.4 shows the combined spectra of Fig 4.2 and Fig 4.3. For later use with the mid-IR Raman lidar calculation, Fig. 4.5 shows the same combined spectrum but out to 2.5m. Another modification to the Hitran -PC is the use of PNNL database26 in order to generate IR transmission sp ectra for a large number of ch emicals not available in the current HITRAN edition.

PAGE 48

37 Figure 4.3 Composite transmission spec trum of the atmosphere in the UV region calculated from the HITRAN2004 database UV cross-sections data for typical urban concentrations of BrO, N2 O, NO2, NO3, O3 using the modified HITRAN-PC program (a – w ithout Rayleigh scattering, b – with Rayleigh scattering). Path: 2000 m. Partial pressures: BrO (1-2ppt), N2O (0.3 ppm), NO2 (0.8 ppm), NO3 (0.35 ppb), O3 (0.5 ppm) (a) (b)

PAGE 49

38 Figure 4.4 Combined spectra of Fig 4.2 and 4.3 for a 250 – 850 nm spectral range (2 km path).

PAGE 50

39 Figure 4.5 Combined UV cross sec tion and line-by-line spectra obtai ned with the Hitran-PC program for a 250 nm – 2.3 m region (2 km path).

PAGE 51

40 CHAPTER 5. MODIFICATION OF LI DAR EQUATION SUITABLE FOR LIBS AND RAMAN LIDAR This chapter explains the changes which were made to the current lidar equation for it to be suitable for LIBS and Raman LIDAR simulations. 5.1 Current LIDAR equations The lidar equation is used to predict the returned lidar signal if the initial parameters of the system used for the LIDAR measurements as well as the absorption coefficients of the atmosphere at the m easurement wavelengths are know. A typical lidar setup is presented in Fig. 5.1. As can be seen typical com ponents of most lidar systems consist of a laser, beam forming optic s or transmit telescope, illumination of the distant target, a telescope used to collect th e backscattered return signal, and an optical detector. Figure 5.2 shows a schematic of the lidar beam transmission to a remote target and the detection of the back scattered radiation. Based upon this simple schematic, the lidar equation can be given as28 Pr = Pt A e -2 R / R2 (1) where Pr is the returned lidar signal in Watts, Pt is the transmitted laser power, is the reflectivity of the target area defined for scattering into steradians, A is the receiving

PAGE 52

41 Figure 5.1 Schematic of a typical lidar or differential-absorption lidar system (From Killinger28; reprinted with permission).

PAGE 53

42 Figure 5.2 Pictorial schematic of lidar equation geometry showing backscatter of the li dar signal and collecti on by a telescope (From Killinger28; reprinted with permission).

PAGE 54

43 telescope area, K is the optical effici ency of the overall lidar system, is the attenuation of the atmosphere at the laser wavelength, and R is the range to the target area. If the parameters of the system, such as laser power, pulse width, telescope diameter etc. are known, the lidar equation may be used for the prediction of the signal detected by the system. There is a number of modified lidar equa tions introduced to suit various applications and lidar types, the general form of th e equation being the same for all of them but with the value of the effective reflectivity being different for Raman lidar, fluorescence lidar, aerosol lidar and hard target lidar. 5.2 Modified LIDAR equation suit able for LIBS and Raman LIDAR In order to model LIBS and Raman signa l attenuation as a function of range a LIDAR equation needs to be modified in orde r to take into account that the signal only travels one way from the excitation site to th e detector as opposed to the typical two-way propagation ( ie. e -2 R ) of the lidar equation. One of the most simple ways to modify the lidar equation suitable for remote LIBS or Raman lidar, is to use only a one-w ay propagation term in the Beer-Lambert term, and to modify the emission from the remote hard target area as Ps instead of Pt In this case, the LIBS/Raman lidar equation is Pr = Ps A e R / R2 (2) The above equation is suitable for horizontal path calculations and makes it possible to predict the spectrum at the detection site from the spectrum at the excitation site (or vise

PAGE 55

44 versa) by applying the modifi ed lidar equation with the atmospheric attenuation coefficients calculated using the modified HI TRAN-PC software package. As such, the value of the LIBS or Raman lidar signal, Pr can be calculated as a function of wavelength over the complete sp ectrum of the LIBS/Raman em ission and as a function of range. Then the S/N values can be determined from the ratio of Pr / NEP, where NEP is the Noise-Equivalent-Power of the detector.

PAGE 56

45 CHAPTER 6. ESTIMATED LIBS AND RAMAN LIDAR EMISSION SPECTRUM AND INTENSITY An estimate of the value for the remo te LIBS and Raman emission spectral intensity, Ps at the excitation site is made in this chapter. The calculation procedure for this value is explained. The spectra intensity obtained will later be used in the modified lidar equation for the S/N versus range lidar simulations. 6.1 Estimated LIBS power spectrum at the excitation site An experimental value for the power spectrum or spectral intensity at the excitation site, i.e. Ps as a function of wavelength, has not been reported as far as we know. As such, a reasonable prediction for su ch a value was obtained by using one of the reported LIBS spectral plots and calculating a reasonable intensity level for the entire spectrum. Toward this goal, the LIBS sp ectrum for an aluminum sample presented in Fig. 2.5 was converted to a BMP image usi ng Adobe Photoshop and digitized using the “ Get Data Graph Digitizer ” software.27 The obtained spectrum was then normalized by calculating the area under the curve and equalizi ng the area to the total power of the laser pulse on target times a multiplication factor which estimates the LIBS emission efficiency factor, Leff or Ps ( ) = Pt Leff (3)

PAGE 57

46 As an estimate, we used a value of Leff of about 0.01 (ie. 1% ). Here, Leff incorporates the target reflectivity term. It should be noted that the LIBS emission power, Ps is approximately the absorbed energy divided by th e pulse length or lifetime of the plasma pulse t. A typical value of the plasma lifetime, t, is about 0.01 to 1s.30 Our estimate of Leff may not take into account a ll the factors contri buting to the dissip ation of the laser radiation because th e initial power Pt of the laser pulse is also absorbed by the target and scattered from the surface. The LIBS emi ssion from the plasma may have a blackbody radiation component and that from atom ic/ionic emission. The latter is the only component measured in a LIBS experiment. As an example, using the LIBS spect rum in Fig. 2.5 and the laser power of (108mJ/1s) 0.01 = 1KW (108mJ, t = 1 s, Leff = 0.01) within the total LIBS spectral output, the resultant LIBS power spectrum calculated for Ps ( ) is presented in Fig. 6.1. As can be seen, the estimated power of the LI BS spectrum at the excitation site is on the order of 20 W change for the line near 400 nm. Additional estimates of the LIBS emission power Ps ( ) were made for other cases and are being used to help produce a LIBS database of such emission spectral power curves but are not given in this thesis.

PAGE 58

47 Figure 6.1 Estimated power spectrum of a LIBS signal at the excitation site, Ps( ), for a peak laser pulse power of 1 kW and 1% LIBS efficiency calculated from the LIBS spectrum shown in Fig 2.5.

PAGE 59

48 6.2 Estimated Raman power spect rum at the excitation site A procedure analogous to the LIBS power spectrum calculation has been carried out in order to obtain an estimate for the sa mple Raman emission power spectrum. As an example, the naphthalene Raman spectru m presented in Fig. 3.4 for 10s integration was digitized. An estimate of the total pow er of the laser pulse was taken to be 25 MW, which is equivalent to a nanosecond pulse with pulse energy of 25mJ as reported in the article.20 The efficiency factor, however, for c onversion of excitation laser power into Raman emission is related to the Raman cro ss section. The Raman cross section for many Raman transitions is several orders of magnitude 10-6 to 10-10 less than that for the absorption of light.28 As such, the equivalent Raman efficiency factor, Reff may be on the order of 10-6 to 10-10 Thus, the power spectrum for the remote Raman emission may be given as Ps ( ) = Pt Reff (4) Using the Raman spectra from Fig. 3.4, laser pulse power of 25 MW, and an efficiency factor of 10-10, the resulting remote Raman emission power spectrum was calculated and the results are presented in Fig 6.2. As can be seen, the estimated power of the remote Raman spectrum at the excitation site is considerably smaller than that of a typical LIBS signal, the remote Raman po wer spectrum being on the order of several Watts.

PAGE 60

49 Figure 6.2 Estimated power spectrum of a Raman signal at the excitation site, Ps( ), for a peak laser power of 25 MW and 10-10Raman efficiency, and laser excitation wavelength of 266 nm Raman spectrum was calculated from Raman spectrum shown in Fig 3.4 (Napht halene), obtained at 532 nm, no 1/ 4correction was used.

PAGE 61

50 Along these lines of thought, we have al so started to compile representative remote Raman emission curves for several other examples and will be compiling a database for eventual use by our Raman lidar simulation program. For example, Fig 6.3 shows the calculated Raman emission lines from the remote hard target using the listed Raman shifts, 1/ 4 dependence, and laser excitation wavelengths of 266 nm, 532 nm, 1064 nm, and 1.57 m. As can be seen, the intensities drop-off toward the mid-IR due to the 1/ 4 dependence and also shows the wide spectral separation of the lines in the Mid-IR. Also, we have used the Raman spectru m of Fig 3.4 to gene rate the shifted spectrum for these four laser wavelengths. The results are shown in Fig. 6.4 (without -4 dependence to better sh ow the shifted lines). Finally, Fig. 6.5 is a composite plot of the Raman emission (Fig 6.4) and the atmospheric transmission (Fig 4.5) to show the overlap of the potential Raman emission spectra with that of the atmosphere.

PAGE 62

51 Figure 6.3 Estimated Raman emission lines for different laser excitation wavelengths and several different gas molecules located at the distant hard target (with -4 dependence taken into account). m

PAGE 63

52 Figure 6.4 Estimated Raman emission for the spectrum in Fig 6.2 (remote hard target of naphthalene; 266 nm excitation) for different laser excitation at 266 nm, 532 nm, 1064 nm, and 1.54 m (without -4 dependence). Wavelength, nm Power mW

PAGE 64

53 Raman emission Atmospheric transmission Figure 6.5 Composite plot of Rama n emission for four different laser excitation wavelengths (Fig 6.4) and transmission spectrum of the atmosphere for a 2 km path (Fig. 4.5).

PAGE 65

54 CHAPTER 7. LIBS AND RAMAN LIDAR SIMULATIONS In this chapter the power versus rang e and S/N versus range calculations for Raman and LIBS lidar signals at different wavelengths are carri ed out. The calculations are based on the atmospheric attenuation co efficients obtained with the modified HITRAN-PC software described in Chapter 4, modified lidar equation described in Chapter 5, as well as the sample calculated pow er spectra of LIBS and Raman signals at the excitation sites described in Chapter 6. 7.1 LIBS Lidar simulations In order to combine the atmospheric a ttenuation spectrum with the sample LIBS spectrum, a program was written in Visual Basic which interpolated the sample LIBS spectrum to match the wavelength points of the atmospheric attenuation file using linear interpolation. After the interpolation proce dure the two spectra were substituted into the modified lidar equation for the calculations of the remote LIBS power spectrum as a function of range. The result of the pow er spectrum range calculation for the LIBS sample spectrum from Fig. 6.1 as a function of wavelength is presented in Fig. 7.1. As can be seen, the intensity of the signal qui ckly decreases with range due to the 1/R2 term in the modified lidar equation, while the spec tral features of th e attenuated signal are preserved. After determining the power versus range dependence, the S/N ratio for the remote LIBS signal was determined at different wavelengths. A typi cal detectivity value

PAGE 66

55 Figure 7.1 Calculated power spectru m of LIBS as a function of range (1 – 25 m to the target) and waveleng th for a peak lase r pulse power of 1 KW, 1% LIBS efficiency, A = 0.01 m2 telescope area, K = 0.01 lidar efficiency. (Based upon Fig. 2.5 and power spectrum in Fig 6.1).

PAGE 67

56 of 1014 cmHz1/2/W (NEP = 10-12 W) for an average detector has been used in order to determine the S/N ratio at different wavele ngths and ranges by di viding the signal power by the NEP of the detector.28 The resulting S/N spectrum is presented in Fig. 7.2. As can be seen, the S/N values of many of the LIBS emission lin es fall below 1 for ranges of about 10000 m or so. Here, we have used NEP = (A B)1/2 / D* = (0.1 cm2 106 Hz)1/2 / (1014 cmHz1/2/W) = 10-12 W (3) for a 1 mm2 detector size and se lected bandwidth of 1/ t = 106 Hz. The computer program used for these calcu lations is listed in Appendix 1. The output of the program is a 3-column file (Range; Wavelength; Power or S/N) which was then imported to an elementary MATLAB rou tine, listed in Appendix 2, to carry out 3D plotting of the data. 7.2 Raman Lidar Simulations The calculations for range dependence of the Raman power spectrum have been carried out using a similar program. For example, the Raman spectra in Fig. 3.4 was used to generate the expected Raman emi ssion for an excitation of 266 nm. The resultant spectral power curve was then used in the modified lidar equation. The remote Raman power spectrum calculated values as a function of range for different wavelengths for a value of K = 0.01 and A = 0.01 m2 is presented in Fig. 7.3. As can be seen the lidar reduction of the signal is identical to

PAGE 68

57 Figure 7.2 Calculated S/N spectrum of a LIBS signal as a function of range (10 – 15 km to the target) a nd wavelength for a peak laser pulse power of 1 KW and 1% LIBS efficiency, A = 0.01 m2 telescope area, K = 0.01 lidar efficiency. (Based upon Fig. 2.5 and power spectrum in Fig 6.1).

PAGE 69

58 Figure 7.3 Calculated power spectrum of a Raman signal as a function of range (1 –25 m to the target) and wavelength for a peak laser pulse power of 25 MW and 10-10Raman efficiency, telescope area A = 0.01 m2, lidar efficiency K = 0.01. (Based upon spectrum in Fig. 3.4 and power spectrum in 6.2 using 266 nm excitation).

PAGE 70

59 that of the LIBS signal-range simulations, but with much smaller power of the Raman signal. The same estimate for the NEP of the detector was used for the Raman S/N versus range value. The results of the S/N versus range calculations for different wavelengths are presented in Fig. 7.4 and i ndicate ranges on the order of 100m for a S/N of 1 for the more prominent lines. As can be seen, the S/N ratio of the Raman signal is considerably smaller than that for the LIBS signal simulation, which is de termined by the initial power of the spectrum at the excitation site. It is worth pointing out, th at the actual NEP or detectivity values of the detectors used fo r remote LIBS and Raman measurements may differ in that, detectors used for Raman signal measurements may have increased effective NEP values, if the use of narrow band rejection optical filter to reduce the strong Rayleigh scattering at the exci tation wavelength is not effective. In general, the final results obtained in th is Chapter indicate that the attenuation of the remote LIBS and Raman signal is primarily due to the 1/R2 term in the lidar equation, atmospheric attenuation only being significant for selected regions in the UVIR spectrum. The contribution of the atmo spheric attenuation will depend on the laser wavelength used for the measurements as well as the compounds being analyzed. 7.3 Comparison of Initial Simulations with Selected LIBS and Raman LIDAR Studies It is interesting to note that recent experiments using 50 to 100 mJ Nd:YAG lasers have yielded LIBS and Raman lidar detecti on ranges on the order of a few 100m. Our

PAGE 71

60 Figure 7.4 Calculated S/N spectrum of a Raman signal as a function of range (100 – 300 m to the target) and wavelength for a peak laser pulse power of 25 MW and 10-10 Raman efficiency, A = 0.01 m2 telescope area, K = 0.01 lidar efficiency.

PAGE 72

61 Raman lidar calculations appear to be in this range, but our LIBS lidar calculation is a factor of 100 greater. It may be that our estimate of the Leff of 1% for LIBS emission may be high, and a value near 0.01% may be more consistent with our simulations and experiments. Additiona l studies are planned. Finally, a recent paper by Palanco31 has suggested that th e range dependence of the LIBS lidar signal may be further modified in that the value for the laser emission intensity at the target, Ps, may depend upon the range. Palanc o’s theory indicated that Ps may be proportional to 1/R3, but in addition the laser beam mode quality and other factors may influence the LIBS emission intensity. We plan to investigate this further both in theory and with experiments.

PAGE 73

62 CHAPTER 8. SUMMARY AN D FUTURE RESEARCH Research studies presented in this thes is have show that stand-off LIBS and Raman lidar can be modeled through use of the appropriate a tmospheric attenuation models, appropriate lidar equation, and use of sp ectral responses of the detection system. Our initial results indicate that LIBS and Ra man lidar should be able to detect chemical compounds at several tens of meters if not km ranges. The computational programs are documented in this thesis and were develope d for expediency and not for general use. Some of our future plans included comp aring these simulation results to 266 nm LIBS and Raman lidar experimental measur ements being conducted under a different program. In addition, we may incorporate so me of these initial lidar simulation studies into a more user friendly program that can th en be used to compare the predicted values to that from experimental LIBS and Raman lidar measurements. Toward this end, we may modify the current LIDAR-PC program developed at USF to incorporate these ideas.

PAGE 74

63 LIST OF REFERENCES 1. D. Cremers, L. Radziemski, Handbook of Laser-Induced breakdown spectroscopy John Wiley & Sons, Ltd, 2006. 2. A. Miziolek, V. Palleschi, I. Schechter Laser-Induced Breakdown Spectroscopy Cambridge University Press, 2006. 3. W. Liu, H Xu, G. Mejean, Y. Kamali, J. Diag le, A. Azarm, P. Simard, P. Mathieu, G. Roy, S. Chin, “Efficient non-gated filame nt-induced breakdown spectroscopy of metallic sample”, Spectrochimica Acta B 62 76-81 (2007). 4. S. K. Sharma, A. K. Misra, B. Shar ma, “Portable remote Raman system for monitoring hydrocarbon, gas hydrates and explosives in the environment”, Spechtrochimica Acta Part A 61 2404 – 2412 (2005). 5. T. Cheng, J. M. J. Madey, F.. M. Prince, S. K. Sharma, B. Lienert, “Remote Raman Spectra of Benzene Obtained from 217 Me ters Using a Single 532 nm Laser Pulse”, Applied Spectroscopy 61 (6), 624 – 629 (2007).

PAGE 75

64 6. E. R. Menzel, Laser Spectroscopy: Techniques and Applications Marcel Dekker, New York, 1995. 7. S. Palanco, S. Conesa, J. Laserna, “Aalytic al control of liquid steel in an induction melting furnace using remote laser-induced plasma spectrometer”, J. Anal. Atom. Spectrom. 19 462 – 467 (2004). 8. H. Xu, W. Liu, S Chin, “Remote time -resolved filament-induced breakdown spectroscopy of biological materials”, Optics Letters 31 (10), 1540 1542 (2006). 9. S. Tzortzakis, D. Anglos, D. Gray, “Ultr aviolet laser filaments for remote laserinduced breakdown spectroscopy (LIBS) analysis : applications in cultural heritage monitoring”, Optics Express 31 (8), 1139 – 1141 (2006). 10. J. Kasparian, M. Rodriguiez, G. Mejean, J. Yu, E. Salmon, H. Wille, R. Bourayou, S. Frey, Y.-B. Andre, A. Mysyrowicz, R. Saue rbrey, J.-P. Wolf, L. Woste, “White light filaments for atmospheric analysis”, Science 301 61 – 64 (2003). 11. R. Gronlund, M. Lundqvist, S. Svanberg, “Remote Imaging Laser-induced Breakdown spectroscopy and Laser-Indu ced Fluorescence Spectroscopy Using Nanosecond Pulses from a Mobile Lidar System”, Applied Spectroscopy 60 (8), 853 – 859 (2006).

PAGE 76

65 12. Ph. Rohwetter, J. Yu, G. Mejean, K. Stel maszczyk, E. Salmon, J, Kasparian, J. Wolf, L. Woste, “Remote LIBS with ultrashort pulses: characteristics in picosecond and femtosecond regimes”, J. Anal. At. Spectrom. 19 437-444 (2004). 13. Ph. Rohwetter, K. Stelmaszczyk, I. Woste, R. Ackermann, G. Mejean, E. Salmon, J. Kasparian, J, Yu, J. Wolf, “Fillament-induced remote surface ab lation for long range laser-induced breakdown spect roscopy operation”, Spectrochimica Acta Part B 60 1025 – 1033 (2005). 14. K. Stelmaszczyk, P. Rohwetter, “Long distance remote la ser-induced breakdown spectroscopy using filamentation in air”, Applied Physics Letters 85 (18), 3977 – 3979 (2004). 15. R. Bogue, “LIBS range extended through the use of a transportable terawatt laser system”, Sensor Review 25 (2), 105 – 108 (2005). 16. H. Xu, J. Bernhardt, “Understanding th e advantage of remote femtosecond laserinduced breakdown spectroscopy of metallic targets”, Journal of Applied Physics 101 033124 (2007). 17. D. Leonard, “Observation of Raman scatte ring from the atmosphere using a pulsed nitrogen ultraviolet laser”, Nature 216 142 – 143 (1967).

PAGE 77

66 18. K. Okamoto, Global environment remote sensing IOS Press, 2001. 19. A. Yariv, Quantum Electronics 3rd edition, John Wiley & Sons, 1989. 20. S. Sharma, “New trends in telescopic remote Raman spectroscopy”, Spectrochimica Acta Part A, (2007), in press. 21. F. Duarte, Tunable Laser Applications Marcel Dekker, Inc., 1995. 22. E. Trakhovsky, A. Ben-Shalom, U. Oppenheim A. Devir, L. Balfour, M. Engel, “Contribution of oxygen to attenuation in the solar blind UV spectral region”, Applied Optics 28 1588 – 1591 (1998). 23. M. Wu, M. Ray, K. Hang Fung, M. Ruckman, D. Harder, A. Sedlacek III, “Stand-off Detection of Chemicals by UV Raman Spectroscopy”, Applied Spectroscopy 54 (6), 800 – 806 (2000). 24. R. Weins, S. Sharma, J. Thompson, A. Misr a, P. Lucey, “Joint analyses by laserinduced breakdown spectroscopy (LIBS) and Raman spectroscopy at stand-off distances”, Spectrochimica Acta Part A 61 2324 – 2334 (2005). 25. L. Rothman, et. al., “The HITRAN 2004 molecular spectr oscopic database”, JQSRT 96 139 – 204 (2005).

PAGE 78

67 26. S. Sharpe, T. Johnson, R.Sams, P. Chu, G. Rhoderick, P. Johnson, “Gas-Phase Databases for Quantitative Infrared Spectroscopy”, Applied Spectroscopy 58 (12), 1452 – 1461 (2004). 27. Get Data Graph Digitizer http://getdata-graph-digitizer.com/ 28. D. K. Killinger, “LIDAR a nd Laser Remote Sensing”, Handbook of vibrational spectroscopy Ed. J. Chalmers and P. Griffiths, John Willey & Sons Ltd, Chichester, 2002. 29. R. Wiens, S. Sharma, J. Thompson, A. Misr a, P. Lucey, “Joint analyses by laserinduced breakdown spectroscopy (LIBS) and Raman spectroscopy at stand-off distances”, Spectrochimica Acta A 61 2324 – 2334 (2005). 30. D. K. Killinger, S. Allen, R. Waterbur y, C. Stefano, E. Dottery, “Enhancement of Nd:YAG LIBS emission of a remote target using a simultaneous CO2 laser pulse” Optics Express 15 (20), 12905 12915 (2007). 31. Santiago Palanco, Cristina Lopez-Moreno, J. Javier Laserna, “D esign, construction and assessment of a field-deployable la ser-induced breakdown spectrometer for remote elemental sensing”, Spectrochimica Acta Part B 61 88-95 (2006).

PAGE 79

68 APPENDIX A. MODIFICATION OF HITRAN-PC PROGRAM TO ADD UV CROSS SECTION ABSORBTION A modification to the HITRAN-PC program (d eveloped at USF) was made to permit the addition of the absorption due to the UV cross section database in HITRAN to that of the line-by-line HITRAN database. Th e related portion (3 pages) of the complete Hitran-PC (300 pages) is given below.

PAGE 80

69 Option Strict On Option Explicit On Module X_Sect_Module '/ Sub adds all x-Section headers (for all records and files) to the array containing the headers information. Public Sub Read_X_Sect_Headers_From_Dir(ByVal directory As String) Dim file_names As System.Collections.ObjectModel.ReadOnlyCollection(Of String) '/ Collection to contain the names of th e ifles in current Dim i As Integer '/ Auxiliary loop variable Dim num_of_files As Integer '/ Number of files in the ReadOnlyCollection Dim str_position As Integer '/ Length of the file name string Dim extension As String '/ Extension of the file to work with file_names = FileIO.FileSystem.GetFiles(directory) '/ Forming the Collection of file names within the cirrent director y num_of_files = file_names.Count '/ Saving the number of files in the directory For i = 0 To num_of_files 1 If X_Sect_Headers.Length() >= 2501 Then Exit For '/ Making sure the number of elements is not higher than 2500 str_position = file_names.Item(i).Length 3 '/ Determinig the starting position to form the su bstring extension = file_names.Item(i).Substring(str_position, 3) '/ Determining the extension of the file (should b e 3 symbols) extension = extension.ToUpper() '/ Making all leters capital in the extension If extension = "XSC" Then Add_X_Sect_Headers_from_file_TWO(file_names.Item(i)) '/ Calling the function below in ord er to fill out the headers Next i End Sub '/ Sub adds the headers of the records (in the file provided) to the array containing the headers information. '/ This procedure works fine but it is slow, the procedure below is 6 times faster, TWO should be used instead Public Sub Add_X_Sect_Headers_from_file_ONE(ByVal file_name As String) Dim file As New Microsoft.VisualBasic.FileIO.TextFieldParser(file_name) '/ File to read the cross section data from Dim header_str As String() '/ Array of strings to obtain the header information from a file Dim i As Integer '/ Auxiliary loop variable Dim new_item_index As Integer '/ Index of the new head (to be added) in the head containing array Dim record_num As Integer '/ Record number in the current file we work with Dim aux_var As Double '/ Used to determine the number of lines beween the header lines (need to be read in) Dim line As String '/ Just an auxiliary variable to read in the line form the cross section file Dim num_lines As Integer '/ Number of lines to read between the header information strings record_num = 1 '/ Reading the header of the file (information about the molecule and cross sections) file.TextFieldType = FileIO.FieldType.FixedWidth file.FieldWidths = New Integer() {20, 10, 10, 7, 7, 6, 10, 5, 15, 4, 3, 3} '/ HITRAN 2004 defiened lengths of the field s While Not file.EndOfData() '/ Determining the index of the array element new_item_index = X_Sect_Headers.Length() '/ Writing information into the Tmain status strip and the label on the X-sect tab of the info panel Write_Info("Reading record & new_item_index.ToString())

PAGE 81

70 If new_item_index >= 2501 Then '/ This is necessary because the formating is impared after 2500 MsgBox("The number of cross section records in the list is limited to 2500") Exit While End If ReDim Preserve X_Sect_Headers(new_item_index) '/ Resizing the array to add new header information Try header_str = file.ReadFields Catch ex As Microsoft.VisualBasic.FileIO.MalformedLineException MsgBox("File: & file_name & ". Line & ex.Message) '/ Even though tere is a mistake beacuse of the type in the database it seem to be reading '/ everything OK anyway End Try X_Sect_Headers(new_item_index).Chem_formula = header_str(0) '/ String Chemical formula of the molecule X_Sect_Headers(new_item_index).wn_min = System.Convert.ToDouble(header_str(1)) '/ Double Starting wavenumb er X_Sect_Headers(new_item_index).WN_max = System.Convert.ToDouble(header_str(2)) '/ Double Ending wavenumber X_Sect_Headers(new_item_index).Num_of_points = System.Convert.ToInt32(header_str(3)) '/ Integer contains the num ber of datapoints X_Sect_Headers(new_item_index).Temp = System.Convert.ToDouble(header_str(4)) '/ Double contains the temp eratur for which the X_Sect_Headers(new_item_index).Press = 1 '/ Double contains the part ial pressure of the X_Sect_Headers(new_item_index).Max_X_sect = System.Convert.ToDouble(header_str(6)) '/ Double contains the maxi mum cross section value X_Sect_Header_Data(new_item_index).Resolution = System.Convert.ToDouble(header_str(7)) '/ Double (sometimes thi s is a string) contains X_Sect_Headers(new_item_index).Name = header_str(8).ToString '/ String Name of the molec ule X_Sect_Header_Data(new_item_index).Not_used() '/ This field is not used at this point X_Sect_Headers(new_item_index).Broadener = header_str(10).ToString() '/ String Type of broadenin g (probably air or self) X_Sect_Headers(new_item_index).Ref_num = System.Convert.ToInt32(header_str(11)) '/ Integer reference number X_Sect_Headers(new_item_index).File_name = file_name '/ The name of the file with the current record X_Sect_Headers(new_item_index).Rec_Num_In_File = record_num '/ The number of the record i n the file X_Sect_Headers(new_item_index).UV_IR = UV_or_IR(file_name) '/ String indicating the type of specrum (UV or IR) '/ This portion is used in order to skip the cross section data and advance to the next '/ record with a different temperature value aux_var = X_Sect_Headers(new_item_index).Num_of_points / 10 aux_var = System.Math.Ceiling(aux_var) num_lines = System.Convert.ToInt32(aux_var) For i = 1 To num_lines Try line = file.ReadLine() '/ Used to be assigned to Line variable but this is not necessary since the string i s not used Catch ex As Microsoft.VisualBasic.FileIO.MalformedLineException MsgBox("File: & file_name & ". Line & ex.Message) End Try Next i record_num = record_num + 1 '/ Increasing the record number in the current file, to write sown for information End While file.Close() '/ Clearing the information in the Tmain window Write_Info("") End Sub

PAGE 82

71 '/ Sub adds the headers of the records (in the file provided) to the array containing the headers information. '/ This procedure is 6 times faster than the one above (ONE), should be used all the time Public Sub Add_X_Sect_Headers_from_file_TWO(ByVal file_name As String) '/ file_name -> name of the file to work with Dim header_aux_str As String '/ String in which a line form the file is read Dim header_str(11) As String '/ Array of strings to obtain the header information from a file Dim new_item_index As Integer '/ Index of the new header (to be added) in the header containing array Dim record_num As Integer '/ Record number in the current file we work with Dim cur_location As Long '/ Current location in the file Dim new_location As Long '/ New location in the file Dim shift_location As Long '/ Shift from the current location in the file to advance to the new record Dim aux_var As Double '/ Variable used for the calculation of the number of lines to skip Dim num_lines As Long '/ Number of lines to be skipped after the header line record_num = 1 FileOpen(10, file_name, OpenMode.Input, OpenAccess.Read) While Not EOF(10) new_item_index = X_Sect_Headers.Length() '/ Determining the index of the array element If new_item_index >= 2501 Then '/ This is necessary because the formating is impared after 2500 MsgBox("The number of cross section records in the list is limited to 2500") Exit While End If ReDim Preserve X_Sect_Headers(new_item_index) '/ Resizing the array to add new header information '/ Writing information into the Tmain status strip and the label on the X-sect tab of the info panel Write_Info("Reading record & new_item_index.ToString()) '/ Reading the header of the file (information about the molecule and cross sections) header_aux_str = InputString(10, 102) '/ Forming the array of strings with the fields formated accordingly header_str(0) = header_aux_str.Substring(0, 20) header_str(1) = header_aux_str.Substring(20, 10) header_str(2) = header_aux_str.Substring(30, 10) header_str(3) = header_aux_str.Substring(40, 7)

PAGE 83

72 APPENDIX B. PROCEDURE TO RE SCALE THE DIGITIZED SPECTUM TO THE RESOLUTION OF THE ATMO SPHERIC ATTENUATION SPECTRUM

PAGE 84

73 '/ Procedure adds a spectrum to the plot, does the resize to the current resolution also Private Sub Add_Spectum_MainMenu_Click( ByVal sender As System.Object, ByVal e As System.EventArgs) Handles Add_Spectum_MainMenu.Click Dim file As Microsoft.VisualBasic.FileIO.TextFieldParser '/ File with the spectral information to parse Dim spectrum_str() As String '/ Array of strings to obtain the header information from a file Dim spectrum_arr(,) As Double '/ Array containing the spectrum which has been read in Dim new_item_index As Integer '/ Index of the next (new) array element Dim freq As Double Dim min_rec As Integer Dim max_rec As Integer Dim i As Integer Dim y_sum As Double '/ Sum of the values along the Y axis '/ Values used for the interpolation procedure Dim x As Double Dim x1 As Double Dim x2 As Double Dim y1 As Double Dim y2 As Double HitranPC_OpenFileDialog.Filter = "Text datafiles (*.txt)|*.txt|All files (*.*)|*.*" If HitranPC_OpenFileDialog.ShowDialog() = Windows.Forms.DialogResult.OK Then file = New Microsoft.VisualBasic.FileIO.TextFieldParser(HitranPC_OpenFileDialog.FileName) '/ Initially resizing the array ReDim spectrum_arr(1, 0) '/ Indicating the file we read in is delimited file.TextFieldType = FileIO.FieldType.Delimited file.Delimiters = New String () {vbTab} While Not file.EndOfData() '/ Determining the index of the last array element new_item_index = System.Convert.ToInt32(spectrum_arr.Length() / 2) ReDim Preserve spectrum_arr(1, new_item_index) Try spectrum_str = file.ReadFields Catch ex As Microsoft.VisualBasic.FileIO.MalformedLineException MsgBox( "Error importing file: & ex.Message) End Try '/ Reading in the wavenumber and the intensity

PAGE 85

74 spectrum_arr(0, new_item_index) = System.Convert.ToDouble(spectrum_str(0)) spectrum_arr(1, new_item_index) = System.Convert.ToDouble(spectrum_str(1)) End While file.Close() ReDim Read_In_Data(1, new_item_index) For i = 1 To new_item_index Read_In_Data(0, i) = spectrum_arr(0, 1) + i ((spectrum_arr(0, new_item_index) spectrum_arr(0, 1)) / (new_item_index 1)) x = Read_In_Data(0, i) Find_record_number(x, spectrum_arr, new_item_index, min_rec, max_rec) x1 = spectrum_arr(0, min_rec) x2 = spectrum_arr(0, max_rec) y1 = spectrum_arr(1, min_rec) y2 = spectrum_arr(1, max_rec) Read_In_Data(1, i) = Linear_Interpolate(x, x1, x2, y1, y2) Next i y_sum = 0 FileOpen(23, "C:\output.txt" OpenMode.Output) For i = 1 To new_item_index WriteLine(23, Read_In_Data(0, i), Read_In_Data(1, i)) y_sum = y_sum + Read_In_Data(1, i) Next i WriteLine(23, y_sum) FileClose(23) End If End Sub

PAGE 86

75 APPENDIX C. MATLAB PROGRAM TO LINK THE ATMOSPHERIC ATTENUATION DATA, POWER SPECTR UM DATA AND PLOT 3D S/N VS RANGE VS. WAVELENGTH FOR LIBS/RAMAN LIDAR fid = fopen('atten.hpa', 'r'); libs = fopen('output.txt', 'r'); a = fscanf(fid, '%g,%g', [2 1875]); % It has two rows now. qs = fscanf(libs, '%g,%g', [2 1875]); % It has two rows now. l = qs(2,:); for i = 1:1875 r(i) = 1 + (i 1) (24 / (1875 1)); end w = 10000 ./ a(1,:); k = a(2,:); [X, S] = meshgrid(r,l); [Y, Y] = meshgrid(w,w); [M, M] = meshgrid(k,k); fprintf('%f\n', a(2,:)); % Ratio taken into account Z = ((25 .* 10e-4 .* 10e-4 ./ 91672797.2371662)) .* S .* exp(-X .* M ./ 1000) ./ (pi .* X .* X) ; surf(X,Y,Z, 'FaceColor','red', 'EdgeColor','none'); camlight hsv; lighting gouraud; fclose(fid); fclose(libs);

PAGE 87

76 The MATLAB program interacted with the output of the HITRAN-PC program as shown in the following flowchart: Initial BMP spectrum Rescaling the power spectrum to the resolution of Hitran-PC generated atmospheric transmission file ( A pp endix B ) Line-by-line + UV cross section data generated atmospheric transmission (Appendix A)Hitran-PCMATLAB routine (Applying the modified LIDAR equation and plotting) (Appendix C) Power and S/N vs. range vs. wavelength plots Digitizing the spectrum (Get Data Digitizer) and normalizing according to the total power