X-ray spectral Analysis

Size: px
Start display at page:

Download "X-ray spectral Analysis"

Transcription

1 X-ray spectral Analysis xamples Problems Algorithm Model types Statistic Fitting procedure Practical problems with Chandra CCD data Disk N(N) Disk N(S) Nucleus Plume Based partly on the lecture given by Keith Arnaud (NASA/GSFC & UMCP) Pietsch et al. (2001) XMM-Newton PC PN Spectra XMM/Newton RGS spectrum of the nucleus region + plume (ULX) Disk: thermal temperature: kev, assuming solar abundance A2125 Z=0.25 Red: kev Green: 2-4 kev Blue: 4-8 kev

2 ACS Spectrum of A2125 Unfold spectrum kt ~ 4 kev N H ~ 2x cm -2 Practical X-ray spectroscopy Can we start with this Most X-ray spectra are of moderate or low resolution (eg Chandra ACS or XMM-Newton PC). However, the spectra generally cover a bandpass of more than 1.5 decades in energy. Moreover, the continuum shape often provides important physical information. Therefore, most uses of X-ray spectra have to rely on a simultaneous analysis of the entire spectrum. and deduce this The Basic Problem Suppose we observe D() counts in channel (of N) from some source. Then : D() = T R(,) A() S() d T is the observation length (in seconds) R(,) is the spectral matrices --- the probability of an incoming photon of energy being registered in channel (dimensionless) A() is the energy-dependent effective area of the telescope and detector system (in cm 2 ) S() is the source flux at the front of the telescope (in photons/cm 2 /s/kev)

3 An example R(,) The Basic Problem escape escape photopeak photopeak fluorescence fluourescenc e D() = T R(,) A() S() d T, A() and R(,) are known. We need to solve this integral equation for S(). We can divide the energy range of interest into M bins and turn this into a matrix equation: D i = T R ij A j S j where S j is now the flux in photons/cm 2 /s in energy bin J. We want to find S j. The Basic Problem This has a solution: D i = T R ij A j S j S j =1/(T A j ) (R ij ) - But this does not work! The S j derived in this way are very sensitive to slight changes in the data D i. This is a great method for amplifying noise. 1 D i A (brief) Mathematical Digression This is known as the remote sensing problem and arises in many areas of astronomy as well as eg geophysics and medical imaging. n mathematics the integral is known as a Fredholm equation of the first kind. Tikhonov showed that such equations can be solved using regularization - applying a priori knowledge to damp the noise. A familiar example is maximum entropy but there are a host of others. Some of these have been tried on X-ray spectra - none have had any impact on the field. Forward-fitting arf rmf χ 2 This standard method of analyzing X-ray spectra comprises the following step: Calculate a model spectrum. Multiply this spectrum by the effective area. Further convolve it with an instrumental response matrix. Compare the result with the actual observed data by calculating some statistic. Modify the model spectrum and repeat till the best value of a chosen statistic is obtained. Forward-fitting (cont.) Works if the number of fitting parameters is reasonably small so that the model can be varied in some sensible fashion. The aim of the forward-fitting is then to obtain the best-fit and confidence ranges of these parameters.

4 Spectral fitting programs XSPC - part of HAsoft. General spectral fitting program with many models available. Sherpa-part of CAO. Multi-dimensional fitting program which includes the XSPC model library and can be used for spectral fitting. SPX - from SRON in the Netherlands. Spectral fitting program specializing in collisional plasmas and high resolution spectroscopy. SS - from the MT Chandra HTG group. Mainly intended for the analysis of grating data. Models All models are wrong, but some are useful! X-ray spectroscopic models are usually built up from individual components. These can be thought of as two basic types: Additive (an emission component e.g. blackbody, line, ) Multiplicative (something which modifies the spectrum e.g. absorption). Model = M 1 * M 2 * (A 1 + A 2 + M 3 *A 3 ) + A 4 Additive Models Basic additive (emission) models include : Blackbody and its variations (e.g., neutron stars) Power-law (AGNs, X-ray binaries) Collisional plasma (e.g, raymond, mekal; SNRs, the hot SM/GM) Gaussian lines Many more models: Accretion disks Comptonized plasmas Non-equilibrium ionization plasmas Multi-temperature collisional plasmas Multiplicative (and Other) Models Photoelectric absorption due to our Galaxy. Photoelectric absorption due to ionized material. High energy exponential roll-off. Other types of XSPC model components (e.g., convolution) which are used like a multiplicative model but perform more complicated operations on the current model. Roll Your Own Models Statistic. A simple XSPC model interface is available. You can write new models and fit them to your data. Specifically, ither write a subroutine that takes in the energies on which to calculate the model and writes out the fluxes (in photons/cm2/s). Or a table model - files of a standard format that contain model spectra, which can then be fit to data. χ 2 =Σ(d i -m i ) 2 /σ 2 i Convenient and goodness of fit To use the χ 2 statistic, the rule-of thumb is that each data point contains > 6 counts Thus rebinning of the data is often necessary (e.g., using the FTOOLS routine grppha). Cash Statistic C=-2log(L), where L=Π i p(d i m i ) Poisson likelihood is an unbiased estimator, whereas χ 2 can be biased, depending on the choice of No goodness of fit criterion C-statistic cannot be used to analyze backgroundsubtracted data, since the difference of two Poissondistributed variables is not Poisson distributed. The background and the source data can be fit simultaneously.

