The μ Tau Association of Stars

I’m a research advisor at the Planétarium Rio Tinto Alcan of Espace pour la vie, and I’ll use this page to keep you updated with some of our original research. I recently completed the study of a young association of stars named μ Tau. This study was just published in The Astrophysical Journal. I will take this opportunity to tell you about this association of stars, how we were able to determine its age, and what new avenues of research this discovery is opening up.

This blog post is also available in French here.

It is about two years ago that we started up this project, when a collaborator brought a very interesting set of stars to my attention during an email exchange. This collaborator is Eric Mamajek, and he is a researcher in astrophysics at the Jet Propulsion Laboratory, in California. He is probably one of the experts that best knows the scientific literature about young stars in the solar neighborhood. Several years ago, he had noticed that a bright star named μ (the Greek letter “mu”) of the Taurus constellation (Latin for “the Bull”)—this star is also called μ Tau for short—seemed to move on the sky in a way that is similar to a small number of other stars. He showed this to me with this figure that he made:

A figure that shows the sky position of the stars in Eric’s list (open circles). The straight lines show the direction in which the stars are moving on the sky. Each star would move to the tip of its respective line if we waited for 350,000 years. Credit: Eric Mamajek.

In the figure above, each circle represents the position of a star on the sky, and each line represents the direction in which the star is moving over time. This lateral motion of stars on the sky is called proper motion. The horizontal location of each star on this chart indicates its right ascension, which we designate with the Greek letter alpha (α). The vertical location represents its declination, designated with the Greek letter delta (𝛿). The right ascension and declination of a star represent its coordinates on the sky, and they are a bit analogous to how we use longitude and latitude to represent the position of an object on the surface of the Earth. Both systems of coordinates are even aligned with each other; right ascension and declination are a lot like longitude and latitude, but projected outwards on the sky. There is just a slight difference: because the Earth rotates, the grid of longitudes rotates a full turn every day with respect to the grid of right ascensions.

We can therefore think of the figure above as a small fragment of the full map of the sky, showing a group of stars that move together. The stars displayed on it move very slowly: we would need to wait about 350,000 years for the stars to reach the tip of their respective straight lines. However, high-precision instruments allow us to measure such small displacements after a few years only!

The size of the map shown above is about as large as the Big Dipper asterism. The size that this ensemble of stars spans on the sky becomes more obvious if we display them on a true map of the sky: they cover a large part of the Taurus constellation, and they even cross over parts of Aries (to the right) and Eridanus (at the bottom).

A representation of the stars in Eric’s list (full red circles) on a map of the sky. The white circles indicate other stars in the vicinity. Larger circles indicate brighter stars, and full lines mark the constellations. The star μ Tau is shown with an orange arrow. A high-resolution version of this map in vector format is available here. Credit: Marc Jobin.

If you have already seen my first blog post, you might have an intuition of why these stars are all moving together. They may represent the tip of the iceberg of a young stellar association, with hundreds of fainter stars that also move in the same way. Stellar associations are groups of stars that were born together from the same molecular cloud. They are similar to open clusters, but some of them contain a much smaller number of stars, and may be fragments that detached from an open cluster, or simply a small cluster that is slowly dissipating. This tends to happen over time as some of its stars get torn apart by the gravitational attraction of other clusters or massive stars on its path around the Milky Way.

The open cluster NGC 265, as seen by the NASA Hubble telescope. Credit: HST/NASA/ESA.

These young associations are of great value to astrophysicists, because they provide us with a set of stars that are all of the same age. This makes it much easier to characterize these stars, because it is extremely hard to measure the age of an individual star in isolation. Most of the stars in a young association only give us subtle hints on how old they may be, and it is only when we combine all of these hints together, over hundreds of stars, that we can determine a precise age for the whole association. If we are lucky, we might even find rare types of stars in the association that allow us to determine their age precisely. When this happens, the improved accuracy of our age estimation is reflected on all other stars in the association!

At the time when Eric found this potential new association of stars, we did not have enough data to reconstruct the velocity of each individual star in three dimensions—this requires a measurement of each star’s distance, and how fast they are moving towards or away from us, which we call the radial velocity. Reconstructing these 3D velocities would thus have required that researchers observe each star individually with specialized instruments at the telescope, during dozens of nights spread out over several years. This set of co-moving stars therefore remained quite obscure for years, without anyone else noticing them.

A 3D visualization of the co-moving stars listed by Eric. We initially get closer to the group of stars, and then the camera stops moving and we go forward in time through millions of years, to illustrate the common motion of these stars. The green circles are centered on the Sun with radii between 10 and 100 light years. A high-resolution version of this video is available here.

