¿Alguien puede decirme cómo puedo encontrar los datos de brillo de estrellas para estudiar exoplanetas usando el método de tránsito? El archivo debe estar en formato de valores separados por comas (CSV) o cualquier otro formato que pueda convertirse posteriormente a ese formato.
Los siguientes datos de muestra muestran el brillo de la estrella junto con el tiempo. Usando estos datos podemos trazar la curva de tránsito. Quiero tales datos para la estrella real:
No conozco una fuente para los datos CSV directamente, pero si está de acuerdo con un poco de Python, esto se puede hacer con el Archivo de exoplanetas de la NASA . Mirando un ejemplo famoso (HD 189733b), si realiza una búsqueda de este objeto en la página principal, debería llevarlo a una página de información detallada y enlaces a conjuntos de datos. Al expandir la sección 'Información auxiliar' se mostrarán los enlaces para los datos (este enlace debería ser equivalente)
Los archivos de datos están en formato de tabla IPAC, pero la clase de AstroPy puede leerlos fácilmente Table
. A continuación se muestra un breve ejemplo de esto:
from astropy.table import Table
import matplotlib.pyplot as plt
# Read photometric table
phot_table = Table.read("https://exoplanetarchive.ipac.caltech.edu/data/ExoData/0098/0098505/data/UID_0098505_PLC_025.tbl", format="ipac")
# Subtract off integer part of first JD to make plotting easier
t0 = int(phot_table['HJD'][0])
plt.figure()
plt.errorbar(phot_table['HJD']-t0, phot_table['Relative_Flux'], yerr=phot_table['Relative_Flux_Uncertainty'], color='r', fmt="+", capsize=3)
plt.minorticks_on()
plt.xlabel('HJD-{:.1f} [days]'.format(t0))
plt.ylabel("Relative Flux")
plt.title("HD 189733b Transit Light Curve")
plt.savefig("transit_lc.png")
# Read radial velocity data
rv_table = Table.read("https://exoplanetarchive.ipac.caltech.edu/data/ExoData/0098/0098505/data/UID_0098505_RVC_001.tbl", format="ipac")
t0 = int(rv_table['JD'][0])
plt.clf()
plt.errorbar(rv_table['JD']-t0, rv_table['Radial_Velocity'], yerr=rv_table['Radial_Velocity_Uncertainty'], color='r', fmt='+')
# Zoom in on one of the nights where Rossiter_McLaughlin effect was being measured
plt.xlim(395.45, 395.70)
plt.minorticks_on()
plt.xlabel('JD-{:.1f} [days]'.format(t0))
plt.ylabel('Radial Velocity [m/s]')
plt.title("HD 189733b Radial Velocity Curve")
plt.savefig("transit_rv.png")
# Example export to CSV format
rv_table.write("transit_RV.csv", format='csv')
La parte superior del archivo CSV exportado se ve así:
JD,Radial_Velocity,Radial_Velocity_Uncertainty
2453946.612661,-2167.93,0.89
2453946.620219,-2155.63,0.85
2453946.627291,-2142.88,0.78
si esto es más fácil para un análisis posterior.
b--rian
b--rian
Usuario123
priyash mistry
fabricante de planetas
b--rian
fabricante de planetas
fabricante de planetas
priyash mistry
b--rian