Spectral Line Shapes in Plasmas and GasesView this Special Issue
Research Article | Open Access
Time-Resolved Spectroscopy Diagnostic of Laser-Induced Optical Breakdown
Transient laser plasma is generated in laser-induced optical breakdown (LIOB). Here we report experiments conducted with 10.6-micron laser radiation, and with 1.064-micron fundamental, 0.532-micron frequency-doubled, 0.355-micron frequency-tripled Nd:YAG laser radiation. Characterization of laser induced plasma utilizes laser-induced breakdown spectroscopy (LIBS) techniques. Atomic hydrogen Balmer series emissions show electron number density of measured approximately 10 s and 1 s after optical breakdown for and Nd:YAG laser radiation, respectively. Recorded molecular recombination emission spectra of CN and Swan bands indicate an equilibrium temperature in excess of 7000 Kelvin, inferred for these diatomic molecules. Reported are also graphite ablation experiments where we use unfocused laser radiation that is favorable for observation of neutral emission due to reduced cation formation. Our analysis is based on computation of diatomic molecular spectra that includes accurate determination of rotational line strengths, or Hönl-London factors.
Laser-induced breakdown spectroscopy (LIBS) is a valuable technique for determining elemental composition with the ability to analyze solids, liquids, and gases with little or no sample preparation suitable for onsite analyses . The success of LIBS is in part due to the ease of availability of nominal 10-nanosecond Nd:YAG laser radiation from compact devices. Several other laser sources however have been historically applied for generation of microplasma with subsequent measurement methods based on use of atomic emission spectroscopy . Recent interest includes applying dual- and multipulse excitation for the purpose of increasing sensitivity of LIBS. Examples of advantages of multipulse excitation include an increase of the plasma volume or plasma reheating by the second pulse, in turn, enhancement of detection limits for LIBS 10- to 100-fold and/or a decrease in relative standard deviation when comparing single- with double- pulse bursts . Following short-pulse UV-excitation a CO2 laser may be used to enhance detection from a distance of atomic and molecular species [4–6]. Some applications are also designed with eye-safety in mind, required for remote and/or field-safe LIBS systems [7, 8]. Applications of CO2 lasers also include aerosol measurements [9, 10].
In this work, we utilize a CO2 transversely excited atmospheric-pressure (TEA) laser at the 10.6-micron wavelength to generate the LIBS plasma. The CO2 laser-induced plasma is expected to show a larger homogenous core than the Nd:YAG plasma. This is primarily due to prolonged, intensive, and spatially extended excitation that masks the temporal and spatial variation of the CO2 laser beam from pulse to pulse. Recorded atomic and molecular emission spectra are evaluated and compared with diagnostics that we applied previously for characterization of Nd:YAG laser-induced microplasma .
For laser ablation studies, we use the Nd:YAG laser operated in the second and third harmonic radiation modes. Graphite laser-ablation is one possible process that can lead to carbon nanoparticle formation [12, 13]. The laser-plasma conditions during the nominal nanosecond laser exposure during ablation can be measured by recording temporally and spatially averaged spectra. In this work, time-resolved methods are applied to investigate the evolution of C2 and C3 formation. In addition, the carbon layer is studied which is deposited at the walls of our Laser-Induced Breakdown Spectroscopy (LIBS) cell. The application of visible absorption spectroscopy and transmission electron microcopy allows us to reach conclusions about the graphite ablation product.
2. Experimental Details
The experimental arrangement is identical to a standard LIBS setup, that is, laser, spectrometer, detector-system, and a sample. The CO2 laser experiment details are as follows: the CO2 laser is the model TEA-820 Lumonics, the spectrometer is the 1/2 meter model 500 SpectraPro Acton Research Corporation, and the detector is the intensified model 1460 EG&G Princeton Applied Research detector/controller Optical Multichannel Analyzer (OMA). The laser plasma is first generated for alignment purposes and initial diagnostic studies in laboratory air using focusing. The plasma is then generated in CO2 gas expanding into the focal volume at a rate of 10 L/m into the laboratory. Survey spectra from laser-induced optical breakdown (LIOB) in CO gas are reported here.
Typical operating parameters of our CO2 TEA laser are (i) pulse width of 70–100 ns full-width at half-maximum, or than Nd:YAG pulses, (ii) pulse energy of 2000–3000 mJoule, or typically 10– than nominal Nd:YAG laser devices, and (iii) wavelength of 10.6 , or 10x longer than the 1.064 of the fundamental Nd:YAG laser radiation.
The Nd:YAG laser experiment details are as follows: the Nd:YAG laser is the model 10 Hz Quantel Brilliant operated in the second and third harmonic modes at 532 and 355 nm, respectively. The laser-induced plasma was created using unfocussed radiation impinging upon a 10 mm diameter spectrally pure graphite pellet that was cut from a graphite rod. The target was shifted from measurement set to measurement set, but was kept stationary during individual experiments.
The laser pulse energy and pulse-width of the fundamental mode (1064 nm) was adjusted electronically by changing the delay between the flash-lamp and the Q-switch trigger. At an optimum delay (190 ), the laser pulse energy amounted to 120 mJ at 532 nm, and 60 mJ at 355 nm. Electronic adjustments of the Q-switch trigger allowed us to reduce the energy per pulse by a factor of 3. An argon background gas was used at various pressures between 10 and 200 Torr.
Time-resolved emission spectra were taken using the model Andor Technology IStar ICCD and the associated Andor Shamrock spectrograph. The detection system shows a temporal resolution of 10 ns. The overall spectral resolution amounted to 8 nm and 1 nm for the 150 grooves/mm and 1200 g/mm grating, respectively. No entrance slit was used for the spectrometer, instead we employed a 50-micron optical fiber.
3. Experimental Results
Temporal characteristics of CO2 TEA laser radiation include a reasonably steep rise-time followed by an extended tail. Figure 1 shows a recorded pulse profile.
The emissions from the plasma are strong enough to allow us to monitor spectra from single breakdown events; however, we elected to average over 50 individual laser-plasma events. Typical recorded atomic and molecular spectra are illustrated in what follows.
Figure 2 shows the recorded and fitted spectrum of the CN sequence of the CN Violet System. The spectral resolution amounted to 0.12 nm when using the 2400 grooves/mm grating. Particularly strong CN emissions occur and dominate the molecular emissions in the visible region when using CO2 flow expanding into laboratory air. We used a gate width of 20 and a time delay of 20 . The inferred temperature amounts to K, using our fitting routines for diatomic spectra [14–16]. Figure 2 also shows an insert near the wavelength region of 4-4 band of CN that we excluded in the fitting routine. The additional spectroscopic feature indicates presence of the C2 Deslandres D’Azambuja system .
Figure 3 shows the recorded spectrum of the C2 sequence of the Swan bands. The spectral resolution amounts to 0.4 nm for the 1200 g/mm grating. Also indicated are atomic carbon lines , recorded for a delay time of 10 after optical breakdown, using an intensifier gate width of 100 . After subtraction of the atomic signatures and/or exclusion of regions with atomic lines, fitting of the molecular spectra may be accomplished using our standard fitting procedures [19, 20]. Results of C2 Swan bands fitting (not illustrated here) for 20 time delay and 20 gate-width (compare Figure 2) indicate a temperature of K.
In our fitting procedures of diatomic molecular spectra , we usually employ local thermodynamic equilibrium (LTE) considerations early in the plasma decay. Focused laser radiation causes a reasonably high electron density, and increased particle-electron collision frequency, which facilitates LTE. One can estimate  the electron density, , above which local equilibrium conditions prevail:
Here, T is the electron temperature and E is the transition energy in units of electron volts. For the Swan bands of C2 typically at 500 nm, and with an electron temperature of 10 000 K, LTE is reached whenever is of greater magnitude than . This limit is easily surpassed for focused laser irradiation, as can also be seen from plasma diagnostic using atomic line profiles (see Figure 4).
In our analysis , the rotational line-strength factor, also called the Hönl-London factor (HLF), is one of the more important terms in the equations for diatomic line intensities. The relative intensities of diatomic spectral lines are primarily controlled by the state population density and the HLF. Computed diatomic spectra are of great value because a comparison of a synthetic spectrum with a recorded spectrum allows one to compare experiment with theory. A comprehensive discussion on spectroscopy of carbon containing diatomic molecules addresses aspects of our analysis of measured diatomic molecular spectra .
Figure 4 shows the recorded spectrum of the H-beta line profile at a time delay of 10 , using a 10 gate. The source of hydrogen is water content in the flowing gas sample. The Stark width and the separation of the double peak structure is a measure for the free electron number density [23–26].
In our Nd:YAG laser experiments, repeated laser ablation of the stationary target reduced the plasma intensity and influenced appearance of kinetic curves for emission. Subsequent laser pulses usually harden the graphite surface, consequently, kinetic spectra were taken from a well-exposed target. Photographs of the plasma were captured using a standard Kodak digital camera. The photographs indicate that light emission occurred from three distinct sections: (i) from the plasma plume, (ii) from an illuminated laser beam, and (iii) from the ground glass vacuum joint of the LIBS cuvette. Spectra were recorded only from the laser plume area with the aid of a small fused silica focusing lens at the entry point of the fiber.
The kinetic spectral mode of the ICCD spectrograph indicated a strong time dependence of the carbon plasma emission. Figures 5(a)–5(e) show typical results of C3 broad 400 nm feature Swings bands and C2 Swan bands 450–575 nm , 0, progression. These emission spectra were recorded following 355 nm Nd:YAG laser ablation of graphite in 70 Torr argon. A gate width of ns was used at the indicated time delays of 250, 370, 490, 610, and 730 ns.
The decay of the Swan emission was further investigated versus time delay . Figure 6 shows the 0-0 band peaks/detector-counts above baseline, and it also shows an exponential fit. The fitted decay time amounts to ns.
Vibration-rotation temperatures for the Swan band lines for the C2 radical were inferred using an exact spectral simulation method [19, 20] that employs computed line strength values [15, 16]. Figure 7 shows results for the temperature versus time delay from optical breakdown, and it also shows error bars. The vertical bars were estimated from fitting the , , progression. The horizontal bars indicate the intensifier gate width.
4. Discussion of Experiments
The spectra recorded in our CO2 laser experiments allow us to infer an electron number density of from the hydrogen-beta Stark-broadened profiles that were measured 10 after plasma generation. Earlier in the plasma decay for time delays of a few microseconds, the profiles indicate higher electron number density in excess of . In comparison with nominal Nd:YAG laser-induced optical breakdown , we find that such free-electron number densities occur at time delay when using nominal nanosecond Nd:YAG laser radiation for optical breakdown.
The spatial extent of the focal spot diameter for 10.6 CO2 radiation is at least than for 1.064 Nd:YAG radiation, as can be computed from Gaussian beam propagation: 2. The focal volume established with the laser will also be at least than the one established with the Nd:YAG laser, as can be expected from the Gaussian beam result: 19.4. Obviously, the actual focal spot diameter for 10.6 radiation differs from Gaussian focusing due to typical spatial variations of the irradiance found for CO2 laser beams. For a 200 focal radius the focal volume is estimated to be approximately 10−3 cm3. Also, for a focal radius of 100 and 100 ns pulses the peak irradiance amounts to 1 GW/cm2 for the 10.6 CO2 laser radiation. In comparison, a peak irradiance of typically 100–1000 GW/cm2 is achieved for optical breakdown in air with nominal nanosecond 1.064 Nd:YAG laser radiation.
The recombination emission spectra from diatomic species, for example, CN or C2, clearly indicate temperatures in excess of 7000 Kelvin. Studies of the CO2 TEA laser-induced microplasmas show highly excited, high-temperature molecular transitions several tens of microseconds after plasma generation, mixed with signatures of Stark-broadened atomic lines. Our initial investigations reveal that the spectra recorded from the 10.6-micron laser-radiation-induced plasma show in principle similar plasma decay features as for the 1.064-micron Nd:YAG laser-generated plasma. Noteworthy is the spatially larger CO2 laser generated plasma, with a homogenous appearance due to extended, approx 100 ns laser excitation/pumping. Spectroscopic fitting of superposed emission spectra will be of future interest, such as fitting of the C2 Deslandres D’Azambuja system near the 4-4 band of the CN sequence of the CN Violet System.
We explored the possibility of exponentially fitting the time dependence recorded in our Nd:YAG laser experiment, see Figure 7. We found that the results are influenced by the ablation alterations of the graphite target. Moreover, graphite laser ablation can show double-peaked C2 emission , due to recombinational production of C2 molecules and due to dissociative laser-ablation of fragments that may be excited C2 molecules. The recombinational and dissociative mechanism is expected to favorably occur for higher and lower laser irradiance, respectively. The irradiance levels used in our time-resolved experiments appear to induce both processes. Revival of measured C2 emission signals seems to be indicated in Figure 5(d). This minute revival may lead to elevated temperature for that time delay as illustrated in Figure 7. Increased C2 emissions can also be caused by expansion dynamics in laser ablation . It will be of interest in future work to address time-dependent revival of C2 signals. However, the time decay fit illustrated in Figure 6 appears to contraindicate double peaking in temporal behavior. Moreover, double peaking may favorably occur only for longer time scale. It will also be of interest to investigate effects of the ablation modifications of the graphite target on the plasma emissions during our measurements.
The deposited ablation products on the walls of the quartz LIBS cell were further studied with visible absorption spectroscopy. For these measurements, we used the model Ocean Optics SP2000 miniature spectrometer. Separate deposits on small copper grids were studied employing transmission electron microscopy (TEM).
The absorption spectrum of the deposit is primarily due to amorphous carbon. The spectra exhibited no distinct features, only a gradual decrease of light absorption from the blue to the red regions. TEM pictures reinforced this conclusion as they show a uniform noncrystalline deposit devoid of special features.
The weakness of C3 emissions in our spectra is somewhat surprising as it has been known for a long time that the C3 radical is a dominant component of thermal equilibrium carbon vapors . Some likely explanations for the weak emissions can be found by considering nonequilibrium conditions , electron collision cross-sections [31, 32] and C3 cation formation  and slow intersystem crossing from triplet C3 to singlet C3  in laser ablation plasma.
The time-resolved LIB spectra of laser-ablated graphite show a sharp temporal decay of the C2 Swan bands, as well as under some circumstances, the appearance of broad features near 400 nm. The 400 nm feature in the laser-induced plasma spectra is likely due to occurrence of the Swings bands of the C3 radical. Theoretical calculations in recent literature suggest [31, 32] that the electron collision excitation cross-section of the C3 Swings band upper singlet state is smaller by a factor of about 35 than the excitation cross-section of the upper triplet state of the C2 Swan bands, therefore, the plasma conditions favor C2 spontaneous emission over C3 spontaneous emission.
In addition, regenerative sooting discharge studies  indicate that in electric discharges, similar in nature to laser-induced plasmas, the C3 cation () dominates, so it is likely that in our spectra the neutral C3 emission is weakened also by electron collision ionization of C3. The use of lower laser irradiance, by unfocusing the beam and by reducing the energy per pulse, results in diminished electron number density and causes at least a three-fold effect: firstly, local thermodynamic equilibrium (LTE) conditions are not achieved; secondly, electronic excitation of the C3 Swings bands is lessened; thirdly, and most importantly, C3 cation formation is lessened, thus it becomes more favorable to observe neutral C3 emission in our plasma. Another factor that may contribute to the weakness of C3 Swings emission is that C3 arises in the plasma in triplet state and the conversion from this triplet to the upper state singlet state in Swings emission is a very slow, millisecond process .
Measurement and analysis of time-resolved emission spectra following laser-induced optical breakdown (LIOB) allows us to characterize microplasmas that contain carbon. Presence of hydrogen in the sample will result in Balmer series emissions, and particularly the emission will be useful for determining plasma parameters. Detailed analysis of selected diatomic spectra, such as the Swan spectra, is accomplished by using line-strength files. Applications of our methods will be useful in the study of laser-induced breakdown phenomena, for example, laser-induced combustion and/or combustion augmentation, or laser-initiated nanoparticle generation.
This work is in part supported by UTSI’s Center for Laser Applications and by the Hungarian Basic Research Foundation (OTKA), under contract no. T046271. The time-resolved spectroscopic equipment in the Research Institute for Chemistry of the Hungarian Academy of Sciences was supported by the European Infrastructure Project KMA-KPI GVOP-3.2.1.- 2004-04-0059/3.0.
- A. W. Miziolek, V. Palleschi, and I. Schechter, Eds., Laser Induced Breakdown Spectroscopy, Cambridge University Press, New York, NY, USA, 2006.
- D. E. Cremers and L. J. Radziemski, Handbook of Laser-Induced Breakdown Spectroscopy, John Wiley & Sons, New York, NY, USA, 2006.
- J. Pender, B. Pearman, J. Scaffidi, S. R. Goode, and S. M. Angel, “Laser-induced breakdown spectroscopy using sequential laser pulses,” in Laser Induced Breakdown Spectroscopy, A. W. Miziolek, V. Palleschi, and I. Schechter, Eds., chapter 15, Cambridge University Press, New York, NY, USA, 2006.
- D. K. Killinger, S. D. Allen, R. D. Waterbury, C. Stefano, and E. L. Dottery, “Enhancement of Nd: YAG LIBS emission of a remote target using a simultaneous laser pulse,” Optics Express, vol. 15, no. 20, pp. 12905–12915, 2007.
- R. D. Waterbury, A. Pal, D. K. Killinger, J. Rose, E. L. Dottery, and G. Ontai, “Standoff LIBS measurements of energetic materials using a 266nm excitation laser,” in Chemical, Biological, Radiological, Nuclear, and Explosives Sensing IX (CBRNE '08), vol. 6954 of Proceedings of SPIE, The Optical Society of America, Orlando, Fla, USA, March 2008, CD-ROM, paper LThC4.
- M. Weidman, S. Palanco, and M. Richardson, “Thermodynamic properties of a Nd: YAG-double pulse laser-induced plasma and its applications under stand-off conditions,” in Proceedings of the Applications to Chemical, Security and Environmental Analysis (LACSEA '08), The Optical Society of America, St. Petersburg, Fla, USA, March 2008, CD-ROM, paper LThC5.
- C. López-Moreno, S. Palanco, J. J. Laserna et al., “Test of a stand-off laser-induced breakdown spectroscopy sensor for the detection of explosive residues on solid surfaces,” Journal of Analytical Atomic Spectrometry, vol. 21, no. 1, pp. 55–60, 2006.
- A. W. Miziolek, F. C. DeLucia, C. A. Munson, and J. L. Gottfried, “Recent progress in LIBS-based technologies for security applications,” in Proceedings of the Applications to Chemical, Security and Environmental Analysis (LACSEA '08), The Optical Society of America, St. Petersburg, Fla, USA, March 2008, CD-ROM, paper LThC1.
- M. Z. Martin, M.-D. Cheng, and R. C. Martin, “Aerosol measurement by laser-induced plasma technique: a review,” Aerosol Science and Technology, vol. 31, no. 6, pp. 409–421, 1999.
- P. I. Singh, “Experiments on Aerosol Induced Breakdown by Pulsed Lasers,” in Proceedings of the 17th AIAA Aerospace Sciences Meeting, New Orleans, La, USA, January 1979, AIAA paper 1979-0249.
- C. G. Parigger, “Laser-induced breakdown in gases: experiments and simulation,” in Laser Induced Breakdown Spectroscopy, A. W. Miziolek, V. Palleschi, and I. Schechter, Eds., chapter 4, Cambridge University Press, New York, NY, USA, 2006.
- C. G. Parigger, J. O. Hornkohl, A. M. Keszler, and L. Nemes, “Measurement and analysis of atomic and diatomic carbon spectra from laser ablation of graphite,” Applied Optics, vol. 42, no. 30, pp. 6192–6198, 2003.
- L. Nemes, A. M. Keszler, J. O. Hornkohl, and C. G. Parigger, “Laser-induced carbon plasma emission spectroscopic measurements on solid targets and in gas-phase optical breakdown,” Applied Optics, vol. 44, no. 18, pp. 3661–3667, 2005.
- J. O. Hornkohl, C. Parigger, and J. W. L. Lewis, “Temperature measurements from CN spectra in a laser-induced plasma,” Journal of Quantitative Spectroscopy and Radiative Transfer, vol. 46, no. 5, pp. 405–411, 1991.
- J. O. Hornkohl and C. G. Parigger, Boltzmann Equilibrium Spectrum Program(BESP), The University of Tennessee Space Institute, Tullahoma, Tenn, USA, 2002, http://view.utsi.edu/besp/.
- J. O. Hornkohl, C. G. Parigger, and L. Nemes, “Diatomic Hönl-London factor computer program,” Applied Optics, vol. 44, no. 18, pp. 3686–3695, 2005.
- H. S. Park, S. H. Nam, and S. M. Park, “Time-resolved optical emission studies on the laser ablation of a graphite target: The effects of ambient gases,” Journal of Applied Physics, Article ID 113103, 2005.
- “Atomic emission line data are obtained using the NIST Atomic Spectra Database Lines Form,” http://physics.nist.gov/PhysRefData/ASD/.
- C. Parigger, D. H. Plemmons, J. O. Hornkohl, and J. W. L. Lewis, “Spectroscopic temperature measurements in a decaying laser-induced plasma using the Swan system,” Journal of Quantitative Spectroscopy and Radiative Transfer, vol. 52, no. 6, pp. 707–711, 1994.
- L. Nemes, A. M. Keszler, J. O. Hornkohl, and C. G. Parigger, “Laser-induced carbon plasma emission spectroscopic measurements on solid targets and in gas-phase optical breakdown,” Applied Optics, vol. 44, no. 18, pp. 3661–3667, 2005.
- A. Thorne, U. Litzén, and S. Johansson, Spectrophysics, Principles and Applications, Springer, Berlin, Germany, 1999.
- J. O. Hornkohl, L. Nemes, and C. G. Parigger, “Spectroscopy, dynamics and molecular theory of carbon plasmas and vapors,” in Advances in the Understanding of the Most Complex High-Temperature Elemental System, L. Nemes and S. Irle, Eds., World Scientific, 2009.
- R. Zĭkić, M. A. Gigosos, M. Ivković, M. Á. González, and N. Konjević, “A program for the evaluation of electron number density from experimental hydrogen balmer beta line profiles,” Spectrochimica Acta Part B, vol. 57, no. 5, pp. 987–998, 2002.
- C. G. Parigger, D. H. Plemmons, and E. Oks, “Balmer series measurements in a laser-induced hydrogen plasma,” Applied Optics, vol. 42, no. 30, pp. 5992–6000, 2003.
- C. G. Parigger, M. Dackman, and J. O. Hornkohl, “Time-resolved spectroscopy measurements of hydrogen-alpha, -beta, and -gamma emissions,” Applied Optics, vol. 47, no. 31, pp. G1–G6, 2008.
- S. Djurovic, M. Cirišan, A. V. Demura et al., “Measurements of Stark central asymmetry and its analysis through standard theory and computer simulations,” Physical Review E, vol. 79, no. 4, Article ID 046402, 14 pages, 2009.
- S. S. Harilal, “Expansion dynamics of laser ablated carbon plasma plume in helium ambient,” Applied Surface Science, vol. 172, no. 1-2, pp. 103–109, 2001.
- S. S. Harilal, C. V. Bindhu, M. S. Tillack, F. Najmabadi, and A. C. Gaeris, “Internal structure and expansion dynamics of laser ablation plumes into ambient gases,” Journal of Applied Physics, vol. 93, no. 5, pp. 2380–2388, 2003.
- A. V. Orden and R. J. Saykally, “Small Carbon Clusters: Spectroscopy, Structure and Energetics,” Chemical Reviews, vol. 98, p. 2313, 1998.
- L. Nemes, A. M. Keszler, C. G. Parigger, J. O. Hornkohl, H. A. Michelsen, and V. Stakhursky, “Spontaneous emission from the radical in carbon plasma,” Applied Optics, vol. 46, no. 19, pp. 4032–4040, 2007.
- G. Halmová, J. D. Gorfinkiel, and J. Tennyson, “Low-energy electron collisions with using the R-matrix method,” Journal of Physics B, vol. 39, no. 12, pp. 2849–2860, 2006.
- H. Munjal and K. L. Baluja, “Elastic and excitation processes of electron impact on using the R-matrix method,” Journal of Physics B, vol. 39, no. 16, pp. 3185–3198, 2006.
- S. A. Janjua, M. Ahmad, S. Khan, R. Khalid, and S. Ahmad, “ as the dominant carbon cluster in high pressure discharges in graphite hollow cathodes,” Journal of Physics D, vol. 40, pp. 1416–1421, 2007.
- G. Larsson, The C3 molecule: A Literature Study and Spectroscopic Investigations in Flames and on Graphite, Lund Reports on Combustion Physics, LRCP-82, Lund Institute of Technology, Lund, Sweden, 2002.
Copyright © 2010 Christian G. Parigger 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.