Topic: Variable Star--RR Lyrae
Team TYC-cepheids
Name : | Roll no: | Contribution : |
Amrit Gupta | 195120038 | Observation Plan, Photometry, Distance to the RR-Lyrae |
Anish Chaurasiya | 180260007 | Target selection, Photometry |
Sudarshan Das | 195120033 | Target selection, Photometry |
Koushik Das | 195120031 | Photometry, Period Verification |
Goals:
Our target is ASAS J094809+0000.1 star which is a variable star (RRab- Lyrae)
Our first goal is to plot the light curve for this star
Second goal is to use the light curve to find the distance of star from earth.
Third goal is to estimate the period of our star using Lomb-Scargle Periodogram algorithm.
Target selection:
Pre data analysis:
-After selecting our target star, based on the period of star we uniformly spread our time of observation to get all the variability in star
-As asked we got our data, but due to some bad weather condition we missed our one day data which had 2 observations. In order to compensate we have been provided 2 observation on the next day.
-To check whether we covered the entire cycle we plotted them to see the phase coverage
-Since it was uniformly spread across its phase we were good to go.
Softwares used:
Origin
Python- Jupyter Notebook
Aperture Photometry Tool (APT)
SAOImageDS9*
Google Sheets
Data analysis: Major Steps (for single image)
-Reference stars - Zero point
-Checking: variable stars
-Eliminate Background sky noise
-Magnitude and uncertainty of ref stars
-Calculate final mag and error of target star
Data analysis: SAOImageDS9
Data Analysis: APT (Aperture Photometry Tool )
-Now the next step is to find the magnitude of all the star using APT
- Open the same image in APT
- The first thing is to subtract the background sky noise
- For this we have to select a sky annulus median subtraction (sky algorithm from ) more setting option
- Then we have to set the aperture overlay and sky annulus overlay , generally the radius of the inner circle of sky annulus should be 3 times the radius of aperture overlay
- The image below shows the image in APT with aperture overlay in red and inner and outer sky overlay annulus
-Then you have to also see the Radial profile of your image, in this you have to minimize the chi square value by choosing a good sigma value. From the value of FWHM you have to set the radius of red circle equal to FWHM and for green and yellow circles accordingly as mentioned
-Note that the above analysis we will do only for our target star, we will use the same parameters for all the reference stars
-Now as we know the x,y coordinate of stars we will read the magnitude and the errors for all the selected image
Data analysis:
Now we have to calculate Zero point
Zero point mag = absolute r mag - calculated mag in APT
Apparent mag of target star = calculated mag + Zero point mag
Error in magnitude,
Data Analysis:Magnitudes and errors of all the images we got:
Plots:
Light Curve:
Simple Mag plot including error bars:
Light curve with 2 phase :
Plots:
Light curve:
Magnitude plot with smooth curve:
Magnitude plot with smooth curve for two phase:
Distance to star:
Absolute magnitude of star
To calculate the absolute magnitude we use the formula
Where M is absolute magnitude of star
MSUN is absolute magnitude of Sun. MSUN = + 4.74
L is luminosity of star
LO is solar luminosity
We taken the luminosity of our star from the GAIA DR2 catalog.That is 56.749 LO .
Using all these values, we have calculated the absolute magnitude of our star to be 0.355104 .
Distance to star:
We calculate the distance to the star from the Earth using the formula,
Where m is apparent magnitude
M is absolute magnitude
d is distance from earth in parsec
Apparent magnitude is found to be 11.82008 magnitude from the area under the light curve.
Calculated value of distance of star from earth is 1963.3886 parsec.
Results: Light curve
Actual Light curve:
Taken from ASAS-SN variable stars database
Plotted Light curve:
Using Python notebook
Results : Period Estimation
Actual period:
Actual Period of this star is 0.833903 days (20.013672 hours), as taken from ASAS-SN variable star database.
Estimated Period:
Estimated period of the star using the Lomb-Scargle Periodogram algorithm is found to be 0.81314 days.
Results: Distance
Actual Distance:
Actual Distance of this star is 1847.84 parsec , as taken from ASAS-SN variable star database.
Calculated Distance:
The distance of the star is found to be 1963.39 parsec.
Here we are not able to calculate the uncertainty in our calculation because the method we used to find the area of light curve (average apparent magnitude) does not give the uncertainty in the area.
Error:
Period of star
The calculated period deviate from the value given in ASAS-SN variable star database by 2.49%.
Cause: It may be due to less number of observations available to us.
Distance of star
Our calculated distance deviates from the distance given in ASAS-SN by 6.25%.
Cause:It may be due to the error in estimation of area under the light curve.Also there are uncertainties in the calculated magnitude.
Link of Python Notebook for period verification-
https://colab.research.google.com/drive/1QuBIDFusSo1_CnzWhPprsVByHFzmpj9z?usp=sharing
Thank you !! hope you liked it.