Wednesday, March 9, 2016

Space and Time - Revisited


Well that was a long hiatus. Between pushing out a paper, TA-ing two classes, and getting distracted with fiction writing, this blog went way on the back burner for a couple of months.

I was frankly surprised and flabbergasted by the degree of interest my first post got. I just play an astrobiology type on the internet, I'm really just a molecular biologist with a hobby. I got a number of interesting pointers to recent research, a lot of interesting questions, and spoke to a number of very interesting people actually working in the field of astrobiology! As a result, I'm putting the next topic (our own solar system and its history) on hold for a little bit and going a little more in depth into our solar system's cosmological position in space and time. The more I look into this, the more clear it becomes to me that we are utterly typical in terms of where and when you'd expect a complex biosphere to form – on a coarse scale, at least. Some interesting questions still arise about the sort of star we find ourselves around.

Anyway, in my last post I mentioned a recent study (Sobral 2012, http://arxiv.org/abs/1202.3436) which looked at star formation rates across cosmological time using proxy measurements, finding a precipitous decline in star formation rates for many gigayears. I can do better than this – I think I can put some rough numbers on our star, planet, and biosphere's relative position in time and how typical it is drawing from multiple sources.

Numerous publications for decades at this point have performed similar measurements of star formation rate over time, and all those with sufficient detail reveal one, simple, important fact: the rate of star formation per unit mass in the universe as a whole is decaying exponentially, and has been for most of the history of the universe. I could point to a number of publications for this, but for the purposes of this post I will use Yuskel, 2008: http://arxiv.org/abs/0804.4008. This group collated a number of datasets and fit an equation to the star formation rate going all the way back to when the universe was less than a billion years old. I've reproduced a figure from Horuchi, 2010 (http://arxiv.org/abs/1006.5751) which while not actually being ABOUT star formation rates takes the fit from Yuskel and compares it to multiple data sets, supporting its validity:



There's a couple of things going on here. Firstly, note that the star formation rate is given in terms of something called Z – the redshift. Astronomers love talking about their observations in terms of redshift because it is directly observable and related to both distance and age of what you are looking at. You look and you see something that should be radiating at wavelength X instead radiating at wavelength X*(Z+1). That observation is generally incontestable. But to turn that observation into a distance or time you need to start messing around with cosmological models which, although they have been converging to a pretty tight focus in recent years are still open to refinement and change. It's generally better to keep your data in the raw observational form as a result.

Secondly, note the FORM of the equation of best fit. Since around Z = 1, corresponding to more than eight billion years ago, the star formation rate has been cratering. As the universe has expanded over that timeframe, star formation rate has very closely tracked [size of the universe]^-3.4th. The plateau of steady high star formation rates lasted only about 3 or 4 gigayears, long before our Sun formed. This represents huge bursts of star formation from young galaxies and giant ellipticals that formed and burned out young.

If I want to get an idea of the total number of stars that have been born and will ever be born, I need to turn this function of star formation rate at a given Z (expressed in solar masses formed per comoving cubic megaparsec per year so as to normalize for the expansion of the universe) into a function of star formation rate at a given absolute time, project it forward with reasonable assumptions, and integrate. Here is a graph of Z with respect to T from the best modern cosmological models:





  Redshift changes with time very differently when the universe is matter-dominated (blue line), has a roughly even amount of matter and dark energy (red line), and when it is dark-energy-dominated (green line). Whatever the case, I can use the relationship between redshift and time to turn the graph of star formation rate across redshift into star formation rate over time for the history of the universe so far:

The red dot represents the formation of the Sun, and the end of the blue line the present day at t = 13.82 gigayears. As you can see, the vast vast majority of star formation occurred before our Sun was born.

How to project this plot forward into the future in order to determine how many stars will ever live? I have made two executive decisions. Firstly, I have decided to simply project the exponential decay that has been occurring since Z=1 forward into the future, for reasons that I will get to later. Secondly, I have decided to break with a slavish devotion to the redshift numbers of the far future, once the universe has turned completely dark-energy dominated. Instead, I will keep the same function of Z with T and therefore the same function of star formation over time that has held since approximately Z=1, because I fail to see how simply increasing the distance between galaxy clusters more and more rapidly as will happen once the universe is dark energy dominated will change star formation rates within the clusters. It doesn't make a huge difference, going one way or the other changes the final numbers by less than five percent.

Anyways, after much obsessive faffing around with these equations in matlab to deal with the above set of assumptions, I projected star formation rates forward into the future to get the following graph. The red dot continues to represent the formation of the Sun, and the green dot the present day:


The decline is striking.

By integrating this curve forward, we can get an idea of where the Sun lays in the final star-order of the universe. I have normalized the graph to the (finite) number of stars that will ever exist when you integrate out to infinity:


In short, we find that the Sun was born when ~79% of stars that will ever exist already existed, and at the present moment ~90% of all stars that will ever exist already exist. Thus, the sun is a relative but by no means extreme latecomer to the universe, and despite existing near the beginning of an apparently open-ended universe its time of formation is not terribly special.

The universe is full to bursting with hydrogen and helium, only the tiniest fraction having been converted into heavy elements by being consumed in stars. Naively one might assume that all this gas would eventually condense down into stars one day. Recent results in astrophysics are suggesting reasons, however, that this probably won't happen – reasons for my continued use of the exponential decay of the universal star formation rate in my above analysis, as processes happening today continue. In my last post, I mentioned the life histories of various types of galaxies and how they suggested that star formation might be closer to finishing than starting. I feel compelled to go into a little bit more detail here, speaking more in terms of what happens to all the gas that could form stars but doesn't using recent astrophysical results.

When you look out into the universe, the vast vast majority of elliptical type galaxies are very red due to the age of their stars and are not forming stars, whether they have an internal reservoir of gas within their dark matter halo or not – see http://www.dailygalaxy.com/my_weblog/2014/02/giant-elliptical-galaxies-why-are-they-red-and-dead.html for a discussion on this. Big spirals are mostly forming stars at a steady clip, with only a few tapering down and turning 'green' or eventually 'red' from their initial 'blue' status. Recently, a project called GalaxyZoo which has automated and crowdsourced the analysis of huge numbers of new galaxies observed in the Sloan Digital Sky Survey has taken a very quantitative look at star formation across galaxy types in the universe, and come up with some striking conclusions:

These studies were able to get more information than the instantaneous rate of star formation, and look back along the history of the galaxies by looking at light of different frequencies – huge stars that don’t live long make lots of ultraviolet, stars like our Sun peak in the green light, while long lived stars peak in the red. They were able to see that among elliptical galaxies, the tiny fraction that are star-forming mostly show evidence of recently being involved in mergers, and that all those that are red and green colored show spectral patterns indicative of very rapid shutdown of star formation, faster than can be accounted for by star formation eating up available gas. They call this fast star-formation shutdown 'quenching'. Something about their formation, either primordially or via mergers of spirals, puts their gas into forms that cannot form stars. The prime suspect is the initiation of regular energetic outbursts from their large central black holes, heating the gas and rendering it too turbulent.

This actually dovetails interestingly with another problem in astrophysics: the 'cooling paradox'. As I mentioned, about 90% of the baryonic mass of the universe is in the form of X-ray hot gas clouds blanketing entire galaxy clusters (largely outside the dark matter halos of individual galaxies). This gas is ridiculously thin and immensely hot, and radiating energy rapidly in the X-rays. It turns out that when you figure out how much mass is in these gas clouds and how much energy they are radiating in the X-rays, they should cool and sink down to the centers of the clusters on a timescale of gigayears, probably turning into cool gas flows onto the large galaxies at the centers of these clusters. But they don't. Looking back in time across the universe they are at more or less the same temperature now as they always have been and never seem to cool despite the fact that they are radiating energy. In recent years, for various reasons (images of turbulence in the gas, calculations of the available energy) the prime suspect for the energy source keeping these gas clouds energized has become supermassive black hole jets.

Anyways, as for spiral galaxies, they were able to model the distribution they saw (most of which are forming stars at a steady rate, some of which are tapering off, and some of which are red and dead) as a mixture of populations. One population is forming stars at a steady slowly decreasing rate, much like ours. Another is quenching on a much slower timescale than ellipticals, indicative of a cut-off of gas inflow into their star-forming discs and star formation then slowly depleting their reservoir of gas over a 1-2 gigayear timescale, likely caused by events in their immediate galactic neighborhood disrupting the inflows of cool gas within their dark matter halos onto their star forming discs.

All of this suggests to me that the decline in star formation rate seen across the cosmological past represents larger and larger fractions of galaxies quenching along with slow decreases in star formation rates within individual galaxies, and that the numbers I produced above have at least some semblance of validity.  I can't, of course, rule out the possibility that galaxy quenching is a temporary or cyclical phenomenon  on very long timescales, or that there's some special subset of galaxies that will never quench and will use up all their gas.  But these numbers are a good start.  

****

These numbers are, however, numbers about STARS and the Sun's position in star-order.  And as vital as stars are for life, we don't live on a star.  We live on a terrestrial planet.  And this makes a difference.

The huge bolus of star formation early in the universe consisted of many low-metallicity stars at the starts of spiral galaxies and the fast star formation bursts of elliptical galaxies.  Many of these stars are probably not suitable for the creation of terrestrial planets and thus biospheres-as-we-know-them.  In order to get a handle on Earth's position in planet-order we need to normalize this.  I am utterly unprepared to do this rigorously since the astronomical community as a whole hasn't got a good handle on planet formation - if there's one incontrovertible takeaway from the Kepler mission, this is it.  I can, however, pull up some numbers  that are better than nothing and layer them on top of my star formation numbers and see what comes out.

I will be taking data from a paper (Behroozi & Peeples, 2015: http://arxiv.org/pdf/1508.01202v1.pdf) that made the rounds last year suggesting that Earth came in the 8th percentile of terrestrial planets (that is, 92% come after us).  Their conclusion is suspect because it comes from the assumption that ALL gas within galactic dark matter halos will eventually form stars.  However, the authors include a nice set of metallicity normalizations that they apply to the early universe that I will take in its entirety, as it is far too complicated for someone like me who doesn't study cosmology professionally to critique.

In this work, all stars are assumed to form terrestrial planets with a power-law dependence on metallicity and a sharp metallicity cutoff for 'gas planet' formation is assumed.  I will use both these functions, on the chance that either of them is relevant - again, nobody
really understands planet formation.  I find myself suspecting that the latter is more relevant due to some talks I've seen on solar system structure and planet formation modeling, but I really don't know.

This work finds that using a power law metallicity dependence, the Eartth is younger than ~83% of currently existing planets, and that with a sharp metallicity cutoff our planetary system is younger than ~64% of currently existing systems (both of these numbers shamelessly scraped from graphs).  Now I apply an almost certainly oversimplified and wrong assumption:  that the fraction of star systems that have been planet-forming has been constant since the formation of the Sun, on the theory that the pollution of heavy elements is more or less complete.

Applying these numbers, I estimate the following positions of the Earth and by extension our biosphere in planet-order:

Under the power-law metallicity assumption,
Earth shows up as younger than 72% of planetary systems that will exist.

Under a sharp metallicity-cutoff assumption,
Earth shows up as younger than 51% of planetary systems that will exist.  

In the words of the first person I showed these calculations to, these numbers are very interesting in that they are very boring.  We find ourselves in an unremarkable position in terms of planet-formation order. We are not early and other explanations for the so-called Fermi pardox must be invoked.

Wow this post has gotten long.  I still have some thoughts on other large-scale considerations about our position in space and time - namely why we find ourselves around a star as large and relatively short-lived as the Sun (larger than most stars) and not, say, 30 billion years into the lifespan of a small star born at the same time as the Sun but with a 40 billion year lifetime.  I will save that for another revisit to this topic.  Stay tuned for my original intent, talking about our own solar system.