When Eric sent me this list in 2018, the situation was completely different. The first data set of the Gaia mission had recently been released by the European Space Agency (ESA). This data set contained extremely precise measurements of distance and proper motions for more than a billion stars, on top of radial velocity measurements for the million brighter stars. These data are revolutionizing our understanding of the galactic neighborhood, by giving us a 3D view of the positions and velocities of many nearby stars. The list of stars assembled by Eric, when viewed in 3D with the help of Gaia data, seemed consistent with his starting hypothesis that we may have been looking at a new association of stars.

At that point, I decided to start looking for additional faint members in this putative new association. To do this, I built a 3D model to represent how the positions and velocities of Eric’s stars are distributed in the Galaxy. During my postdoctoral fellowship at the Carnegie Institution for Science, I built a set of computer codes, named BANYAN Σ (Σ is the capitalized Greek letter ”sigma”). This computer tool uses a model of average position and velocity to identify new stars with similar motions. The power of this tool resides in its ability to identify even the stars for which we do not have all measurements required to build their 3D positions and velocities. This allows us, for example, to recover the members of a stellar association that are too faint for Gaia to have measured their radial velocity. We may know their lateral motion on the sky (proper motion), but not their velocity towards us, preventing us from reconstructing their 3D motion. BANYAN Σ uses the magic of statistics to predict which stars will have the right 3D velocity to match a stellar association!

A 3D visualization of the model that I built for the galactic position of the stars in Eric’s list, which are represented here with black spheres. The model is an ellipsoid, the 3D equivalent of an ellipse, and it is shown with a shaded orange region. Shadows of the stars are shown on each 3D plane with blue spheres, and shadows of the model are shown with orange lines. Thin green lines are also shown connecting the 3D position of each star with its shadow on the XY plane, to help visualize the 3D perspective. The X, Y and Z positions are indicated in units of “parsecs” in this figure, often used in scientific publications. One parsec is the same length as 3.26 light years.

The figure above shows the position of stars in Eric’s list along three Cartesian directions , X, Y and Z, which we also call galactic positions. In this coordinate system, the Sun is located at X = 0, Y = 0 and Z = 0 by definition. The positive values of X point towards the Galactic center, those in Y point towards the direction of rotation of stars around the Galaxy (clockwise). The Z coordinate takes positive values above the Galactic plane, and vice versa.

A visualization of the Sun’s position (red star) in the Galaxy, with a representation of the three axes for the XYZ coordinate system. The X axis points towards the center of the Galaxy, Y towards the direction of rotation of stars around the Galaxy, and Z points upwards, outside of the Galactic plane.

The 3D velocity of a star can also be described in a similar coordinate system. We call this the space velocity of a star; its components U, V and W indicate the rate at which its respective X, Y and Z Galactic positions are changing over time.

A visualization of the Sun’s position in the Galaxy (red star), with an example of an arbitrary velocity vector (blue array) located inside the
plane of the Galaxy (W = 0), with a small U component and a large V component.

Eric’s list of stars can be represented in this system of XYZ Galactic positions, and their 3D velocities can also be calculated as UVW space velocities. To calculate XYZ, we need to measure their position on the sky and their distance, but to calculate UVW we also need to know their proper motion and radial velocity.

Galactic positions (red dots) and space velocities (blue lines) for stars in Eric’s list. We can see that they are physically grouped together, and move in the same direction. The Sun is shown here as a black star.

As you can see in the drawing above, Eric’s stars have similar XYZ positions and UVW velocities. This is one of the most obvious characteristics of stars that belong in a young association: they all formed from the same molecular cloud, and they will all keep moving like their parent cloud for a few hundred million years (I explain this more in this post).

Stars that are very faint are usually the most difficult members of a stellar association to identify. This by how difficult it is to measure their radial velocity: we need to observe them with long camera exposures, using a very large telescope. Because of this, can’t simply calculate their UVW velocities and compare them to those of the known stars that make up a stellar association. We usually know two components of their velocity from their proper motion (which has two components), tangential to the sky. However, these two tangential directions are usually not aligned with the UVW axes. Even worse, the misalignment between proper motions and UVW velocities changes depending on the direction where a star is located. This makes it much harder to compare the velocities of two stars that are far apart in the sky. This comparison is the main goal of BANYAN Σ, the mathematical tool which I built. It uses a statistical approach to determine the probability that a star has the same UVW velocities as known stellar associations, even if we know only its proper motion.

