Ricky Leon Murphy The Appendices: Project – HET606 Appendix 1 – Known Semester 2 – 2004 Appendix 2 – Location of Tau Bootis Appendix 3 – Location of HD 209458 Appendix 4 – Image Reduction

Search for Other Worlds Introduction As of September 2004, there are 136 known outside our (http://exoplanets.org). These extra-solar planets, or exoplanets, are one of the most current and highly studied subjects in today and it is one of the very few subjects that involve both amateur and professional astronomers. The huge telescopes perched atop Mauna Kea in Hawaii are pointed at these objects, as are 8” telescopes purchased from the local shopping mall – and many others around the world. Why are we finding these planets now if only 8” telescopes can detect them? Simple; we know what we are looking for and we have better tools to get the job done. While telescope size does not seem to matter with the search and study of exoplanets, it’s what you do not hear about that really matters: improved sensitivity in CCD cameras, improved resolution in spectroscopy, and fast computers to perform the mathematics. The goal of gathering repeatable data is very important when studying exoplanets. The rewards of such study carry implications across the board in Astronomy: we can learn about our own solar system and test the theories of solar system formation and evolution, improve the sensitivity to detect small -like planets, and possibly provide targets for the spaced based telescopes and SETI projects; however, the most important implication is that perhaps for the first time in history, amateurs and professionals from around the world are engaged in this subject and working together to share the data. The greatest reason for this collaboration is telescope time: there are only a finite number of professional telescopes with tightly guarded schedules that limit prolonged data collection. Amateurs have all the time they can spend and the equipment to help. There are several methods of detecting these exoplanets; however, amateur astronomers are only capable of performing only two – the measurement of and the method, both will be discussed in detail. The Formation of a Solar System

1 The foundation is the base in which ideas are built upon. In this case, the foundation of this project will be to briefly visit the current theory of how a solar system is formed. This is important because this gives clues for us to know what we are looking for when we study other systems. The formation of our Solar System can be traced all the way back to the Big Bang. With hydrogen being the most abundant element in the Universe, clouds of hydrogen began to form under its own gravity. By gravity and rotation, these clouds compressed to form the first (and ) of the Universe – Population II stars. Population II is the designation for stars that do not contain heavy elements – that is heavier than helium. The natural cycle of these stars resulted in supernova explosions that introduced the heavier elements into the interstellar spaces due to intense heat and pressure of the explosion. While hydrogen is still the most abundant element, other heavier elements are now present – elements like iron, carbon, silicone and many others. Supernova stimulates nearby hydrogen clouds and introduces the heavy elements. As the stimulated clouds collapse, they form something a little different, a proto-planetary disk, as well as the central proto-star (figure 1). When these metal-rich stars – called Population I – are formed, they may host a ring of molecular material that begins to collide with one another resulting in the sweeping up of material.

Figure 1: The molecular cloud begins to collapse, eventually forming a dense area at the heart of the cloud. The continual contraction raises the of the cloud causing rotation. Eventually, there will be enough heat and density that a T Tauri star will form at the heart of this cloud with the remaining disk material possibly forming planets. A: the slowly rotating molecular cloud B: Faster rotation with denser and hotter central region C: Faster rotation with T Tauri star at the center of rotation

By retaining orbital momentum, the shape of the disk is formed. The consequence of this is the formation of larger objects that also contain gravity and also spin as a result of the orbiting momentum and the further collection of material as it “sweeps” through the ring.

2 These objects are called planetesimals, and continue to collect more material as they continue to the host star. Proof to support this Solar theory has been discovered by the study of ancient meteorites - called chondrites - found on Earth (Beatty et al, 1999). The material makeup of the chondrites, which originate from space, is found to contain the same elements found on Earth – one of which is an isotope of hydrogen called deuterium. Additional proof to this theory comes from remarkable images taken by the . Looking deep into the heart of the Great Nebula (figure 2), the Hubble Telescope was able to spy tiny solar systems in the making (figure 3). These images show the T Tauri star in the center of a protoplanetary disks, or proplyds. This offers us a remarkable look into the very earliest history into the formation of a stellar system (Figures 4, 5, and 6).

Figure 2: This beautiful image of the Orion nebula was captured using special filters by amateur astronomer Russell Crowman using a 14.5 inch telescope and a Santa Barbara Instrument Group ST-11000 CCD camera.

Figure 3: This close-up image of the center of Orion nebula – care of the Hubble Space Telescope – shows several knotted looking objects. These are the proplyds.

3 While these proplyds look very impressive, the host T Tauri stars are stars that have not initiated hydrogen fusing – and must be observed in the near infrared. Much of the material that surrounds the proto-star will be blown away by the shockwave resulting from the initiation of hydrogen fusion at the heart of the star (Ostlie and Carroll, 1996).

Figures 4, 5, and 6 (clockwise): These are Hubble Telescope close-up images of three of these knots. They show disks of dusty material surrounding their host T-Tauri stars. Because of the density of the dust, these images are photographed in the near-infrared so the host star is visible.

Evidence of possible planetary formation Have you ever looked into a telescope at a star, only to find that the star does not appear any larger than with the unaided eye? This same phenomenon is familiar with even the largest telescopes. As a matter of fact, the only star that has been able to be resolved into a disk on a consistent basis is - the star forming the upper left shoulder of the of Orion (Burnham, 1978). This red is only 520 light away, and with a diameter between 550 to 920 times our (Betelgeuse is a variable star, so its size fluctuates) it can be easily resolved by our largest optical telescopes (figure 7). With this exception, the majority of stars cannot be resolved as a disk.

4 Figure 7: This image of Betelgeuse was taken by the 50cm COAST telescope. It shows the surface of Betelgeuse, which from Earth is only 0.1 arc-seconds across. For comparison, the Pluto is also only 0.1 arc-seconds across. On average, the planet is around 3 arc-seconds across.

With this fact alone, it seems impossible to detect a small planet orbiting a star – after all, if we cannot even resolve the star, how can we possibly detect a planet which is much smaller and does not give off any light? There is one method of imaging that will allow us a view of the end stage of the protoplanetary disk. By masking out the bright central star, it is possible to image the residual disk material called the circumstellar disk (figure 8).

Figure 8: By masking out this central star (star designation HD 141569A), imaging of the circumstellar disk is possible. All the stars in the field are overexposed, so look much larger than they would be under normal imaging circumstances.

Masking of the image is important as the host object is a main-sequence star – that is, a star that has already initiated hydrogen fusing. While much of the protoplanetary nebula can be swept away during the initiation of hydrogen fusion, these images of the circumstellar disks around a main-sequence star tells us that planetary formation is even more likely than with evidence of the protoplanetary disks.

5 While the proplyds and circumstellar disks offer evidence to material that can result in planetary formation, another type of object can also be used to look for evidence of circumstellar disk formation. Herbig-Haro objects (figure 9) are T Tauri stars with an active circumstellar accretion disk (Ostlie and Carroll, 1996). The rotation of this disk is shown by massive lobes of gas that appear perpendicular to the rotating disk. These tell us that the protoplanetary material does in fact rotate about their host star, which can result in planet formation.

Figure 9: These are four Herbig-Haro objects photographed by the Hubble Space Telescope. The green jets are the expulsion of gas from the perpendicular circumstellar disk. These jets are present as a result of circumstellar disk rotation.

Perhaps the most remarkable image of a circumstellar disk is that of Beta Pictoris. This Hubble Space Telescope image (figure 10) is using the masking technique to block the light of the over-exposed star to reveal what looks like several distinct of dense material. While the orbits appear to be very elliptical, it is remarkable that such orbits can be resolved at all.

6

Figure 10: The first image of what appears to be four distinct rings orbiting the star Beta Pictoris. While no planet has been detected in either of these orbits, this supports our theory that as the disk material continues to rotate; they begin to form individual planetesimals. This image may show the early stages of such evolution.

Detecting Exoplanets Radial Velocity - Now that we have identified the presence of planet making debris around other stars, let’s focus on how formed planets are detected around stars. Until direct detections are made, we must rely on indirect methods. The two main techniques in detecting these planets are radial velocity measurements and the transit method. Both are used by professional and amateur astronomers. The first exoplanets were discovered by what is called the “wobble.” This sounds low tech, but this is very significant. When two objects orbiting each other contain any , they will have an affect on each other resulting in an inertial center - called the barycenter (Mayor and Frei, 2003). Using our own solar system as an example, and provide enough mass that the effect is a wobble of the Sun (Marcy and Butler, 1997). The net affect of both planets produce a wobble of around 13 m/s (and if Jupiter were the only

7 planet, a wobble of around 12 m/s would be present). Studies of pulsars and binary stars show that both stars rotate about a common point and not around each other (Mayor and Frei, 2003). The implication is the smaller the companion, the less dramatic this rotation will be. Because of such small variances in stellar wobble, detection is only possible via measurement of the Doppler shift (Butler et al, 1996). Either way, you can think of this effect being analogous to a very heavy person and a very light person on a see-saw with a moveable focus. The Doppler shift is a method used in almost every area of Astronomy, from mapping out the rotation of our to study of the expanding Universe. When an object like a star (or anything that give off or reflects light) is in a stationary position, it can be viewed with a spectrograph to reveal a color spectrum with a specific footprint. Because various chemicals exist, these result in missing spaces within a spectrum called absorption lines. An example would be a star that contains all hydrogen1. Light given off by this star would reveal the familiar color spectrum, but because it contains hydrogen, a certain frequency will be absorbed because of interactions with the hydrogen atom. The electron in the atom absorbs some of the light energy and moves to a higher orbit around the nucleus. This absorbed energy - called a Balmer line when referring to the hydrogen atom - results in a missing portion of the spectrum as shown by example in figure 11 (Freedman and Kaufmann, 2001). Figure 11: This shows an example of a single absorption line. The image at the top shows the object (a galaxy in this case) that is not moving, with the hash marks indicating where that line should be. The middle image shows the object moving away, causing the line to be shifted towards the red – called the redshift. The bottom image shows the object moving towards us causing the line to shift towards the blue – called the blue shift.

1 For analogy purpose only: Stars produce helium through hydrogen fusion, and can contain many other elements is the atmosphere.

8 This shift can be measured. By using simple mathematics, comparing a reference object that contains identical spectra to the shifted object, it is possible to determine the velocity of the object. Equation 1 Equation 1 is the basis of determining the Determining Doppler redshift and blue orbital velocity of the object orbiting the shift: affected star or determining the radial velocity ∆λ λ − λο v z = = = of the affected star. Figure 12 shows how this λ λο c works. Once the orbital velocity is determined, z = redshift simple usage of Kepler’s Third Law will ∆λ = shift in wavelength determine the distance the planet is to the host λ = wavelength of stationary object λο = stationary wavelength – reference star. spectra Kepler’s Third Law: 2 3 v = velocity P = a c = speed of light (300,000 km/s) P = object’s in years It is important to state that: a = object’s distance to star in Astronomical v z = Units c

Figure 12: The unseen orbiting planet creates an inertial center of this system resulting in a wobble that can be detected by the shifting of the stellar spectra. Using equation 1, it is possible to determine the speed of this wobble thereby determining the of the unseen planet. This is also called the radial velocity.

9 The tool used to measure the spectra of an object is called a spectrograph. There are many different flavors of spectrographs, but all of them work using the same principle: separate the visible light into their fundamental wavelengths. Figure 13 shows a very basic spectrograph that is using a special diffracting prism, which is a prism with at least one 60 degree angle. In general, prisms cannot be rotated to adjust resolution.

Figure 13: This is an example of a standard spectrograph. A slit is used to block unwanted light while the lens magnifies the image onto the screen – or our in our case, the CCD sensor.

The sensitivity of the spectrograph is a very important consideration and there are several factors that determine the sensitivity of the spectrograph: how much light can it see, the angle in which the spectra is being observed, quality of design, design of light beam travel within the device, as well as slit dimensions. The most important factor in spectrograph sensitivity is the use of diffraction gratings2 versus a diffracting prism for one reason only: the diffraction grating can be adjusted by the turning of a knob to improve resolution (Tonkin, 2004). Discussing the theory and various flavors of spectrographs is beyond the scope of this paper, so we will focus on the spectrographs used by both amateur and professional planet hunters. I highly recommend Stephen Tonkin’s book Practical Amateur Spectroscopy for additional study. Additional precision of the Doppler shift measurements are possible by extending the focal length of the light within the spectrograph. This is performed by using mirrors to fold the light while collimators (lenses) are used to keep the light focused. An example of this type of spectrograph is the Echelle (figure 14), which happens to be the design used by the California & Carnegie Planet Search group (http://exoplanets.org/).

2 A diffraction grating is a piece of glass that contains hundreds of evenly spaced grooves cut at an angle. These provide the same affect of viewing the spectra as a prism, but can be rotated to increase the viewing angle for improved resolution. The larger surface area also improves resolution over a prism.

10 Figure 14: The Echelle spectrograph extends the focal length of the light by using mirrors and collimators. The collimators keep the light focused. The diffraction grating is at the bottom of the image – labeled Echelle.

The spectrum of a star can contain a large number of absorption lines as a result of the many elements present in the atmospheres of stars. Because of this, the spectra must be extended to allow viewing of all the absorption lines. This is the main reason why the light path within the spectrograph must be extended (Kitchin, 1998). The spacing between the absorption lines allows the Doppler shift to be determined (equation 1). While the Echelle is capable of Doppler precision measurements of around 15 m/s, a method of introducing iodine gas near the slit entrance has allowed for precision measurements of up to 3 m/s (Butler et al, 1996). The iodine is use to create a composite spectra to overlay the analyzed star that enhances our view of the absorption lines while acting like a ruler. By eliminating any uncertainty between stellar absorption lines with a laboratory standard, precision measurements are attained. While the use of iodine has enhanced our abilities to take accurate measurements, the choice of targets also play an important role. The majority of exoplanets discovered have been around metal rich main sequence stars (with a few exceptions) – specifically stars with a spectral class of F, G and H (Butler et al, 2000). The reasons are twofold: 1. F, G, and K type stars are “normal” sized stars like our Sun, and will more than likely exist long enough for planets to form. Larger, hotter burning stars end their lives much sooner so the possibility of mature planets to form is highly unlikely; although planets have been found to orbit stellar remnants such as pulsars.

11 2. Metal rich stars contain heavier elements in their atmospheres as a result of enriched molecular clouds from which they have formed. This results in more absorption lines to be examined. There is a disadvantage to using the Doppler shift to measure radial velocities: the star must be as close to the host star are possible. For example, one of the first exoplanets discovered – the companion to 51 Pegasi – is only 0.05 AU’s3 (Mayor and Frei, 2003). That means the planet, which is 0.5 times the mass of Jupiter, is much closer than the orbit of is about our Sun. For planets that orbit at a larger distance from the star, more precise measurements are desired – mostly because larger orbits require many years to study versus days of a closer orbiting planet. Both radial velocity and astrometric observations by professionals have revealed a number of exoplanets; however, astrometry will not be covered here as this is beyond the current capabilities of the amateur astronomer (for now!). The Transit Method - Another important method of detection is by measuring the transit of a planet over the face of the host star. By performing careful photometric plots of the host star, drops in stellar brightness as a planet moves across the face of the star can be measured (figure 15). The first ever measurement of a stellar transit was made by the Elodie group (discoverers of the companion to 51 Pegasi) and the David Latham group (Mayor and Frei, 2003). The results were shared (as with almost all things in Astronomy) to David Charbonneau and Timothy Brown who are the directors of project STARE (STellar Astrophysics & Research on Exoplanets). What is remarkable about project STARE is their equipment: A 12” Schmidt telescope and a 2K by 2K CCC camera mounted on a Meade LX200 computer controlled mount (http://www.hao.ucar.edu/public/research/stare/stare.html).

3 1 AU is the Earth-Sun distance, or 93 million miles.

12 Figure 15: As a planet moves across the face of a star, the brightness curve of the star drops and can be measured using sensitive CCD camera and computer software.

The transit method can only be used for planet systems that face us edge on (as the orientation of figure 15 indicates). An additional limitation is the size of the planet. As we will see later, a planet 0.64 times the mass of Jupiter decreased the brightness of its host star by only 0.0011 magnitudes. Determining the transit of an is not as difficult as performing radial velocity measurements. While determining the radial velocity requires carefully calibrated equipment, specialized spectroscopes, and lots of patience, determining the transit only requires the skill of photometry and a personal computer. Quite simply, photometry is the study of stellar brightness. Stars with a particular brightness (called by astronomers) have associated relations to size and spectral class. Using online databases and star charts, we can determine the accepted value of brightness for any given star. However, there are a group of stars which can fluctuate in brightness called variable stars. As a star leaves the main sequence and begins to burn up what little hydrogen is left near the core, the outer layers of the star expand and contract. This is very convenient because this “breathing” of the star can be studied using same iodine infused high resolution spectroscopy to determine the speed and duration of this breathing. By gathering a large sample of variable stars that inhabit the Cepheid variable strip (the numerous and most common type of variable star) and evaluating them over a long period of time, it has been concluded that such stars are photometrically stable, and demonstrate peaks of radial pulsation anywhere between 50 to 80 days (Butler, 1998). While variable stars have a particular sequence in their variations in brightness, a handful of normal main sequence stars have been shown to oscillate. Once again, high resolution iodine induced spectrography has revealed a definite pattern of oscillation that can also induce slight

13 variations in stellar brightness (Bedding et al, 2001) when measured photometrically. While this seems to carry serious implications to the accuracy of transit measurements, variations in luminosity as a result of a transit fall in between oscillations – which can occur rapidly, and variations due to the expanding shell of a older star, the stability of which is shown by a much greater delay (table 1). Table 1: Approximate time variations as a result of competing causes. Stellar Oscillations Several times a second Transit of an exoplanet A few days Photometrically stable variable star 50 to 80 days (Butler, 1998)(Bedding et al, 2001)

Radial Velocity – Tau Bootis Star Name Tau Bootis Distance 15.6 4.5 Spectral Class F7 0.28 Planet Mass 4.13 time the mass of Jupiter Orbital Distance 0.05 AU Orbital duration 3.313 days Measured Radial Velocity 15 km/s (Fischer et al, 2001) A group of dedicated amateur astronomers at the Spectrashift.com extrasolar planet search project have successfully measure the radial velocity of a known extrasolar planet around Tau Bootis. To capture the Doppler shift using store bought equipment of such a small target is no easy task, however many of the hurdles that would have blocked any attempt to get accurate readings have been successfully avoided by designing a device specific to the measure of radial velocity. Issues to overcome: 1. Telescope aperture 2. Spectral resolution 3. Signal to noise 4. Light loss 5. Spectroscope stability

14 Telescope aperture: Since project STARE has in its program a 12” telescope, telescope aperture does not play a vital role in photometric evaluations; however, the use of a spectrometer will require a greater aperture as the light entering the spectrograph is already reduced by the entrance slit. The Spectrashift.com group has selected a Meade 16” SCT telescope (figure 16) for spectroscopy; however a 1.1 meter telescope is currently under construction.

Figure 16: While their 1.1 meter telescope is still under construction, the Spectrashift.com team successfully recorded the radial velocity values of Tau Bootis with this Meade 16” SCT.

Spectral Resolution: This can be a major setback for amateurs as the design and stability of the spectrograph play an important role in resolution. This group designed and constructed their own spectrograph, so many of the inherent design limitations from commercial models have been eliminated – mostly because commercial models are designed for the study of a wide variety of spectra and compromise resolution as a result4. None the less, the best resolution of this particular piece of equipment can only measure radial velocities of 200 m/s (Tonkin, 2003). As a result, only large exoplanets can be studies. The spectrograph uses a Czerny-Turner design that uses internal beam folding similar to the Echelle (figure 17).

4 Not all objects in the Universe benefit from accurate measure of radial velocity. For example, such specific designs cannot be used to identify all of the variety of elements in a stars atmosphere.

15 Figure 17: The Czerny-Turner design of the spectrograph allows for good resolution similar to the Echelle (figure 14) while at the same time allows for the swapping of internal parts to adjust resolution.

Signal to Noise: The sources of noise when using a spectrograph are the equipment itself as well as the CCD camera use to capture the images. By using a stable table for the equipment, noise introduced by vibration can be eliminated. In addition, internal light reflections are eliminated by blackening the entire inner structures of the spectrograph. To eliminate noise from the CCD camera, additional cooling is required5. This group used a novel approach and used an office water cooler to feed cooled water through the cooling tubes of the CCD camera.

Light Loss: Although directly related to telescope aperture and slit dimensions, the light loss of concern is in the fabrication of the fiber. Because of the large size of the spectrometer (as well as the desire to be stable), it resides on a table some distance from the telescope, and a fiber optic cable is fed from an eye piece adapter in the telescope to the entrance slit of the spectrograph. To best illustrate the design requirements of the fiber optic cable, I used Adobe Illustrator to show the fiber orientation:

5 Noise is the result of the CCD interpreting heat as a signal. While image reduction can filter out the noise, it’s best to start with an already cooled CCD camera.

16 Using small strands of fiber instead of a solid glass fiber is ideal. The fiber exits the eyepiece in its traditional round format, while the other end is terminated fiber over fiber so all the light can enter the spectrograph slit.

In addition to fiber design, careful polishing of the fiber ends is very important to maximize light transmission. Careful alignment on both ends is possible by using a reference light source as a guide. It will be necessary to have at least two of these fibers: one for the eyepiece and one for the reference light source – in spectroscopy, a reference light source (Argon was chosen by this group) is provided to calibrate the spectrometer (Tonkin, 2003). While fiber optics can be used on just about any telescope, telescopes with a large focal ration are desirable. The longer the focal length of the image, the more narrow the image cone. As a result, the fiber will be able to use all available light (Kitchin, 1998). Spectroscope stability: The stability of the spectroscope is just as important as maximizing available light. Any vibrations induced on the instrument will prevent accurate imaging of spectra, and may prevent imaging altogether as the already small target can shift as a result of cable movement. This issue is the reason telescope mounted spectroscopes – especially in the amateur world – are avoided. A very sturdy work bench, preferably isolated from any walkways – is highly desired to eliminate any induced vibrations.

Looking at the design requirements, the Spectrashift.com group managed to construct a very nice setup using a Meade 16” SCT telescope on a computerized mount, a table- mounted hand made Czerny-Turner spectroscope, an Argon reference lamp, hand-made fiber optic cables, and an Apogee CCD (www.ccd.com) with a 512 x 512 pixel array with

17 24 micron square pixels6. The software of choice is MaxImDL (http://www.cyanogen.com/) for camera control and intermediate image processing, IRAF (http://iraf.noao.edu/) which is the standard for astronomical image analysis, and Microsoft Excel to create the plots. The process of gathering data, image reduction, and analysis is a very time consuming endeavor and will not be covered here. Instead, here is the sequence of events to serve as an overview. 1. Inventory equipment and decide what will be used 2. Test the equipment to ensure working condition 3. Turn on computers, CCD cameras, spectroscopes, reference light sources, and any other pieces of equipment to ensure temperature equilibrium 4. Connect the fiber optic cable from the Argon reference source to the spectroscope 5. Ensure the argon spectrum will overlay the target spectrum – try this with a test star 6. Using computer control, slew the telescope to the desired object (in this case, Tau Bootis) 7. Confirm target in the eyepiece 8. Remove the eyepiece and place the fiber optic cable 9. Begin image capture – 45 minutes per exposure is typical, possible via computer tracking 10. Periodically capture other stars in the same field of view for reference 11. Once enough images are captured – the more the better – image calibration can commence. Please see the attached Appendix: Image Reduction Step by Step 12. Use of the IRAF software is to be used at this point (http://iraf.noao.edu/) which will create any data points required 13. Data points are entered into an Excel spreadsheet. Scatter plots are preferred. The results of all this hard work is an Excel scatter plot that graphs out the radial velocities of the orbiting object (figure 18). The positive numbers on the y-axis indicate radial velocity toward us, and the negative numbers on the y-axis indicate radial velocity away from us. The x-axis is time. Notice the : measured 3.41 days here which is in good agreement with the published results. Once the data has been reduced and analyzed, it is always a good idea to have an independent group evaluate the data. This data was sent to NOAO (National Optical Astronomy Observatory: http://www.noao.edu/) for analysis, with the results equal to the Spectrashift.com group (figure 19).

6 The size of the pixels in analogous to film speed; 24 microns is very sensitive, but a lower number means higher resolution. In this case, we want the most light possible so a more sensitive CCD is required.

18

Figure 19: This graph, the result of independent Figure 18: The result of numerous images of spectra analysis from NOAO, shows almost identical from Tau Bootis. results.

Final Test: The final test as to the accuracy attained by an amateur group is to compare the results with the published data (figure 20): The California & Carnegie Planet Search The Spectrashift.com Team: Team:

Figure 20: Orbital Period: 3.312 days; Orbital Period: 3.41 days; Amplitude: 420 Amplitude: 471 km/s km/s

While these results are not exact, it shows that a dedicated group of amateurs can yield results very similar to the professionals. So why are the numbers not exact? According to the Spectrashift.com website, the initial series of data was out of phase 180%, and the mathematics was not tested accurate 99 times out of 100. Basically the errors were in the image processing and not in the techniques used to capture the data.

19 Transit Method – HD209458 Star Name HD 209458 or SAO 107623 Distance 47 parsecs Apparent Magnitude 7.65 Spectral Class G0 Metallicity 0.04 Planet Mass 0.62 time the mass of Jupiter Orbital Distance 0.046 AU Orbital duration 3.5239 days Differential Magnitude7 0.0011 (Henry et al, 2000) Greg Laughlin and Tim Castellano – founders of http://transitsearch.org - have demonstrated that photometry to measure a stellar transit can be obtained with a 8” telescope, an entry-level CCD camera, and over the counter Astronomy software (not to mention of course clear skies). Specifically, an Meade 8” LX200 telescope (www.meade.com) armed with a Santa Barbara Instrument Group (www.sbig.com) ST-7 CCD camera and CCDSoft (www.bisque.com) software will yield very impressive results (figure 21).

Figure 21: This Meade 8” telescope and SBIG ST- 7 camera can be purchased for around $7500.00. With it, it is possible to obtain professional quality light curves of a transiting exoplanet.

CCD cameras are very sensitive to changes in brightness of a star. With a properly reduced image, changes in luminosity of as little as 0.011 magnitudes are possible. The procedure of gathering photometric data is simple, but requires patience and skill with a telescope and CCD imaging software:

7 The differential magnitude is the level of magnitude change when the planet is at maximum transit

20 1. Gather a target list. In this case, the target is HD 209458 2. The Meade LX200 is computer controlled. Software Bisque makes a wonderful software package called TheSky – which is a computer planetarium and offers telescope control. Click on the desired star, and tell the software to move the telescope in position. 3. Use CCDSoft – the included CCD camera control software when purchasing an SBIG camera– to begin capturing a series of 30 second images with the ST-7. 4. Perform image reduction of all the images within CCDSoft. 5. Use the imbedded photometry tools within CCDSoft to gather photometric points of the target star, as well as a few other stars that are in the same field. 6. Input the target points into an Excel spreadsheet – a scatter plot is preferred. Since the interest is to evaluate changes in brightness, photometric calibration of accurate stellar magnitudes is not required. However, it is a good idea to make sure the surrounding stars do not exhibit the same changes in brightness. When compared to a plot generated by professional equipment (figure 23), it is clear the amateur has much to offer (figure 22).

Figure 23: This plot, released in the Astrophysical Figure 22: This data plot is the result of a Meade Journal Letters by a professional Astronomer shows 8” telescope, an SBIG ST-7 CCD camera, and it’s an identical photometric curve of star HD 209458-b provided CCD imaging software – CCDSoft. The (Charbonneau et al, 1999). scatter plot was created in Microsoft Excel.

Both professional and amateur plots reveal the orbital period of HD 209458 to be 3.52 days.

21 Other methods

The radial velocity and transit methods of planetary detection are the most common used techniques in the search of exoplanets, and the only methods used by amateurs; however, there are many more techniques available to the professional, as well as future space missions designed by NASA and the ESA for the sole purpose of improving our resolution capabilities. • Astrometry – this method is used for long term accurate measure of the star apparent motion in the sky, and is used to detect planets greater than around 3 AU’s from the host star. This method does not seem very popular and is passed over in favor of the accurate measurements of the host stars radial velocity. • Microlensing – this technique is used to in attempt to locate Dark Matter and black holes, but has been very successful in looking for orbiting planets. During a microlensing effect on a star, a brief but noticeable deviation of the light curve, as shown by figure 24 can be imaged. • Optical and Infrared Interferometry – an interferometer is a device that is used to combine the wave sources from two or more instruments and combing them to produce an image of much higher resolution. The Keck 1 and 2 in Hawaii, and the VLT Interferometer (VLTI) in Chile are the two current interferometers used to help detect exoplanets. The VLTI is in operation, but is also a work in progress. The resolution capabilities of this system hope to reach 10 micro-arcseconds (Mayor and Frei, 2003). Figure 24: This photometric microlens record shows an Earth sized planet orbiting pulsar PSR B1257+12. The light source for this lens is a distant galactic bulge, and the deviation is a result of the orbiting planet being on either the front or backside of the pular resulting in a net increase of mass thereby producing a more powerful lens. A lens effect is the result of the deviation of light as a result of a massive object placed in between the source of light and the observer.

22 Here is a nice summary of detection methods – both present methods and proposed methods:

A Brief Window into the Future Telescopes in space offer tremendous benefits: there is no atmosphere to affect the quality of the images, and the already low temperature will ensure better noise control and infrared images require as low as possible. Three major space-based projects are in the design stage: NASA’s Kepler project, Terrestrial Planet Finder – or TPF and ESO’s DARWIN project (also known as the IRSI). The Kepler (figure 25) mission will utilize a very sensitive photometer to examine the transits of Earth-sized planets. With a planned launch date of 2007, 100,000 stars will be evaluated with the goal of obtaining a list of targets for the following Terrestrial Planet Finder mission.

Figure 25: The Kepler is still in the design stage. With a proposed launch date around 2007, its 37” photometry lens will study 100,000 stars to look for transits of Earth-sized planets.

23 The TPF will use two space-based, infrared sensitive telescopes in concert to create an infrared interferometer. The use of spectroscopy in the infrared will allow the study of cooler objects that orbit the stars instead of the star itself. The sensitivity of the TPF has a goal to view Earth-like atmospheres around Earth sized planets that orbit within the habitable zone. While this zone will be difficult to determine due to differences in and temperature, the idea is this zone is approximately the Earth-Sun distance. While obtaining spectroscopic data on the various minor gaseous elements in an atmosphere will prove to be difficult, the goal is to at least obtain spectra of water, ozone, and carbon dioxide. Such elements in an atmosphere would be considered Earth-like (Mayor and Frei, 2003). The has their own project to look for exoplanets as well, called DARWIN. The goal of DARWIN is very much the same as the TPF: to look for Earth-like atmospheres around planets orbiting within the habitable zone. DARWIN has a target launch date around 2015, so we have some waiting to do. In the meantime, have a look at the proposed design of the DAWRIN interferometer (figure 26):

Figure 26: The DARWIN spaced-based interferometer will hopefully launch by around 2015. Its mission is to study the characteristics of nearby exoplanets, search for Earth-like atmospheres, and perform some “general” Astronomy. NASA’s TPF will use a very similar design.

More information on Kepler can be found here: http://discovery.nasa.gov/kepler.html

24 More information on the TPF can be found here: http://planetquest.jpl.nasa.gov/TPF/tpf_index.html More information on DARWIN can be found here: http://ast.star.rl.ac.uk/darwin/ Conclusion The study of exoplanets is ongoing. Continued advances in professional astronomy are allowing for increased sensitivity. The two most common and successful tools used by both amateurs and professionals are the measure of transit brightness and radial velocity. By collaborating with amateur astronomers, professional telescope time is preserved. The gathering of orbital data is vital to ensure repeatability with software analysis, and the amateur is poised to provide this important data. Most of all, we have shown that these two methods of planetary detection is possible, and that amateurs can also join the hunt. Greg Laughlin and Tim Castellano of http://transitsearch.org are actively recruiting amateur astronomers that wish to participate. This website contains wonderful information on how to put together a telescope ready to capture transit data, and also coordinates target stars with participating members to avoid any overlap or missing data. With careful planning, it is possible to duplicate the methods used by the Spectrashift.com group to gather data on radial velocity of exoplanets. While no professional group is seeking amateurs to provide spectroscopic data, this is sure to change as the contribution of the amateur have proven valuable for those involved in transit searches. The hunt is on……… Recommended Internet Resources: California & Carnegie Planet Search: http://exoplanets.org Anglo-Australian Planet Search: http://www.aao.gov.au/local/www/cgt/planet/aat.html The European Southern Observatory VLTI: http://www.eso.org/projects/vlti/ NASA Origins of Solar Systems amateur project: http://origins.jpl.nasa.gov/index1.html The Geneva Extrasolar Planet Search: http://obswww.unige.ch/~udry/planet/planet.html NASA Terrestrial Planet Finder: http://planetquest.jpl.nasa.gov/TPF/tpf_index.html Advanced Fiber Optic Echelle Program: http://cfa-www.harvard.edu/afoe/index.html Project STARE: http://www.hao.ucar.edu/public/research/stare/overview.html Planet Homepage - Microlensing: http://planet.iap.fr/

25 References: Beatty, J. Kelly. Carolyn C. Petersen and Andrew Chaikin. The New Solar System 4th Edition. Cambridge University Press, 1999. Bedding, Thomas et al. “Evidence for Solar-Like Oscillations in ß Hydri.” The Astrophysical Journal, 549: L105-L108, March 1, 2001. Burnham, Robert Jr. Burnham’s Celestial Handbook – Volume Two. Dover Publications, Inc., New York, 1978. Butler, Paul. “A precision Velocity Study of Photometrically Stable Stars in the Cepheid Instability Strip.” The Astrophysical Journal, 494: 342-365, February 10, 1998. Butler, Paul, et al. “Attaining Doppler Precision of 3 m s -1.” Publications of the Astronomical Society of the Pacific, v 108: 500-509, June 1996. Butler, Paul et al. “Planetary Companions to the Metal-Rich Stars BD -10o 3166 and HD 52265.” The Astrophysical Journal, 545: 504 – 511, December 10, 2000. Charbonneau, David et al. “Detection of Planetary Transits Across a Sun-like Star.” Astrophysical Journal Letters, 23 November, 1999. Fischer, Debra et al. “Planetary Companions to HD 12661, HD92788, and HD 38529 And Variations in Keplerian Residuals of Extrasolar Planets.” The Astrophysical Journal, 551: 1107 – 1118, April 20, 2001. Freedman, Rodger and William Kaufmann. Universe 6th Edition. W.H. Freeman and Company, New York 2001. Henry, Gregory et al. “A Transiting “51 Peg-Like” Planet.” The Astrophysical Journal, 529: L41 – L44, January 20, 2000. Kitchin, C. R. Astrophysical Techniques 3rd Edition. Institute of Physics Publishing. Bristol, 1998. Marcy, G. W. and Paul Butler. “Characteristics of Observed Extrasolar Planets.” The Tenth Cambridge Workshop on Cool Stars, Stellar Systems and the Sun. Cambridge, Massachusetts. July 16-20, 1997. Mayor, Michael and Pierre-Yves Frei. New Worlds in the Cosmos. The Discovery of Exoplanets. Cambridge University Press, 2003. Ostlie, Dale. Bradley Carroll. An Introduction to Modern Stellar Astrophysics. Addison- Wesley Publishing Company, Inc. Reading, Massachusetts, 1996.

26 Tonkin, Stephen. Practical Amateur Spectroscopy. Springer. London, 2003.

Image Credits: Figure 1: http://www.star.ucl.ac.uk/~rhdt/diploma/lecture_1/contraction.jpg Figure 2: http://www.rc-astro.com/nebulae/m42_2004-01-27.htm Figure 3: http://hubblesite.org/newscenter/newsdesk/archive/releases/1995/49/image/b Figure 4: http://hubblesite.org/newscenter/newsdesk/archive/releases/1995/45/image/b Figure 5: http://hubblesite.org/newscenter/newsdesk/archive/releases/1995/45/image/b Figure 6: http://hubblesite.org/newscenter/newsdesk/archive/releases/1995/45/image/b Figure 7: http://www.mrao.cam.ac.uk/telescopes/coast/betel.html Figure 8: http://hubblesite.org/newscenter/newsdesk/archive/releases/2003/02/image/b Figure 9: http://hubblesite.org/newscenter/newsdesk/archive/releases/1999/05/image/c Figure 10: http://hubblesite.org/newscenter/newsdesk/archive/releases/2000/02/ Figure 11: http://www.gb.nrao.edu/~rmaddale/Education/Wvsta'98/200c.gif Figure 12: http://cfa-www.harvard.edu/afoe/doppler-shift.gif Figure 13: http://www.astrophys-assist.com/educate/solarobs/ses01p16.htm Figure 14: http://msowww.anu.edu.au/observing/74in/Echelle/ech_get_go_quick.html Figure 15: http://www.hao.ucar.edu/public/research/stare/overview.html Figure 16: http://www.spectrashift.com/meade.jpg Figure 17: http://www.spectrashift.com/spectro.html Figure 18: http://www.spectrashift.com/tauboo.html Figure 19: http://www.spectrashift.com/tauboo.html Figure 20: http://exoplanets.org Figure 21: http://transitsearch.org Figure 22: http://transitsearch.org Figure 23: Charbonneau et al, 1999. Figure 24: http://www.nd.edu/~srhie/MPS/ Figure 25: http://discovery.nasa.gov/kepler.html Figure 26: http://ast.star.rl.ac.uk/darwin/pics/alcatel_ff_jul99.jpg

27 The Appendices – Candidate Exoplanets as of October 4, 2004. Source: The Extra-Solar Planet Encyclopedia (http://www.obspm.fr/encycl/encycl.html). Appendix 1A: Candidate Planets around Main Sequence Stars

NAME M[.SINI] SEM-MAJ. PERIOD ECC. INCL. Update click Jup. mass AXIS (AU) days (deg) for more (note 1) OGLE-TR-56 b 1.45 0.0225 1.2 - 81.0 07/10/04 * OGLE-TR-113 b 1.35 0.0228 1.43 - - 18/08/04* OGLE-TR-132 b 1.01 0.0306 1.69 - - 18/08/04* HD 73256 b 1.85 0.037 2.54863 0.038 - 14/04/03* GJ 436 b 0.067 0.0278 2.6441 0.12 - 31/08/04* 55 Cnc e 0.045 0.038 2.81 0.174 25? 01/10/04 * b 0.84 0.11 14.65 0.02 25? - * c 0.21? 0.24? 44.28? 0.34? 25? - * d 4.05 5.9 5360 0.16 25? - * HD 83443 b 0.41 0.04 2.985 0.08 - 03/10/04 * HD 46375 b 0.249 0.041 3.024 0.04 - 15/08/04 * TrES-1 0.75 0.0393 3.030065 0.0 - 30/09/04 * HD 179949 b 0.84 0.045 3.093 0.05 - 31/08/04 * HD 187123 b 0.52 0.042 3.097 0.03 - 24/09/02 * Tau Boo b 3.87 0.0462 3.3128 0.018 - 21/02/04 * HD 330075 b 0.76 0.043 3.369 0.0 - 15/08/04 * HD 88133 b 0.29 0.046 3.415 0.11 - 07/09/04 BD-10_3166 b 0.48 0.046 3.487 0. <84.3 27/04/00 * HD 75289 b 0.42 0.046 3.51 0.054 - 20/10/03 * HD 209458 b 0.69 0.045 3.524738 0.0 86.1 17/08/04 * HD 76700 b 0.197 0.049 3.971 0.0 - 16/07/02 * OGLE-TR-111 b 0.53 0.047 4.02 - - 17/09/04 * 51 Peg b 0.468 0.052 4.23077 0.0 - 29/02/04 * Ups And b 0.69 0.059 4.6170 0.012 - 17/09/04 * c 1.19 0.829 241.5 0.28 - - * d 3.75 2.53 1284. 0.27 - - * HD 49674 b 0.12 0.0568 4.948 0. - 13/06/02 * HD 68988 b 1.90 0.071 6.276 0.14 - 07/02/03 * HD 168746 b 0.23 0.065 6.403 0.081 - 26/02/02 * HD 217107 b 1.28 0.07 7.11 0.14 - 31/08/04 *

28 HD 162020 b 13.75 0.072 8.428198 0.277 - 26/02/02 * HD 160691 d 0.042 0.09 9.55 0. - 16/09/04 * b 1.7 1.5 638 0.31 - - * c 3.1 4.17 2986 0.8? - - * HD 130322 b 1.08 0.088 10.724 0.048 - 30/03/00 * HD 108147 b 0.41 0.104 10.901 0.498 - 26/02/02 * HD 38529 b 0.78 0.129 14.309 0.29 - 01/10/04 * c 12.70 3.68 2174.3 0.36 - - * Gl 86 b 4. 0.11 15.78 0.046 - 17/09/04 * HD 195019 b 3.43 0.14 18.3 0.05 - 20/04/03 * HD 102117 b 0.18 ~0.15 20.8 0.08 - 16/09/04 HD 6434 b 0.48 0.15 22.09 0.30 - 07/02/03* HD 192263 b 0.72 0.15 24.348 0.0 - 23/05/03 HD 117618 b 0.16 ~0.15 25.8 0.31 - 16/09/04 Gliese 876 c 0.56 0.13 30.1 0.12 ? 08/04/04 * b 1.98 0.21 61.02 0.27 84 - * rho CrB b 1.04 0.22 39.845 0.04 - 21/02/04 * HD 74156 b 1.86 0.294 51.643 0.636 - 01/10/04* c >6.17 3.40 2025 0.583 - - HD 37605 b 2.85 0.26 55.2 0.736 - 08/07/04 * HD 168443 b 7.7 0.29 58.116 0.529 - 30/04/04 * c 16.9 2.85 1739.50 0.228 - - * HD 3651 b 0.2 0.284 62.23 0.63 - 16/01/03 * HD 121504 b 0.89 0.32 64.6 0.13 - 21/08/00* HD 178911 B b 6.292 0.32 71.487 0.1243 - 20/04/03 * HD 16141 b 0.23 0.35 75.560 0.28 - 10/07/03 * HD 114762 b 11. 0.3 84.03 0.334 low ? 20/04/03 * HD 80606 b 3.41 0.439 111.78 0.927 - 20/04/03 * 70 Vir b 7.44 0.48 116.689 0.4 - 21/02/04 * HD 216770 b 0.65 0.46 118.45 0.37 - 02/09/03 * HD 52265 b 1.13 0.49 118.96 0.29 - 31/08/04 * HD 208487 b 0.41 ~0.5 129 0.47 - 16/09/04 GJ 3021 b 3.21 0.49 133.82 0.505 - 31/08/04 * HD 37124 b 0.75 0.54 152.4 0.10 - 01/10/04 * c 1.2 2.5 1495 0.69 - - * HD 219449 b 2.9 ~ 0.3 182 - - 10/10/03 HD 73526 b 3.0 0.66 190.5 0.34 - 05/07/02 *

29 HD 104985 b 6.3 0.78 198.2 0.03 - 26/06/03 * HD 82943 b 0.88 0.73 221.6 0.54 - 31/08/04* c 1.63 1.16 444.6 0.41 - -* HD 169830 b 2.88 0.81 225.62 0.31 - 30/04/04* c 4.04 3.60 2102 0.33 - -* HD 8574 b 2.23 0.76 228.8 0.40 - 01/09/03 HD 89744 b 7.99 0.89 256.6 0.67 - 19/04/04 * HD 134987 b 1.58 0.78 260 0.25 - 31/08/04 * HD 40979 b 3.32 0.811 267.2 0.25 - 25/05/03 * HD 12661 b 2.30 0.83 263.6 0.096 - 01/09/03 * c 1.57 2.56 1444.5 <0.1 - - * HD 150706 b 1.0 0.82 264.9 0.38 - 19/06/02 HD 59686 b 6.5 ~ 0.8 303 - - 30/10/03 HR 810 b 2.26 0.925 320.1 0.161 - 11/07/03 * HD 142 b 1.36 0.980 338.0 0.37 - 31/08/04 * HD 92788 b 3.8 0.94 340 0.36 - 24/09/02 * HD 28185 b 5.6 1.0 385 0.06 - 29/10/01 * HD 142415 b 1.62 1.05 386.3 0.5 - 02/09/03* HD 177830 b 1.28 1.00 391 0.43 - 22/11/00 * HD 154857 b 1.80 1.11 398 0.51 - 16/09/04 * HD 108874 b 1.65 1.07 401 0.20 - 15/06/02 * HD 4203 b 1.65 1.09 400.944 0.46 - 10/07/03 * HD 128311 b 2.63 1.06 414 0.21 - 13/06/02 * HD 27442 b 1.28 1.18 423.841 0.07 - 11/07/03 * HD 210277 b 1.28 1.097 437 0.45 - 14/06/01 * HD 19994 b 2.0 1.3 454 0.2 - 20/04/03* HD 20367 b 1.07 1.25 500 0.23 - 19/06/02 HD 114783 b 0.9 1.20 501.0 0.1 - 15/08/04 * HD 147513 b 1. 1.26 540.4 0.52 - 27/01/04* HIP 75458 b 8.64 1.34 550.651 0.71 - 29/03/02 * HD 222582 b 5.11 1.35 572.0 0.76 - 07/05/03 * HD 65216 b 1.21 1.37 613.1 0.41 - 02/09/03 * HD 141937 b 9.7 1.52 653.22 0.41 - 26/02/02 * HD 41004A b 2.3 1.31 655 0.39 - 15/08/04 HD 47536 b 4.96 - 9.67 1.61 - 2.25 712.13 0.20 - 18/12/02 * HD 23079 b 2.61 1.65 738.459 0.10 - 31/08/04 *

30 16 CygB b 1.69 1.67 798.938 0.67 - 21/02/04 * HD 4208 b 0.80 1.67 812.197 0.05 - 31/08/04 * HD 114386 b 0.99 1.62 872 0.28 - 19/06/02* gamma Cephei b 1.59 2.03 902.96 0.2 - 17/09/04 HD 213240 b 4.5 2.03 951 0.45 - 29/10/01 * HD 10647 b 0.91 2.10 1040 0.18 - 02/09/03 HD 10697 b 6.12 2.13 1077.906 0.11 - 10/07/03 * 47 Uma b 2.41 2.10 1095 0.096 - 01/04/04 * c 0.76 3.73 2594 <0.1 - - * HD 190228 b 4.99 2.31 1127 0.43 - 24/09/02 HD 114729 b 0.82 2.08 1131.478 0.31 - 07/05/03 * HD 111232 b 6.8 1.97 1143 0.20 - 02/09/03 * HD 2039 b 4.85 2.19 1192.582 0.68 - 10/07/03 * HD 136118 b 11.9 2.335 1209.6 0.366 - 11/02/02 * HD 50554 b 4.9 2.38 1279.0 0.42 - 01/09/03 * HD 196050 b 3.0 2.5 1289 0.28 - 07/10/02 * HD 216437 b 2.1 2.7 1294 0.34 - 27/01/04 * HD 216435 b 1.49 2.7 1442.919 0.34 - 22/05/03 * HD 106252 b 6.81 2.61 1500 0.54 - 01/09/03 * HD 23596 b 7.19 2.72 1558 0.314 - 01/09/03 14 Her b 4.74 2.80 1796.4 0.338 - 08/03/04 * OGLE-235/MOA-53 1.5-2.5 2.8 - 3 ? ? - 15/04/04 * HD 39091 b 10.35 3.29 2063.818 0.62 - 27/01/04 * HD 72659 b 2.55 3.24 2185 0.18 - 08/03/04 * HD 70642 b 2.0 3.3 2231 0.1 - 27/01/04 * HD 33636 b 9.28 3.56 2447.292 0.53 - 31/08/04 * Epsilon Eridanib 0.86 3.3 2502.1 0.608 46? 17/09/04 * c 0.1?? 40?? 280 y?? 0.3?? ?? - HD 30177 b 9.17 3.86 2819.654 0.30 - 11/07/03 * Gl 777A b 1.33 4.8 2902 0.48 - 08/03/04* b ~ 5 > 55 > 2450 yr - - 10/09/04*

An * means that the planet is published (or accepted) in a refereed journal Appendix 1B: Pulsar Planets

Name M[.sini] Semi-maj. Period Ecc. Inclin. Update

31 Jup. mass:(J) Axis (AU) years (y) (deg) Earth mass:(E) days (d) (note 1) PSR 1257+12 0.020 (E) 0.19 25.262(d) 0.0 - 11/09/03 * 4.3 (E) 0.36 66.5419(d) 0.0186 53 or 127 - * 3.9 (E) 0.46 98.2114(d) 0.0252 47 or 133 - * cometary mass ? ? ~ 3 (y) ? ? - - ~ 100 (E) 40 ~ 170 (y) - - - * PSR B1620-26 2.5 (J) 23 ~ 100 (y) - 55 01/10/04 *

[*: Planet published or accepted in a refereed journal]. Appendix 1C: Cluster and Free Floating Planets

NAME MASS TEMPERATURE (K) RADIUS DISTANCE Update click Jup. mass RJ pc [ *: published] for more 3 1,100 1.6 440 14/11/03

An * means that the planet is published (or accepted) in a refereed journal

Appendix 1D: Unconfirmed, controversial and/or retracted planets

Name M[.sini] Semi-maj. Period Ecc. Inclin. Update Jup. Axis (AU) years(y) (deg) mass:(J) days(d) Earth mass:(E) (note 1) OGLE-TR 3 0.5 ? 0.025 1.1899 - 90 ? 12/11/03 GM Aurigae 1.7 ? 2.5 ? 4 yrs ? - 54 ? 12/02/03 HD 83443 0.17 0.174 29.83 0.42 - 26/11/02 HD 13507 3.46 2.39 1318 0.13 - 09/10/02 HD 223084 1.21 (J) 0.44 104.090 0.48 - 01/10/02 HR7875 0.69(J) - 42.5 (d) 0.429 - 15/07/98 Q0957+561 ~ a few (E) - - - - - Lal. 21185 0.9 (J) - 5.8 (y) - - 05/12/97 1.6 (J) - 30 ? (y) - - - CM Dra - - - - - 25/03/00 Alpha Tau 11 (J) 1.3 - 1.4 653.8 (d) 0.182 #177 - 04/12/97

32 0.065 TMR-1C 2- 3 (J) ? > 1,500 > 40,000 ? ? 03/04/01 (y) OGLE-2002- ~ 1 - 10 (J) ? ? ? ? 03/12/02 BLG-055 98-BLG-35 ~ 0.4 - 1.5 1.5 or 2.3 ? ? ? 07/03/01 (E) 97-BLG-41 ~ 3 (J) ~ 7 ? ? ? 05/04/02 95-BLG-3 ~ 2 (J) > 5 - 10 ? - - - 23/12/96 94-BLG-4 ~ 5 (J) ~ 1 - - - 23/12/96 Beta Pic - > 6 2000 (d) - < 1 o 02/09/99 BD +31 o 643 disk rad. disk: - - < 10 o 10/03/97 6,600 W 3 (OH) 10-4 ? (E) 2,000? - - - 28/10/99

Appendix 2

The location of Tau Bootis: Image from Starry Night 4.5

33 Appendix 3

The location of HD 209458: Image from http://transitsearch.org

Appendix 4: Image Reduction – Step by Step:

This appendix discusses the method of image reduction of a cluster of star – a previous project of mine. The methods outlined are standard in any type of imaging in Astronomy.

34 Our software of choice will be MaxImDL. This software is capable of controlling a CCD camera, and performs some very powerful image edit using a wide variety of tools. With this software, we are able to remove any bad or damaged pixels, calibrate images using the bias and flat frames, create a color composite, and perform photometry – which is the only method used to create the data points on a CMD. In addition to MaxImDL, we will also use Microsoft Excel to store the photometric data as well as perform image calibration and create the actual diagram. The image reduction process itself is a fairly simple concept, but time consuming. While the attached appendix documents every step of the reduction process, a guided tour is provided: Each image has, in addition to the actual image (or flat or bias), an area of unexposed pixels used to store bias information. This area is called ‘overscan.’ Because we are calibrating various images using various filters, the overscan area is of no use and must be removed. There is a thin column of nothing on the far right of the image on the left. This thin area is the overscan and is present in every image provided by the McDonald Observatory. This area can be mapped out within MaxImDL and applied to every image. In addition to the overscan, the border of the entire image, which is one pixel in width, must also be removed. This has been included in the map.

35 The remove bad pixels tool under the process menu is able to remember the selected pixels in one image. The unfortunate is that every pixel has to be selected by the mouse, or entered by hand – if you know the exact pixel location. In this case, I renamed my map1 to Remove Overscan – SAO project. To apply this map to other images, I open this tool and click the process button.

The result of the overscan removal is seen here. Notice the dark area on the far right is no longer present. This image of M67 is oriented properly versus the mirrored image above. Once the overscan has been removed from all cluster images, standard fields, bias and flat images, the headers of each image must be fixed.

36

Looking at the view menu within MaxImDL, there is a tool to view the fits header. Every image captured on professional CCD cameras store important information about the image captured. For example, the size of the chip, date of image capture, duration of image capture, the airmass value, and filter used is recorded in this header. Also, demographic, telescope type, and astrometry information can also be a part of this header. Within this header, the areas of overscan have been removed, and the image size has been updated to reflect the actual image size (details are available in the image reduction appendix). In the image above, the value ‘airmass’ is highlighted. This value is an important one as it is used to determine the airmass value of the atmosphere. Every night, the quality of the atmosphere (called the “seeing”) is given a numerical value. This value also changes as the area photographed is closer to the horizon. The airmass values included in the fits header will probably be incorrect. The good news is there are airmass calculators available so the correct airmass value can be found. The bias and flat field images do not require a correct airmass value as no object is being photographed through the atmosphere. Once the overscan has been removed and the fits headers correct, image calibration can begin. As above, there is a detailed report of image reduction in the image reduction diary; however, a brief tour will be give. Bias and flat images are more effective of there are more of them. These files are combined into a master file. The first file to create is the bias.

37 As you can see, there is nothing fancy about a bias image. The sole purpose is to create a calibration as to the levels of brightness for every image the bias has been applied. On close inspection, there are very small white dots on this image. These are not a part of the bias measurement and will be averaged out when combined with other bias frames.

Combining files is a two step process. The first step is to locate the images you wish to combine using the combine tool under the file menu. The more images the better, but in our case we have 5 (only for are selected here, but all five images were selected to create our bias frame.

38 When the combine button is clicked, the second step of the process is available. Auto – star matching is selected automatically, but will be ignored in this case since there are no stars to match. It is important to select median under the output option. This averages the information in all the images and combines them to a single master image. This is the process that removed the white dot artifacts on the individual bias images. This process of combining files to a master file will work for both the bias and the flat images. Due to the purpose of the images, the science images are not to be combined!

The result of the combine is an equally boring, but important image. In addition to combing the bias frames, the flat frames will also be combined. Since we are using BRI filters for our analysis, flat frames for filters BRI are required; however, before any flat frames are combined, the bias image must be applied to the flat fields. The calibration tool, under the process menu has two variants: set calibration or calibration wizard.

39 Version 4 of MaxImDL has a very capable calibration wizard that I highly suggest. The wizard walks you through the entire process. The image on the left shows the calibration setting for the science images, but the calibration tool can also use only the bias image to calibrate the flat fields. Once the master bias file is selected, simply open all of the flat frames and select calibrate all under the process menu. This method will also be used for the science images.

The image on the left shows a single flat image (from the I flat field) and the image on the right is a median combine of 5 flat images. This image will be used to calibrate all of the science images using the I filter. The red and the blue flat images look similar to the images above, and will not be demonstrated to save space. Using the same calibration above, apply the calibration to the appropriate images. Every image will use the master bias frame, but the filter specific images will require the master file counterpart – i.e. images through the B filter must be calibrated with the master bias and the master B flat.

40

The result is a nicely calibrated science image with no gross defects. Now that we have calibrated images, we can perform photometry on all the images. For the purpose of calibration to the Landolt system, we will use images of globular cluster NGC4147, Landolt Standard Area 104 (SA104), and Landolt Standard Area 107 (SA107). As mentioned above, airmass plays a role in a telescopes ability to ‘see’ a star. This affect can also interfere with the color term. As a result, the fields NGC4147, SA104 and SA107 will be images at various times through the night, so the position of these areas will cause a difference in airmass value. These values, shown later, will affect the outcome of the color term. The apparent magnitude of selected stars in each of the

41 calibration fields will need to be documented. The best method is to use the photometry tool within MaxImDL.

The photometry and information windows work in tandem. In order to create a photometry plot, a reference star is to be selected. A single reference star from the Landolt standards can be selected, and the magnitude of the reference star is required for the Ref Mag field. The remaining stars are selected as objects. When the plot is viewed, the results can be saved as a CSV (comma separated values) file. To the left is a portion of an image being analyzed for photometry. The Ref1 is the reference star, and the Obj1, Obj2…... are the object stars. Once the photometry plots are saved, the data can be entered in an Excel spreadsheet. The attached spreadsheet contains all of the instrument magnitudes of each filter (the photometric data) as well as the Landolt standards.

42