5 Finding the best-fit Finding the best-fit means minimizing the statistic value (e.g., χ 2 or C-statistic). Many algorithms available to do this in a computationally efficient fashion (see Numerical Recipes). Most methods used to find the best-fit are local i.e. they use some information around the current parameters to guess a new set of parameters. All these methods are liable to get stuck in a local minimum. Watch out for this! The more complicated your model and the more highly correlated the parameters then the more likely that the algorithm will not find the absolute best-fit. Global Minimization There are global minimization methods available - simulated annealing, genetic algorithms, - but they require many function evaluations (so are slow) and are still not guaranteed to find the true minimum. A new technique called Markov Chain Monte Carlo, which provides an intelligent sampling of parameter space, looks promising but it is not yet widely available (not yet available to XSPC). Confidence ntervals 64%, 90%, 99% contours corresponding to χ 2 = χ 2 min + (1, 2.7, and?) Final Advice and Admonitions The purpose of spectral fitting is to attain understanding, not fill up tables of numbers. Try to test whether you really have found the best-fit. X-ray data analysis is half science and half art. A fit only means that the model is consistent with the data: There are always alternative models that can fit the data equally well or better. Astrophysical intuition is absolutely important! Don t misuse the F-test. nergy resolution of the ACS CT correction for ACS- data The energy resolution of the S3 and 3 as a function of row number. These data were taken at -120C. Radiation damage of F chips the need for the CT correction before spectral analysis

1 Statistics Aneta Siemiginowska a chapter for X-ray Astronomy Handbook October 2008

1 Statistics Aneta Siemiginowska a chapter for X-ray Astronomy Handbook October 2008 1 Statistics Aneta Siemiginowska a chapter for X-ray Astronomy Handbook October 2008 1.1 Introduction Why do we need statistic? Wall and Jenkins (2003) give a good description of the scientific analysis

More information

Spectral Analysis of High Resolution X-ray Binary Data

Spectral Analysis of High Resolution X-ray Binary Data Spectral Analysis of High Resolution X-ray Binary Data Michael Nowak, mnowak@space.mit.edu X-ray Astronomy School; Aug. 1-5, 2011 Introduction This exercise takes a look at X-ray binary observations using

More information

Andy Ptak NASA/GSFC Active Galactic Nuclei

Andy Ptak NASA/GSFC Active Galactic Nuclei Andy Ptak NASA/GSFC andrew.f.ptak@nasa.gov Active Galactic Nuclei Outline AGN X-ray spectral features overview CCD resolution spectral fitting Examples Active Galactic Nuclei Powered by accretion onto

More information

NuSTAR observation of the Arches cluster: X-ray spectrum extraction from a 2D image

NuSTAR observation of the Arches cluster: X-ray spectrum extraction from a 2D image NuSTAR observation of the Arches cluster: X-ray spectrum extraction from a 2D image Space Research Institute (IKI), Moscow, Russia E-mail: krivonos@iki.rssi.ru The NuSTAR mission performed a long (200

More information

Analysis of Off-Nuclear X-Ray Sources in Galaxy NGC Sarah M. Harrison

Analysis of Off-Nuclear X-Ray Sources in Galaxy NGC Sarah M. Harrison SLAC-TN-6-19 August 26 Analysis of Off-Nuclear X-Ray Sources in Galaxy NGC 4945 Sarah M. Harrison Office of Science, Science Undergraduate Laboratory Internship (SULI) Massachusetts Institute of Technology

More information

XSPEC. Keith Arnaud. Center for Research Excellence in Space Science and Technology

XSPEC. Keith Arnaud. Center for Research Excellence in Space Science and Technology XSPEC Keith Arnaud Center for Research Excellence in Space Science and Technology University of Maryland College Park and NASA s Goddard Space Flight Center 9,525 refereed papers containing the word xspec

More information

RELATIVISTIC SPECTROSCOPY OF BLACK HOLES

RELATIVISTIC SPECTROSCOPY OF BLACK HOLES RELATIVISTIC SPECTROSCOPY OF BLACK HOLES Michael Parker ESAC science seminar 24/5/18 BLACK HOLES 101 For an object to just escape a massive body, it needs the sum: Kinetic energy + gravitational binding

More information

X-ray data analysis. Andrea Marinucci. Università degli Studi Roma Tre

X-ray data analysis. Andrea Marinucci. Università degli Studi Roma Tre X-ray data analysis Andrea Marinucci Università degli Studi Roma Tre marinucci@fis.uniroma3.it Goal of these lectures X-ray data analysis why? what? how? How? This is intentended to be a general reference

More information

Monte Carlo Simulator to Study High Mass X-ray Binary System

Monte Carlo Simulator to Study High Mass X-ray Binary System SLAC-PUB-11350 Monte Carlo Simulator to Study High Mass X-ray Binary System S. Watanabe, F. Nagase, T. Takahashi ISAS/JAXA, Sagamihara, Kanagawa 229-8510, Japan M. Sako, S.M. Kahn KIPAC/Stanford, Stanford,

More information

JEREMY DRAKE, PETE RATZLAFF, VINAY KASHYAP AND THE MC CALIBRATION UNCERTAINTIES TEAM MONTE CARLO CONSTRAINTS ON INSTRUMENT CALIBRATION

JEREMY DRAKE, PETE RATZLAFF, VINAY KASHYAP AND THE MC CALIBRATION UNCERTAINTIES TEAM MONTE CARLO CONSTRAINTS ON INSTRUMENT CALIBRATION JEREMY DRAKE, PETE RATZLAFF, VINAY KASHYAP AND THE MC CALIBRATION UNCERTAINTIES TEAM MONTE CARLO CONSTRAINTS ON INSTRUMENT CALIBRATION 11th IACHEC Meeting, Pune March 2016 JEREMY DRAKE, PETE RATZLAFF,

More information

3.1.1 Lightcurve, colour-colour and hardness intensity diagram

3.1.1 Lightcurve, colour-colour and hardness intensity diagram Chapter 3 X ray data analysis methods 3.1 Data Analysis Procedure The analysis and reduction procedure of astronomical data can be broadly classified into two categories - (1) count rate variations as

More information

MONTE CARLO METHODS FOR TREATING AND UNDERSTANDING HIGHLY-CORRELATED INSTRUMENT CALIBRATION UNCERTAINTIES

MONTE CARLO METHODS FOR TREATING AND UNDERSTANDING HIGHLY-CORRELATED INSTRUMENT CALIBRATION UNCERTAINTIES JEREMY DRAKE, PETE RATZLAFF, VINAY KASHYAP AND THE MC CALIBRATION UNCERTAINTIES TEAM MONTE CARLO METHODS FOR TREATING AND UNDERSTANDING HIGHLY-CORRELATED INSTRUMENT CALIBRATION UNCERTAINTIES Harvard Astrostatistics

More information

X-ray Observations of Nearby Galaxies

X-ray Observations of Nearby Galaxies X-ray Observations of Nearby Galaxies Wolfgang Pietsch Max Planck Institut für extraterrestrische Physik XMM-Newton EPIC view of the starburst galaxy NGC 253 Wolfgang Pietsch Max Planck Institut für extraterrestrische

More information

X-ray Radiation, Absorption, and Scattering

X-ray Radiation, Absorption, and Scattering X-ray Radiation, Absorption, and Scattering What we can learn from data depend on our understanding of various X-ray emission, scattering, and absorption processes. We will discuss some basic processes:

More information

Radiation Detection and Measurement

Radiation Detection and Measurement Radiation Detection and Measurement June 2008 Tom Lewellen Tkldog@u.washington.edu Types of radiation relevant to Nuclear Medicine Particle Symbol Mass (MeV/c 2 ) Charge Electron e-,! - 0.511-1 Positron

More information

Broadband X-ray emission from radio-quiet Active Galactic Nuclei

Broadband X-ray emission from radio-quiet Active Galactic Nuclei 29 th ASI Meeting ASI Conference Series, 2011, Vol. 3, pp 19 23 Edited by Pushpa Khare & C. H. Ishwara-Chandra Broadband X-ray emission from radio-quiet Active Galactic Nuclei G. C. Dewangan Inter-University

More information

X-ray Photoionized Models of Emission Spectra Tim Kallman NASA/GSFC. Why Models Emission vs. absorpeon An example ImplicaEons for lab astro

X-ray Photoionized Models of Emission Spectra Tim Kallman NASA/GSFC. Why Models Emission vs. absorpeon An example ImplicaEons for lab astro X-ray Photoionized Models of Emission Spectra Tim Kallman NASA/GSFC Why Models Emission vs. absorpeon An example ImplicaEons for lab astro Many astrophysical objects exhibit line emission associated with

More information

High Resolution X-Ray Spectra A DIAGNOSTIC TOOL OF THE HOT UNIVERSE

High Resolution X-Ray Spectra A DIAGNOSTIC TOOL OF THE HOT UNIVERSE High Resolution X-Ray Spectra A DIAGNOSTIC TOOL OF THE HOT UNIVERSE Katarzyna Bensch Rosario González-Riestra María Santos-Lleó XMM-NEWTON / ESAC The RGS Instruments THE REFLECTION GRATING SPECTROMETER

More information

A Cluster of Galaxies, Abell 496

A Cluster of Galaxies, Abell 496 A Cluster of Galaxies, Abell 496 F. Nobels and W. Mulder MSc Astronomy & Kapteyn Astronomical Institute, University of Groningen (Dated: November 11, 2015) As part of the course High-Energy Astrophysics,

More information

Chapter 0 Introduction X-RAY BINARIES

Chapter 0 Introduction X-RAY BINARIES X-RAY BINARIES 1 Structure of this course 0. Introduction 1. Compact stars: formation and observational appearance. Mass transfer in binaries 3. Observational properties of XRBs 4. Formation and evolution

More information

Brandon C. Kelly (Harvard Smithsonian Center for Astrophysics)

Brandon C. Kelly (Harvard Smithsonian Center for Astrophysics) Brandon C. Kelly (Harvard Smithsonian Center for Astrophysics) Probability quantifies randomness and uncertainty How do I estimate the normalization and logarithmic slope of a X ray continuum, assuming

More information

Dan Evans (Harvard), Julia Lee (Harvard), Jane Turner (UMBC/GSFC), Kim Weaver (GSFC), Herman Marshall (MIT) NGC 2110 Spectroscopy

Dan Evans (Harvard), Julia Lee (Harvard), Jane Turner (UMBC/GSFC), Kim Weaver (GSFC), Herman Marshall (MIT) NGC 2110 Spectroscopy Dan Evans (Harvard), Julia Lee (Harvard), Jane Turner (UMBC/GSFC), Kim Weaver (GSFC), Herman Marshall (MIT) Fe Kα Lines and Reflection: AGN Geometry The Fe Kα line complex in general consists of a narrow

More information

Lecture 20 High-Energy Astronomy. HEA intro X-ray astrophysics a very brief run through. Swift & GRBs 6.4 kev Fe line and the Kerr metric

Lecture 20 High-Energy Astronomy. HEA intro X-ray astrophysics a very brief run through. Swift & GRBs 6.4 kev Fe line and the Kerr metric Lecture 20 High-Energy Astronomy HEA intro X-ray astrophysics a very brief run through. Swift & GRBs 6.4 kev Fe line and the Kerr metric Tut 5 remarks Generally much better. However: Beam area. T inst

More information

Detection of X-Rays. Solid state detectors Proportional counters Microcalorimeters Detector characteristics

Detection of X-Rays. Solid state detectors Proportional counters Microcalorimeters Detector characteristics Detection of X-Rays Solid state detectors Proportional counters Microcalorimeters Detector characteristics Solid State X-ray Detectors X-ray interacts in material to produce photoelectrons which are collected

More information

X-ray Radiation, Absorption, and Scattering

X-ray Radiation, Absorption, and Scattering X-ray Radiation, Absorption, and Scattering What we can learn from data depend on our understanding of various X-ray emission, scattering, and absorption processes. We will discuss some basic processes:

More information

Observational appearance of strong gravity in X-rays

Observational appearance of strong gravity in X-rays Observational appearance of strong gravity in X-rays Compact object Companion star (5000 K) Schwarzschild radius Schwarzschild metric R s 2GM Schwarzschild radius (= 2.95 km for 1М ) c 2 Accretion Black

More information

THE OBSERVATION AND ANALYSIS OF STELLAR PHOTOSPHERES

THE OBSERVATION AND ANALYSIS OF STELLAR PHOTOSPHERES THE OBSERVATION AND ANALYSIS OF STELLAR PHOTOSPHERES DAVID F. GRAY University of Western Ontario, London, Ontario, Canada CAMBRIDGE UNIVERSITY PRESS Contents Preface to the first edition Preface to the

More information

X-ray spectroscopy of nearby galactic nuclear regions

X-ray spectroscopy of nearby galactic nuclear regions X-ray spectroscopy of nearby galactic nuclear regions 1. How intermittent are AGNs? 2. What about the silent majority of the SMBHs? 3. How do SMBHs interplay with their environments? Q. Daniel Wang University

More information

Chandra Analysis of a Possible Cooling Core Galaxy Cluster at z = 1.03

Chandra Analysis of a Possible Cooling Core Galaxy Cluster at z = 1.03 Chandra Analysis of a Possible Cooling Core Galaxy Cluster at z = 1.03 Kyle Dolan Institute for Astronomy, University of Hawaii at Manoa (Dated: 12 August 2005) We present an analysis of Chandra observations

More information

Line Broadening. φ(ν) = Γ/4π 2 (ν ν 0 ) 2 + (Γ/4π) 2, (3) where now Γ = γ +2ν col includes contributions from both natural broadening and collisions.

Line Broadening. φ(ν) = Γ/4π 2 (ν ν 0 ) 2 + (Γ/4π) 2, (3) where now Γ = γ +2ν col includes contributions from both natural broadening and collisions. Line Broadening Spectral lines are not arbitrarily sharp. There are a variety of mechanisms that give them finite width, and some of those mechanisms contain significant information. We ll consider a few

More information

AGN Physics of the Ionized Gas Physical conditions in the NLR Physical conditions in the BLR LINERs Emission-Line Diagnostics High-Energy Effects

AGN Physics of the Ionized Gas Physical conditions in the NLR Physical conditions in the BLR LINERs Emission-Line Diagnostics High-Energy Effects AGN Physics of the Ionized Gas Physical conditions in the NLR Physical conditions in the BLR LINERs Emission-Line Diagnostics High-Energy Effects 1 Evidence for Photoionization - continuum and Hβ luminosity

More information

Problem Solving. radians. 180 radians Stars & Elementary Astrophysics: Introduction Press F1 for Help 41. f s. picture. equation.

Problem Solving. radians. 180 radians Stars & Elementary Astrophysics: Introduction Press F1 for Help 41. f s. picture. equation. Problem Solving picture θ f = 10 m s =1 cm equation rearrange numbers with units θ factors to change units s θ = = f sinθ fθ = s / cm 10 m f 1 m 100 cm check dimensions 1 3 π 180 radians = 10 60 arcmin

More information

Example: model a star using a two layer model: Radiation starts from the inner layer as blackbody radiation at temperature T in. T out.

Example: model a star using a two layer model: Radiation starts from the inner layer as blackbody radiation at temperature T in. T out. Next, consider an optically thick source: Already shown that in the interior, radiation will be described by the Planck function. Radiation escaping from the source will be modified because the temperature

More information

Characteristic temperatures

Characteristic temperatures Characteristic temperatures Effective temperature Most sources are only roughly blackbodies (if that). So we integrate the flux over frequency and define: F = I cosθ dω d = σ T e 4 i.e. a source of effective

More information

Soft X-ray Emission Lines in Active Galactic Nuclei. Mat Page

Soft X-ray Emission Lines in Active Galactic Nuclei. Mat Page Soft X-ray Emission Lines in Active Galactic Nuclei Mat Page MSSL-UCL Observations of soft X-ray line emission allow us to investigate highly ionized plasmas in galaxies and AGN. I ll start with the most

More information

Chemical Engineering 412

Chemical Engineering 412 Chemical Engineering 412 Introductory Nuclear Engineering Lecture 26 Radiation Detection & Measurement II Spiritual Thought 2 I would not hold the position in the Church I hold today had I not followed

More information

Citation for published version (APA): Wang, Y. (2018). Disc reflection in low-mass X-ray binaries. [Groningen]: Rijksuniversiteit Groningen.

Citation for published version (APA): Wang, Y. (2018). Disc reflection in low-mass X-ray binaries. [Groningen]: Rijksuniversiteit Groningen. University of Groningen Disc reflection in low-mass X-ray binaries Wang, Yanan IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from it. Please check

More information

NuSTAR spectral analysis of the two bright Seyfert 1 galaxies: MCG and NGC Alessia Tortosa Università Roma Tre

NuSTAR spectral analysis of the two bright Seyfert 1 galaxies: MCG and NGC Alessia Tortosa Università Roma Tre NuSTAR spectral analysis of the two bright Seyfert 1 galaxies: MCG 8-11-11 and NGC 6814 Alessia Tortosa Università Roma Tre From the Dolomites to the event horizon: sledging down the black hole potential

More information

Absorption models in SPEX. Katrien C. Steenbrugge St John s College, University of Oxford

Absorption models in SPEX. Katrien C. Steenbrugge St John s College, University of Oxford Absorption models in SPEX Katrien C. Steenbrugge St John s College, University of Oxford Overview Introduction Absm model Hot model Slab model Xabs model Warm model Introduction Collisionally ionized absorption

More information

Ultra Luminous X-ray sources ~one of the most curious objects in the universe~

Ultra Luminous X-ray sources ~one of the most curious objects in the universe~ Ultra Luminous X-ray sources ~one of the most curious objects in the universe~ Shogo B. Kobayashi the University of Tokyo ULX workshop@isas 1 The discovery of the enigmatic sources pfabbiano & Trincheri

More information

Sample Spectroscopy System Hardware

Sample Spectroscopy System Hardware Semiconductor Detectors vs. Scintillator+PMT Detectors Semiconductors are emerging technology - Scint.PMT systems relatively unchanged in 50 years. NaI(Tl) excellent for single-photon, new scintillation

More information

Accretion Disks. 1. Accretion Efficiency. 2. Eddington Luminosity. 3. Bondi-Hoyle Accretion. 4. Temperature profile and spectrum of accretion disk

Accretion Disks. 1. Accretion Efficiency. 2. Eddington Luminosity. 3. Bondi-Hoyle Accretion. 4. Temperature profile and spectrum of accretion disk Accretion Disks Accretion Disks 1. Accretion Efficiency 2. Eddington Luminosity 3. Bondi-Hoyle Accretion 4. Temperature profile and spectrum of accretion disk 5. Spectra of AGN 5.1 Continuum 5.2 Line Emission

More information

Where are oxygen synthesized in stars?

Where are oxygen synthesized in stars? The oxygen abundance from X-rays : methods and prospects K. Matsushita Where are oxygen synthesized in stars? Hot intracluster medium (ICM) Warm-hot intergalactic medium? Hot interstellar medium in early-type

More information

The X-Ray Universe. Potsdam University. Dr. Lidia Oskinova Wintersemester 2008/09

The X-Ray Universe. Potsdam University. Dr. Lidia Oskinova Wintersemester 2008/09 The X-Ray Universe The X-Ray Universe Potsdam University Dr. Lidia Oskinova Wintersemester 2008/09 lida@astro.physik.uni-potsdam.de astro.physik.uni-potsdam.de/~lida/x-ray.html Chandra X-ray Observatory

More information

Modeling Compton Thick Obscuration in Accreting Supermassive Black Hole Systems

Modeling Compton Thick Obscuration in Accreting Supermassive Black Hole Systems Modeling Compton Thick Obscuration in Accreting Supermassive Black Hole Systems Kendrah Murphy MIT Kavli Institute for Astrophysics & Space Research Tahir Yaqoob Johns Hopkins University/NASA GSFC Overview

More information

Detecting high energy photons. Interactions of photons with matter Properties of detectors (with examples)

Detecting high energy photons. Interactions of photons with matter Properties of detectors (with examples) Detecting high energy photons Interactions of photons with matter Properties of detectors (with examples) Interactions of high energy photons with matter Cross section/attenution length/optical depth Photoelectric

More information

X- ray surface brightness fluctuations and turbulence in galaxy clusters. Jeremy Sanders. Andy Fabian. Sanders & Fabian 2011, MNRAS, submitted

X- ray surface brightness fluctuations and turbulence in galaxy clusters. Jeremy Sanders. Andy Fabian. Sanders & Fabian 2011, MNRAS, submitted X- ray surface brightness fluctuations and turbulence in galaxy clusters Jeremy Sanders Andy Fabian Sanders & Fabian 2011, MNRAS, submitted Simulations predict that in galaxy clusters turbulent energy

More information

Foundations of Astrophysics

Foundations of Astrophysics Foundations of Astrophysics Barbara Ryden The Ohio State University Bradley M. Peterson The Ohio State University Preface xi 1 Early Astronomy 1 1.1 The Celestial Sphere 1 1.2 Coordinate Systems on a Sphere

More information

High-resolution spectroscopy of the low-mass X-ray binary EXO

High-resolution spectroscopy of the low-mass X-ray binary EXO A&A 365, L277 L281 (2001) DOI: 10.1051/0004-6361:20000053 c ESO 2001 Astronomy & Astrophysics High-resolution spectroscopy of the low-mass X-ray binary EXO 0748 67 J. Cottam 1,S.M.Kahn 1,A.C.Brinkman 2,J.W.denHerder

More information

Ay Fall 2004 Lecture 6 (given by Tony Travouillon)

Ay Fall 2004 Lecture 6 (given by Tony Travouillon) Ay 122 - Fall 2004 Lecture 6 (given by Tony Travouillon) Stellar atmospheres, classification of stellar spectra (Many slides c/o Phil Armitage) Formation of spectral lines: 1.excitation Two key questions:

More information

Astrophysical techniques. Lecture 3: X-ray astronomy

Astrophysical techniques. Lecture 3: X-ray astronomy Astrophysical techniques Lecture 3: X-ray astronomy Peter Wheatley P.J.Wheatley@warwick.ac.uk Outline Introduction X-ray emission processes Instrumentation Practicalities Multi-wavelength astronomy UV

More information

Ultra-fast disk wind from a high accretion rate black hole 1H

Ultra-fast disk wind from a high accretion rate black hole 1H Ultra-fast disk wind from a high accretion rate black hole 1H 0707-495 Kouichi Hagino (ISAS/JAXA) H. Odaka, C. Done, R. Tomaru, S. Watanabe, T. Takahashi K. Hagino et al. 2016, MNRAS, 461, 3954 BREAKING

More information

A relativistically smeared line profile in the spectrum of the bright Z-source GX 340+0

A relativistically smeared line profile in the spectrum of the bright Z-source GX 340+0 A relativistically smeared line profile in the spectrum of the bright Z-source GX 340+0 Dipartimento di Scienze Fisiche ed Astronomiche. Università degli Studi di Palermo, Italy E-mail: dai@fisica.unipa.it

More information

Measurement Error and Linear Regression of Astronomical Data. Brandon Kelly Penn State Summer School in Astrostatistics, June 2007

Measurement Error and Linear Regression of Astronomical Data. Brandon Kelly Penn State Summer School in Astrostatistics, June 2007 Measurement Error and Linear Regression of Astronomical Data Brandon Kelly Penn State Summer School in Astrostatistics, June 2007 Classical Regression Model Collect n data points, denote i th pair as (η

More information

XMM-Newton - Chandra Synergy: Maximizing their Future Science Impact over the Next Decade

XMM-Newton - Chandra Synergy: Maximizing their Future Science Impact over the Next Decade XMM-Newton - Chandra Synergy: Maximizing their Future Science Impact over the Next Decade Norbert Schartel XMM-Newton Project Scientist, ESA, ESAC, Spain 17.8.2016 Instruments Draft N. Schartel XMM-NEWTON

More information

Chandra was launched aboard Space Shuttle Columbia on July 23, 1999!!!

Chandra was launched aboard Space Shuttle Columbia on July 23, 1999!!! Chandra was launched aboard Space Shuttle Columbia on July 23, 1999!!! Crew Lost During Re-Entry Modern X-ray Telescopes and Detectors X-ray Telescopes X-ray Instruments Some early highlights Observations

More information

Hot-subdwarf stars: a new class of X-ray sources

Hot-subdwarf stars: a new class of X-ray sources Hot-subdwarf stars: a new class of X-ray sources Nicola La Palombara 1 S. Mereghetti 1, A. Tiengo 1,2,3, P. Esposito 1 1 INAF - IASF Milano, 2 IUSS Pavia, 3 INFN Pavia Dublin June 16-19, 2014 Hot subdwarf

More information

Chandra Observation of Point Sources in the X-Ray Elliptical Galaxy NGC 1407

Chandra Observation of Point Sources in the X-Ray Elliptical Galaxy NGC 1407 Chin. J. Astron. Astrophys. Vol. 4 (2004), No. 3, 221 230 ( http: /www.chjaa.org or http: /chjaa.bao.ac.cn ) Chinese Journal of Astronomy and Astrophysics Chandra Observation of Point Sources in the X-Ray

More information

Probing Neutron Star Physics using Thermonuclear X-ray Bursts

Probing Neutron Star Physics using Thermonuclear X-ray Bursts Probing Neutron Star Physics using Thermonuclear X-ray Bursts Sudip Bhattacharyya University of Maryland (CRESST) NASA s Goddard Space Flight Center Outline Neutron Stars: why do we care? Thermonuclear

More information

Simultaneous XMM-Newton Radio Observations of the Mode-switching Radio Pulsar PSR B Wim Hermsen 1,2

Simultaneous XMM-Newton Radio Observations of the Mode-switching Radio Pulsar PSR B Wim Hermsen 1,2 Simultaneous XMM-Newton Radio Observations of the Mode-switching Radio Pulsar PSR B1822-09 Wim Hermsen 1,2 Collaborators: J.W.T. Hessels 3,2, L. Kuiper 1, J. van Leeuwen 3,2, D. Mitra 4, J.M. Rankin 2,5,

More information

Randall Smith (CfA) Collaborators: N. Brickhouse, A. Foster, H. Yamaguchi (CfA), J. Wilms (Erlangen), Li Ji (PMO)

Randall Smith (CfA) Collaborators: N. Brickhouse, A. Foster, H. Yamaguchi (CfA), J. Wilms (Erlangen), Li Ji (PMO) Randall Smith (CfA) Collaborators: N. Brickhouse, A. Foster, H. Yamaguchi (CfA), J. Wilms (Erlangen), Li Ji (PMO) X-ray Line Emission from Capella Mg, Si, S, and Fe are unambiguously detected Inconsistent

More information

Principal Component Analysis

Principal Component Analysis Principal Component Analysis Revealing AGN Spectral Variability Michael Parker with Andy Fabian, Giorgio Matt, Erin Kara, Andrea Marinucci, Dom Walton, Karri Koljonen, Will Alston, Laura Brenneman and

More information

Physics and Chemistry of the Interstellar Medium

Physics and Chemistry of the Interstellar Medium Physics and Chemistry of the Interstellar Medium Sun Kwok The University of Hong Kong UNIVERSITY SCIENCE BOOKS Sausalito, California * Preface xi The Interstellar Medium.1.1 States of Matter in the ISM

More information

PoS(INTEGRAL 2012)049

PoS(INTEGRAL 2012)049 Testing reflection features in 4U 175 44 with XMM-Newton, BeppoSAX, and RXTE in the hard and soft states Università di Cagliari - Italy E-mail: emailelise.egron@dsf.unica.it Tiziana Di Salvo Università

More information

How does your eye form an Refraction

How does your eye form an Refraction Astronomical Instruments Eyes and Cameras: Everyday Light Sensors How does your eye form an image? How do we record images? How does your eye form an image? Refraction Refraction is the bending of light

More information

Future X-rayX Spectroscopy Missions. Jan-Willem den Herder

Future X-rayX Spectroscopy Missions. Jan-Willem den Herder Future X-rayX Spectroscopy Missions Jan-Willem den Herder contents Plasma diagnostics in the 0.1 to 10 kev band with resolution > 100 X-ray spectrometers: instrumental promises Future missions (a dream)

More information

Topics Covered in Chapter. Light and Other Electromagnetic Radiation. A Subatomic Interlude II. A Subatomic Interlude. A Subatomic Interlude III

Topics Covered in Chapter. Light and Other Electromagnetic Radiation. A Subatomic Interlude II. A Subatomic Interlude. A Subatomic Interlude III Light and Other Electromagnetic Radiation Topics Covered in Chapter 1.Structure of Atoms 2.Origins of Electromagnetic Radiation 3.Objects with Different Temperature and their Electromagnetic Radiation

More information

Light and Other Electromagnetic Radiation

Light and Other Electromagnetic Radiation Light and Other Electromagnetic Radiation 1 Topics Covered in Chapter 1.Structure of Atoms 2.Origins of Electromagnetic Radiation 3.Objects with Different Temperature and their Electromagnetic Radiation

More information

Gamma-ray Astrophysics with VERITAS: Exploring the violent Universe

Gamma-ray Astrophysics with VERITAS: Exploring the violent Universe Gamma-ray Astrophysics with VERITAS: Exploring the violent Universe K. Ragan McGill University Soup & Science 11-Jan-2008 Soup & Science Jan. '08 1 How do we know about the Universe? Historically, all

More information

Accounting for Calibration Uncertainty in Spectral Analysis. David A. van Dyk

Accounting for Calibration Uncertainty in Spectral Analysis. David A. van Dyk Bayesian Analysis of Accounting for in Spectral Analysis David A van Dyk Statistics Section, Imperial College London Joint work with Vinay Kashyap, Jin Xu, Alanna Connors, and Aneta Siegminowska ISIS May

More information

arxiv:astro-ph/ v1 6 May 2004

arxiv:astro-ph/ v1 6 May 2004 XMM-NEWTON OBSERVATIONS OF THREE HIGH REDSHIFT RADIO GALAXIES arxiv:astro-ph/0405116 v1 6 May 2004 Abstract E. Belsole, D.M. Worrall, M. J. Hardcastle Department of Physics - University of Bristol Tyndall

More information

Tidal disruption events from the first XMM-Newton Slew Survey

Tidal disruption events from the first XMM-Newton Slew Survey Tidal disruption events from the first XMM-Newton Slew Survey Pilar Esquej Centro de Astrobiologia (INTA-CSIC) Richard Saxton, Stefanie Komossa, Andy Read, M. Sanchez-Portal et al. Tidal disruption events

More information

Atomic Structure & Radiative Transitions

Atomic Structure & Radiative Transitions Atomic Structure & Radiative Transitions electron kinetic energy nucleus-electron interaction electron-electron interaction Remember the meaning of spherical harmonics Y l, m (θ, ϕ) n specifies the

More information

A8824: Statistics Notes David Weinberg, Astronomy 8824: Statistics Notes 6 Estimating Errors From Data

A8824: Statistics Notes David Weinberg, Astronomy 8824: Statistics Notes 6 Estimating Errors From Data Where does the error bar go? Astronomy 8824: Statistics otes 6 Estimating Errors From Data Suppose you measure the average depression of flux in a quasar caused by absorption from the Lyman-alpha forest.

More information

Response Function of the BGO and NaI(Tl) Detectors Using Monte Carlo Simulations

Response Function of the BGO and NaI(Tl) Detectors Using Monte Carlo Simulations Response Function of the BGO and NaI(Tl) Detectors Using Monte Carlo Simulations 271 I. ORION 1,2 1 AND L. WIELOPOLSKI 2, 1 St. Luke=s/Roosevelt Hospital, Columbia University, New York, NY 10025,USA 2

More information

arxiv:astro-ph/ v1 17 Dec 2001

arxiv:astro-ph/ v1 17 Dec 2001 Spectral properties of the Narrow-Line Seyfert 1 Quasar PG1211+143 arxiv:astro-ph/0112387v1 17 Dec 2001 A. Janiuk 1, B. Czerny 1, G.M. Madejski 2 1) N. Copernicus Astronomical Centre, Bartycka 18, 00-716,

More information

New Suzaku Results on Active Galaxies. or Reflections on AGN. James Reeves (Keele) and Suzaku team (given by Lance Miller)

New Suzaku Results on Active Galaxies. or Reflections on AGN. James Reeves (Keele) and Suzaku team (given by Lance Miller) New Suzaku Results on Active Galaxies or Reflections on AGN James Reeves (Keele) and Suzaku team (given by Lance Miller) Overview Overview of Suzaku and instrument performance. Why Suzaku is important

More information

Practical Statistics

Practical Statistics Practical Statistics Lecture 1 (Nov. 9): - Correlation - Hypothesis Testing Lecture 2 (Nov. 16): - Error Estimation - Bayesian Analysis - Rejecting Outliers Lecture 3 (Nov. 18) - Monte Carlo Modeling -

More information

Some HI is in reasonably well defined clouds. Motions inside the cloud, and motion of the cloud will broaden and shift the observed lines!

Some HI is in reasonably well defined clouds. Motions inside the cloud, and motion of the cloud will broaden and shift the observed lines! Some HI is in reasonably well defined clouds. Motions inside the cloud, and motion of the cloud will broaden and shift the observed lines Idealized 21cm spectra Example observed 21cm spectra HI densities

More information

v Characteristics v Possible Interpretations L X = erg s -1

v Characteristics v Possible Interpretations L X = erg s -1 Ultra Luminous X-ray sources v Characteristics L X = 39-41 erg s -1 >> M black hole(bh) L Edd = 39 erg s -1. Usually seen in star forming regions. On arms of spiral galaxies (right). In starburst or irregular

More information

Synchrotron Radiation II

Synchrotron Radiation II Synchrotron Radiation II Cyclotron v

More information

arxiv:astro-ph/ v2 12 Jan 2004

arxiv:astro-ph/ v2 12 Jan 2004 Astronomy & Astrophysics manuscript no. rxj0720 December 31, 2013 (DOI: will be inserted by hand later) Long-term variability in the X-ray emission of RX J0720.4-3125 Cor P. de Vries 1, Jacco Vink 1, Mariano

More information

Problem Set 2 Solutions

Problem Set 2 Solutions Problem Set 2 Solutions Problem 1: A A hot blackbody will emit more photons per unit time per unit surface area than a cold blackbody. It does not, however, necessarily need to have a higher luminosity,

More information

Comets observed with XMM-Newton

Comets observed with XMM-Newton A major step in cometary X-ray research Konrad Dennerl Max-Planck-Institut für extraterrestrische Physik Comets the most surprising X-ray objects in space 1996: Comet Hyakutake (C/1996 B2) EPIC Consortium

More information

Mapping the non thermal emission in Coma cluster of galaxies using the FeXXV/FeXXVI line ratio

Mapping the non thermal emission in Coma cluster of galaxies using the FeXXV/FeXXVI line ratio Mapping the non thermal emission in Coma cluster of galaxies using the FeXXV/FeXXVI line ratio J. Nevalainen 1 & D. Eckert 2 1 Helsinki University Observatory, Finland 2 ISDC, Geneva, Switzerland Evidence

More information

Cooling Limits for the

Cooling Limits for the Cooling Limits for the Page et al. 2004 Youngest Neutron Stars Cooling from the Youngest NSs SNR Zone NSs younger than ~50 kyr offer strong constraints on rapid cooling - the associated physical processes

More information

Light and Atoms. ASTR 1120 General Astronomy: Stars & Galaxies. ASTR 1120 General Astronomy: Stars & Galaxies !ATH REVIEW: #AST CLASS: "OMEWORK #1

Light and Atoms. ASTR 1120 General Astronomy: Stars & Galaxies. ASTR 1120 General Astronomy: Stars & Galaxies !ATH REVIEW: #AST CLASS: OMEWORK #1 ASTR 1120 General Astronomy: Stars & Galaxies!ATH REVIEW: Tonight, 5-6pm, in RAMY N1B23 "OMEWORK #1 -Due THU, Sept. 10, by 5pm, on Mastering Astronomy CLASS RECORDED STARTED - INFO WILL BE POSTED on CULEARN

More information

The Magnificent Seven Similarities and Differences

The Magnificent Seven Similarities and Differences The Magnificent Seven Similarities and Differences Frank Haberl Max-Planck-Institut für extraterrestrische Physik (MPE), Garching The discovery of thermal, radio quiet isolated neutron stars New XMM-Newton

More information

arxiv: v1 [astro-ph] 30 Jan 2008

arxiv: v1 [astro-ph] 30 Jan 2008 Astronomy & Astrophysics manuscript no. aa xmm c ESO 2018 May 26, 2018 XMM-Newton observations of the Small Magellanic Cloud: X-ray outburst of the 6.85 s pulsar XTE J0103-728 F. Haberl and W. Pietsch

More information

Emission from Isolated Neutron Stars. Observations and model applications

Emission from Isolated Neutron Stars. Observations and model applications Thermal Emission from Isolated Neutron Stars Theoretical aspects Observations and model applications Slava Zavlin (MPE, Garching) A Short History Chi & Salpeter (1964) and Tsuruta (1964): thermal radiation

More information

Notes on Photoionized Regions Wednesday, January 12, 2011

Notes on Photoionized Regions Wednesday, January 12, 2011 Notes on Photoionized Regions Wednesday, January 12, 2011 CONTENTS: 1. Introduction 2. Hydrogen Nebulae A. Ionization equations B. Recombination coefficients and cross sections C. Structure of the hydrogen

More information

Laura Barragán. Universidad Complutense de Madrid

Laura Barragán. Universidad Complutense de Madrid . Universidad Complutense de Madrid Tutors: Erik Kuulkers & Peter Kretschmar INTEGRAL: International Gamma-Ray Astrophysics Laboratory The payload The mission The project IBIS (Imager on-board the Integral

More information

X-ray spectroscopy of low-mass X-ray binaries

X-ray spectroscopy of low-mass X-ray binaries X-ray spectroscopy of low-mass X-ray binaries Laurence Boirin Observatoire astronomique de Strasbourg (Handout version) Artistic impression of a low-mass X-ray binary An X-ray binary spectrum (from the

More information

What can we learn from the AGN radio quiet continuum with SIMBOL-X?

What can we learn from the AGN radio quiet continuum with SIMBOL-X? What can we learn from the AGN radio quiet continuum with SIMBOL-X? Pierre-Olivier Petrucci LAOG, Grenoble, France Generalities about the RQ AGN X-ray emission - continuum - reflection component Expected

More information

Prentice Hall EARTH SCIENCE

Prentice Hall EARTH SCIENCE Prentice Hall EARTH SCIENCE Tarbuck Lutgens Chapter 24 Studying the Sun 24.1 The Study of Light Electromagnetic Radiation Electromagnetic radiation includes gamma rays, X-rays, ultraviolet light, visible

More information

Applied Nuclear Physics (Fall 2006) Lecture 21 (11/29/06) Detection of Nuclear Radiation: Pulse Height Spectra

Applied Nuclear Physics (Fall 2006) Lecture 21 (11/29/06) Detection of Nuclear Radiation: Pulse Height Spectra 22.101 Applied Nuclear Physics (Fall 2006) Lecture 21 (11/29/06) Detection of Nuclear Radiation: Pulse Height Spectra References: W. E. Meyerhof, Elements of Nuclear Physics (McGraw-Hill, New York, 1967),

More information

X-raying galactic feedback in nearby disk galaxies. Q. Daniel Wang University of Massachusetts

X-raying galactic feedback in nearby disk galaxies. Q. Daniel Wang University of Massachusetts X-raying galactic feedback in nearby disk galaxies Q. Daniel Wang University of Massachusetts Chandra survey of diffuse X-ray emission from 53 edge-on galaxies i > 60 o, D < 30 Mpc (Li, J.-T. & Wang, Q.D.

More information

X-ray variability of AGN

X-ray variability of AGN X-ray variability of AGN Magnus Axelsson October 20, 2006 Abstract X-ray variability has proven to be an effective diagnostic both for Galactic black-hole binaries and active galactic nuclei (AGN). This

More information

IMPLEMENTATION OF THE MONTE CARLO-LIBRARY LEAST- SQUARES APPROACH TO ENERGY DISPERSIVE X-RAY FLUORESCENCE ANALYSIS

IMPLEMENTATION OF THE MONTE CARLO-LIBRARY LEAST- SQUARES APPROACH TO ENERGY DISPERSIVE X-RAY FLUORESCENCE ANALYSIS 227 IMPLEMENTATION OF THE MONTE CARLO-LIBRARY LEAST- SQUARES APPROACH TO ENERGY DISPERSIVE X-RAY FLUORESCENCE ANALYSIS Fusheng Li, Weijun Guo, and Robin P. Gardner Center for Engineering Applications of

More information

NuStar has similar area kev, Astro-H better at E> 80 Kev

NuStar has similar area kev, Astro-H better at E> 80 Kev NuStar has similar area 20-80 kev, Astro-H better at E> 80 Kev 10 3 1 HXI Effective Area [cm 2 ] 10 2 10 1 5 10 30 80 Energy [kev] NuStar and Astro-H HXI should have similar FOV HXI becomes background

More information