Evaluation of XRF Spectra from Basics to Advanced Systems

Total Page:16

File Type:pdf, Size:1020Kb

Evaluation of XRF Spectra from Basics to Advanced Systems Joint ICTP-IAEA School on Novel Experimental Methodologies for Synchrotron Radiation Applications in Nano-science and Environmental Monitoring Evaluation of XRF Spectra from basics to advanced systems Piet Van Espen [email protected] 20 Nov 2014 1 Content 1. Introduction: some concepts 2. Simple peak integration 3. Method of least squares 4. Fitting of x-ray spectra 5. Improvements to the model 6. Final remarks 2 Some basic concepts Nature and science: a personal view Nature = Signal + Noise Science = Model + Statistics 3 Why spectrum evaluation? element concentrations ⇔ net intensity of fluorescence lines interference free continuum corrected But: frequent peak overlap presence of a continuum Especially in energy-dispersive spectra 4 4 Information content of a spectrum Spectrum contains • Information: energy and intensity of x-rays the signal • Amplitude noise: due to Poisson statistics ► fluctuations in the spectrum the noise • Energy noise: finite resolution of the detector ► nearly Gaussian peaks with a width of ~160 eV 5 Amplitude noise Counting events involves Poisson statistics Poisson probability density function: The probability to observe N counts if the true number is µ Property: Poisson distribution Normal distribution ≅ µ µ and σ2 = µ approximation is very good for µ (or N) ≥ 9 Poisson : P(N | µ=3) Normal : P(x | µ=3 σ2 = 3) 6 2 Spectrum analysis Energy Noise Resolution of ED-XRF spectrometers 2.35⇥ F E (9) Full Width at Half Maximum× (FWHM)× of a peak 2 2 2 FWHMPeak = FWHMElec + FWHMDet (10) 2 Spectrum analysis Electronic noise Intrinsic contribution ~100 eV 2.35⇥ F E (9) × × FWHM2 = FWHM2 +Det2 (10) Peakϵ energy to createElec e-h pairPealk 3.85 eV F Fano factor ~0.114 E x-ray energy in eV Mn Kα @ 5.895 keV FWHMDet = 120 eV FWHMElec = 100 eV => FWHMPeak = 156 eV 7 2 2 Energy Noise Resolution of ED-XRF spectrometers Cr – Mn – Fe overlap at ~20 eV Cr – Mn – Fe overlap at ~160 eV 8 Without amplitude noise (counting statistics) there would be NO PROBLEM But it is part of the nature We can only measure longer or with a more efficient system Without energy noise there would be LITTLE PROBLEM The natural line width of x-rays is only a few eV!!! The observed peak width is the result of the detection process with a fundamental limitation imposed by the Fano factor 9 Information content of a spectrum If no energy noise or no amplitude noise ► could determine the “information” unambiguously Need methods to extract information in a optimum way These methods rely on “addition” information (knowledge) to extract the useful information Not the method itself is important (if implemented correctly) but the correctness of the additional information. 10 10 2. Simple peak integration 11 Simple peak integration Estimate Uncertainty As good as is can get We have to make assumptions if the assumptions (model) are correct! integration limits linear background no interference 12 3. Method of Least Squares 13 Method of least squares Need to “estimate” the net peak area with highest possible • correctness (no systematic error) • precision (smallest random error) Least-squares estimation (fitting): • unbiased • minimum variance Limiting factors: • counting statistical fluctuations (precision) • accuracy of the fitting model 14 Method of least squares, straight line 2 SpectrumData: analysis {xi,yi}, i=1, 2, …, N 2 SpectrumModel: y(i analysis) = b + b x noise 2.35p0 1F i E (9) ⇥ ⇥ 2.35p F E (9) Fitting the model: ⇥ ⇥ estimating2 b0 and2 b1 2 FWHMPeak = FWHMElec + FWHMDet (10) 2 2 2 FWHMPeak = FWHMElec + FWHMDet (10) Criterion: SS = [y y(i)]2 = [y b b x ]2 = min (11) i − i − 0 − 1 i i SS = [y y(ii)]2 = [y b b x ]2 = min (11) X i − X i − 0 − 1 i Xi Xi ⇤SS Set of 2 equations in 2 unknowns b and b =0⇤SS yi = b0n + b1 xi (12) 0 1 ⇤b0 ! =0i yi = b0in + b1 xi Normal equations(12) ⇤b X ! X 0 i i ⇤SS X X 2 =0 xiyi = b0 xi + b1 x (13) ⇤b ⇤SS! i Direct2 analytical solution 1 =0i xiyi = b0 xii + b1 xi (13) ⇤b1 X ! X X Xi X Xi 15 2 2 4. Fitting X-ray Spectra 16 Least-squares estimate of x-ray spectrum parameters position area 1 (x x )2 Peak described by a Gaussian y(i)=b + A exp i − p σp2⇡ 2σ2 linear parameter continuum width non-linear n2 2 1 1 2 Criterion, agreement between model and data χ = [y(i) yi] (1) ⌫ wi − iX=n1 y(i)=yCont(i)+ Aj RijP (i, Eij) (2) " # Minimum:Xj Xk No direct analytical solution 2 Gain 2 1 E(i) Ejk P = G(i, Ejk)=Search χ exp for minimum − (3) σ p2⇡ −2 σjk jk " ✓ ◆ # E(i)=Zero + Gain i (4) ⇥ 17 1/2 Noise 2 σ = + ✏ FanoE (5) jk p jk "✓2 2 ln 2◆ # Gain E(i) Ejk E(i) Ejk 1 T (i, Ejk)= exp − erfc − + 2βσ exp 1 βσ σp2 p2β − 2β2 h i Gain E(i) Ejk S(i, Ejk)= erfc − 2Ejk σp2 P (i, Ejk)=G(i, Ejk)+fSS(i, Ejk)+fT T (i, Ejk) (6) fS(Ejk)=µDet(Ejk)(a0 + a1Ejk) fTK↵(Ejk)=b0 + b1µDet(Ejk) fTKβ(Ejk)=c0 + c1µDet(Ejk) β(Ejk)=d0 + d1Ejk Gain E1(i) E ↵ K − jk , L σjkp2⇡ σjkp2 2σjkp2! K(x, y)=Re exp( z2)erfc( iz) z = x + iy (7) − − ⇥ ⇤ 1 We can still apply the concept of least squares minimising the square of the differences between the model and the data 2 2 1 1 2 χ = χ (b, h, x0,w)= [yi y(xi,b,h,x0,w)] ⌫ i yi − P The sum of squares is a function of the values of the parameters and for a given set of values should be minimum In this case SS describes a 4 dimensional hyper-surface in a 5- dimension space We can only “see” in 3-dimensions but mathematically we can search in a higher dimensional space to locate the minimum h Starting from some initial values we can modify the parameter values until the minimum is reached. w 18 General form of such a search algorithm 1. Select starting values for all parameters bj χ = χ(b) and calculate the ch-square 2. Obtain (calculate, guess...) a change (increment or decrement) Δbj such that one χ(b +∆b) <χ(b) moves towards the minimum: 3. Replace the old parameter values with the new ones b ← b + Δb 4. repeat step 2 until the “true” minimum is found Iterative process AXIL = Analysis of X-ray spectra by Iterative Least-squares 19 n2 2 1 1 2 χ = [y(i) yi] (1) ⌫ wi − iX=n1 Analytically important parameters: net peak areas y(i)=yCont(i)+ Aj RijP (i, Eij) (2) " # Xj Xk Statistical optimal estimate: 2 using correct weight (PoissonGain statistics w1 i =E yi()i) Ejk P = G(i, Ejk)= exp − (3) σ p2⇡ −2 σjk jk " ✓ ◆ # In general y(i) is non-linear → Marquardt – Leverberg algorithm Gradient search ↔ linearisationE(i)=Zero + Gain i (4) Reliable error estimated ⇥ But unstable 1/2 Noise 2 σ = + ✏ FanoE (5) jk p jk "✓2 2 ln 2◆ # Gain E(i) Ejk E(i) Ejk 1 T (i, Ejk)= exp − erfc − + 2βσ exp 1 βσ20 σp2 p2β − 2β2 h i Gain E(i) Ejk S(i, Ejk)= erfc − 2Ejk σp2 P (i, Ejk)=G(i, Ejk)+fSS(i, Ejk)+fT T (i, Ejk) (6) fS(Ejk)=µDet(Ejk)(a0 + a1Ejk) fTK↵(Ejk)=b0 + b1µDet(Ejk) fTKβ(Ejk)=c0 + c1µDet(Ejk) β(Ejk)=d0 + d1Ejk Gain E(i) E ↵ K − jk , L σjkp2⇡ σjkp2 2σjkp2! K(x, y)=Re exp( z2)erfc( iz) z = x + iy (7) − − ⇥ ⇤ 1 10 peaks ⇒ > 30 parameters !!!! WANT WORK in practice!!!! Need parsimony!!! But we can do better 2 We known the energies of the x-ray lines (in most cases)1 (xi xp) y(i)=b +1 A (expx x )−2 Where they are depends on the energy calibration σp2⇡ i p 2σ2 y(i)=b + A exp − 2 2 (the same applies to the width of the peaks: resolutionσp2⇡ 1 calibration)2(σxi xp) y(i)=b + A exp − 2 σp2⇡ 2σ 2 Gain (Ei E(i)) ⇒ Add additional information to the modelG(i, E)= exp − 2 Gainσ(E)p2⇡(Ei E(2iσ))2(E) G(i, E)= Gainexp −(E E)2 Gaussian peak shape G(i,σ E(E)=)p2⇡ exp2σ2(Ei −) σ(E)p2⇡ 2σ2(E) E(i)=Zero + Gain i ⇥ Energy relation: E(i)=Zero + Gain i ⇥ 1/2 Noise 2 Resolution relation: σ(E)= 2 + ✏FanoE1/2 Noise" 2p2 ln 2 # σ(E)= ✓ + ✏◆FanoE p Only 4 non-linear parameters"✓ 2 2 ln 2◆ # For 10 peaks only 14 parameters 21 1 1 1 1 (x x )2 y(i)=b + A exp i − p σp2⇡ 2σ2 2 1 (xi xp) 2 y(i)=bGain+ A exp(E E−(i))2 G(i, E)= σpexp2⇡ i − 2σ σ(E)p2⇡ 2σ2(E) Gain (E E(i))2 G(i, E)= exp i − E(i)=σ(ZeroE)p2+⇡Gain i 2σ2(E) ⇥ Already better, but we know more 1/2 E(iNoise)=Zero2 + Gain i σ(E)= + ✏FanoE⇥ p We know (to some extend)"✓2 2 lnthe 2 ◆ratio between# lines of an element 1/2 Noise 2 σ(E)=IK↵2 IKβ + ✏FanoE ,p ... I"K✓↵12 2I lnK↵ 2◆ # I I We can group lines togetherK↵2 (“peakgroup”), Kβ ... with one “area” and fixed intensity ratiosIK↵1 IK↵ y(i)=yCont(i)+ Aj RjkP (i, Ejk) ( ) Xj Xk Continuum Area Line Peak function ratio shape for j elements (or peak groups) 10 elements ⇒ 10 Area’s + 4 calibration parameters 22 1 1 Further refinements: escape peaks Known position (energy) intensity (escape probability) 23 Sum peaks Known position relative intensity 24 And more Different background models polynomial exponential polynomial Bremsstrahlung background filter background Parameter constraining and more..
Recommended publications
  • Field Guide to Continuous Probability Distributions
    Field Guide to Continuous Probability Distributions Gavin E. Crooks v 1.0.0 2019 G. E. Crooks – Field Guide to Probability Distributions v 1.0.0 Copyright © 2010-2019 Gavin E. Crooks ISBN: 978-1-7339381-0-5 http://threeplusone.com/fieldguide Berkeley Institute for Theoretical Sciences (BITS) typeset on 2019-04-10 with XeTeX version 0.99999 fonts: Trump Mediaeval (text), Euler (math) 271828182845904 2 G. E. Crooks – Field Guide to Probability Distributions Preface: The search for GUD A common problem is that of describing the probability distribution of a single, continuous variable. A few distributions, such as the normal and exponential, were discovered in the 1800’s or earlier. But about a century ago the great statistician, Karl Pearson, realized that the known probabil- ity distributions were not sufficient to handle all of the phenomena then under investigation, and set out to create new distributions with useful properties. During the 20th century this process continued with abandon and a vast menagerie of distinct mathematical forms were discovered and invented, investigated, analyzed, rediscovered and renamed, all for the purpose of de- scribing the probability of some interesting variable. There are hundreds of named distributions and synonyms in current usage. The apparent diver- sity is unending and disorienting. Fortunately, the situation is less confused than it might at first appear. Most common, continuous, univariate, unimodal distributions can be orga- nized into a small number of distinct families, which are all special cases of a single Grand Unified Distribution. This compendium details these hun- dred or so simple distributions, their properties and their interrelations.
    [Show full text]
  • Detection of Voigt Spectral Line Profiles of Hydrogen Radio
    Detection of Voigt Spectral Line Profiles of Hydrogen Radio Recombination Lines toward Sagittarius B2(N) Azrael A. von Proch´azka1,2, Anthony J.Remijan3,4, Dana S. Balser3, Robert S. I. Ryans1, Adele H. Marshall5 Fredric R.Schwab3, Jan M. Hollis6, Philip R. Jewell3,4, and Frank J. Lovas4,7 ABSTRACT We report the detection of Voigt spectral line profiles of radio recombination lines (RRLs) toward Sagittarius B2(N) with the 100-m Green Bank Telescope (GBT). At radio wavelengths, astronomical spectra are highly populated with RRLs, which serve as ideal probes of the phys- ical conditions in molecular cloud complexes. An analysis of the Hnα lines presented herein shows that RRLs of higher principal quantum number (n > 90) are generally divergent from their expected Gaussian profiles and, moreover, are well described by their respective Voigt profiles. This is in agreement with the theory that spectral lines experience pressure broadening as a result of electron collisions at lower radio frequencies. Given the inherent technical dif- ficulties regarding the detection and profiling of true RRL wing spans and shapes, it is crucial that the observing instrumentation produce flat baselines as well as high sensitivity, high reso- lution data. The GBT has demonstrated its capabilities regarding all of these aspects, and we believe that future observations of RRL emission via the GBT will be crucial towards advancing our knowledge of the larger-scale extended structures of ionized gas in the interstellar medium (ISM). Subject headings: ISM: H II regions — ISM: radio continuum — ISM: radio recombination lines — radiation mechanisms: nonthermal — radiation mechanisms: pressure broadening — line: profiles 1.
    [Show full text]
  • Integro-Differential Equations to Model Generalized Voigt Profiles
    1 Integro-differential equations to model generalized Voigt profiles: a fractional diffusion approach Gianni Pagnini, Francesco Mainardi CRS4, Polaris Bldg. 1, I-09010 Pula (CA), Italy E-mail: [email protected] Department of Physics, University of Bologna, and INFN, via Irnerio 46, I-40126 Bologna, Italy E-mail: [email protected] URL: www.fracalmo.org Summary. In molecular spectroscopy and atmospheric radiative transfer, the combined effects of Doppler and pressure broadenings lead to the Voigt profile function, which turns out to be the convo- lution of the Gaussian (due to the Doppler broadening) and the Lorentzian (due to the pressure broad- ening) distributions. Here we are interested to study the Voigt profile function when the widths are not constant but depending on a scale-factor with a power law. Since in probability theory the Gaussian and the Lorentzian distributions are known to belong to the class of symmetric L´evy stable distributions, in this framework we propose to generalize the Voigt function by adopting the convolution of two arbitrary symmetric L´evy distributions. Moreover, we provide the integro-differential equations with respect to the scale-factor satisfied by the generalized Voigt profiles. These evolution equations can be interpreted as space-fractional diffusion equations of double order. Key words: Voigt profile function, L´evy stable distributions, integro-differential equations, Riesz space-fractional derivative, Mellin-Barnes integrals 2 Gianni Pagnini, Francesco Mainardi 1 Introduction In the present paper some mathematical aspects of the Voigt profile are discussed. This func- tion emerges as a profile spectrum in molecular spectroscopy and atmospheric radiative trans- fer due to the combined effects of Doppler and pressure broadeningsand it turns out to be the convolution of the Gaussian (due to the Doppler broadening) and the Lorentzian (due to the pressure broadening) probability densities.
    [Show full text]
  • Detecting Damped Ly&Thinsp;Α Absorbers with Gaussian Processes
    MNRAS 472, 1850–1865 (2017) doi:10.1093/mnras/stx1958 Advance Access publication 2017 August 2 Detecting damped Ly α absorbers with Gaussian processes Roman Garnett,1‹ Shirley Ho,2 Simeon Bird3 and Jeff Schneider4 1Department of Computer Science and Engineering, Washington University in St Louis, One Brookings Drive, St Louis, MO 63130, USA 2Department of Physics, Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, PA 15213, USA 3Department of Physics & Astronomy, Johns Hopkins University, 3400 N. Charles Street, Baltimore, MD 21218, USA 4School of Computer Science, Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, PA 15213, USA Accepted 2017 July 28. Received 2017 July 28; in original form 2016 November 19 ABSTRACT We develop an automated technique for detecting damped Ly α absorbers (DLAs) along spectroscopic lines of sight to quasi-stellar objects (QSOs or quasars). The detection of DLAs in large-scale spectroscopic surveys such as SDSS III sheds light on galaxy formation at high redshift, showing the nucleation of galaxies from diffuse gas. We use nearly 50 000 QSO spectra to learn a novel tailored Gaussian process model for quasar emission spectra, which we apply to the DLA detection problem via Bayesian model selection. We propose models for identifying an arbitrary number of DLAs along a given line of sight. We demonstrate our method’s effectiveness using a large-scale validation experiment, with excellent performance. We also provide a catalogue of our results applied to 162 858 spectra from SDSS-III data release 12. Key words: methods: statistical – intergalactic medium – quasars: absorption lines – galaxies: statistics. had some success in detecting large DLA catalogues, their reliance 1 INTRODUCTION entirely on templates made them subject to hard-to-quantify sys- The damped Ly α systems (DLAs; Wolfe et al.
    [Show full text]
  • Mellin-Barnes Integrals for Stable Distributions and Their Convolutions
    MELLIN-BARNES INTEGRALS FOR STABLE DISTRIBUTIONS AND THEIR CONVOLUTIONS Francesco Mainardi 1 and Gianni Pagnini 2 Dedicated to Professor Anatoly A. Kilbas Dept. of Mathematics and Mechanics, University of Minsk, Belarus, on the occasion of his 60th birthday (20 July 2008) Abstract In this expository paper we ¯rst survey the method of Mellin-Barnes integrals to represent the ® stable L¶evydistributions in probability theory (0 < ® · 2). These integrals are known to be useful for obtaining conver- gent and asymptotic series representations of the corresponding probability density functions. The novelty concerns the convolution between two stable probability densities of di®erent L¶evyindex, which turns to be a probabil- ity law of physical interest, even if it is no longer stable and self-similar. A particular but interesting case of convolution is obtained combining the Cauchy-Lorentz density (® = 1) with the Gaussian density (® = 2) that yields the so-called Voigt pro¯le. Our machinery can be applied to derive the fundamental solutions of space-fractional di®usion equations of two or- ders. 2000 Mathematics Subject Classi¯cation: 26A33, 33C60, 42A38, 44A15, 44A35, 60G18, 60G52 Key Words and Phrases: Mellin-Barnes integrals, Mellin transforms, stable distributions, Voigt pro¯le, space-fractional di®usion 444 F. Mainardi, G. Pagnini Introduction Mellin-Barnes integrals are a family of integrals in the complex plane whose integrand is given by the ratio of products of Gamma functions. De- spite of the name, the Mellin-Barnes integrals were initially introduced in 1888 by the Italian mathematician S. Pincherle [27] in a couple of papers on the duality principle between linear di®erential equations and linear dif- ference equations with rational coe±cients, as discussed by the Authors in [20].
    [Show full text]
  • Deconvolution of Spectral Voigt Profiles Using Inverse Methods And
    Deconvolution of Spectral Voigt Profiles Using Inverse Methods and Fourier Transforms Genia Vogman Senior Thesis June 2010 Contents 1 Introduction 1 2 Light Spectrum and Individual Line Profiles 2 2.1 Gaussian Profile and Doppler Broadening . 3 2.2 Lorentzian Profile and Stark Broadening . 3 2.3 Other Profile-Changing Effects . 5 2.4 Convolution of Gaussian and Lorentzian Functions . 5 3 Direct Inverse Method for Resolving Temperature and Density 5 3.1 Recovering Temperature and Density from a Voigt Fit . 5 3.2 Parameter Search to Verify Best Fit . 6 4 Fourier Transform Method 9 4.1 Fourier Transforms Applied To Convolutions . 9 4.2 Fitting Routine in the Fourier Domain . 11 4.3 Comparison to Direct Inverse Method . 13 5 Conclusions 14 A Nomenclature 17 B Matlab Code 18 B.1 Optimization Function and Parameter Space Search . 18 B.2 Convolution Function . 21 B.3 Fitting Function . 21 1 Introduction Emissive spectroscopy, which is a means of measuring light radiation emitted by high-energy par- ticles, is the basis for many experimental studies. Through spectroscopy it is possible to determine particle temperature and density by decomposing emitted light into constituent wavelengths and resolving the shapes of individual spectral profiles. This is particularly useful for the study of stel- lar and terrestrial plasmas, in which particles can reach temperatures that are so high that their properties cannot be measured by any other means. 1 Figure 1: An emission light spectrum [9] of intensity, in arbitrary units, versus wavelength from 200 nm (ultraviolet) to 700 nm (infrared). What makes spectroscopy particularly useful is the ready availability of data – that is, the retrieval of spectroscopy data is much easier than the use of other diagnostics.
    [Show full text]
  • Coupling Nuclear Induced Phonon Propagation with Conversion Electron Mössbauer Spectroscopy
    Air Force Institute of Technology AFIT Scholar Theses and Dissertations Student Graduate Works 6-18-2015 Coupling Nuclear Induced Phonon Propagation with Conversion Electron Mössbauer Spectroscopy Michael J. Parker Follow this and additional works at: https://scholar.afit.edu/etd Part of the Atomic, Molecular and Optical Physics Commons Recommended Citation Parker, Michael J., "Coupling Nuclear Induced Phonon Propagation with Conversion Electron Mössbauer Spectroscopy" (2015). Theses and Dissertations. 194. https://scholar.afit.edu/etd/194 This Thesis is brought to you for free and open access by the Student Graduate Works at AFIT Scholar. It has been accepted for inclusion in Theses and Dissertations by an authorized administrator of AFIT Scholar. For more information, please contact [email protected]. COUPLING NUCLEAR INDUCED PHONON PROPAGATION WITH CONVERSION ELECTRON MÖSSBAUER SPECTROSCOPY THESIS Michael J. Parker, Capt, USAF AFIT-ENP-MS-15-J-054 DEPARTMENT OF THE AIR FORCE AIR UNIVERSITY AIR FORCE INSTITUTE OF TECHNOLOGY Wright-Patterson Air Force Base, Ohio DISTRIBUTION STATEMENT A. APPROVED FOR PUBLIC RELEASE; DISTRIBUTION UNLIMITED. The views expressed in this thesis are those of the author and do not reflect the official policy or position of the United States Air Force, Department of Defense, or the United States Government. This material is declared a work of the United States Government and is not subject to copyright protection in the United States. AFIT-ENP-MS-15-J-054 COUPLING NUCLEAR INDUCED PHONON PROPAGATION WITH CONVERSION ELECTRON MÖSSBAUER SPECTROSCOPY THESIS Presented to the Faculty Department of Engineering Physics Graduate School of Engineering and Management Air Force Institute of Technology Air University Air Education and Training Command In Partial Fulfillment of the Requirements for the Degree of Master of Science in Nuclear Engineering Michael J.
    [Show full text]
  • Bayesvp: a Bayesian Voigt Profile Fitting Package
    Mon. Not. R. Astron. Soc. 000, 1–?? (2015) Printed 30 October 2017 (MN LATEX style file v2.2) BayesVP: a Bayesian Voigt profile fitting package Cameron J. Liang1? and Andrey V. Kravtsov1;2 1Department of Astronomy & Astrophysics, and Kavli Institute for Cosmological Physics, University of Chicago, Chicago IL 60637 2Enrico Fermi Institute, The University of Chicago, Chicago, IL 60637, USA 30 October 2017 ABSTRACT We introduce a Bayesian approach for modeling Voigt profiles in absorption spec- troscopy and its implementation in the python package, BayesVP , publicly available at https://github.com/cameronliang/BayesVP. The code fits the absorption line profiles within specified wavelength ranges and generates posterior distributions for the col- umn density, Doppler parameter, and redshifts of the corresponding absorbers. The code uses publicly available efficient parallel sampling packages to sample posterior and thus can be run on parallel platforms. BayesVP supports simultaneous fitting for multiple absorption com- ponents in high-dimensional parameter space. We provide other useful utilities in the pack- age, such as explicit specification of priors of model parameters, continuum model, Bayesian model comparison criteria, and posterior sampling convergence check. Key words: methods:data analysis – quasars:absorption lines 1 INTRODUCTION Therefore, uncertainties of the best-fit parameters can be under- or overestimate estimated in a coarse grid of the parameter space. Fur- Absorption line spectroscopy is undoubtedly one of the most thermore, the grid search becomes prohibitively computationally important tools for probing the physical properties of absorbing gas expensive or impossible when fitting multiple absorption compo- in a variety of astrophysical environments. Spectral analyses are nents because the number of parameters and thus the number of ubiquitous in astrophysics and have a long and rich history in stud- dimensions of parameter space is large.
    [Show full text]
  • Stochastic Methods – Definitions, Random Variables, Distributions
    Stochastic Methods – Definitions, Random Variables, Distributions Sigma-algebra In mathematics, a σ-algebra (pronounced sigma-algebra) or σ-field over a set X is a family Σ of subsets of X that is closed under countable set operations; σ-algebras are mainly used in order to define measures on X. Formally, Σ is a σ-algebra if and only if it has the following properties: 1. The empty set is in Σ. 2. If E is in Σ then so is the complement X\E of E. 3. If E1, E2, E3, ... is a (countable) sequence in Σ then their (countable) union is also in Σ. From 1 and 2 it follows that X is in Σ; from 2 and 3 it follows that the σ-algebra is also closed under countable intersections (via De Morgan's laws). Elements of the σ-algebra are called measurable sets. An ordered pair (X, Σ), where X is a set and Σ is a σ-algebra over X, is called a measurable space. Measures are defined as certain types of functions from a σ-algebra to [0,∞]. Random variables Some consider the expression random variable a misnomer, as a random variable is not a variable but rather a function that maps events to numbers. Let A be a σ-algebra and Ω the space of events relevant to the experiment being performed. In the die-rolling example, the space of events is just the possible outcomes of a roll, i.e. Ω = { 1, 2, 3, 4, 5, 6 }, and A would be the power set of Ω.
    [Show full text]
  • Signal Processing Manuscript Draft Manuscript Number
    Signal Processing Manuscript Draft Manuscript Number: SIGPRO-D-14-00375 Title: Optimum Linear Regression in Additive Cauchy-Gaussian Noise Article Type: Fast Communication Keywords: Impulsive noise, Cauchy distribution, Gaussian distribution, mixture noise, Voigt profile, maximum likelihood estimator, pseudo-Voigt function, M-estimator Abstract: In this paper, we study the estimation problem of linear regression in the presence of a new impulsive noise model, which is a sum of Cauchy and Gaussian random variables in time domain. The probability density function (PDF) of this mixture noise, referred to as the Voigt profile, is derived from the convolution of the Cauchy and Gaussian PDFs. To determine the linear regression parameters, the maximum likelihood estimator is first developed. Since the Voigt profile suffers from a complicated analytical form, an M-estimator with the pseudo-Voigt function is also derived. In our algorithm development, both scenarios of known and unknown density parameters are considered. In the unknown scenario, density parameters need to be estimated prior to proposals, by utilizing the empirical characteristic function and characteristic function. Simulation results show that the performance of both proposed methods can attain the Cram\'{e}r-Rao lower bound. Highlights (for review) Highlights: An additive mixture noise is studied in this paper and the corresponding noise PDF, i.e., the Voigt function is derived. To determine the parameters of a linear regression model, the maximum likelihood estimator (MLE) is developed, where both the scenarios of known and unknown density parameters, are considered. To reduce the computational complexity of the MLE, an M-estimator with pseudo-Voigt function is presented.
    [Show full text]
  • Theoretical Model of Diffraction Line Profiles As Combinations of Gaussian and Cauchy Distributions
    Journal of Crystallization Process and Technology, 2014, 4, 145-155 Published Online July 2014 in SciRes. http://www.scirp.org/journal/jcpt http://dx.doi.org/10.4236/jcpt.2014.43019 Theoretical Model of Diffraction Line Profiles as Combinations of Gaussian and Cauchy Distributions Girija Bhushan Mitra* Indian Association for the Cultivation of Science, Jadavpur, Kolkata, India Email: [email protected] Received 6 May 2014; revised 29 May 2014; accepted 17 June 2014 Copyright © 2014 by author and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/ Abstract Previously we derived equations determining line broadening in ax-ray diffraction profile due to stacking faults. Here, we will consider line broadening due to particle size and strain which are the other factors affecting line broadening in a diffraction profile. When line broadening in a diffrac- tion profile is due to particle size and strain, the theoretical model of the sample under study is either a Gaussian or a Cauchy function or a combination of these functions, e.g. Voigt and Pseudo- voigt functions. Although the overall nature of these functions can be determined by Mitra’s R(x) test and the Pearson and Hartley χ test, details of a predicted model will be lacking. Development of a mathematical model to predict various parameters before embarking upon the actual experiment would enable correction of significant sources of error prior to calculations. Therefore, in this study, predictors of integral width, Fourier Transform, Second and Fourth Moment and Fourth Cumulant of samples represented by Gauss, Cauchy, Voigt and Pseudovoigt functions have been worked out.
    [Show full text]
  • Etd-11062007-190817.Pdf (3.508 Mb )
    LASER INDUCED BREAKDOWN SPECTROSCOPY: INVESTIGATION OF LINE PROFILES, SLURRIES AND ARTIFICAL NEURAL NETWORK PREDICTION By Seong Yong Oh A Dissertation Submitted to the Faculty of Mississippi State University in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy in Engineering with an Emphasis in Applied Physics in the Department of Physics and Astronomy Mississippi State, Mississippi December 2007 Copyright by Seong Yong Oh 2007 LASER INDUCED BREAKDOWN SPECTROSCOPY: INVESTIGATION OF LINE PROFILES, SLURRIES AND ARTIFICAL NEURAL NETWORK PREDICTION By Seong Yong Oh Approved: _________________________________ _________________________________ Jagdish P. Singh David L. Monts Adjunct Professor of Physics and Professor and Graduate Coordinator of Research Professor of Institute for Clean Physics Energy Technology (Co-Major Professor) (Director of Dissertation and Major Professor) _________________________________ _________________________________ Fang Yu Yueh Chun Fu Su Research Scientist II of Institute for Clean Professor of Physics Energy Technology (Committee Member) (Committee Member) _________________________________ _________________________________ Thomas Philip Chuji Wang Professor of Computer Science and Assistant Professor of Physics Engineering (Committee Member) (Committee Member) ________________________________ W. Glenn Steele, Jr. Interim Dean of the Bagley College of Engineering Name: Seong Yong Oh Date of Degree: December 14, 2007 Institution: Mississippi State University Major Field: Applied Physics Major Professor: Dr. Jagdish P. Singh Title of Study: LASER INDUCED BREAKDOWN SPECTROSCOPY: INVESTIGATION OF LINE PROFILES, SLURRIES AND ARTIFICAL NEURAL NETWORK PREDICTION Pages in Study: 129 Candidate for Degree of Doctor of Philosophy Laser induced breakdown spectroscopy (LIBS) was tested to examine its applicability to remote and in suit analysis in inaccessible situation. Two types of liquid sample (slurry) prepared for simulating vitrification of liquid hazardous wastes was tested.
    [Show full text]