You are here

Kevin Jardine's blog

Mapping TGAS - Part 1

Submitted by Kevin Jardine on 2 November, 2016 - 18:53

Mapping is all about providing physical context for data. On Earth, the question "Where am I?" is answered on a map by showing the user a set of hierarchical context, including

  • local streets
  • nearby buildings
  • elevation
  • parks and rivers
  • transportation systems

and so on.

A good map of the solar neighbourhood would have its own set of structures:

  • dust clouds
  • regions of ionized gas
  • hydrogen concentrations
  • star formation regions
  • supernova remnants

and so on.

I have placed dust clouds, some ionized gas and a few supernova remnants on the Tycho Galaxy interactive TGAS display.

However, the truth is that Gaia has already made this information outdated and much more accurate information will be available in a few years, especially once Gaia scientists publish a catalog containing the distances and spectral types for a billion stars. For example, by comparing a star's real spectral type with its colour index as seen from Earth, we can determine its reddening and therefore the dust that lies between Earth and that star. With reddening data for a billion stars, astronomers will be able to construct an incredibly detailed 3D map of dust and gas in the solar neighbourhood.

But for some things we don't have to wait.

In principle, the TGAS data in Gaia DR1 already allows us to produce detailed maps of the star formation regions within about 600 parsecs. The key is the location of the hot stars.

Hot stars tend to be young and young stars have not drifted far from the sites where they were born.

We can compute a star's temperature from its colour index. A star's colour index can be determined using the BT and VT magnitudes from the Tycho-2 catalog. Specifically,

Ci = 0.85*(BT - VT)

Usually a hot star is considered to be an O and B class star (Ci < 0) or even only O stars and B stars down to B3 (Ci < -0.18). However, we have to consider that many stars in the Tycho-2 catalog are reddened by dust, and so some hot stars might have a positive colour index.

Once we have the temperatures for the hot stars, we can create a temperature density function, T(x,y,z), that essentially tells us how close any point (x,y,z) in space is to hot stars (and how hot these are). It is this "hotness" function that will help us map the structure of star formation regions.

See Bouy, H., and J. Alves. "Cosmography of OB stars in the solar neighbourhood." Astronomy & Astrophysics 584 (2015): A26 for a similar approach using Hipparcos data.

So how can we create T(x,y,z) and how can we use it to map the solar neighbourhood once we have it? Check my next blog post for many more details.

A Void in TGAS

Submitted by Kevin Jardine on 27 October, 2016 - 19:39

In my last blog post, I drew attention to a hot star concentration that I labelled (or rather mislabelled) "Cepheus" in a temperature density image. Here is the image again:

I called the concentration Cepheus because it appears at about 95° galactic longitude and the constellation Cepheus is located around this longitude above the galactic plane.

However, it turns out that the concentration is created by a thin wall of hot stars located at a distance of -300 < z < -150 parsecs below the galactic plane. Here is a temperature density map restricted to -300 < z < -150 pc:

It looks like the "wall" (the brightest part of this image) is part of a larger complex that forms the boundary of an enormous void in the lower half of the first quadrant. Could it be a bubble?

I did a second height map animation that makes the wall and the "bubble" look quite impressive:

I did a preliminary calculation that shows that the centre of the "bubble" is somewhere in the direction of Aquarius.

Then, this morning, Gaia scientist Ronald Drimmel sent out this tweet with his latest TGAS completeness image:

The biggest gap is ... somewhere in the direction of Aquarius.

So there is a void in TGAS in the lower first quadrant, but in the data, not in space! Gaia had simply not scanned that part of the sky much yet when the first data release was prepared and TGAS is missing most stars in that direction.

The wall around the void appears simply because there are two nearby gaps in the TGAS data below the galactic plane and the "wall" is the narrow region in between.

I have been thinking of TGAS as a donut that starts to fade out around 600 pc with a hole around the Sun caused by a lack of bright (as seen from Earth) and high proper motion stars. And to a first approximation it is - but a donut with a few bites taken out of it.

The Mountains of Tycho

Submitted by Kevin Jardine on 26 October, 2016 - 09:50

Suppose that you wanted to make a map of Europe and all you had was a satellite image taken at night. You might start with something like the image below.

If you did a careful analysis of the distribution of the lights, you could extract quite a bit of information from this image, including the location of major cities and most of the coast line.

