Astronomy Lab

Determination of Stellar Magnitudes from Digital Images - Revised to use DSLR images

OBJECTIVES
Note: We will do Procedure I on one day, then return later to Procedure II to calculate magnitudes.

INTRODUCTION. Astronomers describe the apparent brightness of stars in magnitudes (see Comins p. 179f). Hipparchus developed the magnitude scale in ancient Greece. Arbitrarily, he said that a 6th magnitude star was 100 times fainter than a 1st magnitude star. A magnitude difference means a certain ratio in intensities. Let r be the ratio of intensities for stars differing by one magnitude. r x r = r2 is the ratio of intensities for stars differing by two magnitudes, etc. r5 = 100. Thus, 
 The magnitude scale is logarithmic, i. e. the magnitude of a star is related to the exponent of the intensity. Another strange characteristic of the magnitude scale is that it is backwards. A magnitude 6 star is 100 times fainter than a magnitude 1 star. Thus magnitude is a logarithmic measure of "faintness".
 
 

Suppose we have a standard star of magnitude mo which yields an intensity reading of Io from a photometer. Another star reads I on the same photometer. A simple formula for the magnitude m is given by

m = m0 - 2.5*log(I/I0)

Although this formula is rather cryptic involving logarithms, it presents no problem for a calculator or computer spreadsheet.

PROCEDURE I. Measure brightness with digital imaging. The class has previously photographed a number of stars of varying brightness using a digital SLR camera.  Cassiopeia, Fall, 2008.   You will load digital images of these stars on machines in the Spidel Computer Laboratory or Physics Laboratory using Matlab.  Matlab is a sophisticated and powerful analytical package that is programmable.  In future years, a photometry package that runs without the expensive Matlab program may be available. 

Click here to examine a small portion of a digital star image.

We will measure selected stars in the constellation Cassiopeia - photographed Oct. November 6, 2008 by Gracie McCarroll, Chandler Jones, and Bayle Owens.


PART II Use spreadsheet to calculate the magnitudes. After the table of the star distances and brightnesses has been completed, enter all the data into a table in a spreadsheet.  All computers can run the spreadsheet (excel).
   m = mDelta - 2.5*log(I/IDelta) or m = 2.7 - 2.5*log(I/IDelta)
For the "I" in the above formula you point and click to the cell containing the ADU brightness values.  This will be demonstrated in class.
    • Be very careful how you point to IDelta in the above formula.  Remember that you entered another row for IDelta. You must use "absolute" addressing.  When you point (or highlight) the cell coordinates in the formula, you should "absolutize" the coordinates.  This is done using the function key <F4> in Excel.  You will notice that the cell coordinates have a dollar sign ($) placed in front of the row and column numbers.  If the <F4> doesn't work you may explicitly type the $ signs in front of both row number and column number in the formula.  The reason for this will become apparent in the next step of replicating the formula.
where d is the distance from the earth in LY.  Enter the column heading and the formula for the first line of data.  Copy the formula down the column as before.
10^((4.83 - M)/2.5)
Where 4.83 is the intrinsic magnitude of the sun.  In this formula, (4.83 - M)/2.5 is the exponent, or the power of 10.  This is also taking the antilog of a number.  When typing this strange formula in the spread sheet, you have to write it in the following format:  10^((4.83-M)/2.5).  Use the pointing feature of the spreadsheet to point to the intrinsic magnitude cell.  You must be very careful with the parentheses.  Note that first you subtract numbers, then divide the difference by 2.5, then rais 10 to the result!  The brightest star should be larger than 30,000 times brighter than the sun.  The "^" operator means "raise to a power".

Summarize the data.  Discuss the results, and any surprises you observe.  Compare the brightest, second brightest, and faintest star (as viewed from earth) with your "eyeball" predictions made earlier.  Why is the star with the greatest intrinsic magnitude not very bright in the photograph?    For the star(s) that is(are) fainter than the sun where does(do) it(they) rank in the apparent brightness?


Several stars in Cassiopeia are very interesting.  Gamma Cas is a variable star with very interesting properties.  See the link to the AAVSO discussion.  Beta Cas is another variable star.  Bright variable stars are very difficult for professional astronomers to measure because they saturate the detectors.  A possible research project would be for a student to measure the brightness of  some of these stars over an extended period using the digital camera technique.  This would require special attention to detail, correcting carefully for flat fielding, etc., but students should be able to make a research contribution to science.

Cassiopeia


Star

Distance (LY)

alpha
229
beta
54.5
gamma
613
chi
204
delta
99
epsilon
442
eta
19.4
kappa
4129
lambda
355
mu
24.6
theta
137
upsilon1
229
upsilon2
206
zeta
597

Cassiopeia - Nov. 6, 2008.  Sum of 20 exposures at 15 sec. using DSLR Camera "piggy-backed" onto Meade Clock-drive telescope.  Focal length: 55 mm.  Photo made by Gracie McCarroll, Chandler Jones, and Bayle Owens.