Binary Supernova Ejected Runaway Stars Code Dias 2021 [https://vizier.cds.unistra.fr/viz-bin/VizieR-3?-source=J/MNRAS/504/356&-out.max=50&-out.form=HTML%20Table&-out.add=_r&-out.add=_RAJ,_DEJ&-sort=_r&-oc.form=sexa ]
Sample Color-Magnitude Diagram
Part 0: Import necessary libraries Import the functions necessary for the code:
% run runaway_functionsv2 % matplotlib qt This also imports a list of young open clusters from (based on Dias+ 2021, Gaia DR2):
Table length=399 Cluster RA_ICRS DE_ICRS r50 Diameter r50_table2 N pmRA e_pmRA pmDE e_pmDE Plx e_Plx RV e_RV NRV Dist e_Dist logage e_logage __Fe_H_ e__Fe_H_ Av e_Av FileName SimbadName _RA.icrs _DE.icrs deg deg deg mas / yr mas / yr mas / yr mas / yr mas mas km / s km / s pc pc log(yr) log(yr) mag mag deg deg str16 float64 float64 float32 float64 float64 int16 float32 float32 float32 float32 float32 float32 float64 float32 int16 int16 int16 float32 float32 float32 float32 float32 float32 str30 str31 float64 float64 ASCC_107 297.1623 22.0071 0.156 20.88 0.174 59 -0.144 0.135 -5.158 0.141 1.118 0.055 -- -- -- 864 30 7.440 0.121 0.353 0.103 1.372 0.129 clusters1/ASCC_107.dat [KPR2005] 107 297.1623 22.0071 ASCC_114 324.9790 53.9990 0.180 25.92 0.216 149 -3.754 0.210 -3.435 0.145 1.063 0.039 -- -- -- 911 12 7.632 0.271 0.035 0.078 1.216 0.091 clusters1/ASCC_114.dat [KPR2005] 114 324.9790 53.9990 ASCC_127 347.1807 64.9151 0.541 75.24 0.627 113 7.490 0.261 -1.781 0.319 2.618 0.080 -11.267 2.676 16 365 10 7.496 0.131 0.152 0.115 0.668 0.080 clusters1/ASCC_127.dat [KPR2005] 127 347.1806 64.9151 ASCC_13 78.3057 44.4212 0.564 73.08 0.609 110 -0.477 0.111 -1.737 0.108 0.899 0.076 -- -- -- 1066 26 7.615 0.098 -0.075 0.078 0.915 0.027 clusters1/ASCC_13.dat [KPR2005] 13 78.3057 44.4212 ASCC_16 81.2025 1.6256 0.367 45.12 0.376 207 1.363 0.280 0.002 0.274 2.844 0.113 21.308 1.696 12 348 3 7.088 0.061 -0.062 0.069 0.224 0.045 clusters1/ASCC_16.dat [KPR2005] 16 81.2025 1.6256 ASCC_19 82.0035 -1.9617 0.613 72.6 0.605 173 1.112 0.263 -1.303 0.241 2.756 0.088 23.576 2.719 10 356 2 7.139 0.030 0.076 0.077 0.189 0.043 clusters1/ASCC_19.dat [KPR2005] 19 82.0035 -1.9617 ASCC_21 82.1423 3.4771 0.419 49.2 0.41 102 1.381 0.292 -0.610 0.237 2.893 0.132 15.313 3.818 8 343 5 7.102 0.038 -0.008 0.029 0.236 0.048 clusters1/ASCC_21.dat [KPR2005] 21 82.1423 3.4771 ASCC_32 105.7112 -26.5758 0.646 78.72 0.656 255 -3.317 0.232 3.475 0.126 1.240 0.067 34.607 4.624 10 792 11 7.432 0.022 -0.003 0.048 0.220 0.019 clusters1/ASCC_32.dat [KPR2005] 32 105.7112 -26.5758 ASCC_67 175.2892 -60.9906 0.165 21.96 0.183 46 -6.775 0.064 0.925 0.059 0.482 0.026 -- -- -- 1921 89 7.483 0.227 0.215 0.095 0.810 0.044 clusters1/ASCC_67.dat [KPR2005] 67 175.2893 -60.9906 ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... UPK_540 114.5354 -58.4348 0.766 98.64 0.822 48 -4.815 0.212 7.661 0.210 2.663 0.097 14.456 3.188 3 365 4 7.513 0.043 -0.008 0.070 0.448 0.071 clusters1/UPK_540.dat UPK 540 114.5354 -58.4348 UPK_604 224.3164 -59.8095 0.260 42.36 0.353 43 -4.548 0.144 -3.711 0.199 1.307 0.079 -- -- -- 745 9 7.113 0.492 -0.209 0.314 1.773 0.417 clusters1/UPK_604.dat UPK 604 224.3164 -59.8095 UPK_606 216.1298 -46.3629 0.716 93.84 0.782 46 -20.147 0.688 -16.551 0.686 5.882 0.184 10.435 2.725 7 167 2 7.231 0.142 -0.052 0.175 0.133 0.284 clusters1/UPK_606.dat UPK 606 216.1299 -46.3628 UPK_62 289.7268 20.8263 0.110 13.92 0.116 33 -0.452 0.111 -5.418 0.128 1.075 0.056 -- -- -- 885 21 7.039 0.246 0.021 0.254 3.421 0.257 clusters1/UPK_62.dat UPK 62 289.7268 20.8263 UPK_621 237.1990 -54.3853 0.425 56.52 0.471 57 -2.471 0.150 -3.101 0.100 1.126 0.058 -- -- -- 878 32 7.559 0.229 0.150 0.134 0.942 0.206 clusters1/UPK_621.dat UPK 621 237.1990 -54.3853 UPK_640 250.4137 -39.5740 1.231 163.32 1.361 540 -12.014 0.917 -21.350 0.779 5.666 0.239 1.174 2.002 50 173 1 7.379 0.091 0.149 0.102 0.189 0.101 clusters1/UPK_640.dat UPK 640 250.4138 -39.5739 vdBergh_130 304.4624 39.3404 0.049 5.88 nan 62 -3.609 0.308 -5.075 0.292 0.521 0.154 -- -- -- 1456 240 6.974 0.091 -0.029 0.163 2.356 0.042 clusters2/vdBergh_130.dat Cl VDB 130 304.4624 39.3404 vdBergh_80 97.7471 -9.6215 0.151 17.16 0.143 60 -3.285 0.430 0.481 0.361 1.026 0.112 -- -- -- 947 2 6.790 0.046 -0.148 0.091 1.726 0.219 clusters1/vdBergh_80.dat Cl VDB 80 97.7471 -9.6215 vdBergh_85 101.7288 1.3329 0.045 4.8 0.04 29 -0.973 0.147 0.345 0.164 0.550 0.049 -- -- -- 1720 167 7.104 0.125 -0.055 0.124 1.206 0.270 clusters1/vdBergh_85.dat Cl VDB 85 101.7288 1.3329 vdBergh_92 106.0426 -11.4884 0.114 13.44 0.112 154 -4.539 0.219 1.607 0.211 0.834 0.091 27.580 6.680 2 1114 42 6.749 0.074 0.025 0.087 0.984 0.062 clusters1/vdBergh_92.dat Cl VDB 92 106.0426 -11.4884
Part 1: Obtain the cluster class object using theget_cluster function from therunaway_functions.py, with thecluster_name as the input, we obtain the parameters of the cluster.
example usage:
cluster_name = 'ASCC_21' import os from astropy .table import Table ,Column from runaway_functions import get_cluster cluster = get_cluster (cluster_name )display (cluster )This imports all the details for the cluster.Various parameters of the cluster can be accessed:
Name Diameter (r50 Diameter: Diameter within which 50% of the cluster members lie) Number of Cluster members etc.All parameters can be accessed together by: cluster_name = 'ASCC_21' cluster = Cluster (cluster_name )cluster .all Row index=0
Cluster RA_ICRS DE_ICRS r50 Diameter r50_table2 N pmRA e_pmRA pmDE e_pmDE Plx e_Plx RV e_RV NRV Dist e_Dist logage e_logage __Fe_H_ e__Fe_H_ Av e_Av FileName SimbadName _RA.icrs _DE.icrs deg deg deg mas / yr mas / yr mas / yr mas / yr mas mas km / s km / s pc pc log(yr) log(yr) mag mag deg deg str16 float64 float64 float32 float64 float64 int16 float32 float32 float32 float32 float32 float32 float64 float32 int16 int16 int16 float32 float32 float32 float32 float32 float32 str30 str31 float64 float64 ASCC_21 82.1423 3.4771 0.419 49.2 0.41 102 1.381 0.292 -0.610 0.237 2.893 0.132 15.313 3.818 8 343 5 7.102 0.038 -0.008 0.029 0.236 0.048 clusters1/ASCC_21.dat [KPR2005] 21 82.1423 3.4771
Part 2: Search stars in a region around the cluster Using thecalculate_search_arcmin function from runaway_functions, calculate the region to be searched around the cluster. by default it is$10\ \mathrm{pc}$ around the clusters (from the edge of the cluster). This returns an astropy quantity object and prints its value.
cluster .calculate_search_arcmin ()$124.771 \mathrm{{}^{\prime}}$
We can also visualize this search region using:
cluster .plot_search_region ()
Using this as the search radius for a conic search around the cluster center coordinates, we find a table of all the stars in the cone.
Getting runaways (all functions necessary included): cluster = Cluster ('Ruprecht_170' )cluster .generate_tables ()theoretical_data = theoretical_isochrone (cluster ,output = "table" ,printing = False )fs = cluster .read_table ('fs' )runaways = get_runaways (cluster ,fs ,theoretical_data )display (runaways )