We have a similar situation with the TGAS data set. The distribution of the stars, especially the hotter stars, is by no means random. Using some mathematical tools, we can extract quite a bit of information about the solar neighbourhood out to about 800 parsecs (beyond this distance, the limited accuracy of the parallax measurements for even the brightest stars makes them impossible to place on a map).

One key tool is temperature density. The Tycho-2 catalog provides B and V magnitudes for almost all the stars. The difference B-V is called the colour index and it can be used to estimate the temperature of a star.

We are more likely to find structures to map using the hotter stars because these tend to be younger and younger stars are located close to the star formation regions within which they were born. (We can think of a star formation region as analogous to a city in a map of Earth.) Older, cooler stars often drift in random directions from their origin over time and so are less useful for mapping purposes.

Astronomers usually use the hottest O and B class stars to map star formation regions. These correspond to B - V < 0. However, I've been a bit more generous in my analysis because stars embedded in dust clouds can be reddened, increasing their colour index. So I've selected all the Tycho-2 stars with B-V < 0.1 to include some of the reddened B-class stars. In some cases this pulls in some hotter A-class stars but that should make little difference for the analysis.

I've interpolated Eric Mamajek's very useful table to convert colour index to effective temperature.

As usual, I am starting with the approximately 1 million stars in the TGAS data set with err/parallax < 0.2 for the reasons explained in my previous blog post on TGAS limitations.

In order to find structures, you have to have a way to aggregate individual star data. I've done this in two steps:

  • Bin the data
  • Smooth the data

In my first experiment, I calculated the x, y and z values in parsecs relative to the Sun. I defined my bins as all the stars with the same integer x and y values. For this first experiment, I ignored the z value, so this adds together all the stars with the same x and y parsec values above and below the galactic plane regardless of their z-height. I then added together the temperatures for all the stars in each bin with B-V < 0.1.

To smooth the data, I started by taking the square roots of the temperature sums to reduce the spikiness of regions with a lot of hot stars. I then used gaussian smoothing with a sigma (standard deviation) of 15 parsecs. The result of my first experiment is below. I have added the position of the sun at the centre, an arrow pointing in the direction of the galactic nucleus, and names for each of the four identified hot star concentrations. The full image (right to the edge of the rectangle) is 800x800 pc. You can see that the hot star density drops well before 800 pc.

It is much easier to visualise these density distributions as height maps, so I created and animated one in the 3D graphics application Blender. You can see the result on Youtube:

(I suggest going to full screen and right-clicking on the video to set the loop option as the animation is fairly fast.)

There are some surprising structures visible in these images, especially in the hot star concentration that I labelled Cepheus. I'll discuss some of them in my next blog post.

TGAS Limitations

Submitted by Kevin Jardine on 16 October, 2016 - 18:02

The Tycho-Gaia Astrometric Solution (TGAS) star parallax catalog, released as part of Gaia DR1 on September 14, 2016, was created by combining star position data from the Tycho-2 catalog (produced in the late 1990s) with observations from the first few months of Gaia observations. Because of the short scanning period and the dependence on older observations, it has a number of limitations.

The most obvious of these limitations is the accuracy of the data.

An important paper published in 2015, Bailer-Jones, Coryn AL. "Estimating distances from parallaxes." Publications of the Astronomical Society of the Pacific 127.956 (2015): 994 (read in arXiv), concludes that converting parallax measurements with errors to distances is not straightforward unless the estimated error/parallax ratio is < 0.2. If the error/parallax ratio is higher, not only are the results less reliable, but the formula depends upon a model of the distribution of the stars in the Milky Way. In other words, to place these higher error stars on a map, you essentially already have to have a map!

Is the lack of bright stars near the Sun real or a TGAS limitation? (View in Tycho Galaxy)

Half the TGAS stars have a parallax error of 0.32 mas according to the Gaia DR1 documentation. Plugging this error into Bailer-Jones's formula shows that most of the TGAS results can only be reliably placed on a map for distances less than 625 parsecs (about 2000 light years). The green disk in the image above shows this distance superimposed on an artist's model of the Milky Way. As you can see, 625 parsecs only covers the solar neighbourhood and does not even reach any of the galaxy's major spiral arms.

