Mettre tous les spectres (archive de ces spectres) sur lesquels faire les mesures dans un répertoire.
Définir ce répertoire comme répertoire de travail : dans Audace, menu Configuration/Répertoires.
La commande est disponible depuis le menu "SpcAudace/Astrophysique/Construit le spectre dynamique 2D d'une série de profils".
Synthaxe : spc_dynagraph lambda_deb lambda_fin lambda_reference interpolation(o/n) RA_d RA_m RA_s DEC_h DEC_m DEC_s
Prendre les coordonnées de l'astre sur Simbad, format FK5, ici :16 00 20.005 -22 37 18.14
Création du spectre dynamique : spc_dynagraph 6545 6585 6562.82 o 16 00 20.005 -22 37 18.14.
Le résultat sont les fichiers image del_Sco_dynagraph.png et del_Sco_dynagraph.ps.
Remarques : on voit bien les variations du pic bleu de la raie d'émission, les tirets noirs le long de l'axe vertical de droite indiquent la position des profils utilisés.
La fonction spc_dynagraph est un outil tout en un qui :
met les spectres à l'échelle au continuum sur une longueur d'onde commune ;
corrige de la vitesse héliocentrique ;
découpe sur la zone demandée ;
trace les axes gradués et légendés autour du spectre 2D dynamique ;
produit trois versions du spectre dynamique : en niveaux de gris, colorisé et postscript (utile pour les publications).
Un exemple de spectre dynamique (cliquer dessus pour agrandir) obtenu par cette fonction appliquée sur plus de 250 spectres de eps Aur à l'occasion de la rédaction d'un article traitant de la compagne d'observations 2010-2011 :
Application d'une commande à tous les spectres :
Mettre tous les spectres (archive de ces spectres) sur lesquels faire les mesures dans un répertoire.
Définir ce répertoire comme répertoire de travail : dans Audace, menu Configuration/Répertoires.
La commande est disponible depuis le menu "SpcAudace/Astrophysique/Construit le spectre dynamique 2D d'une série de profils".
Synthaxe : bm_cmd "la commande audela a executer sur les fichiers notes %s"
Découpage sur une plage de longueur d'onde d'une série de spectre : bm_cmd "spc_select %s 6550 6570"
Normalisation du continuum d'une série de spectre : bm_cmd "spc_autonorma %s"
Ces actions sont appliquées d'une traite à tous les profils de raies présents dans le répertoire de travail. Cela fait gagner du temps !
Détemination de la température de surface d'une étoile :
Définir ce répertoire comme répertoire de travail : dans Audace, menu Configuration/Répertoires.
La commande est disponible depuis le menu "SpcAudace/Astrophysique/Déterminationde la température stellaire par ajustement de Planck".
Synthaxe : spc_ajustplanck nom_profil_calibré ?pas du calcul (1000)?
Réalisation de la mesure de la température par pas de 500 K : spc_ajustplanck vega--profil-1c.fit 500.
Conformité du résultat : Véga est une étoile de type A0V, donc de température comprise entre 7500 et 10000 K. Dans notre exemple, la mesure de la température de surface est correcte (9602 ± 180 K, Astronomy and Astrophysics 391 (3): 1039-1052). Cependant, la méthode par ajustement d'un profil de Planck est plus délicate pour les étoiles froides et nécessite des spectres corrigés de l'absorption atmosphérique pour parvenir à un résultat plus fiable.
Analyse des spectres - mesures de largeur équivalente (EW) :
Mettre le spectre à étudier (del_sco_20100429.fit) dans le répertoire de travail.
La commande est disponible depuis le menu "SpcAudace/Astrophysique/Largeur équivalente d'une raie d'un spectre complexe".
Synthaxe : spc_ew nom_profil_raies_calibré lambda_debut lambda_fin ?taux_doucissage_continuum (0-[6]-15)? ?efface_continuum(o)? ?degré_polynomes_continuum_methode_passebas(2)?
Réalisation de la mesure de EW sur l'intervalle 6550-6578 A contenant la raie en émission : spc_ew del_sco_20100429.fit 6550 6578.
Résultat de la mesure dans la console :
# Date: 2010-04-29T00:39:53.00
# JD: 2455315.5277
# EW(28.00=6550.00-6578.00)=-8.344 A.
# Sigma(EW)=0.093 A.
# SNR=261.
Mesure avec affichage du continuum utilisé pour la mesure : spc_ew del_sco_20100429.fit 6550 6578 6 n
Évidemment, la valeur de EW mesurée est la même que précédemment car 6 est la valeur utilisée par défaut pour le doucissage du continuum.
Mesure avec un continuum adpaté pour la mesure : spc_ew del_sco_20100429.fit 6550 6578 4 n
Résultats de la mesure dans console :
# Date: 2010-04-29T00:39:53.00
# JD: 2455315.5277
# EW(28.00=6550.00-6578.00)=-6.755 A.
# Sigma(EW)=0.089 A.
# SNR=261.
On constate que le continuum utilisé a une influence non négligeable sur la valeur de EW mesurée. Je vous renvoie à la lecture de la page décrivant les paramètres influençant la valeur mesurée de EW.
Il est clair que le continuum utilisé pour cette deuxième mesure épouse mieux le continuum de l'étoile : la mesure de EW n'en sera que meilleure.
Mettre tous les spectres (archive de ces spectres) sur lesquels faire les mesures dans un répertoire vide.
Définir ce répertoire comme répertoire de travail : dans Audace, menu Configuration/Répertoires.
La commande est disponible depuis le menu "SpcAudace/Astrophysique/Courbe EW=f(t) d'une série".
Synthaxe : spc_ewcourbe lambda_raie/lambda_deb lambda_fin ?degré polynome du continuum(1)?
Réalisation des mesures de EW entre 6540 A et 6585 A sur tous les spectres du répertoire de travail (comme décrit précédemment) : spc_ewcourbe 6540 6585
Le fichier ewcourbe.dat contenant les mesures de EW de la série de spectres est crée ainsi que le graphique ew_courbe.png au format PNG :
Pour chaque mesure, une incertitude est évaluée et est reportée dans le graphique.
Détermination de la période de variation par périodogramme à partir de mesures obtenues :
Mettre tous les spectres (archive de ces spectres) sur lesquels faire les mesures dans un répertoire.
Définir ce répertoire comme répertoire de travail : dans Audace, menu Configuration/Répertoires.
Réalisation des mesures de EW entre 6540 A et 6585 A sur tous les spectres du répertoire de travail (comme décrit précédemment) : spc_ewcourbe 6540 6585
Le fichier ewcourbe.dat contenant les mesures de EW de la série de spectres est crée.
La commande est disponible depuis le menu "SpcAudace/Astrophysique/Construit le periodogramme et affiche la sinusoïde des mesures issues d'un fichier ascii".
Synthaxe : spc_periodogram data_filename.dat time_unit measured_quantity ?nb_periodes_plausibles (10)? ?period_min (0.00001)? ?period_max (=duree enregistrement des mesures)? ?valeur minimum autorisee pour le pas d'echantillonage du periodogramme?
Recherche de périodes sans a priori : spc_periodogram ewcourbe_20100421.dat "julian days" "EW (A)"
Remarques :
Une fenêtre affiche le périodogramme et une deuxième l'ajustement sinusoïdal obtenu avec la période trouvée ayant l'intensité la plus importante.
On constate dans le graphique du périodogramme qu'un maximum vers 30 jours est plus problable. On doit donc affiner la recherche.
Affinons la recherche de la période en supposant que 10 périodes sont couvertes et que la période vaut au plus 40 jours : spc_periodogram ewcourbe_20100421.dat "julian days" "EW (A)" 10 0.01 40
Les fenêtres des deux graphiques colorés en rouge sont superposées. Il faut en déplacer une pour voir le périodogramme et l'ajustement sinusoïdal confortablement.
Remarques :
Dans la console sont affichés les maxima trouvés dans le périodogramme.
Une fenêtre affiche le périodogramme et une deuxième l'ajustement sinusoïdal obtenu avec la période trouvée ayant l'intensité la plus importante.
Le périodogramme est sauvé au format ascii DAT et PNG pour usages ultérieurs.
Affichage des résultats dans la console :
# Graphique sauvé sous periodogram.png
# Le nombre de maxima 2 trouves est inferieur au nombre demande
# Maximum N°1 : période valant 32.9705078125 et d'intensité 3.91682490474
# Maximum N°2 : période valant 19.5363671875 et d'intensité 2.12947752888
#
Ajustement sinusoidal :
# Sinefit : RMS= 1.30594942471
Estimated amplitude : 0.89
Phase : 0.19
Form of the fitting function:y(t)=0.84+0.887013362985*sin(2pi*t/32.97-0.19)
# periodogram.dat