2 Properties of light, magnitudes, errors and error analysis
|
|
- Job Park
- 6 years ago
- Views:
Transcription
1 Astr 535 Class Notes Fall Properties of light, magnitudes, errors and error analysis 2.1 Light (References: Birney chapter 5, Sutton , Rieke 1.1) astronomy - learn about objects outside of our atmosphere. need information from them to do so: light! light, by quantum mechanics, is photons, has characteristics of both waves and particles. Wavelength/frequency corresponds to energy: E = hν = hc λ electromagnetic spectrum: gamma rays - X rays - UV - optical - IR - mm - radio. Different units often used for wavelength in different parts of spectrum: 1Å = m (used in UV, optical), 1nm = m (used in UV, optical), 1µ = m (used in IR), 1mm = m Numerical wavelengths of different parts of spectrum (roughly, there is no established string vocabulary!): far-uv ( µ, Å), near-uv ( µ, Å), optical (0.35 1µ, Å), near-ir (1 10µ), mid-ir (10 100µ), far-ir ( µ). Of course, some people used frequency instead of wavelength! And others, especially for high energy radiation, use energy! We can describe the amount of light an object emits or that we receive by three fundamental quantities: intensity or surface brightness, I, flux F, or luminosity, L. The surface brightness is defined as the amount of energy received in a unit surface element per unit time per unit frequency (or wavelength) from a unit solid angle in the direction (θ, φ), where θ is the angle away from the normal to the surface element, and φ the azimuthal angle. Note units of solid angle: sterradian, square degree, square arcsec. The flux is the amount of energy passing through a unit surface element in all directions, defined by F ν = I ν cos θdω
2 Astr 535 Class Notes Fall where dω is the solid angle element, and the integration is over the entire solid angle. Usually, the angle subtended by an object is very small, so the cos θ term is well approximated by unity. The luminosity is the intrinsic energy emitted by the source per second. For an isotropically emitting source, where d = distance to source. L = 4πd 2 F What do we measure for sources? For resolved sources we can directly measure their surface brightness (intensity) distribution on the sky, usually over some bandpass or wavelength interval; for unresolved sources, we measure the flux. We can only calculate the luminosity to a source if we know the distance. Questions What are the dimensions of the three quantities: luminosity, surface brightness (intensity), flux? How do the three quantities, L, I, and F, depend on distance to the source? (Note cosmological effects at large distances.) To what quantity is apparent magnitude of a star related? To what quantity the absolute magnitude? Note, there are variations in terminology, especially between disciplines (e.g. astronomy and engineering) and you should make sure you understand what is adopted by someone presenting data. While astronomers are reasonably consistent in the use of these terms, you may run into things like radiance( surface brightness), irradiance (flux), radiant flux (luminosity), spectral intensity, and others, which are all related to intensity or flux. Amount of light emitted is a function of wavelength, so we actually are often interested in estimates of the monochromatic flux/intensity/luminosity, sometimes known as flux/intensity/luminosity density = flux / unit wavelength (or unit frequency), also sometimes known as specific flux/intensity/luminosity. F ν : flux per unit frequency. F λ : flux per unit wavelength. nu2 nu1 lambda2 F ν dν = F λ dλ lambda1
3 Astr 535 Class Notes Fall F ν = F λ dλ/dν c = F λ ν = F λ 2 2 λ c Similarly, intensity and luminosity can be given per unit wavelength (or frequency). Note that a constant F λ implies a non-constant F ν and vice versa! Integral of flux/brightness over all wavelengths/frequencies gives the bolometric flux/brightness. Units: astronomers often (not always) work in CGS units, although, as discussed below, they most often work in a dimensionless unit... magnitudes. You may also run into the Jansky, a flux density unit corresponding to W/m 2 /Hz In practice, with most detectors, we measure photon flux (photon counting devices), rather than energy flux (bolometers). The photon flux is given by Terminology of measurements: photonflux = flux/energyperphoton = F λ λ hc dλ photometry (broad-band flux measurement), spectroscopy (relative measurement of fluxes at different wavelengths), spectrophotometry (absolute measurement of fluxes at different wavelengths). astrometry: concerned with positions of observed flux; morphology: intensity as a function of position; often, absolute measurements are unimportant Have a complete understanding of the difference between intensity, flux, and luminosity and their units. Recognize and understand that these can be specified per unit wavelength or per unit frequency and how to convert between the two.
4 Astr 535 Class Notes Fall Magnitudes and photometric systems In astronomy, however, magnitude units are often used instead of measuring the basic quantities in energy or photon flux. Magnitudes are a dimensionless quantities, and are related to flux (same holds for surface brightness or luminosity) by: or m = 2.5 log F F 0 m = 2.5 log F log F 0 where the coefficient of proportionality, F 0, depends on the definition of photometric system; the quantity 2.5 log F 0 may be referred to as the photometric system zeropoint. Inverting, one gets: F = F m Just as fluxes can be represented in magnitude units, flux densities can be specified by monochromatic magnitudes: F λ = F 0 (λ)10 0.4m(λ) although spectra are more often given in flux units than in magnitude units. Note that it is possible that F 0 is a function of wavelength! Note that since magnitudes are logarithmic, the difference between magnitudes corresponds to a ratio of fluxes; ratios of magnitudes are generally unphysical! If one is just doing relative measurements of brightness between objects, this can be done without knowledge of F 0 (or, equivalently, the system zeropoint); objects that differ in brightness by M mag have the same ratio of brightness ( M ) regardless of what photometric system they are in. The photometric system definitions and zeropoints are only needed when converting between calibrated magnitudes and fluxes. Note that this means that if one references the brightness of one object relative to that of another, a magnitude system can be set up relative the brightness of the reference source. However, the utility of a system when doing astrophysics generally requires an understanding of the actual fluxes. Know how magnitudes are defined, and that relative fluxes can be represented as magnitudes independent of the magnitude system. There are three main types of magnitude systems in use in astronomy. We start by describing the two simpler ones:the STMAG and the ABNU mag system. In these simple system, the reference flux is just a constant value in F λ or F ν. However, these
5 Astr 535 Class Notes Fall are not always the most widely used systems in astronomy, because no natural source exists with a flat spectrum. In the STMAG system, F 0,λ = 3.63E 9ergs/cm 2 /s/å, which is the flux of Vega at 5500Å; hence a star of Vega s brightness at 5500Å is defined to have m=0. Alternatively, we can write m ST MAG = 2.5 log F λ 21.1 (for F λ in cgs units). In the ABNU system, things are defined for F ν instead of F λ, and we have or F 0,ν = erg/cm 2 /s/hz10 0.4m ν m ABNU = 2.5 log F ν 48.6 (for F ν in cgs units). Again, the constant comes from the flux of Vega. Usually, when using magnitudes, people are talking about flux integrated over a spectral bandpass. In this case, F and F 0 refer to fluxes integrated over the bandpass. The STMAG and ABMAG integrated systems are defined relative to sources of constant F λ and F ν systems, respectively. m ST MAG = 2.5 log Fλ λdλ λdλ (the factor of λ comes in for photon counting detectors). m ABNU = 2.5 log (Fν /ν)dν ( /ν)dν (where the units are implicitly cgs with these numerical fluxes for Vega). Note that these systems differ by more than a constant, because one is defined by units of F λ and the other by F ν, so the difference between the systems is a function of wavelength. They are defined to be the same at 5500Å. (Question: what s the relation between m ST MAG and m ABNU?) Note also that, using magnitudes, the measured magnitude is nearly independent of bandpass width (a broader bandpass does not imply a brighter (smaller) magnitude), which is not the case for fluxes! The standard UBVRI broadband photometric system, as well as several other magnitude systems, however, are not defined for a constant F λ or F ν spectrum; rather, they are defined relative to the spectrum of an A0V star. Most systems are defined (or at least were originally) to have the magnitude of Vega be zero in all bandpasses
6 Astr 535 Class Notes Fall Figure 1: Differences between different magnitude systems (VEGAMAGS); if you ever get into this in detail, note that this is not exactly true for the UBVRI system. For the broadband UBVRI system, we have m UBV RI 2.5 log UBV RI F λ(object)λdλ UBV RI F λ(v ega)λdλ (as above, the factor of λ comes in for photon counting detectors). Figure 1 shows a plot to demonstrate the difference between the different systems. Why do the different systems exist? While it seems that STMAG and ABNU systems are more straightforward, in practice it is difficult to measure absolute fluxes, and much easier to measure relative fluxes between objects. Hence, historically observations were tied to observations of Vega (or to stars which themselves were tied to Vega), so VEGAMAGs made sense, and the issue of determining physical fluxes boiled down to measuring the physical flux of Vega. Today, in some cases, it may be more accurate to measure the absolute throughput of an instrumental system, and using STMAG or ABNU makes more sense.
7 Astr 535 Class Notes Fall Know that there are several different magnitude systems in use, and understand how they differ. Know when it is important to know what the magnitude system is, and when it isn t Colors Working in magnitudes, the difference in magnitudes between different bandpasses (called the color index, or simply, color) is related to the flux ratio between the bandpasses, i.e., the color. In the UBVRI system, the difference between magnitudes gives the ratio of the fluxes in different bandpasses relative to the ratio of the fluxes of an A0V star in the different bandpasses (for VEGAMAG). Note the typical colors of astronomical objects which are different for the different photometric systems! Which is closer to the UBVRI system, STMAG or ABNU? What would typical colors be in an STMAG or ABNU system? Understand how colors are represented by a difference in magnitude, and recognize how colors expressed in magnitudes are related to the shape of the underlying spectrum, with differences for different magnitude systems Magnitude-flux conversion How would one go about converting Vega-based magnitudes to fluxes? Roughly, just look up the flux of Vega at the center of the passband ( e.g., here (from Bessell et al 1998 or here (see references within), or here; note, however, if the spectrum of the object differs from that of Vega, this won t be perfectly accurate. Given UBVRI magnitudes of an object in the desired band, filter profiles (e.g. Bessell 1990, PASP 102,1181), and absolute spectrophotometry of Vega (e.g., Bohlin & Gilliland 2004, AJ 127, 3508, one can determine the flux. If one wanted to estimate the flux of some object in arbitrary bandpass given just the V magnitude of an object (a common situation used when trying to predict exposures times, see below), this can be done if an estimate of the spectral energy distribution (SED) can be made (e.g., from the spectral type, or more generally, the stellar parameters T eff, log g, and metallicity). Given the filter profiles, one can compute the integral of the SED over the V bandpass, determine the scaling by comparing with the integral of the Vega spectrum over the same bandpass, then use the normalized SED to compute the flux in any desired bandpass. Some possibly useful references for SEDs are: Pickles atlas, MILES library, Bruzual, Persson, Gunn, & Stryker; Hunter, Christian, & Jacoby; Kurucz).
8 Astr 535 Class Notes Fall Things are certainly simpler in the ABNU or STMAG system, and there has been some movement in this direction: the STScI gives STMAG calibrations for HST instruments, and the SDSS photometric system is close to an ABNU system. Note, however, that even when the systems are conceptually well defined, determining the absolute calibration of any photometric system is very difficult in reality, and determining absolute fluxes to the 1% level is very challenging. As a separate note on magnitudes themselves, note that some people, in particular, the SDSS imaging survey, have adopted a modified type of magnitudes, called asinh magnitudes, which behave like normal (also known as Pogson) magnitude for brighter objects, but have different behavior for very faint objects (near the detection threshold); see Lupton, Gunn, & Szalay 1999 AJ 118, 1406 for details. 2.3 Observed fluxes and the count equation What if you are measuring flux with an actual instrument, i.e. counting photons? The intrinsic photon flux from the source is not trivial to determine from the observed photon flux, i.e., the number of photons that you count. The observed flux depends on the area of your photon collector (telescope), photon losses and gains from the Earth s atmosphere (which changes with conditions), and the efficiency of your collection/detection apparatus (which can change with time). Generally, the astronomical signal (which might be a flux or a surface brightness, depending on whether the object is resolved) can be written S = T t Fλ hc λ a λ tel λ inst λ filt λ det λ dλ T ts where S is the observed photon flux (the signal), T is the telescope collecting area, t is the integration time, a λ is the atmospheric transmission (more later) and the other terms refer to the efficiency of various components of the system (telescope, instrument, filter, detector). S is an observed flux rate, i.e. with all of the real details of the observing system included. I refer to this as the count equation. Usually, however, one doesn t use this information to go backward from S to F λ because it is very difficult to measure all of the terms precisely, and some of them (e.g. a, and perhaps some of the system efficiencies) are time-variable; a is also spatially variable. Instead, most observations are performed differentially to a set of other stars of known brightness. If the stars of known brightness are observed in the same observation, then the atmospheric term is (approximately) the same for all stars; this is known as differential photometry. From the photon flux of the object with known brightness, one could determine an exposure efficiency or an effective area for this exposure. Equivalently, and more commonly, one can calculate an instrumental
9 Astr 535 Class Notes Fall magnitude: m = 2.5 log S/t (i.e., normalize by the exposure time to get counts/sec, although this is not strictly necessary) and then determine the zeropoint that needs to be added to give the calibrated magnitude (M, make sure you recognize that this is still an apparent magnitude!): M = m + z Note that in the real world, one has to also consider sensitivity differences (e.g., slightly different filter profiles) between a given experimental setup and the setup used to measure the reference brightnesses, so this is only a first approximation (i.e., the zeropoint may be different for different stars with different spectral properties)! If using instrumental mags that include exposure time normalization (as above), the zeropoint gives the magnitude of a star that will give 1 count/second. If there are no stars of known brightness in the same observation, then calibration must be done against stars in other observations. This then requires that the different effects of the Earth s atmosphere in different locations in the sky be accounted for. This is known as all-sky, or absolute, photometry. To do this requires that the sky is well-behaved, i.e. one can accurately predict the atmospheric throughput as a function of position. This requires that there be no clouds, i.e. photometric weather. Differential photometry can be done in non-photometric weather, hence it is much simpler! Of course, it is always possible to obtain differential photometry and then go back later and obtain absolute photometry of the reference stars. Of course, at some point, someone needs to figure out what the fluxes of the calibrating stars really are, and this requires understanding all of the terms in the count equation. It is challenging, and often, absolute calibration of a system is uncertain to a couple of percent! It is also common to stop with differential photometry, even if there are no stars of known brightness in your field, if you are studying variable objects, i.e. where you are just interested in the change in brightness of an object, not the absolute flux level. In this case, one only has to reference the brightness of the target object relative some other object (or ensemble of objects) in the field that are non-variable. While the count equation isn t usually used for calibration, it is very commonly used for computing the approximate number of photons you will receive from a given source in a given amount of time for a given observational setup. This number is critical to know in order to estimate your expected errors and exposure times in observing proposals, observing runs, etc. Understanding errors in absolutely critical in all sciences, and maybe even more so in astronomy, where objects are faint, photons are scarce, and errors are not at all insignificant. The count equation provides the basis for exposure time calculator (ETC) programs, because it gives an expectation
10 Astr 535 Class Notes Fall of the number of photons that will be received by a given instrument as a function of exposure time. As we will see shortly, this provides the information we need to calculate the uncertainty in the measurement as a function of exposure time. Understand the count equation and the terms in it. Understand the distinction between estimating count rates from an understanding of all of the terms in the count equation vs. measuring the overall throughput (zeropoint) by observing stars of known brightness. Know what instrumental magnitudes and zeropoints are. 2.4 Uncertainties in photon rates For a given rate of emitted photons, there s a probability function which gives the number of photons we detect, even assuming 100% detection efficiency, because of statistical uncertainties. In addition, there may also be instrumental uncertainties. Consequently, we now turn to the concepts of probability distributions, with particular interest in the distribution which applies to the detection of photons. Distributions and characteristics thereof concept of a distribution : define p(x)dx as probability of event occuring in x + dx: p(x)dx = 1 Some definitions relating to values which characterize a distribution: mean µ = xp(x)dx variance σ 2 = (x µ) 2 p(x)dx standarddeviation σ = variance median : mid-point value. xmedian p(x)dx p(x)dx = 1 2 mode : most probable value Note that the geometric interpretation of above quantities depends on the nature of the distribution; although we all carry around the picture of the mean and the variance for a Gaussian distribution, these pictures are not applicable to other distributions, but the quantities are still well-defined.
11 Astr 535 Class Notes Fall Also, note that there is a difference between the sample mean, variance, etc. and the population quantities. The latter apply to the true distribution, while the former are estimates of the latter from some finite sample (N measurements) of the population. The sample quantities are derived from: xi sample mean : x N (xi x) 2 x 2 sample variance = i ( x i ) 2 /N N 1 N 1 The sample mean and variance approach the true mean and variance as N approaches infinity. But note, especially for small samples, your estimate of the mean and variance may differ from their true (population) values (more below)! Now we consider what distribution is appropriate for the detection of photons. The photon distribution can be derived from the binomial distribution, which gives the probability of observing the number, x, of some possible event, given a total number of events n, and the probability of observing the particular event (among all other possibilities) in any single event, p, under the assumption that all events are independent of each other: P (x, n, p) = n!px (1 p) n x x!(n x)! For the binomial distribution, one can derive: mean xp(x)dx = np variance σ 2 (x µ) 2 p(x)dx = np(1 p) Understand the concept of probability distribution functions and basic quantities used to describe them: mean, variance, standard deviation, median, and mode. Understand the difference between population quantities and sample quantities. The Poisson distribution In the case of detecting photons, n is the total number of photons emitted, and p is the probability of detecting a photon during our observation out of the total emitted. We don t know either of these numbers! However, we do know that p << 1 and we know, or at least we can estimate, the mean number detected:. µ = np
12 Astr 535 Class Notes Fall In this limit, the binomial distribution asymtotically approaches the Poisson distribution: P (x, µ) = µx exp µ x! From the expressions for the binomial distribution in this limit, the mean of the distribution is µ, and the variance is variance = x [(x µ) 2 p(x, µ)] variance = np = µ This is an important result. Note that the Poisson distribution is generally the appropriate distribution not only for counting photons, but for any sort of counting experiment where a series of events occurs with a known average rate, and are independent of time since the last event. Figure 2 shows some plots of the Poisson distribution for µ = 2, µ = 10, and µ = Understand the Poisson distribution and when it applies. Know how the variance/standard deviation of the Poisson distribution is related to the mean of the distribution. The normal, or Gaussian, distribution Note, for large µ, the Poisson distribution is well-approximated around the peak by a Gaussian, or normal distribution: P (x, µ, σ) = 1 e (x µ)2 2σ 2 2πσ This is important because it allows us to use simple least squares techniques to fit observational data, because these generally assume normally distributed data. However, beware that in the tails of the distribution, and at low mean rates, the Poisson distribution can differ significantly from a Gaussian distribution. In these cases, least-squares may not be appropriate to model observational data; instead, one might need to consider maximum likelihood techniques instead. The normal distribution is also important because many physical variables seem to be distributed accordingly. This may not be an accident because of the central
13 Astr 535 Class Notes Fall Figure 2: The Poisson distribution for three values of µ
14 Astr 535 Class Notes Fall limit theorem: if a quantity depends on a number of independent random variables with ANY distribution, the quantity itself will be distributed normally (see statistics texts). In observational techniques, we encounter the normal distribution because one important source of instrumental noise, readout noise, is distributed normally. Know what a Gaussian (normal) distribution is, including the full functional form of it. Understand under what circumstances the Poisson distribution is similar to a normal distribution. Importance of error distribution analysis You need to understand the expected uncertainties in your observations in order to: predict the amount of observing time you ll need to get uncertainties as small as you need them to do your science, answer the question: is scatter in observed data consistent with expected uncertainties? If the answer is no, they you know you ve either learned some astrophysics or you don t understand something about your observations. This is especially important in astronomy where objects are faint and many projects are pushing down into the noise as far as possible. Really we can usually only answer this probabilistically. Generally, tests compute the probability that the observations are consistent with an expected distribution (the null hypothesis). You can then look to see if this probability is low, and if so, you can reject the null hypothesis. interpret your results in the context of a scientific prediction Confidence levels For example, say we want to know whether some single point is consistent with expectations, e.g., we see a bright point in multiple measurements of a star, and want to know whether the star flared. Say we have a time sequence with known mean and variance, and we obtain a new point, and want to know whether it is consistent with known distribution? If the form of the probability distribution is known, then you can calculate the probability of getting a measurement more than some observed distance from the mean. In the case where the observed distribution is Gaussian (or approximately so), this is done using the error function (sometimes called erf(x)), which is the integral of a gaussian from some starting value. Some simple guidelines to keep in mind follow (the actual situation often requires more sophisticated statistical tests). First, for Gaussian distributions, you can calculate that 68% of the points should fall within plus or minus one sigma from the
15 Astr 535 Class Notes Fall mean, and 95.3% between plus or minus two sigma from the mean. Thus, if you have a time line of photon fluxes for a star, with N observed points, and a photon noise σ on each measurement, you can test whether the number of points deviating more than 2σ from the mean is much larger than expected. To decide whether any single point is really significantly different, you might want to use more stringent criterion, e.g., 5σ rather than a 2σ criterion; a 5σ has much higher level of significance. On the other hand, if you have far more points in the range 2 4σ brighter or fainter than you would expect, you may also have a significant detection of intensity variations (provided you really understand your uncertainties on the measurements!). Also, note that your observed distribution should be consistent with your uncertainty estimates given the above guidelines. If you have a whole set of points, that all fall within 1σ of each other, something is wrong with your uncertainty estimates (or perhaps your measurements are correlated with each other)! For a series of measurements, one can calculate the χ 2 statistic, and determine how probable this value is, given the number of points. χ 2 = (observed(i) model(i)) 2 σ 2 i For a given value of χ 2 and number of measurements, one can calculate the probability that the measurements are consistent with the model (and the uncertainties are correctly predicted). A quick estimate of the consistency of the model with the observed data points can be made using reduced χ 2, defined as χ 2 divided by the degrees of freedom (number of data points minus number of parameters), which should be near unity of the measurements are consistent wth the model Noise equation: how do we predict expected uncertainties? Signal-to-noise Astronomers often describe uncertainties in terms of the fractional error, e.g. the amplitude of the uncertainty divided by the amplitude of the quantity being measured; often, the inverse of this, referred to as the signal-to-noise ratio is used. Given an estimate the number of photons expected from an object in an observation, we can calulate the signal-to-noise ratio: S N = S σ 2 which is the inverse of the predicted fractional error (N/S). Consider an object with observed photon flux (per unit area and time, e.g. from the count equation above), S, leading to a signal, S = S T t where T is the telescope
16 Astr 535 Class Notes Fall area and t is the exposure time. In the simplest case, the only noise source is Poisson statistics from the source, in which case: σ 2 = S = S T t S N = S = S T t In other words, the S/N increases as the square root of the object brightness, telescope area, efficiency, or exposure time. Note that S is directly observable, so one can calculate the S/N for an observation without knowing the telescope area or exposure time! We ve just broken S down so that you can specifically see the dependence on telescope area and/or exposure time. Understand the concept of S/N and fractional error. Know how S/N depends on the signal for the Poisson-limited case. Background noise A more realistic case includes the noise contributed from Poisson statistics of background light (more on the physical nature of this later), B, which has units of flux per area on the sky (i.e. a surface brightness); note that this is also usually given in magnitudes. B = Bλ hc λ q λ dλ The amount of background in our measurement depends on how we choose to make the measurement (how much sky area we observe). Say we just use an aperture with area, A, so the total observed background counts is AB = AB T t Again, B T t is the directly observable quantity, but we split it into the quantities on which it depends to understand what factors are important in determining S/N. The total number of photons observed, O, is O = S + AB = (S + AB )T t The variance of the total observed counts, from Poisson statistics, is: σ 2 O = O = S + AB = (S + AB )T t To get the desired signal from the object only, we will need to measure separately the total signal and the background signal to estimate: S S T t = O A < B >
17 Astr 535 Class Notes Fall where < B > is some estimate we have obtained of the background surface brightness. The noise in the measurement is σ 2 S σ 2 O = S + AB = (S + AB )T t where the approximation is accurate if the background is determined to high accuracy, which one can do if one measures the background over a large area, thus getting a large number of background counts (with correspondingly small fractional error in the measurement). This leads to a common form of the noise equation: S N = S S + AB Breaking out the dependence on exposure time and telescope area, this is: S N = S S + AB T t In the signal-limited case, S >> B A, we get S N = S = S tt In the background limited case, B A >> S, and S N = S AB = S B A tt As one goes to fainter objects, the S/N drops, and it drops faster when you re background limited. This illustrates the importance of dark-sky sites, and also the importance of good image quality. Consider two telescopes of collecting area, T 1 and T 2. If we observe for the same exposure time on each and want to know how much fainter we can see with the larger telescope at a given S/N, we find: for the signal-limited case, but for the background-limited case. S 2 = T 1 T 2 S 1 S 2 = T1 T 2 S 1
18 Astr 535 Class Notes Fall Understand how Poisson uncertainties in the background contributes to the S/N of an object, and how the background contribution depends on both the brightness of the background but also on the image quality because the amount of background included in the measurement. Instrumental noise In addition to the uncertainties from Poisson statistics (statistical noise), there may be additional terms from instrumental uncertainties. A common example of this that is applicable for CCD detectors is readout noise, which is additive noise (with zero mean!) that comes from the detector and is independent of signal level. For a detector whose readout noise is characterized by σ rn, S N = S S + BApix + σrn 2 if a measurement is made in a single pixel. If an object is spread over N pix pixels, then S N = S S + BA + Npix σrn 2 For large σ rn, the behavior is the same as the background limited case. This makes it clear that if you have readout noise, image quality (and/or proper optics to keep an object from covering too many pixels) is important for maximizing S/N. It is also clear that it is critical to have minimum read-noise for low background applications (e.g., spectroscopy). There are other possible additional terms in the noise equation, arising from things like dark current, digitization noise, uncertainties in sky determination, uncertainties from photometric technique, etc. (we ll discuss some of these later on), but in most applications, the three sources discussed so far signal noise, background noise, and readout noise are the dominant noise sources. Note the applications where one is likely to be signal dominated, background dominated, and readout noise dominated. Understand readout noise and how it represented by a normal distribution with zero mean. Know under what circumstances readout noise is an important contributor to the total noise Error propagation Why are the three uncertainty terms in the noise equation added in quadrature? The measured quantity (S) is a sum of S + B < B > + < R >, where < R >= 0 since
19 Astr 535 Class Notes Fall readout noise has zero mean. The uncertainty in a summed series is computed by addding the individual uncertainties in quadature; in the equation above, we have neglected the uncertainty in < B >. To understand why they add in quadrature, let s consider general error propagation. Now that we know how to estimate uncertainties of observed count rates, let s say we want to make some calculations (e.g., calibration, unit conversion, averaging, conversion to magnitudes, calculation of colors, etc.) using these observations: we need to be able to estimate the uncertainties in the calculated quantities that depend on our measured quantities. Consider what happens if you have several known quantities with known error distributions and you combine these into some new quantity: we wish to know what the uncertainty is in the new quantity. x = f(u, v,...) The question is what is σ x if we know σ u, σ v, etc.? As long as uncertainties are small: x i < x > (u i < u >) ( ) ( ) x x + (v i < v >) +... u v σ 2 x = lim(n > ) 1 N (xi < x >) 2 = σ 2 u ( x u ) 2 + σ 2 v ( ) 2 x + 2σ 2 v uv x u x v +... The last term is the covariance, which relates to whether uncertainties are correlated. σ 2 uv = lim(n > ) 1 N (ui < u >)(v i < v >) If uncertainties are uncorrelated, then σ uv = 0 because there s equal chance of getting opposite signs on v i for any given u i. When working out uncertainties, make sure to consider whether there are correlated errors! If there are, you may be able to reformulate quantities so that they have independent errors: this can be very useful! Examples for uncorrelated errors: addition/subtraction: x = u + v σ 2 x = σ 2 u + σ 2 v In this case, errors are said to add in quadrature.
20 Astr 535 Class Notes Fall multiplication/division: x = uv σ 2 x = v 2 σ 2 u + u 2 σ 2 v logs: x = ln u σ 2 x = σ 2 u/u 2 Note that when dealing with logarithmic quantities, uncertainties in the log correspond to fractional uncertainties in the raw quantity. Distribution of resultant uncertainties When propagating errors, even though you can calculate the variances in the new variables, the distribution of uncertainties in the new variables is not, in general, the same as the distribution of uncertainties in the original variables, e.g. if uncertainties in individual variables are normally distributed, uncertainties in output variable are not necessarily. When two variables are added, however, the output is normally distributed. Know the error propagation formula and how to apply it Determining sample parameters: averaging measurements We ve covered errors in single measurements. Next we turn to averaging measurements. Say we have multiple observations and want the best estimate of the mean and variance of the population, e.g. multiple measurements of stellar brightness. Here we ll define the best estimate of the mean as the value which maximizes the likelihood that our estimate equals the true parent population mean. For equal uncertainties, this estimate just gives our normal expression for the sample mean: xi x = N Using error propagation, the estimate of the error in the sample mean is given by: σ 2 x = σi 2 N = σ2 2 N But what if errors on each observation aren t equal, say for example we have observations made with several different exposure times? Then the optimal determination of the mean is using a: xi σi weightedmean = 2 1 σi 2
21 Astr 535 Class Notes Fall and the estimated error in this value is given by: σ 2 µ = σ2 i σ 4 i ( 1 σ 2 i ) 2 = 1 1 σ 2 i where the σ i s are individual weights/errors (people often talk about the weight of an 1 observation, i.e. : large weight means small error. σi 2 This is a standard result for determining sample means from a set of observations with different weights. However, there can sometimes be a subtlety in applying this formula, which has to do with the question: how do we go about choosing the weights/errors, σ i? We know we can estimate σ using Poisson statistics for a given count rate, but remember that this is a sample variance (which may be based on a single observation!) not the true population variance. This can lead to biases. Consider observations of a star made on three nights, with measurements of 40, 50, and 60 photons. It s clear that the mean observation is 50 photons. However, beware of the being trapped by your error estimates. From each observation alone, you would estimate errors of 40, 50, and 60. If you plug these error estimates into a computation of the weighted mean, you ll get a mean rate of 48.64! Using the individual estimates of the variances, we ll bias values to lower rates, since these will have estimated higher S/N. Note that it s pretty obvious from this example that you should just weight all observations equally. However, note that this certainly isn t always the right thing to do. For example, consider the situation in which you have three exposures of different exposure times, and you are calculating the photon rate (counts/s). Here you probably want to give the longer exposures higher weight (at least, if they are signal or background limited). In this case, you again don t want to use the individual error estimates or you ll introduce a bias. There is a simple solution here also: just weight the observations by the exposure time. However, while this works fine for Poisson errors (variances proportional to count rate), it isn t strictly correct if there are instrumental errors as well which don t scale with exposure time. For example, the presence of readout noise can have this effect; if all exposures are readout noise dominated, then one would want to weight them equally, if readout noise dominates the shorter but not the longer exposures, once would want to weight the longer exposures even higher than expected for the exposure time ratios! The only way to properly average measurements in this case is to estimate a sample mean, then use this value scaled to the appropriate exposure times as the input for the Poisson errors. Another subtlety: averaging counts and converting to magnitudes is not the same as averaging magnitudes!
22 Astr 535 Class Notes Fall Understand how the distinction between sample variance and true variance can lead to biases when calculating a weighted mean, and how to overcome this. Can you split exposures? Although from S/N considerations, one can determine the required number of counts you need (exposure time) to do your science, when observing, one must also consider the question of whether this time should be collected in single or in multiple exposures, i.e. how long individual exposures should be. There are several reasons why one might imagine that it is nicer to have a sequence of shorter exposures rather than one single longer exposure (e.g., tracking, monitoring of photometric conditions, cosmic ray rejection, saturation issues), so we need to consider under what circumstances doing this results in poorer S/N. Consider the object with photon flux S, background surface brightness B, and detector with readout noise σ rn. A single short exposure of time t has a variance: σ 2 S = S T t + B AT t + N pix σ 2 rn If N exposures are summed, the resulting variance will be σ 2 N = Nσ 2 S If a single long exposure of length Nt is taken, we get σ 2 L = S T Nt + B AT Nt + N pix σ 2 rn. The ratio of the noises, or the inverse ratio of the S/N (since the total signal measured is the same in both cases), is σ N = S T Nt + B AT Nt + NN pix σrn 2 σ L S T Nt + B AT Nt + N pix σrn 2 The only difference is in the readout noise term! In the signal- or backgroundlimited regimes, exposures can be added with no loss of S/N. However, if readout noise is significant, then splitting exposures leads to reduced S/N. Understand under what circumstances you can split exposures without a noise penalty, and under what circumstances you cannot. Understand some of the reasons why you might want to split an exposure into shorter pieces.
23 Astr 535 Class Notes Fall Random errors vs systematic errors So far, we ve been discussing random errors. There is an additional, usually more troublesome, type of errors known as systematic errors. These don t occur randomly but rather are correlated with some, possibly unknown, variable relating to your observations, and can have the effect of not just adding spread around the true value that you are trying to measure, but actually measuring the wrong mean. EXAMPLE : flat fielding EXAMPLE : WFPC2 CTE Note also that in some cases, systematic errors can masquerade as random errors in your test observations (or be missing altogether if you don t take data in exactly the same way), but actually be systematic in your science observations. EXAMPLE: flat fielding, subpixel QE variations. Note that error analysis from expected random errors may be the only clue you get to discovering systematic errors. To discover systematic errors, plot residuals vs. everything! Understand the distinction between random and systematic errors.
EXPOSURE TIME ESTIMATION
ASTR 511/O Connell Lec 12 1 EXPOSURE TIME ESTIMATION An essential part of planning any observation is to estimate the total exposure time needed to satisfy your scientific goal. General considerations
Ay Fall 2012 Imaging and Photometry Part I
Ay 122 - Fall 2012 Imaging and Photometry Part I (Many slides today c/o Mike Bolte, UCSC) Imaging and Photometry Now essentially always done with imaging arrays (e.g., CCDs); it used to be with single-channel
The magnitude system. ASTR320 Wednesday January 30, 2019
The magnitude system ASTR320 Wednesday January 30, 2019 What we measure: apparent brightness How bright a star appears to be in the sky depends on: How bright it actually is Luminosity and its distance
Astr 323: Extragalactic Astronomy and Cosmology. Spring Quarter 2014, University of Washington, Željko Ivezić. Lecture 1:
Astr 323: Extragalactic Astronomy and Cosmology Spring Quarter 2014, University of Washington, Željko Ivezić Lecture 1: Review of Stellar Astrophysics 1 Understanding Galaxy Properties and Cosmology The
Photometric Systems. Neil Phillips 08 April Introduction 1. 2 Magnitudes and Colours 2
Photometric Systems Neil Phillips 8 April 28 Contents 1 Introduction 1 2 Magnitudes and Colours 2 3 Performing photometry 2 3.1 The System Response R λ....................................
Photons. Observational Astronomy 2018 Part 1 Prof. S.C. Trager
Photons Observational Astronomy 2018 Part 1 Prof. S.C. Trager Wavelengths, frequencies, and energies of photons Recall that λν=c, where λ is the wavelength of a photon, ν is its frequency, and c is the
Detection ASTR ASTR509 Jasper Wall Fall term. William Sealey Gosset
ASTR509-14 Detection William Sealey Gosset 1876-1937 Best known for his Student s t-test, devised for handling small samples for quality control in brewing. To many in the statistical world "Student" was
STATISTICS OF OBSERVATIONS & SAMPLING THEORY. Parent Distributions
ASTR 511/O Connell Lec 6 1 STATISTICS OF OBSERVATIONS & SAMPLING THEORY References: Bevington Data Reduction & Error Analysis for the Physical Sciences LLM: Appendix B Warning: the introductory literature
WFC3/UVIS Photometric Transformations
Instrument Science Report WFC3 2014-016 WFC3/UVIS Photometric Transformations Kailash Sahu, Susana Deustua and Elena Sabbi January 05, 2017 ABSTRACT We provide photometric transformation coefficients for
point, corresponding to the area it cuts out: θ = (arc length s) / (radius of the circle r) in radians Babylonians:
Astronomische Waarneemtechnieken (Astronomical Observing Techniques) 1 st Lecture: 1 September 11 This lecture: Radiometry Radiative transfer Black body radiation Astronomical magnitudes Preface: The Solid
Astr 511: Galactic Astronomy. Winter Quarter 2015, University of Washington, Željko Ivezić. Lecture 1:
Astr 511: Galactic Astronomy Winter Quarter 2015, University of Washington, Željko Ivezić Lecture 1: Review of Stellar Astrophysics (and other useful stuff) 1 Understanding Galaxy Properties and the Milky
Fluxes. 1 March 2016
Fluxes 1 March 2016 intrinsic luminosity L = energy emi5ed per second (erg/s) intrinsic luminosity L = energy emi5ed per second (erg/s) flux F = energy received per unit area F = L / 4π d 2 (erg/s/cm 2
Ay 20 Basic Astronomy and the Galaxy Problem Set 2
Ay 20 Basic Astronomy and the Galaxy Problem Set 2 October 19, 2008 1 Angular resolutions of radio and other telescopes Angular resolution for a circular aperture is given by the formula, θ min = 1.22λ
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
Measurement of the stellar irradiance
Measurement of the stellar irradiance Definitions Specific Intensity : (monochromatic) per unit area normal to the direction of radiation per unit solid angle per unit wavelength unit (or frequency) per
summary of last lecture
radiation specific intensity flux density bolometric flux summary of last lecture Js 1 m 2 Hz 1 sr 1 Js 1 m 2 Hz 1 Js 1 m 2 blackbody radiation Planck function(s) Wien s Law λ max T = 2898 µm K Js 1 m
Astronomical "color"
Astronomical "color" What color is the star Betelgeuse? It's the bright star at upper left in this picture of Orion taken by a student at the RIT Observatory. Orange? Red? Yellow? These are all reasonable
Photometric Products. Robert Lupton, Princeton University LSST Pipeline/Calibration Scientist PST/SAC PST/SAC,
Photometric Products Robert Lupton, Princeton University LSST Pipeline/Calibration Scientist PST/SAC 2018-02-27 PST/SAC, 2018-02-27 LSST2017 Tucson, AZ August 14-18, 2017 1 Outline Photometry Point Sources
Astronomical imagers. ASTR320 Monday February 18, 2019
Astronomical imagers ASTR320 Monday February 18, 2019 Astronomical imaging Telescopes gather light and focus onto a focal plane, but don t make perfect images Use a camera to improve quality of images
AstroBITS: Open Cluster Project
AstroBITS: Open Cluster Project I. Introduction The observational data that astronomers have gathered over many years indicate that all stars form in clusters. In a cloud of hydrogen gas, laced with helium
Lecture Outline: Chapter 5: Telescopes
Lecture Outline: Chapter 5: Telescopes You don t have to know the different types of optical reflecting and refracting telescopes. It is important to understand the difference between imaging, photometry,
Theoretical quantities: blackbody radiation
Theoretical quantities: blackbody radiation Magnitudes are observed quantities; that is, in practice, optical astronomers typically 1. take pictures of stars 2. measure the apparent brightness of each
Atmospheric Extinction
Atmospheric Extinction Calibrating stellar photometry requires correction for loss of light passing through the atmosphere. Atmospheric Rayleigh and aerosol scattering preferentially redirects blue light
OPTICAL PHOTOMETRY. Observational Astronomy (2011) 1
OPTICAL PHOTOMETRY Observational Astronomy (2011) 1 The optical photons coming from an astronomical object (star, galaxy, quasar, etc) can be registered in the pixels of a frame (or image). Using a ground-based
Astronomy 15 - Problem Set Number 4
Astronomy 15 - Problem Set Number 4 1) Suppose one were to observe a star for parallax with a telescope of 15 meters focal length. Suppose the pixels of the CCD are 24 microns in size (1 micron = 10 6
Open Cluster Research Project
Open Cluster Research Project I. Introduction The observational data indicate that all stars form in clusters. In a cloud of hydrogen gas, laced with helium and a trace of other elements, something triggers
1. Why photons? 2. Photons in a vacuum
Photons and Other Messengers 1. Why photons? Ask class: most of our information about the universe comes from photons. What are the reasons for this? Let s compare them with other possible messengers,
You, too, can make useful and beautiful astronomical images at Mees: Lesson 3
You, too, can make useful and beautiful astronomical images at Mees: Lesson 3 Calibration and data reduction Useful references, besides Lessons 1 and 2: The AST 142 Projects manual: http://www.pas.rochester.edu/~dmw/ast142/projects/project.pdf
Data Processing in DES
Data Processing in DES Brian Yanny Oct 28, 2016 http://data.darkenergysurvey.org/fnalmisc/talk/detrend.p Basic Signal-to-Noise calculation in astronomy: Assuming a perfect atmosphere (fixed PSF of p arcsec
Chapter 10: Unresolved Stellar Populations
Chapter 10: Unresolved Stellar Populations We now consider the case when individual stars are not resolved. So we need to use photometric and spectroscopic observations of integrated magnitudes, colors
arxiv:astro-ph/ v3 7 Apr 2005 Space Telescope Science Institute 2, 3700 San Martin Drive, Baltimore, MD 21218, U.S.A.
Accepted at Publications of the Astronomical Society of the Pacific (June 2005 issue) A cross-calibration between Tycho-2 photometry and HST spectrophotometry arxiv:astro-ph/0504085v3 7 Apr 2005 J. Maíz
A Calibration Method for Wide Field Multicolor. Photometric System 1
A Calibration Method for Wide Field Multicolor Photometric System 1 Xu Zhou,Jiansheng Chen, Wen Xu, Mei Zhang... Beijing Astronomical Observatory,Chinese Academy of Sciences, Beijing 100080, China Beijing
9. Evolution with redshift - z > 1.5. Selection in the rest-frame UV
11-5-10see http://www.strw.leidenuniv.nl/ franx/college/galaxies10 10-c09-1 11-5-10see http://www.strw.leidenuniv.nl/ franx/college/galaxies10 10-c09-2 9. Evolution with redshift - z > 1.5 Selection in
Introduction to Error Analysis
Introduction to Error Analysis Part 1: the Basics Andrei Gritsan based on lectures by Petar Maksimović February 1, 2010 Overview Definitions Reporting results and rounding Accuracy vs precision systematic
X Rays must be viewed from space used for detecting exotic objects such as neutron stars and black holes also observing the Sun.
6/25 How do we get information from the telescope? 1. Galileo drew pictures. 2. With the invention of photography, we began taking pictures of the view in the telescope. With telescopes that would rotate
Galaxy formation and evolution. Astro 850
Galaxy formation and evolution Astro 850 Introduction What are galaxies? Systems containing many galaxies, e.g. 10 11 stars in the Milky Way. But galaxies have different properties. Properties of individual
The Basics of Light. Sunrise from the Space Shuttle, STS-47 mission. The Basics of Light
The Basics of Light The sun as it appears in X-ray light (left) and extreme ultraviolet light (right). Light as energy Light is remarkable. It is something we take for granted every day, but it's not something
Lecture 3: Specific Intensity, Flux and Optical Depth
Lecture 3: Specific Intensity, Flux and Optical Depth We begin a more detailed look at stellar atmospheres by defining the fundamental variable, which is called the Specific Intensity. It may be specified
5.6 Spectrophotometry and Magnitudes
5.6. SPECTROPHOTOMETRY AND MAGNITUDES 87 a constant optical depth, τ = χ L, across the beam cross section (so that the absorption is independent of the point of origin of a beam element emitted from the
AST111 PROBLEM SET 4 SOLUTIONS. Ordinarily the binary has a magnitude of 10 and this is due to the brightness of both stars.
AST111 PROBLEM SET 4 SOLUTIONS Homework problems 1. On Astronomical Magnitudes You observe a binary star. Ordinarily the binary has a magnitude of 10 and this is due to the brightness of both stars. The
AS750 Observational Astronomy
Lecture 9 0) Poisson! (quantum limitation) 1) Diffraction limit 2) Detection (aperture) limit a)simple case b)more realistic case 3) Atmosphere 2) Aperture limit (More realistic case) Aperture has m pixels
arxiv:astro-ph/ v1 17 Oct 2002
The K correction David W. Hogg 1,2, Ivan K. Baldry 3, Michael R. Blanton 1, and Daniel J. Eisenstein 4 ABSTRACT arxiv:astro-ph/0210394 v1 17 Oct 2002 The K correction corrects for the fact that sources
Flux Units and Line Lists
APPENDIX 2 Flux Units and Line Lists In This Appendix... Infrared Flux Units / 255 Formulae / 258 Look-up Tables / 259 Examples / 266 Infrared Line Lists / 266 In this chapter we provide a variety of background
Point Spread Functions. Aperture Photometry. Obs Tech Obs Tech 26 Sep 2017
Point Spread Functions & Aperture Photometry Obs Tech Obs Tech 26 Sep 2017 Point Spread Functions (PSFs) Stars are far away (typically d >> 3x10 15 km) Stars are relatively small (D ~ 10 6 km) Angular
ASTRO 114 Lecture Okay. What we re going to discuss today are what we call radiation laws. We ve
ASTRO 114 Lecture 15 1 Okay. What we re going to discuss today are what we call radiation laws. We ve been spending a lot of time talking about laws. We ve talked about gravitational laws, we ve talked
Fitting a Straight Line to Data
Fitting a Straight Line to Data Thanks for your patience. Finally we ll take a shot at real data! The data set in question is baryonic Tully-Fisher data from http://astroweb.cwru.edu/sparc/btfr Lelli2016a.mrt,
Physics Oct A Quantum Harmonic Oscillator
Physics 301 5-Oct-2005 9-1 A Quantum Harmonic Oscillator The quantum harmonic oscillator (the only kind there is, really) has energy levels given by E n = (n + 1/2) hω, where n 0 is an integer and the
Light and Matter(LC)
Light and Matter(LC) Every astronomy book that I ve seen has at least one chapter dedicated to the physics of light. Why are astronomers so interested in light? Everything* that we know about Astronomical
The Accuracy of WFPC2 Photometric Zeropoints
The Accuracy of WFPC2 Photometric Zeropoints Inge Heyer, Marin Richardson, Brad Whitmore, Lori Lubin July 23, 2004 ABSTRACT The accuracy of WFPC2 photometric zeropoints is examined using two methods. The
Properties of Electromagnetic Radiation Chapter 5. What is light? What is a wave? Radiation carries information
Concepts: Properties of Electromagnetic Radiation Chapter 5 Electromagnetic waves Types of spectra Temperature Blackbody radiation Dual nature of radiation Atomic structure Interaction of light and matter
Active Galactic Nuclei OIII
Active Galactic Nuclei In 1908, Edward Fath (1880-1959) observed NGC 1068 with his spectroscope, which displayed odd (and very strong) emission lines. In 1926 Hubble recorded emission lines of this and
Absolute Flux Calibration for STIS First-Order, Low-Resolution Modes
Instrument Science Report STIS 97-14 Absolute Flux Calibration for STIS First-Order, Low-Resolution Modes Ralph Bohlin, Space Telescope Science Institute Nicholas Collins, Hughes STX/LASP/GSFC Anne Gonnella,
Direct imaging of extra-solar planets
Chapter 6 Direct imaging of extra-solar planets Direct imaging for extra-solar planets means that emission from the planet can be spatially resolved from the emission of the bright central star The two
TECHNICAL REPORT. Doc #: Date: Rev: JWST-STScI , SM-12 August 31, Authors: Karl Gordon, Ralph Bohlin. Phone:
When there is a discrepancy between the information in this technical report and information in JDox, assume JDox is correct. TECHNICAL REPORT Title: Title: JWST Absolute Flux Calibration II: Expanded
Detecting Cosmic Rays in Infrared Data
The 2010 STScI Calibration Workshop Space Telescope Science Institute, 2010 Susana Deustua and Cristina Oliveira, eds. Detecting Cosmic Rays in Infrared Data Rachel E. Anderson and Karl Gordon Space Telescope
Chapter 26: Cosmology
Chapter 26: Cosmology Cosmology means the study of the structure and evolution of the entire universe as a whole. First of all, we need to know whether the universe has changed with time, or if it has
Originality in the Arts and Sciences: Lecture 2: Probability and Statistics
Originality in the Arts and Sciences: Lecture 2: Probability and Statistics Let s face it. Statistics has a really bad reputation. Why? 1. It is boring. 2. It doesn t make a lot of sense. Actually, the
Lecture 12: Distances to stars. Astronomy 111
Lecture 12: Distances to stars Astronomy 111 Why are distances important? Distances are necessary for estimating: Total energy released by an object (Luminosity) Masses of objects from orbital motions
optical / IR: photon counting flux density or magnitude corresponds to number of electrons per second (mean rate)
optical / IR: photon counting flux density or magnitude corresponds to number of electrons per second (mean rate) N electrons/sec = ɛ F λ λa hc/λ 0 efficiency factor flux density x bandpass x collecting
BRIDGE CIRCUITS EXPERIMENT 5: DC AND AC BRIDGE CIRCUITS 10/2/13
EXPERIMENT 5: DC AND AC BRIDGE CIRCUITS 0//3 This experiment demonstrates the use of the Wheatstone Bridge for precise resistance measurements and the use of error propagation to determine the uncertainty
First On-orbit Measurements of the WFC3-IR Count-rate Non-Linearity
First On-orbit Measurements of the WFC3-IR Count-rate Non-Linearity A. G. Riess March 10, 2010 ABSTRACT Previous HgCdTe detectors on HST have suffered from a count-rate dependent non-linearity, motivating
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
Part II: Shining a Light on Visual Magnitude
Part II: Shining a Light on Visual Magnitude Patrick North, Novarah Kazmi Who doesn t like staring up at the stars? One of the most humbling experiences is looking up at the night sky and staring at the
Astronomical image reduction using the Tractor
the Tractor DECaLS Fin Astronomical image reduction using the Tractor Dustin Lang McWilliams Postdoc Fellow Carnegie Mellon University visiting University of Waterloo UW / 2015-03-31 1 Astronomical image
Midterm Exam. IT Posting scores Finding out about missed questions Reminder about dropping lowest of 3
Midterm Exam #&%?@)#$! IT Posting scores Finding out about missed questions Reminder about dropping lowest of 3 Nature of Light 10/3 Apparent versus Actual Brightness 10/6 Electromagnetic Spectrum of Light
Gemini Integration Time Calculators
Gemini Integration Time Calculators Phil Puxley SPE-C-G0076 Version 1.1; June 1998 Version History: 1.0 June 1998; first draft. 1.1 June 1998; following comments from Ted von Hippel and Joe Jensen. 1.
Statistical Tools and Techniques for Solar Astronomers
Statistical Tools and Techniques for Solar Astronomers Alexander W Blocker Nathan Stein SolarStat 2012 Outline Outline 1 Introduction & Objectives 2 Statistical issues with astronomical data 3 Example:
A Cross-Calibration between Tycho-2 Photometry and Hubble Space Telescope Spectrophotometry
Publications of the Astronomical Society of the Pacific, 117:615 619, 2005 June 2005. he Astronomical Society of the Pacific. All rights reserved. Printed in U.S.A. A Cross-Calibration between ycho-2 Photometry
1. Give short answers to the following questions. a. What limits the size of a corrected field of view in AO?
Astronomy 418/518 final practice exam 1. Give short answers to the following questions. a. What limits the size of a corrected field of view in AO? b. Describe the visibility vs. baseline for a two element,
Today in Astronomy 142: observations of stars
Today in Astronomy 142: observations of stars What do we know about individual stars?! Determination of stellar luminosity from measured flux and distance Magnitudes! Determination of stellar surface temperature
Statistics of Accumulating Signal
Instrument Science Report NICMOS 98-008 Statistics of Accumulating Signal W.B. Sparks, STScI May 1998 ABSTRACT NICMOS detectors accumulate charge which can be non-destructively. The charge accumulation
de = j ν dvdωdtdν. (1)
Transfer Equation and Blackbodies Initial questions: There are sources in the centers of some galaxies that are extraordinarily bright in microwaves. What s going on? The brightest galaxies in the universe
The Hertzsprung-Russell Diagram
The Hertzsprung-Russell Diagram Name: Date: 1 Introduction As you may have learned in class, the Hertzsprung-Russell Diagram, or the HR diagram, is one of the most important tools used by astronomers:
Chapter 6: Transforming your data
Why is transformation necessary? Chapter 6: Transforming your data The AAVSO International Database is composed of data collected from many different observers, at different times, from around the globe.
Astro 1050 Wed. Feb. 18, 2015
Astro 1050 Wed. Feb. 18, 2015 Today: Begin Chapter 5: Light the Cosmic Messenger For Friday: Study for Test #1 Be sure to bring green bubble sheet, #2 pencil and a calculator. 1 Chapter 5: Light, the Cosmic
We start with a reminder of a few basic concepts in probability. Let x be a discrete random variable with some probability function p(x).
1 Probability We start with a reminder of a few basic concepts in probability. 1.1 discrete random variables Let x be a discrete random variable with some probability function p(x). 1. The Expectation
How to calibrate interferometric data
National Aeronautics and Space Administration Jet Propulsion Laboratory California Institute of Technology How to calibrate interferometric data Rachel Akeson 28 July 2006 Outline Calibration basics System
ACS CCDs UV and narrowband filters red leak check
Instrument Science Report ACS 2007-03 ACS CCDs UV and narrowband filters red leak check Marco Chiaberge and Marco Sirianni May 01, 2007 ABSTRACT We present results of the observations of the star 15 Mon,
Counting Photons to Calibrate a Photometer for Stellar Intensity Interferometry
Counting Photons to Calibrate a Photometer for Stellar Intensity Interferometry A Senior Project Presented to the Department of Physics California Polytechnic State University, San Luis Obispo In Partial
Project for Observational Astronomy 2018/2019: Colour-magnitude diagram of an open cluster
Project for Observational Astronomy 018/019: Colour-magnitude diagram of an open cluster Søren S. Larsen November 9, 018 1 1 Colour-magnitude diagram for an open cluster 1.1 Background The colour-magnitude
Cross-Talk in the ACS WFC Detectors. I: Description of the Effect
Cross-Talk in the ACS WFC Detectors. I: Description of the Effect Mauro Giavalisco August 10, 2004 ABSTRACT Images acquired with the Wide Field Channel (WFC) of the Advanced Camera for Surveys (ACS) are
Transiting Exoplanet in the Near Infra-red for the XO-3 System
Transiting Exoplanet in the Near Infra-red for the XO-3 System Nathaniel Rodriguez August 26, 2009 Abstract Our research this summer focused on determining if sufficient precision could be gained from
High Signal-to-Noise, Differential NICMOS Spectrophotometry
Instrument Science Report NICMOS 2003-001 High Signal-to-Noise, Differential NICMOS Spectrophotometry R.L. Gilliland, S. Arribas January 7, 2003 ABSTRACT We report analysis for NICMOS CAL/9642, High S/N
A Random Walk Through Astrometry
A Random Walk Through Astrometry Astrometry: The Second Oldest Profession George H. Kaplan Astronomical Applications Department Astrometry Department U.S. Naval Observatory Random Topics to be Covered
Prof. Jeff Kenney Class 4 May 31, 2018
Prof. Jeff Kenney Class 4 May 31, 2018 Which stellar property can you estimate simply by looking at a star on a clear night? A. distance B. diameter C. luminosity D. surface temperature E. mass you can
FLAT FIELDS FROM THE MOONLIT EARTH
Instrument Science Report WFPC2 2008-01 FLAT FIELDS FROM THE MOONLIT EARTH R. C. Bohlin, J. Mack, and J. Biretta 2008 February 4 ABSTRACT The Earth illuminated by light from the full Moon was observed
Galaxies 626. Lecture 10 The history of star formation from far infrared and radio observations
Galaxies 626 Lecture 10 The history of star formation from far infrared and radio observations Cosmic Star Formation History Various probes of the global SF rate: ρ* (z) M yr 1 comoving Mpc 3 UV continuum
Chapter 5: Telescopes
Chapter 5: Telescopes You don t have to know different types of reflecting and refracting telescopes. Why build bigger and bigger telescopes? There are a few reasons. The first is: Light-gathering power:
An Algorithm for Correcting CTE Loss in Spectrophotometry of Point Sources with the STIS CCD
An Algorithm for Correcting CTE Loss in Spectrophotometry of Point Sources with the STIS CCD Ralph Bohlin and Paul Goudfrooij August 8, 2003 ABSTRACT The correction for the change in sensitivity with time
Stellar Photometry: II. Transforming. Ast 401/Phy 580 Fall 2015
Stellar Photometry: II. Transforming Ast 401/Phy 580 Fall 2015 Summary We ve learned how to measure instrumental magnitudes for a star using aperture photometry or PSF-fitting: (a) Add up all of the counts
Photometric Techniques II Data analysis, errors, completeness
Photometric Techniques II Data analysis, errors, completeness Sergio Ortolani Dipartimento di Astronomia Universita di Padova, Italy. The fitting technique assumes the linearity of the intensity values
Astrophysics I ASTR:3771 Fall 2014
Astrophysics I ASTR:3771 Fall 2014 Prof. Kaaret 702 Van Allen Hall 335-1985 philip-kaaret@uiowa.edu Office hours: Tuesday 1:30-2:30 pm, Wednesday 9-11 am, or by appointment, or drop by. Topics We will
Measuring the Properties of Stars (ch. 17) [Material in smaller font on this page will not be present on the exam]
Measuring the Properties of Stars (ch. 17) [Material in smaller font on this page will not be present on the exam] Although we can be certain that other stars are as complex as the Sun, we will try to
4. Energy, Power, and Photons
4. Energy, Power, and Photons Energy in a light wave Why we can often neglect the magnetic field Poynting vector and irradiance The quantum nature of light Photon energy and photon momentum An electromagnetic
3 Effects of the earth s atmosphere
Astr 535 Class Notes Fall 2017 29 3 Effects of the earth s atmosphere The earth s atmosphere has several different effects: it emits light, it absorbs light, it shifts the apparent direction of incoming
Hertzprung-Russel and colormagnitude. ASTR320 Wednesday January 31, 2018
Hertzprung-Russel and colormagnitude diagrams ASTR320 Wednesday January 31, 2018 H-R diagram vs. Color- Magnitude Diagram (CMD) H-R diagram: Plot of Luminosity vs. Temperature CMD: Plot of magnitude vs.
37-6 Watching the electrons (matter waves)
37-6 Watching the electrons (matter waves) 1 testing our proposition: the electrons go either through hole 1 or hole 2 add a very strong light source behind walls between two holes, electrons will scatter
1 The Preliminary Processing
AY 257 Modern Observational Techniques...23 1 The Preliminary Processing Frames must be corrected for a bias level and quantum efficiency variations on all scales. For a minority of CCDs and most near-ir
XMM-Newton Optical-UV Monitor: introduction and calibration status OM instrument and calibration
XMM-Newton Optical-UV Monitor: introduction and calibration status OM instrument and calibration Antonio Talavera XMM-Newton Science Operation Centre, ESAC, ESA OM: Instrument Description 30 cm Ritchey-Chretien
Telescopes. Optical Telescope Design. Reflecting Telescope
Telescopes The science of astronomy was revolutionized after the invention of the telescope in the early 17th century Telescopes and detectors have been constantly improved over time in order to look at