June 15, 2009


Today’s lecture was given by Dr. Dick Canefield,
his lecture focused on radiation in space. One of the
only ways scientists are able to study astrophysics is
 through radiation. He discussed continuous radiation,
line radiation, thermal/ non- thermal radiation,
emission/ absorptions, and photon mean free path. 



After the lecture I met up with Angela and she was able to help me resolve some of my problems with my program but there seems to be a problem with Dana’s programs that will allow me to rotate the topology to match the RHESSI data. A quick fix was to rotate the RHESSI data to match the topology but a minor problem occurred, the topology was not fitted to the magnetogram properly. 


I started to break down the spectrum into smaller intervals. 





Time start-1:54:40
Time end- 1:55






Time start-1:57:40
Time end- 1:58:30




Time start- 1:58:30

Time end-  1:59:40









Time start-1:59:40
Time end- 2:00:50







Time start- 2:00:50
Time end- 2:02:00






Time start- 2:03
Time end- 2:06








I then created a light curve  for the time period 2:03:07 to 2:06 which allowed me to see if I had an accurate amount of counts to make a RHESSI image.







I have enough counts to make images for energy intervals 6-12 Kev, 12-25 Kev, and barley enough for 25-50 but I made an image for it anyways.







Time start-2:03:07
Time end-2:06:00
Energy intervals-6-12






Time start-2:03:07
Time end-2:06:00
Energy intervals-12-25



Time start-2:03:07
Time end-2:06:00
Energy intervals-25-50






I mimicked Angela’s dissertation and did a RHESSI image for the beginning of the flare. 


Time start- 1:54:40
Time end- 1:55
Energy level- 30-100 keV





















June 16,2009




I continued to make corrected images and graphs of my flare.

I made an overlay of lower energy intervals over the foot point to show the movement of X ray flux.

Before the overlay I had to make RHESSI Images of each energy interval from 6-50 keV.
    




                Energy range: 6-12 keV



                          

             Energy range: 12-25 keV

                           

            Energy range: 25-50 keV

         
                           






The final product with the overlay of the four different energy levels including 30-100 keV.

The contour intervals are 90 and 70 percent of the max.


                                   








I made another spectrum graph from time 2:03:07 to 2:06 because
the original spectrum graph start at 2:03 and it was beginning an
attenuator change. The graph did not look as clear as this one below.

                         



                                    







I continued to work on the second half of the
 OSPEX tutorial.
As I was reading I realized how important the
spectral images are going to be for my research
and how this is going to be the challenging part of my project.
I learned that I needed to make another spectral graph
that included a wider range of time plus night time
up to 50-100 keV for background radiation. 







June 17, 2009






I continued to do the OSPEX tutorial and making notes of certain areas that were not clear to me so I was able to pick Angela’s brain.   She made the “Fit” part of the tutorial a little clear. Basically I will be playing with four different parameters to make a line that will best fit the spectral in 20-sec intervals.  Apparently the way you can tell if each component fits the line correctly is through a residual map which is a graph of the fit and the chi-square helps as well.
The different components are as follow:



Variable Thermal- optically thin thermal bremsstrahlung radiation function as differential spectrum seen at Earth units of photon per cm^2 s keV

2 Lines- Single Gaussian function (high) quality width in sigma



bpow-   Broken Power-law function with/without discontinuities in the derivatives.







I also realized I needed to make another spectral image with a higher range of energy.

Energy range- 3-300 keV

                              






I produced an imaged using the various components and a residual map but I’m not too sure if it is correct. Then I made the program loop through the time intervals starting from the one I was previously looking at.  By the end of the working day the program was still looping so I left it running over night.





                                             
                 


June 18, 2009




I was interested to see if the loop was finished since I never left a program running over night and to my surprise it was not finished. It finally finished and I was able to finish the tutorial. 

After I spoke to Angela I learned that I didn’t have to do any looping initially but now I know the procedure to do it.  We went over more details on how get a good fit of the spectrum. Unfortunately, the program was cutting of my spectrum curve in the smaller time intervals but after restarting the OSPEX GUI it proceeded to work and I was able to ask questions. I learned how to make smaller adjustment on the parameters of each component to get a decent fit for my curve.


For the time being I am to work on one time interval at a time to figure out the parameters. I know if I have a good fit by looking at the residual graph and the chi-square. The graph calculates the total components and the actual data to calculate the best fit.  When there is a good fit the line on the residual graph would have a maximum and minimum of 5 but it is centered at origin and goes across the x axis.  The chi-square value should be around one for a good fit. The residual map below had a chi-square value of 79 and I managed to lower its value to 64.39 (this mean it is a horrible fit). The fit map only has the thermal, line 1, line 2, and it is missing the broken power law component.










                         


















June 19, 2009








I went to an amazing lecture presented by Dana Longcope about the radiation zone and the convection zone. The lecture gave me more insight of how the sun is operating and creating a magnetic field.

I met with Angela and went over the spectrum graph I produced with the fit line and discussed the next step in the process.

  

         -Loop other time intervals:

            * with the thermal component only
            *with thermal and broken power law component
  
     - From the thermal component find the relationship between the energy and temperature using the
        Equipartition theorem- E= (3/2) *Kboltzman*T.

    -Use two thermal components without the broken power law after 2:05



Then Angela and I discussed my agenda for next week in her absence.
               -continue working on the spectrum graph with fit lines
               -continue working on the topology overlay
              - Look at the TRACE tutorial