Monday, July 15, 2013

Week of 7/8 - 7/12, 2013

On my plot for all cluster members in the EDisCS survey, I separated galaxies into "low" redshift z∼0.4-0.6 and "high" redshift z∼0.6-0.8. I say "high" and "low" because in reality, these galaxies are really far away. The "low" redshift galaxies are just a bit closer. Here is my improved graph on EDisCS cluster members:


In addition, I sub-plotted each EDisCS cluster by number of galaxies vs. log10 infrared luminosity. This is the breakdown:

Next week I plan to "fix up" these plots to make them look better as well as create a fits file for infrared luminosity, raw counts, corrected counts, and the error for both counts. I will send this data to Alissa Earle who will fit the luminosity function.

Monday, July 8, 2013

Week of 7/1 - 7/3, 2013

This week was cut sort by Independence Day, but I still accomplished a lot. On Monday, my summer research did a ninety degree turn. I switched from analyzing data from the WISE Survey to picking up Brendan Gallagher's project. My new task is to look at the EDisCS Clusters (redshift 0.4-0.8) and follow mostly the same instructions as I did for WISE. In other words, I have to fit the luminosity functions for galaxies in the field of EDisCS and in the clusters of EDisCS.

In doing this I have to look at the EDisCS catalog filled with photometric and spectroscopic information regarding the galaxies. There are flags built into the catalog which help me separate/choose the galaxies specific to each luminosity distribution.

Tuesday, July 2, 2013

Week of 6/24 - 6/28, 2013

This week I fit the function for a second sample of luminosity data using chi-squared least squares fitting. The plot looks like this:

On Wednesday, we visited Union College in Schenectady, New York for the June 2013 UAT Summer Research Telecon. Undergraduate students came from institutions including Skidmore, Lafayette, and West Texas A & M to present their research in astronomy. In addition, students phoned in from San Francisco, Hong Kong, and other places to do the same. I learned more about ALFALFA at the conference and a little bit about how the Arecibo Observatory in Puerto Rico works. 

When it was Siena's turn to speak at the Telecon, we presented our work concerning luminosity functions for the WISE, Local Clusters, and EDisCS Surveys. At this point we only have luminosity distributions for the surveys but will be fitting them next week. Debbie Johnson gave a preview of what  the real fitted data could look like in her breakdown of a sample data set like the one pictured above.

On Friday, I worked on adjusting my code for the luminosity distribution for WISE. To clean up the data to be read in, I flagged (in the NSATLAS catalog) only galaxies observed in the Sloan Digital Sky Survey (SDSS). Here is my python code to do this:


And here is the before/after transformation of my plot:



Monday, June 24, 2013

Week of 6/17 - 6/21, 2013

This week was dedicated to fitting sample luminosity functions. To do this I plotted a sample set of data (number of galaxies vs luminosities). I used the Schecter Luminosity Function:

This function can be fitted using chi-squared. Many fitting problems in python can be expressed as least-squares problems. Chi-squared looks like this:


The numerator represents the physical difference in distance of a data point to the best-fit curve. It is squared to account for points above and below the curve. The denominator, written as sigma squared, represents the poisson error of the point. The whole expression is summed to account for all data points.  In a perfect world, chi-squared should be as small a value as possible. The lowest value I was able to achieve was 50.5. This sounds high but it is a great improvement over my original value (6545.3) that I had before debugging my program.

In this exercise, I plotted two sets of sample data and fitted both using chi-squared. The programs for each are nearly identical except for the file that the program reads in and the number of bins used.
Here is my program:



The program turns up these graphs:







In addition to creating this program I researched the Hydrogen alpha line. The emission line occurs when a hydrogen electron falls from the 3rd to 2nd orbital. Since hydrogen is the most abundant element in the universe, the hydrogen alpha line is the brightest emission line in interstellar space.

Next week I will work more with fitting our luminosity functions as well as construct a website that will display our team's work thus far during Astro Research at Siena. On Wednesday June, 26th, we will travel to Union College to present our work and learn more about undergraduate research in the local area.

Friday, June 14, 2013

Week of 6/10 - 6/14, 2013

This week began with the continued addition of summaries towards our group document. The document holds summaries from Debbie Johnson, Brendan Gallagher, and myself on various topics in astronomy. For example, one paper I reviewed was by Sara Ellison concerning merging galaxies and ultra-luminous infrared galaxies (ULIRGs). ULIRGs are galaxies are galaxies with greater than or equal to 10^12 solar luminosities (these galaxies make our sun look very faint in comparison).

I also learned more about The Galaxy Luminosity Function. This function describes relative number of galaxies at certain luminosities. Using sample data, I was able to plot the Luminosity Function.

Schecter Luminosity Function: (a way to fit the Luminosity Function)
Sample set of data (number of galaxies vs. luminosities):




























































Friday, June 7, 2013

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


      1. 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).
      2. 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?
      3. 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













Friday, May 31, 2013

Week of 5/28 - 5/31, 2013

During the introductory week of AstroPhysics Research at Siena College, I was given the task of plotting, graphing, and analyzing data from the WISE (Wide-Field Infrared Survey Explorer) Survey. While I haven't begun working with its data, I have been developing some tools to help me in the future. This week, I have made significant progress on CodeAcademy: Python edition and completed UNIX tutorials 1, 2, and 4. Knowing the basics of Python and UNIX will allow me to work extensively with the WISE Survey Data and collaborate with Debbie Johnson (Local Cluster Survey) and Brendan Gallagher (ESO Distant Cluster Survey).

I learned valuable information about WISE from looking at the December 2010 version of The Astronomical Journal: The Wide-Field Infrared Survey Explorer (WISE): Mission Description and Initial On-Orbit Performance. In this paper, the design, mission, orbit, and scientific goals of the WISE Survey are outlined in detail.

On Friday, we used Python to create a program that examines the redshift of local cluster MKW11. The   histogram I created looks like this:

Next week, I will continue working with UNIX and Python to create more figures like the one displayed above.

- Mike Englert