Week of 6/3 - 6/7, 2013

This week brought on a whole new set of skills and tasks for my Astro Research. For starters, I was given a specific set of instructions to guide me for the rest of the summer. That list looks like this:

1. The Local Field (WISE)
1. overview of the survey
1. WISE
2. NASA-Sloan Atlas
2. tasks
1. identify environment of galaxies - are they in field/group/cluster
2. write a program that will fit a luminosity function to data - there should be some existing python programs that we can use, but this will require some research
3. match the NASA-Sloan Atlas with the galaxy zoo catalogs

3. plots
1. number of galaxies versus Mg
2. number of galaxies versus g-band luminosity
3. number of galaxies per Mpc^3 versus g-band luminosity
4. fit a schechter luminosity function to the above plot
5. now repeat using 22-micron luminosity
6. plot number of galaxies per Mpc^3 versus 22-micron luminosity for elliptical, spiral, and other (this requires matching the NASA-Sloan Atlas to the Galaxy Zoo catalog).
4. issues to think about
1. completeness - how well do we recover sources of a given flux, and how to we account for this when fitting a luminosity function?  This is going to be a little complicated because the NASA-Sloan Atlas is a mix of many different surveys.  For simplicity, we might first limit ourselves to using only the galaxies that are in the sloan digital sky survey.
2. how do we convert from 22-micron luminosity to total infrared luminosity, and why do we want to do this?
3. how sensitive is WISE compared to the Local Cluster Survey 24-micron observations?
5. compare with previous results and with results from local cluster survey

Just as promised last week, I have begun using actual data from WISE to construct plots based on the above guidelines. Struggling through it, I have completed the first two plots and cleaned them up to make them appear more professional. Here are my plots of number of galaxies versus g-band luminosity:

You may be wondering why I have created two histograms that use identical data. I wanted to explain what "bins" are and how they can be used for coding. The first histogram is constructed using 500 bins. That means that the data (luminosities observed in the g-band) is split into 500 separate containers (bins) . This explains why the maximum value on the y-axis (number of galaxies) on this histogram is much lower (5000 versus 25000) than on the lower histogram. On the lower plot,  the data is put into just 100 bins. Because there are less bins to put data into, more galaxies wil be put into the same bins. Hence, the maximum number of galaxies (y-axis) in a bin on this histogram will be much higher. In other words, there will be more galaxies per bin if there are less bins to put galaxies into.

In addition to building code for our data, I was busy summarizing papers written by astronomers working with the same kind of data that I am. I learned new information about luminous infrared galaxies (LIRGs) as it pertains to star formation rate (SFR), environment of a galaxy (field, group, or cluster), the "flatness problem," critical density, and what the virial radius means. Understanding the science definitely helps me work through all of the scientific terms and units thrown at me during research. 

As far as UNIX goes, I now have a better understanding of git commands and "jobs" commands. With git, I am more familiar with checking statuses, adding files, committing them, and finally pushing them to the cloud repository on github. I also learned how to work with the multiple jobs that can all be running at the same time and how to terminate them. Having less jobs running makes it easier to navigate through the important ones. 

Next week I plan to finish m article summaries and push forward in creating more plots. 

--Mike Englert

P.S. definition: foobar: a universal variable understood to represent whatever is being discussed