--
AriannaSaba - 06 Mar 2019
https://www.eso.org/public/outreach/eduoff/cas/cas2002/cas-projects/germany_ngc1960_1/
https://arxiv.org/pdf/astro-ph/0003327.pdf
- in a log L - log T plot, pp chain has a gradient of 5.6, CNO 8.4
- New HR diagram computed. Fluxes were wrong last time. Data calibrated with the accepted values by adding a constant (the published B and V of the brightest star) to the magnitude V and B measured. The brightest star in the cluster is the brightest star also in our data. Now the published and the measured data agree better.
- No errors shown because too big. Need to fix them.
- turn-off is visible here in a plot of absolute magnitude vs Temperature
- data converted according to these equations:
absolute magnitude: M = m_V - (5 * np.log10(D)) + 5 , D = 1330 pc
temperature: T = 4600 * ((1/((0.92*i)+0.17)) + (1/((0.92*i)+0.62))) , i = m_(B-V)
- isochrones derived from these options and an age of 25 million years:
- tried to fit isochrone with a Z=0.02 to the HR diagram of V vs B-V, using data of absolute magnitude and B-V.
- The data downloaded had to be corrected to match the diagram, (B-V) + 0.225 and M_v + 11.3
-
- Without any modification of the data, other isochrones were plotted. This time the data plotted is in absolute magnitude vs B-V as well as the isochrones.
- with the correction in B-V only, adding a value of 0.225, probably given by the reddening, as also explained by the paper, the data seems to match very well
- tried to fit isochrones in a plot of absolute magnitude vs temperature but the curves were many orders of magnitude away from the data.
- Temperatures in the isochrones range between 4 and 3, Mv goes from 9 to -5.
# Isochrone by Pietrinferni - Cassisi - Salaris - Castelli 2006 ::: Standard Model
#===========================================================================================================================
# Scaled solar model & transformations (Castelli 1999) + BaSeL Library (2002)
#===========================================================================================================================
# Np= 2000 [M/H]= 0.058 Z= 0.0198 Y=0.273 t(Gyr)= 0.0250 normalized: 18 2 2006
#===========================================================================================================================
# (M/Mo)in (M/Mo) log(L/Lo) logTe Mv (U-B) (B-V) (V-I) (V-R) (V-J) (V-K) (V-L) (H-K)
#===========================================================================================================================
0.5000000000 0.4999948559 -1.39726 3.58803 9.382 1.044 1.452 2.059 0.970 2.795 3.571 3.659 0.137
0.5392574727 0.5392574727 -1.29476 3.59741 9.031 1.076 1.440 1.937 0.957 2.666 3.466 3.542 0.129
0.5785149454 0.5785149454 -1.19031 3.60711 8.675 1.098 1.415 1.805 0.928 2.536 3.356 3.422 0.121
0.6177724180 0.6177724180 -1.08264 3.61834 8.304 1.101 1.370 1.654 0.873 2.399 3.226 3.287 0.110
0.6570298907 0.6570298907 -0.97368 3.63101 7.924 1.076 1.304 1.499 0.802 2.256 3.070 3.129 0.098
0.6962873634 0.6962873634 -0.87105 3.64296 7.570 1.030 1.235 1.378 0.743 2.126 2.911 2.966 0.088
JC:
- Was forced to update star detecton algorithm. Used histogram of raw data to re-influence threshold decision:
- Much more successful at detecting stars now:

- Calculated separation between highest flux star in both frames and plotted histogram of minimum separation between this star in clear frame and all corrected star coordinates in the blue frame:
- One can see the corrected coordinates of the highest flux star at d=0 and a subsequent spread for all the other stars
- Next, I plotted both the clear frame centroids and blue frame corrected centroids on same plot to see if stars had matched well:
- Following this, I plotted a HR diagram for these matched stars: