Sign in to follow this  
Vincent_L

Wolf-Rayet

Recommended Posts

Posted (edited)

bonjour,

 

j'y vais aussi de mes spectres d'étoiles de Wolf-Rayet ... .

 

Ces étoiles massives, extrêmement chaudes, expulsent leur couches externes à des vitesses de plusieurs milliers de km/s avant d’exploser en supernova. Les spectres de ces étoiles sont dominés par l’hélium, le carbone ou l’oxygène. Dans les spectres ci-dessous, le continuum ne permet pas de déterminer la température de la zone émettant les photons, le maximum n’étant pas dans le visible. Pour WR133, la partie bleue du continuum est probablement liée à la présence d’une étoile de type O dans ce système binaire.

Un calcul « à la louche » des vitesses des vents stellaires sur WR133 à partir de la largeur des raies de l’hélium ionisé (fit des pics avec une gaussienne), donne des vitesses de l’ordre de 1300-1900 km/s (20% de la vitesse des vents éjectés par une supernova !)…

Spectres obtenus avec  un Alpy600 sur un télescope Newton 250. Traitement des spectres avec ISIS, calculs et tracé en Python.

vincent

 

 

 

 

 

Image1.jpg

Edited by Vincent_L
  • Like 4

Share this post


Link to post
Share on other sites
Publicité
En vous inscrivant sur Astrosurf,
ce type d'annonce ne sera plus affiché.
Photographier la Lune
Guide complet pour la photographier de la Lune.
Information et commande sur www.photographierlalune.com
Il y a 5 heures, Vincent_L a dit :

Un calcul « à la louche » des vitesses des vents stellaires sur WR133 à partir de la largeur des raies de l’hélium ionisé (fit des pics avec une gaussienne), donne des vitesses de l’ordre de 1300-1900 km/s

Très intéressant :)

as-tu un lien vers la manip détaillée? il faut que je regarde dans mes docs

où prends-tu la mesure? à mi-hauteur?

il y en a d'autres des raies de l'hélium non? pourquoi choisir ces deux là?

tu t'es fait une feuille de calcul? 

bon j’arrête 

Share this post


Link to post
Share on other sites

Super intéressant !

J'ai aussi un faible pour ces super chaudes !

Un petit encart sur le calcul en Python et le beau graphe des étoiles données, voir un lien sur le calcul à faire...

 

J'ai quelques spectres de WR qui trainent, j'aimerai tester la manip.

Dans les gaz, tu peux ajouter l'azote pour les WN;)

Share this post


Link to post
Share on other sites

bonjour,

 

Gérard, Mizar:

pour le calcul des vitesses de vent stellaires, il s'agit d'utiliser la formule: v=(FWHM/lambda)*c , dérivé de l'équation Doppler,(voir par exemple l'excellent bouquin de Walker:Spectral atlas for amateur astronomers"); je n'ai pas de lien particulier, à part cette référence ...

 

la mesure de la FWHM peut être faite par ISIS ou VisualSpec je pense. Pour ma part, j'utilise un petit script en python qui vient fitter une gaussienne sur un pic (j'utilise le module lmfit pour ça) (cf dessous) et renvoie la position du centre de la raie et sa largeur à mi-hauteur (en fait j'utilise un modele gaussien + linéaire pour tenir compte du continuum si besoin).

par contre, il faut faire cette mesure sur une raie bien isolée, pas mélangée avec d'autres éléments, et c'est là où c'est plus compliqué ... . et aussi corriger de l'élargissement instrumental (pas fait, d'où le "à la louche" :) ). Mais bon, c'est pas fait pour une publi scientifique et on retrouve les bons ordres de grandeur !

 

pour tracer les courbes, tout est fait avec le module matplotlib en python (https://matplotlib.org/),

 

vincent

 

 

 

def fit_H(name, center, num_fig):
    idx=find_nearest(data,center)
    x=data[idx[0]:idx[1],0]
    y=data[idx[0]:idx[1],1]
    
    def gaussian(x, amp, cen, wid):
        return (amp / (sqrt(2*pi) * wid)) * exp(-(x-cen)**2 / (2*wid**2))

    def line(x, slope, intercept):
        return slope*x + intercept

    mod = Model(gaussian) + Model(line)
    pars = mod.make_params(amp=2, cen=center, wid=2, slope=0, intercept=1)
    result = mod.fit(y, pars, x=x)
    plt.figure(num_fig)
    plt.title(name)
    plt.plot(x, y, 'bo', markersize=4)
    plt.plot(x, result.best_fit, 'r-', linewidth=3)
#    print(result.fit_report())
    A=result.fit_report()
    A=A.split()
    print('position de ',name, '= ', round(float(A[52]),1), 'A ','(',center,' A)\t', 'FWHM= ', abs(round(float(A[60])*2.3548,0)))

 

  • Love 1

Share this post


Link to post
Share on other sites

Merci Vincent 

Je vais regarder ça de plus près 👍🏿

Share this post


Link to post
Share on other sites

un grand merci Vincent, à regarder à tête reposée ! xD

Share this post


Link to post
Share on other sites

Beaux spectres, les Wolf Rayets sont décidément impressionnantes !

Share this post


Link to post
Share on other sites

Create an account or sign in to comment

You need to be a member in order to leave a comment

Create an account

Sign up for a new account in our community. It's easy!

Register a new account

Sign in

Already have an account? Sign in here.

Sign In Now
Sign in to follow this