There are more limitations than this. According to the DR1 release notes:

  • Many bright stars at G≲7 are missing from Gaia DR1;
  • Sources close to bright objects are sometimes missing;
  • High proper motion stars (μ>3.5 arcsec yr-1) are missing;
  • Extremely blue and red sources are missing;

The net effect is that no naked eye stars and few stars close to the Sun are included in TGAS (and therefore appear on Tycho Galaxy). This might explain the odd dearth of bright stars near the Sun in the Tycho Galaxy map of the solar neighbourhood.

There is even more bad news. In addition to the 0.32 mas parallax measurement error for most of the TGAS stars, the Gaia DR1 release notes warn that the Gaia data may have a systematic error and that this error might be as high as 0.3 mas. If we add this to the 0.32 mas measurement error, Bailer-Jones's formula gives us a usable distance of 323 parsecs. This is not even as far as the Orion nebula and not that much further than the Hipparcos results from the 1990s.

I have ignored the possible systematic error in producing Tycho Galaxy but it does place a question mark over much of the map.

Another limitation is caused by Gaia's incomplete sky scans during the DR1 period. As the Gaia data that went into TGAS was gathered only over the first few months of the mission, some of the sky was not completely scanned. Gaia scientist Ronald Drimmel tweated the incompleteness maps:

So am I disappointed with all this? Not really. The scans will be completed and the naked eye and high proper motion stars will be added to a future release. Moreover, the science goal of the Gaia mission is to produce parallaxes with an error of 0.0067 mas for the brighter stars. This is smaller than the size of a euro coin on the moon as seen from the Earth, and is fifty times more accurate than achieved for the TGAS release.

The Gaia scientists say that with more observation and calibration, the mission is still on track to achieve this high accuracy over the next few years. Plugging this reduced error into Bailer-Jones's formula gives us a usable distance that is larger than the entire galaxy. Clearly, the real limitation is the brightness of the stars as seen from Gaia. For those stars not embedded in thick dust clouds, the effective range of Gaia will include about a billion stars distributed through most of the Milky Way on this side of the galactic nucleus.

That is a big map.

The Stars, Like Dust

Submitted by Kevin Jardine on 14 October, 2016 - 13:15

The stars, like dust, encircle me
In living mists of light;
And all of space I seem to see
In one vast burst of sight.

Isaac Asimov

September 14, 2016 was the day when everything changed in astronomy with the release of parallax estimates for more than 2 million stars. My last few blog posts have covered the momentous release of the Tycho-Gaia Astrometric Solution (TGAS) as part of Gaia DR1.

Today I have put up the Tycho Explorer, a zoomable, pannable, clickable interface to the Tycho-2 catalog, both as it appears in the sky (Tycho Sky) and as it appears as a face-on map of the galactic plane (Tycho Galaxy).

Tycho Sky includes 2.5 million stars from the Tycho-2 catalog with provided B and V magnitudes (which is almost all of them), and Tycho Galaxy includes more than a million of these stars with reasonably low parallax errors (measured error/parallax ratio < 0.2) stretching out to a distance of about 700 parsecs (2300 light years). At the highest zoom level you can click on each star to get more information and use links to see the star in Tycho Sky and Tycho Galaxy. Search boxes enable you to zoom into a star in the maps by entering an identifier.

Much more information about both the sky and face-on Tycho viewers can be found by clicking the Help links at the upper right of the viewers. The Tycho Explorer is still very much in beta. The infrared background to Tycho Sky has a number of artifacts and limitations, and not all the functions work yet on mobile devices. TGAS itself does not include any naked eye stars or many stars close to the Sun. I'll mention the limitations in more detail in my next blog post.

The Tycho Explorer will continue to improve and the next Gaia data release at the end of 2017 will have more and more accurate data.

Despite the limitations, the results are already amazing. We are on the verge of the greatest mapping project in human history - one that dwarfs that of the New World. The Milky Way is appearing in all its glory before our eyes.

You can download a 4K wallpaper image of some of Tycho Galaxy here (warning, 10 Mb file).

Update: I think I have fixed the issues for mobile and the Tycho Explorer should now work for Android and iOS devices. The interface is not very convenient for small screens (I recommend a large monitor - it looks spectacular on my 4K Philips monitor). If there is demand I might work on a tailored interface for phones.

The code is now in GitHub: Contributions are welcome.


Subscribe to RSS - Kevin Jardine's blog