# Detecting planets via transitss

You are familiar with the basic idea behind the transit method, right?

The idea is pretty simple: if a star has a planetary companion, AND if a distant observer lies in the orbital plane of that planet, then ever now and then the planet will pass in front of the star. As the planet makes its transit, it will block some of the star's light, causing a brief, small dip in its brightness. Find a series of little dips, and *boom* there's the planet.

So, what's so hard about that?

```

Q:  Really -- why is this so difficult?

```

Well, let's begin by seeing just how large these dips in brightness might be, and how long they last. Let's take as an example the solar system, and two of the most interesting of the planets:

```

(m)                     (m/s)
-----------------------------------------------------------------

Sun                   6.96 x 108               ------

Earth                 6.37 x 106               29,800

Jupiter               7.15 x 107               13,100

------------------------------------------------------------------

Q:   What fraction of the Sun's disk will Earth block?
How long will it take Earth to travel across the Sun?

Q:   What fraction of the Sun's disk will Jupiter block?
How long will it take Jupiter to travel across the Sun?

```

You should find results somewhat like this:

Yes, these dips in brightness are pretty darn small: Jupiter blocks about 1% of the Sun's disk, while the Earth blocks about 0.01%. For this reason, astronomers use the following terms when discussing transits:

• millimags: one millimag is one-thousandth of a magnitude, or 0.001 magnitudes. Jupiter's transit is about 10 millimags deep, while Earth's is about 0.1 millimag deep.
• micromags: one micromag is one-millionth of a magnitude, or 0.000 001 magnitudes. People who look for transits by terrestrial planets sometimes use this unit: Earth produces a dip about 100 micromags deep.

So, if we want to find planets via the transit method, we need high precision photometry.

We also need patience. Consider a poor ground-based astronomer who is looking for new planets with this technique.

```

Q:  How long does a transit from an Earth-like planet last?

Q:  How often will a transit occur?

Q:  What fraction of the time will one particular system
be in transit?

(Compare this to the fraction of time one particular

Q:  How many transits must one observe in order to be sure
that the dip is due to a planet, and not to some
other effect, or due to noise?

```

For example, suppose that you and your friend decide to study a star which may host an exoplanet.

• You decide to use the transit technique

Each person makes one measurement per night. Let's see what each of you can learn as time passes.

Day 1 - 80 Click on the graph to see what happens.

```

Q:  After 80 days, what do the radial velocity measurements
reveal to us?

Q:  After 80 days, what do the transit measurements reveal to us?

```

Day 81 - 130 Click on the graph to see what happens.

```

Q:  After 130 days, what do the radial velocity measurements
reveal to us?

Q:  After 130 days, what do the transit measurements reveal to us?

```

Day 131 - 400 Click on the graph to see what happens.

```

Q:  After 400 days, what do the radial velocity measurements
reveal to us?

Q:  After 400 days, what do the transit measurements reveal to us?

```

Day 401 - 500 Click on the graph to see what happens.

```

Q:  After 500 days, what do the radial velocity measurements
reveal to us?

Q:  After 500 days, what do the transit measurements reveal to us?

Q:  What would happen if the weather was bad, so that
it was only possible to make one measurement
every week or so?

```

#### Complications in transit detection

Well, okay, the signals are small -- maybe very small -- and they are rare. Fine. We have marvelous detectors, such as CCDs, and we have lots of time. Surely it can't be so hard to find these transits?

In one sense, no, it isn't. After all, the first transit ever observed was measured with the 4-inch STARE telescope:

Image courtesy of the STARE project

But, in practice, there are a number of complications which make the detection of transits a lot harder than one might think at first.

extinction
The light from a star must pass through the Earth's atmosphere to reach our ground-based telescopes. As we watch a star over the course of the night, it moves from horizon, to zenith, to horizon; the amount of air (airmass) traversed by its light changes ... quite a bit.

For example, during the course of a long run on WZ Sge on July 27, 2001 at the RIT Observatory, I measured the following raw instrumental magnitudes for a reference star:

```

Q:  How large are the variations due to extinction?

Q:  How large are the variations due to transits?

```

You can read more about extinction in a lecture from the Observational Astronomy course.

There are ways to remove most of its effects. I recommend very, very strongly that one use multiple comparison stars in the same field as the target and perform differential photometry. Even better is to go whole-hog and apply inhomogeneous ensemble photometry. You might even be able to grab some some free software to do the work for you.

Oh, and all this becomes a lot worse if there are any clouds, of course.

```

```

flatfielding
The response of a typical telescope + filters + CCD to light is, alas, not uniform. When one takes a picture of a white card, or the twilight sky, one does not see a uniformly bright image. Instead, one sees something like this:

Ugh! What an ugly mess of large-scale gradients, big donuts, small donuts, and tiny little spots.

```

Q:  What causes the big and small donuts?

```

It might not be as bad as you think; the image above was created with a relatively high contrast. To get a quantitative feel for the size of these variations in sensitivity, examine this slice through the image.

```

Q:  What is the size the sensitivity variation
across one of the small donuts?

```

Now, clearly, these variations cause errors in ordinary photometry of stars across the field of view: those near the center would seem brighter than they ought, those near the edges or under donuts would be fainter.

```

Q:  But does flatfielding matter for transit observations?

```

Given the performance of the RIT 12-inch telescope in guiding measured on May 6, 2013 (and shown below), yes, a failure to correct for flatfield variations WOULD be a big problem.

```

```

photon noise
There's no such thing as a free lunch. Even the signal from the target star itself can cause problems. Since the arrival of photons from the star is a Poisson process, the actual number which reach our detector varies from one second to the next. The statistics are pretty simple:
```

RMS of number of photons  =   sqrt(avg number of photons)

```

One can derive an expression for the fractional error in each measurement of a star due to this random variation in the number of photons measured each time:

```
sqrt(avg number of photons)
fractional uncertainty    =   -----------------------------
(avg number of photons)

1
fractional uncertainty    =   -----------------------------
sqrt(avg number of photons)

```

And, as long as this fractional uncertainty is small, say, less than 0.1, there's a neat little numerical coincidence which ends up giving us an easy-to-remember result:

```
1
uncertain in mag   (approx) =   -----------------------------
sqrt(avg number of photons)

```
```

Q:  So, in order to reach an uncertainty of 0.001 mag,
how many photons must we collect in each image?

Q:  So, in order to reach an uncertainty of 0.0001 mag,
how many photons must we collect in each image?

```

That sounds like a lot of photons ... and it IS a lot of photons. Fortunately for us OPTICAL astronomers, each photon has so little energy that typical stars send us a boatload of them in even a short exposure.

You can find many CCD signal-to-noise calculators which will do the work for you (such as this general-purpose one or this one for NOAO ). Or, you can do it yourself, using the rough starting point

```

star of mag 0 produces  106 photons / sq. cm. / sec

```

Well, this doesn't seem so bad: just exposure for a minute, or two minutes, or ten minutes, and you can beat down the photon noise as much as you'd like. Right?

Wrong!

CCDs are nice, but they aren't perfect. You can read about the "dark side" of CCD detectors for more details, but one particular feature of a CCD is its full well capacity. If too many photons hit one pixel, and liberate too many electrons within that tiny volume of silicon, the electrons will start to leak out into the surrounding pixels.

We call this saturation, or, when it gets worse, bleeding. If you see long trails

or stars with "flat tops" in a radial profile, you should be on the lookout for errors in photometry.

A typical CCD might have a full-well capacity of 100,000 electrons per pixel. Take a look at the closeup of a star in an image below.

```

Q:  Estimate the number of electrons you could collect
from a star in a single image, without saturating
the CCD.

Q:  What is the resulting uncertainty in the measurement?

```

Hmmmm. That might be good enough for giant planets, but it's not precise enough to detect small ones. Is there anything we can do to improve the dynamic range of our CCD measurements?

Yes! We can simply de-focus the telescope! The MOST satellite turned its stars into big donuts:

The COROT satellite did, too, though not to the same extreme.

```

```

scintillation noise
As starlight fights through the Earth's atmosphere, it interacts with little pockets of air. Since these pockets have different temperatures and densities, they refract the light by different amounts; and since the pockets are moving, the light is constantly being tossed this way and that.

The result is scintillation, which is a fancy word for "twinkling."

The appearance of this twinkling depends on the size of one's telescope.

• small telescopes see only 1 or 2 bundles of light rays at a time; the image appears to swoop across the detector, moving this way and that
• large telescopes collect photons from many bundles of rays all at once; adding together all these little speckles yields an image which doesn't move around, but seems to "boil".

The image below was taken by a very large telescope.

Astronomers have studied the nature of the scintillation noise for many years. Some good references are

This scintillation causes an uncertainty in each individual measurement which depends on several factors. A rough estimate can be calculated like so:

```

-0.67           1.75                       -0.5
sig  =  0.09  * (D)      *  (sec Z)      *  exp(-h/h0)  * (2T)

where

sig      is the RMS deviation from the mean intensity (fractional)

D        is the diameter of the telescope aperture in cm

sec(Z)   is the airmass

h        is the altitude of the site

h0       is the scale height of the atmosphere (about 8000 m)

T        is the exposure time, in seconds

```

For the RIT Observatory, which is only 168 meters above sea level, scintillation can be important.

If we pick some particular site, telescope and exposure time, then the scintillation noise is constant. Let's compare it to photon noise.

#### Putting it all together: high-precision photometry is hard!

If you aren't a regular observer, working frequently with the real data, you might over-estimate the quality of your measurements.

A number of years ago, while I was working on the SDSS, I checked the precision of photometric measurements of stars in the SDSS catalogs.

Has this inconvenient truth been inserted into the SDSS catalogs? Well, I recently (Mar 19, 2015) made a query using the SDSS website. Take a look at what I found:

Wierd things often appear at the bright end of stellar distribution. For example, when I examined some images taken with the HDI camera on the WIYN 0.9-m telescope, I found -- well, you take a look: