Friday, July 28, 2006

Gamma Doradus Star


Gamma Doradus Star
Originally uploaded by Astro Guy.

Another day, another neat light curve. Here is 44 days of continuous photometry of a gamma Doradus star. The top panel shows the DFT (amplitude, not power!) and the bottom panel shows the light curve binned by 40 minute intervals. The gamma Doradus stars are a relativity new class of variable stars found in the classical instability strip. They commonly have low amplitudes (4 mmag for this star!) and periods around 1 day, which makes them difficult to pick up from the ground. The low frequency comment (~1 c/d) are thought to be high-n, low-l g-mode (gravity is the restoring force). This star also shows delta-Scuti like pulsations around 9 c/d. Quite and interesting star and will serve well to test current ideas, such as mode identification schemes!

Thursday, July 27, 2006

Spherical harmonic animation.


Spherical harmonic animation.
Originally uploaded by Astro Guy.

So I'm still working on models to explain some extraneous frequencies in the echelle diagram for AQ Leo. The application of non-radial modes has been suggested - although I'm not really a fan. I think the science is more interesting just in explaining the existence of combination frequencies.

Well, to investigate the non-radial suggestion, one has to work with spherical harmonics, which are a set of solutions to Laplace's equation and are handy for decomposing pulsation modes in stars into a series of numbers. If you click on the photo above you'll be linked to my flickr account. Click on the "all sizes" button above the picture, and you'll get to see the animation of a pulsation mode described by the spherical harmonic of l=5, m=2. Of course a star does not pulsate with such large variations, but only a fraction of a percent of the radius, but it's much cooler to see the deformations amplified.

So now I'm inserting non-radial modes by brute force to see what kind of light curve is produced with the addition of rotation and radial pulsations. I don't expect to see anything dramatic, because the important physics is missing.. namely, boundaries and assymetrics within the star. Since an RR Lyrae is a helium burning horizontal branch star the pulsation should travel through convective and radiative zones. You dump a non-radial mode at a transition zone and the radial pulsation becomes preterbed, possibly setting up a standing wave from interference. In photometry this will show up as combination frequencies.

This is definitely work in progress.

Tuesday, July 25, 2006

AQ Leo Echelle Diagram.


AQ Leo Echelle Diagram.
Originally uploaded by Astro Guy.

A beautiful weekend in Vancouver! Temperatures went above 30C! I've got a nice tan now. Let just hope this lasts a while!

I started fiddling around with echelle diagrams on the weekend. An echelle diagram is a method of displaying evenly spaced frequencies. For instance, with the double mode RR Lyrae star AQ Leo you can extract about 60 significant frequencies from a discrete fourier transform. The y-axis on the diagram plots the frequency, the xaxis shows the phase of the frequency which depends on the folding frequency. For the case of AQ Leo I folded at the frequency of the fundamental pulsation mode (marked as F1 in the diagram). The cyan and magenta lines show where conbination frequencies are expected to occur. Almost all frequencies can be accounted for, except the magenta lines are not quite satifactory. I'm going to investigate if a second overtone pulsation will help matters.

Thursday, July 20, 2006

M51 - Spiral Galaxy


M51 - Spiral Galaxy
Originally uploaded by Astro Guy.

I made this photo a couple years back now. I got the observations at the Dominion Astrophysical Observatory (DAO) in Victoria British Columbia with 60 second exposure times in B,V and R filters. I'm in a rush out the door, so I'll comment more on this tomorrow.

Tuesday, July 18, 2006

Jupiter on July 3rd 2006


Jupiter on July 3rd 2006
Originally uploaded by Astro Guy.

Just over two weeks ago I sat down and figured out how to couple by digital camera to my little telescope. I have two cameras. A 35mm canon rebel (takes great shots, uses film) and a digital elph 3D300 also by Canon. This discussion involves the later.

I've seen some neat images of planets created by webcams to a lesser extent by digital cameras. The idea is to mount the camera in front of the eyepiece and project the image on to the detector. This can work quite well for bright objects. For an object like the moon, you can simply point the telescope and then hold the camera with your hands in front of the eyepiece. Once you find the right angle, the moon will show up on the LCD screen. So click away. You can get some cool results!

Planets can be a bit more tricky. I have succeeded in imaging Saturn and Jupiter by hand holding the camera, but this is not ideal as even the steadiest hand will shake. The trick to getting a good image is to combine multiple images. So to help make this easier I pulled out the toolbox. I used an L-bracket and a clamp to hold the camera in front of the eyepiece. I have a f4.5 3inch refractor and I used a 2.5mm long relief eyepiece - the long relief really helps.

At this point I aim the telescope at Jupiter and put the camera into video mode. I recorded at 30 frames/sec for 30 seconds. Next I transfer the AVI on to my computer using gtkcam on Linux. I can convert the AVI into individual PNG images with the mplayer software (just use -vo png). These png images can then be converted into FITS images for the 3 seperate colours with the image magick software "convert". I use FITS format as I do astronomy CCD processing with IRAF. In IRAF you can use the immatch.xregister program to remove any x-y offsets and allign all the frames. The key parameters are to set dxlag and dylag to INDEF to use the x-y shifts from the previous image and to set regions parameter to bracket planetary image.

Once you have a set of alligned images, I compute the average of the all the frames in each colour with imcombine. This gives you a single image for each bandpass (blue, green and red). To make the final colour image, I loaded each fits file with The Gimp and made a colour composite. Then I played with the brightness levels and so forth to bring out some detail. You can see my final product above.

To really improve the image, I like to get my hands on a video capture device that outputs RAW data. No MPEG/JPEG compression. The compression routines kill fine detail. I think I could do a much better job with raw images. I've been investigating webcams, but I don't know much about right now (and I'd like one that works with Linux).

Sunday, July 16, 2006

M55 Blue Stragglers.


M55 CMD
Originally uploaded by Astro Guy.

This is a colour brightness plot (more commonly known as a colour-magnitude diagram or CMD) of stars in the globular cluster M55.
The black dots show the measurements for each star. There is probably about 4000 data points. Overlayed are various models that describe the age of the cluster.

A globular cluster is a great object to understand the evolution of stars. All stars in the globular cluster were born at the same time - around 13 giga-years (Gyr) ago (that's 13 000 000 000 years ago!). The more massive the star, the faster it burns hydrogen. As the hydrogen fuel source runs out, the structure of the star changes. This can be observed as change in colour and brightness. So the 0.1 Gyr line on the plot shows what the diagram would look like 100 million years have the birth of all the stars. The 13 Gyr line does a good job describing what the cluster looks like to us now. For dim stars (V less than 20) you can see that all the model tracks match. These stars,
because of their low mass, evolve very slowly, where as all the bright
massive stars are gone. By locating what is known as the cluster turn off (the bend in the data points at V=18, V-R=0.4) we can estimate the age of the cluster.

Now all this is well known, but what interests me are the data points that have squares and triangles around them. The triangles are RR Lyra stars, like AQ Leo I talked briefly about in a previous post. The squares are pulsating blue stragglers. They are metal poor delta scuti stars in the cluster. The problem, is that these don't match the age of the cluster. The isochrones (the model lines) show ages less than 4 Gyr. These stars are massive enough that if they were the age of the cluster they would of ceased core hydrogen burning a long time ago. They seem to be evolving too slow, hense the same Blue Stragglers.
(Side note: Originally abbrivated as BS --whoops-- are now seen as BSs in literature)
There are two theories to explain the existance of BSs:
1. Captured from the Halo/Disc of the Milky Way as the globular cluster oribits the galaxy.
2. Two stars collided and merged. (most like a binary system through three body interaction)
Each may be pauslible. My idea is to match the pulsation frequencies of the variable stars to different models with different characteristics - such as composition and rotation and see if the results are consistent with either of the two theories. There is also the possibility that both are correct. Should be exciting.

Saturday, July 15, 2006

Interesting...


Interesting...
Originally uploaded by Astro Guy.

Finished another dataset today and this interesting variable turned up! (The x-axis is in days and and the y-axes is in magnitudes).
The white points are the data and the red line is a fit to the 20 strongest frequencies in the fourier transform.
The the "blips" that occur every 5 and bit days is the interesting part! There is also a significant frequency at 3 c/d that has an amplitude of about 500 part per million (yeap,, you heard that right) which gives it a S/N greater than 5 in fourier space.

Another star to put on my "What is the stellar type?" list.

I personally like to reduce datasets blindly. So I don't introduce any personal bias, but I always hate the 1-2 day wait to get the HD number and see what this baby is!

The little town of Verrieres is one of the prettiest in Franche-Comte.


Book for blog.
Originally uploaded by Astro Guy.

Well, not every entry is going to be complete science related. Afterall, a person needs to have a life! I take the bus onto campus everyday and it takes about 45 minutes one-way. This gives me 1.5 h of "quiet" reading time everyday. I've developed a habit of always reading on the bus and I rarely pay attention to what is going on around me. Every now and then I'll forgot by book in the office or at home. In these cases I usually stare out the window and I'm always amazed at the changes that take place. Most of this change is due to construction, and I'll remark to myself, "Hey, that building wasn't there before...".

I just finished reading Scarlet and Black by Stendhal (at least that's the name he assumed for writing the book). It's always commented on as being good representation of France after the revolution (when the story takes place), but I was more capitivated by the hero of the book, Julien Sorel and his story. I think I pitied him more that critics due and as such, enjoyed it more (but what do critics know anyways, this is opinion). This is one of these books that I will hold on to for now to read again later. I usually give to give a book away after I've read it. I'll list this one in my top ten.

Friday, July 14, 2006

MOST satellite detects planetary transit.


MOST satellite detects planetary transit.
Originally uploaded by Astro Guy.

I'm in the middle of finishing a paper on the albedo measurements of an extrasolar planet; HD 209458b. I've been playing with the data for about 6 months. I'm reaching the point of blind hatred. I'll feel much better once I finish the draft of my paper.

I gave talk in Calgary about the work I've been doing on this system, which includes full photometric reductions of about 500 000 CCD frames and fitting for the albedo. The albedo of the planet can be detected by watching the contributed light from the planet change during the orbital cycle. So I watch the planet change phases. It is brightest when fully illuminated by the star (like the full moon) and dimmest when looking at the night side (like a new moon). The planet as transits (goes in front) and is eclipsed (goes behind) its parent star. The transit causes a 2% dip on the light curve as seen in the figure above. The eclipse depth is equal to the amplitude light variations caused by phase changes, which is about 1/20000. So we need stable photometry better than 0.05 mmag over 40+ days of photometry. This can only be done from space.

As I switch from data-reduction monkey mode into paper writing mode, my blogs will entail my writing experience and how progress is going.

Wednesday, July 12, 2006

Subpixel Sensitivity


Subpixel Sensitivity
Originally uploaded by Astro Guy.

I spent most of the day work on better routines to handle partial pixels in aperture photometry. Usually, a simple approximation is used to calculate the contribution of pixels at the aperture boundary. These approximations can add systematic errors to extracted photometry. So I decided to actually integrate each pixel with exact formulae. This ended up being a logic brain bomb. Lots of if statements, but in the end it works. Time to test it out. I also want to include a bilinear pixel interpolation scheme to take into account the slope of the PSF when using small radii. I'll attack that tomorrow.

Now some real issues when dealing with CCD images. The figure shows the dependance of the measured instrumental magntitude as a function of the PSF centroid. There is a clear dependance. Since the PSF is unsampled (FWHM < 2 pixels) the CCD filling factors come into play and have to be corrected for. That'll be an exercise for next week.

Tuesday, July 11, 2006

V471 Tau


Eclipsing System.
Originally uploaded by Astro Guy.

Here is a beautiful light curve for the binary system V471 Tau. It consists of a main sequence star (K2V) and a hot white dwarf (DA). The dip at phase 0.25 is when the white dwarf passes in front of the main sequence star. The overall pattern is due to tidal distortions of the main sequence star (gravitational interaction with the white dwarf) and the occurance of star spots.

Sunday, July 09, 2006

Figured it out...


Crazy Light Curve.
Originally uploaded by Astro Guy.

I sat down today and figured out a scheme to fix the correlation between the PSF size and the measured instrumental magnitude. This was very difficult for star shows large variability (like the one in the light curve pictured). So I extracted the low frequency information from the lightcurve and then examined the correlations, fit them and applied them to the original light curve. I iterated the low frequency fitted solution a few times to make sure any instrumental effects did not survive.

The lightcurve above shows a really interesting variable. I don't know the stellar type yet, but my first guess is a B or O class star that frequently show this kind of variability. I'll inquire more about this star tomorrow when everyone is in the office.

Saturday, July 08, 2006

Another problem...


Another problem...
Originally uploaded by Astro Guy.

The new observing scheme with the MOST satellite is already cause problems for me. This plot shows the correlation between the PSF size and the instrumental magnitude. Time to start hitting the source code and see if I can figure out what is going on.

AQ Leo with the MOST Satellite


AQ Leo with the MOST Satellite
Originally uploaded by Astro Guy.

This has been a good week. I've managed to figure out how to decorrelate the stray light signal from large amplitude pulsators. The attached image shows a 2 day blow up of a 35 day observing run on the double mode RR Lyra AQ Leo. It's an impressive target. The fourier transform shows over 30 frequencies which all be related to harmonics or combination frequencies of the main two modes (fundamental and first overtone). Now I'm racking my brain over a physical explanation for the combination frequencies - shock waves come to mind, but a difficult to model. I'll have to keep thinking about it.