Once I have built a 3D model of how Eric’s stars are distributed in XYZ and UVW, I used the kinematic properties of all 1 billion stars in Gaia—including stars without radial velocity measurements—and I calculated the probability that they share similar UVW velocities to Eric’s list of stars. This allowed me to build a list of 500 new stars that had a high probability of moving in a similar way to Eric’s stars. I then dug up dozens of scientific publications and databases to verify their properties in detail.

A 3D visualization of the final list of members in the μ Tau association. We initially move towards the group of stars, and then the camera stabilizes and we move forward in time, millions of years in the future, to illustrate their common motion. The green circles are centered on the Sun with radii between 10 and 100 light years. A high-resolution version of this video is available here.

This allowed our team to identify a set of 392 additional stars with properties that are consistent with a single population of the same age, and with similar XYZ positions and UVW velocities. We had thus recovered what really seems to be a new young association of stars! We named it the “μ Tau association” after one of its brightest members, as we usually do when we identify a new association.

If we place all of these stars in the μ Tau association on a map of the sky, it becomes more obvious how Eric’s list of stars was really just the tip of the iceberg!

A visualization of all stars in the μ Tau association (full red circles) on a map of the sky. Larger circles indicate brighter stars, but we used a different symbol size scaling in this figure to bring up even the faintest members of the μ Tau association, some of which have magnitudes fainter than 20. A high-resolution version of this map in vector format is available here. Credit: Marc Jobin.

One important tool to determine the age of a young association is what we call the color-magnitude diagram. The horizontal axis of this diagram shows the color of a star, and its vertical axis shows its intrinsic brightness. The intrinsic brightness of a star can be expressed as an absolute magnitude, which corresponds to the magnitude of a star if we placed it at 32.6 light years from us. Here, I will use the terms color-magnitude diagram and brightness specifically in reference to this intrinsic brightness.

A color-magnitude diagram showing all stars in our neighborhood of the Galaxy (in black), observed by the Gaia telescope of the ESA. Each small black and semitransparent circle represents the properties of a single star. Stars that are in the μ Tau association are displayed with red circles. The stars that would look redder to our eyes are located to the right of this figure, and those that would look bluer are located to the left. The brighter stars are located in the upper parts of the figure, and fainter stars are located lower down. Regions of this figure that are black and opaque indicate that a large number of stars overlap at this position. Most stars form a diagonal sequence in the middle of the figure (bright stars tend to be bluer). This feature is called the “main sequence”. A second sequence is visible to the lower left of the figure: it is made of white dwarfs only, which are much fainter than ordinary stars. We can see that one member of the μ Tau association is actually a white dwarf!

 The bluer stars are located to the left in the figure above, and stars that look red are located to the right. The brighter stars are located higher up, and vice versa. Black dots show the properties of all stars in the solar neighborhood, and red circles represent stars that belong to the μ Tau association. The bluer stars (to the left) tend to be brighter (higher up on the diagram); this is due to the fact that massive, bright stars appear bluer because of their hotter temperatures. Stars with lower masses are fainter, but also cooler; this gives them a redder color, and they are thus located to the lower right of the color-magnitude diagram. The horizontals position of stars in this diagram is thus an indication of their temperatures: cooler stars are located further to the right. If you are paying close attention, you might have noticed a few differences between the black dots (all stars in our neighborhood) and red circles (stars of the μ Tau association). For example, the low-mass stars of the μ Tau association seem brighter than those of the neighborhood that have the same horizontal position in the color-magnitude diagram.

This is a well-known property of young stars with low masses: during the few hundred million years that follow their birth, they slowly contract until they reach their final size, which they will keep for billions of years. When we catch them during their young age, they will appear brighter to us because they are still physically larger!

In the case of massive stars (at the upper left of the color-magnitude diagram), this trend reverses. These stars contract very fast after their birth, and they live much shorter lives overall. It is for that reason that there are very few stars of our neighborhood (black dots) at the upper left of the figure. In fact, the more massive a star is, the shorter its life will be. For example, the stars with a mass above 8 times that of the Sun only live for a few million years before they burst into a supernova, and give birth to a neutron star. Stars with about 3 times the mass of the Sun will live for about 300 million years, and then become red giants. When a star becomes a red giant, it puffs up and becomes much larger and cooler. You might have guessed it: when a star reaches its red giant phase, it moves up and to the right in the color-magnitude diagram. After a few more hundred million years, the outer layers of the red giant star will be completely ejected into space, and become a planetary nebula. At the center of this nebula, the core of the former red giant star will still be intact: we call this remnant of the star’s core a white dwarf.

The Messier 57 planetary nebula, also called the “Ring Nebula”. We can see gas in blue, orange and red that comes from the outermost layers of a red giant star that died about 7,000 years ago. The small, bright dot at the center of the nebula is the white dwarf born from the core of the former star. This image was taken by the NASA Hubble telescope. Credit: NASA, ESA and the Hubble Heritage (STScI/AURA)-ESA/Hubble Collaboration. A full-resolution image is available here, it is worth the click!

If we go back once again to the color-magnitude diagram that shows all stars of the μ Tau association, you might notice that it includes several stars which are very bright, which we only rarely find in our galactic neighborhood. The most massive of those is called 29 Tau, and it is about 5.4 times more massive than the Sun. We know that these stars only live for 80 million years, and this is a clear indication that the μ Tau association has to be quite young.

Two examples of typical life trajectories that will be taken by a low-mass star (in red) and a massive star (in blue) in a color-magnitude diagram. The red and blue circles show the respective positions of these two stars at their birth. The star will then follow the thick line as it slowly ages. Every perpendicular graduation mark along this line indicates that 100 million years have passed. The two small arrows (one red and one blue) indicate the respective locations where the low-mass and the massive star will reach their position on the “main sequence”, where the low-mass star will stay for billions of years and the massive star will stay for a few hundred million years. This is the “normal” phase of life of a star. I have not included graduation marks at the main-sequence locations because many of them would overlap. We can see that after having left its main-sequence position, the massive star (in blue) will take about 100 million years to go through the first half of its trajectory along the red giant branch as it puffs up and cools down. This will be followed by the rapid planetary nebula phase, giving birth to a white dwarf. This white dwarf, initially bright (for a white dwarf) and hot, will keep moving down and right in the color-magnitude diagram as it cools down and becomes fainter. As indicated by the large number of graduation marks, this cooling down of the white dwarf is a very slow process.

By comparing the color-magnitude position of all stars in the μ Tau association with those of other associations with known ages, we were able to estimate that the age of the μ Tau association is between 55 and 70 million years. When we do this, it is important to remember that the massive stars that belong to μ Tau are currently in a phase of their stellar life that is a lot more advanced than its stars of lower masses.

One very interesting discovery we have made while studying this young association is that one of its stars has even already died! A white dwarf named WD 0340+103 is moving in space very similarly to the μ Tau association. It is a very unique white dwarf, because it is still extremely hot, and its properties indicate that it comes from a massive star (6.7 times more massive than the Sun), and that it went through its planetary nebula phase only 270,000 years ago! This might seem like a large number, but it represents a time scale so short that we rarely ever encounter it in astrophysics. This discovery was very exciting for us, because physics-based computer models built right here at the University of Montreal allowed us to determine that the star which died and gave birth to WD 0340+103 was itself born about 60 million years ago—right within the age range we had estimated for the μ Tau association!

Even if the planetary nebula associated with the white dwarf WD 0340+103 has already dissipated, the white dwarf is still so hot that it is heating up the interstellar gas in its surroundings, making it bright in H-alpha maps (shown in red). 

This new association of stars will be of great use to researchers in astrophysics, especially so because we currently know very few examples of stellar populations with ages between 50 and 100 million years near us! Searching for exoplanets around the stars that belong in the μ Tau association will be very interesting, because this will allow us to see what the properties of exoplanets are at this specific age after their birth. This is something we can only rarely do, because we do not know the precise age of most exoplanet systems. By comparing the properties of exoplanets in the μ Tau association to those of others young associations—for example the β Pictoris association with its 24 million years of age, or the 110-million-years-old Pleiades association—, we will better understand how their atmospheres change over time, and how they might migrate onto different orbits around their star.

The study of the magnetic activity of stars as they age will also allow us to better understand the habitability of rocky exoplanets in orbit around them. We think that low-mass stars may have rocky exoplanets with the right surface temperatures for liquid water when they are in tight orbit around their star, which could allow life as we know it to exist. However, this close separation to their host star also means that these exoplanets would be withstanding incredible amounts of UV and X-ray radiation, especially when the star is young, that life as we know it would undergo mutations too quickly to survive for long.

There is also a whole other project that I am excited to begin working on, now that we found this new association. We can use it to look for “free-floating” exoplanets, which are wandering alone in space instead of in orbit around a star! The reasons why we need young associations to discover these rare objects are complex, and will be the subject of a future blog post!

I would like to thank Marc Jobin and André Grandchamps for comments and figure contributions.

Published by jgagneastro

I’m a researcher in astrophysics at the Rio Tinto Alcan planetarium of Espace pour la Vie, in Montreal.

2 thoughts on “The μ Tau Association of Stars

Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s

%d bloggers like this: