# Structure stellaire

Les modèles de structure stellaire décrivent la structure interne des étoiles de différentes masses et âges, ainsi que permettent de faire des prédictions sur la luminosité, la couleur et l'évolution future de ces étoiles.

## Transport d'énergie

À partir des différentes grandeurs mesurées et de simulations issues de différents modèles, il est possible de construire une image de l’intérieur d’une étoile, bien qu’il nous soit presque inaccessible — l’astérosismologie permettant littéralement de sonder les étoiles. En l’état actuel de nos connaissances, une étoile est structurée en différentes régions concentriques, décrites ci-après à partir du centre.

### Types de transport d'énergie

Les différentes couches d'une étoile transportent l'énergie produite vers son extérieur. Le transport d'énergie est accompli principalement par convection et transfert de rayonnement, mais la conduction thermique est importante pour les naines blanches.
La convection est le mode dominant de transport d'énergie lorsque le gradient de température est suffisamment grand pour qu'une portion du gaz dans l'étoile monte graduellement via un processus adiabatique. Quand c'est le cas, cette portion de gaz flotte et continue de monter si elle est plus chaude que le gaz environnant ; si, au contraire, la portion de gaz est plus froide, elle retombe à sa hauteur originale.
Dans les régions des étoiles avec de basses températures et une opacité assez basse pour permettre le transport par transfert de rayonnement, c'est justement le rayonnement qui est le mode dominant de transport.

### Types d'étoiles

La structure d'une étoile de la séquence principale dépend uniquement de sa masse. Pour cette section, soient des étoiles de ce type.
Pour les étoiles de masse solaire (0,3–1,5 masse solaire), la fusion nucléaire de l'hydrogène vers l'hélium se fait principalement via la chaîne proton-proton qui n'entraîne pas un gradient de température élevé. Le rayonnement domine donc dans la portion la plus centrale des étoiles comparables au Soleil. La zone plus à l'extérieur de ces étoiles est assez froide pour que l'hydrogène soit sous forme atomique, laquelle est opaque aux photons ultraviolets, ce qui fait que la convection y domine. Somme toute, les étoiles de type solaire ont des noyaux radiatifs et des enveloppes convectives.
Pour les étoiles massives (dont la masse est de plus de 1,5 masse solaire), la température du noyau dépasse les 1,8 × 10$^7$ K, ce qui fait que la fusion nucléaire de l'hydrogène vers l'hélium qui domine est celle du cycle CNO. Dans ce processus, l'énergie générée est proportionnelle à la température à la puissance 17, tandis que pour le processus proton-proton, l'exposant n'était que de 4. À cause de cette grande sensibilité à la température dans le noyau, le gradient de température y est suffisamment élevé pour le rendre convectif. Dans la portion extérieure de l'étoile, le gradient de température est plus bas, mais la température est assez élevée pour que l'hydrogène soit presque totalement ionisé, ce qui fait que cette portion de l'étoile est transparente aux photons ultraviolets. Somme toute, les étoiles massives ont des noyaux convectifs et des enveloppes radiatives.
Les étoiles de plus petites masses n'ont pas de zone radiative, ce qui fait que le transport d'énergie dominant est la convection, sur toute leur profondeur. Les étoiles géantes sont également tout à fait convectives.

## Équations de structure stellaire

Le modèle le plus simple de structure stellaire est un modèle quasi statique qui suppose que l'étoile est dans un état stationnaire, c’est-à-dire qui ne change pas avec le temps, et qui est symétrique de façon sphérique. Ce modèle consiste en quatre équations différentielles du premier degré : deux décrivent comment la matière et la pression varient par rapport au rayon ; les deux autres décrivent comment la température et la luminosité varient, toujours avec le rayon.

### Hypothèse d'équilibre thermodynamique

Dans une étoile, les modèles stellaires supposent souvent qu'il y a équilibre thermodynamique local (ETL). L'ETL n'est pas strictement atteint parce que pour une coquille (une couche de l'étoile, à un rayon r donné à partir du centre) ne « voit » pas la même température au-dessous et au-dessus d'elle. Toutefois, c'est dans la plupart des cas une bonne approximation parce que le libre parcours moyen d'un photon, $\lambda$, est beaucoup plus petit que la distance sur laquelle la température varie de façon considérable, c.-à-d. $\lambda \ll T/|\nabla T|$.

### Équations détaillées

Soient la densité de matière $\rho (r)$, la température $T(r)$, la pression totale (matière plus rayonnement) $P(r)$, la luminosité $l(r)$ et la génération d'énergie par unité de masse $\epsilon (r)$. Ceci, dans une coquille sphérique d'épaisseur ${\mbox{d}}r$ à une distance $r$ du centre de l'étoile.

#### Équilibre hydrostatique

La force vers l'extérieur de l'étoile, due au gradient de pression dans celle-ci, est exactement balancée par la force vers l'intérieur, due à la gravité.
    ${{\mbox{d}}P \over {\mbox{d}}r}=-{Gm\rho \over r^{2}}{\text{ ou }}{{\mbox{d}}P \over {\mbox{d}}m}=-{Gm \over 4\pi r^{4}}$

#### Équation de continuité de masse

L'intégrale de l'équation de continuité de masse, à partir du centre de l'étoile ($r=0$) jusqu'au rayon de l'étoile ($r=R$) donne la masse totale de l'étoile.
    ${{\mbox{d}}m \over {\mbox{d}}r}=4\pi r^{2}\rho$

#### Équation d'énergie

La luminosité quittant une coquille sphérique à un rayon r donne l'équation d'énergie :
    ${{\mbox{d}}l \over {\mbox{d}}r}=4\pi r^{2}\rho (\epsilon -\epsilon _{\nu })$
où $\epsilon _{\nu }$ est la luminosité produite sous forme de neutrinos (qui s'échappent en quasi-totalité de l'étoile sans interagir avec la matière) par unité de masse.
Hors du noyau, il n'y a pas de production d'énergie par fusion nucléaire. Donc, hors du noyau, aucune énergie n'est générée et la luminosité est constante.

#### Équation de transport d'énergie

Pour le transport d'énergie par conduction :
    ${{\mbox{d}}T \over {\mbox{d}}r}=-{1 \over k}{l \over 4\pi r^{2}}$
où k est la conductivité thermique.
Pour le transport d'énergie par rayonnement :
    ${{\mbox{d}}T \over {\mbox{d}}r}=-{3\kappa \rho l \over 64\pi r^{2}\sigma T^{3}}$
où $\kappa$ est l'opacité de la matière, $\sigma$ est la constante de Stefan-Boltzmann et la constante de Boltzmann est posée à 1.
Le cas du transport d'énergie par convection est plus compliqué. L'approche phénoménologique du modèle de longueur de mélange (en) donne :
    ${{\mbox{d}}T \over {\mbox{d}}r}=\left(1-{1 \over \gamma }\right){T \over P}{{\mbox{d}}P \over {\mbox{d}}r}$
où $\gamma =c_{p}/c_{v}$ est l'indice adiabatique, c'est-à-dire le ratio des capacités thermiques du gaz.

### Modèle polytropique

Le rayon et la masse se déduisent de la solution de l'équation de Lane-Emden, qui décrit la structure d'un objet dont l'équation d'état comme celle d'un polytrope, et qui est soumis à l'influence de son propre champ gravitationnel. C'est Jonathan Lane qui le premier propose cette équation, en 1870, dans un travail qui est présenté comme étant le premier visant à étudier la structure interne d'une étoile.

#### Polytropes

Un fluide polytropique, est un fluide dont la pression P est reliée à la masse volumique μ par l'équation d'état polytropique :
    $P=\kappa \mu ^{\gamma }$,
où κ est une constante et γ un nombre appelé indice adiabatique. En pratique, on définit l'indice polytropique n par :
    $\gamma =1+{\frac {1}{n}}$.
La raison du choix du terme d'indice polytropique plutôt que d'indice adiabatique tient à ce qu'en pratique, il peut y avoir des échanges de chaleur au sein du fluide considéré (par exemple du fait du flux d'énergie issu des réactions de fusion nucléaire au cœur d'une étoile). Au lieu de μ, certains ouvrages utilisent ρ.

#### Solutions pour la masse et le rayon

Le rayon et la masse se déduisent de la solution de l'équation de Lane-Emden :
    $R={\sqrt {\frac {(n+1)\kappa }{4\pi G}}}\mu _{\rm {c}}^{\frac {1-n}{2n}}\;\xi _{1}$.
    $M=-4\pi \left({\frac {(n+1)\kappa }{4\pi G}}\right)^{\frac {3}{2}}\mu _{\rm {c}}^{\frac {3-n}{2n}}\xi _{1}^{2}\left.{\frac {{\rm {d}}\theta }{{\rm {d}}\xi }}\right|_{\xi _{1}}$.
avec la quantité θ définie par
    $\mu =\mu _{\rm {c}}\theta ^{n}$,
et la longueur d'échelle ξ par
    $\xi =r{\sqrt {\frac {4\pi G\mu _{\rm {c}}^{2}}{(n+1)P_{\rm {c}}}}}$,
r étant la coordonnée radiale en coordonnées sphériques et l'indice c désignant les quantités évaluées au centre de la configuration. Le point où $\theta$ s'annule est noté $\xi _{1}$.
La longueur d'échelle est une variable de distance sans dimension qui provient d'avoir posé $r=\alpha \xi$, où les constantes sont rassemblées telles que :
    $\alpha ^{2}=(n+1)\kappa \mu _{c}^{{\frac {1}{n}}-1}/4\pi G$

#### Solutions particulières
## Physique du gaz et du rayonnement

Une étoile est constituée d'un mélange d'ions, électrons et photons. Donc la physique stellaire doit couvrir les propriétés du gaz d'ions et d'électrons, les propriétés du rayonnement, ainsi que celle de leurs interactions.

### Équations d'état

En thermodynamique, une équation d'état d'un système à l'équilibre thermodynamique est une relation entre différents paramètres physiques (ou variables d'état) qui déterminent son état. Dans un gaz parfait, il peut s'agir de la relation entre sa température, sa pression, sa densité et son volume. Dans une étoile, les modèles stellaires supposent souvent qu'il y a équilibre thermodynamique local (ETL) et une équation d'état relie la pression, la composition stellaire, la température et la densité. À partir de l'équation d'état caractéristique d'un système physique, en particulier pour une étoile, il est possible de déterminer la totalité des quantités thermodynamiques décrivant ce système et par suite de prédire ses propriétés.

#### Hypothèse de l'ETL

Pour le cas des étoiles de type solaire, les interactions entre particules peuvent être négligées et l'équilibre thermodynamique local (ETL) peut être supposé. Ceci vient de que le ratio de l'énergie électrique et l'énergie cinétique est de l'ordre de 1 %.

#### Calcul

La pression totale est la somme de trois pressions : les pressions d'ions, d'électrons et de rayonnement.
    $\ P=P_{I}+P_{e}+P_{rad}$
Les ions et électrons sont considérés comme formant un gaz, donc on appelle pression du gaz la somme des deux premiers termes :
    $\ P_{gaz}=P_{I}+P_{e}$

##### Pression d'ions

Suivant l'équation d'état à l'ETL :
    $\ P_{I}=n_{I}k_{\mathrm {B} }T$

##### Pression d'électrons

Pour un gaz parfait et classique d'électrons :
    $\ P_{e}=n_{e}k_{\mathrm {B} }T$
Si les électrons sont dégénérés, alors il se produit de la pression de dégénérescence :
    $\ P_{e,deg}={8\pi \over 15m_{e}h^{3}}p_{0}$
où $p_{0}=(3h^{3}n_{e}/8\pi )^{1/3}$.

##### Pression de rayonnement

La pression de rayonnement est due aux photons qui transfèrent leur quantité de mouvement aux particules de gaz lorsque ceux-là sont absorbés ou dispersés :
    $\ P_{rad}={1 \over 3}aT^{4}$
où $a={4\sigma \over c}$ est la constante de rayonnement. Ceci se dérive de la fonction de pression générale pour des collisions élastiques :
    $\ P={1 \over 3}\int _{0}^{\infty }pvn(\nu )d\nu$
avec, pour la lumière, que $p=h\nu /c$ est la quantité de mouvement, que $v=c$ est la vitesse des particules (c.-à-d. de la lumière) et que $n(\nu )d\nu$ est donné par la fonction de Planck.

### Processus adiabatique

Selon la première loi de la thermodynamique :
    $dU=dW+dQ$
et pour un processus adiabatique, il n'y a pas d'échange de d'énergie, donc dQ = 0 et finalement :
    $dU=dW=-PdV$

### Transfert radiatif

Les modèles simples d'étoiles supposent un objet stationnaire (indépendant du temps), à symétrie sphérique et en équilibre thermodynamique local. Cette dernière hypothèse suppose que l'émission est $\epsilon _{\nu }=\kappa _{\nu }B_{\nu }$ où $B_{\nu }$ est donnée par la loi de Planck.
De l'hypothèse de stationnarité on peut en déduire immédiatement que le flux d'énergie sortant de la sphère définie par le rayon r est constant
    $4\pi r^{2}H=C^{ste}$
Cette expression est valide à l'intérieur de l'étoile tout comme à l'extérieur.

#### Assombrissement du bord d'une étoile

L'observation d'une étoile montre que le bord de celle-ci est plus sombre que le centre. Ce phénomène est explicable à partir du profil d'émission dans la partie externe. Cette émission notée S et pouvant contenir un terme de diffusion varie avec la profondeur. On peut l'approcher localement par une variation linéaire en fonction de la profondeur optique $\tau$ comptée suivant la ligne de visée (problème de Milne)
    $S(\tau )=S_{0}(1+\alpha \tau )$ avec $S_{0}=S(\tau =0)$
L'intégration de l'équation de transfert fournit pour l'intensité émergente de l'étoile la valeur
    $I=\int _{0}^{+\infty }Se^{-\tau }d\tau =S_{0}\int _{0}^{+\infty }(1+\alpha \tau )e^{-\tau }d\tau =S_{0}(1+\alpha )\equiv S(\tau =1)$
L'intensité I le long d'une ligne de visée est égale à la valeur de la fonction source au point de cette ligne de visée situé à la profondeur unité. Ce résultat constitue l'approximation connue sous le nom de relation de Barbier-Eddington.
Si on se place en un point de la surface, il est facile de réaliser (figure ci-contre) que puisque l'on « voit » à τ = 1 on sonde dans la direction radiale des couches plus profondes que dans les directions proches du bord, cela à cause du facteur μ = cosθ reliant profondeur radiale et profondeur le long de la ligne de visée.
Au bord d'une étoile, comme le Soleil, on voit donc des couches plus superficielles qu'au centre. Si la température décroît vers l'extérieur, ce qui est la situation la plus générale (et en tout cas ce que prévoient tous les modèles ETL), le bord apparaîtra plus sombre que le centre. Ainsi s'explique le phénomène de l'assombrissement centre-bord sur le disque solaire.

#### Température effective d'une étoile

La température effective (ou apparente) $T_{\it {eff}}$ est la température donnant pour un corps opaque le même flux que la valeur $H$ de l'étoile
    $H=\sigma T_{\it {eff}}^{4}$
où $\sigma =5.670473\times 10^{-8}~W\cdot m^{-2}\cdot K^{-4}$ est la constante de Stefan-Boltzmann.
Cette quantité qui vaut 5770 K pour le soleil n'est pas la température de surface de celui-ci. En effet si on calcule l'intensité émergente à l'incidence $\theta =arcos~\mu$, celle-ci vaut, compte tenu de l'hypothèse ci-dessus pour le terme source
    $I_{\nu }(\mu )=\int _{0}^{\infty }S(\nu ,t)e^{-t/\mu }{\frac {dt}{\mu }}=\int _{0}^{\infty }S_{0}(\nu )(1+\alpha t)e^{-t/\mu }{\frac {dt}{\mu }}=S_{0}(\nu )(1+\alpha \mu )$
Soit un flux
    $H_{\nu }=2\pi \int _{0}^{\infty }S_{0}(\nu )[1+\alpha \mu ]\mu d\mu =2\pi S_{0}(\nu )$
Si on fait l'hypothèse de l'ETL, la fonction source en surface vaut
    $S_{0}=B_{\nu }(T_{0})$ où $T_{0}$ est la température vraie à la surface.
Le flux correspondant est
    $H=\int _{0}^{\infty }H_{\nu }d\nu =2\pi \int _{0}^{\infty }B_{\nu }(T_{0})d\nu =2\sigma T_{0}^{4}$
En comparant cette expression à celle définissant la température apparente, on voit que
    $T_{0}/T_{\rm {eff}}=(1/2)^{1/4}\simeq 0,841$
Par comparaison la résolution exacte du cas gris ($\kappa _{\nu }=C^{ste}$) donne
    $T_{0}/T_{\rm {eff}}=({\sqrt {3}}/4)^{1/4}\simeq 0,811$
ce qui montre que l'approximation d'Eddington est assez bonne. On voit que le paramètre de température effective possède une certaine signification physique mais cela n'implique pas que le rayonnement émergent suive la loi de Planck.
Pour reconnaître une étoile dans le ciel et en déduire ses paramètres physiques, les astronomes utilisent sa température effective et sa luminosité totale. Lorsque l'on porte les étoiles dans un diagramme basé sur ces deux paramètres (le justement célèbre diagramme de Hertzsprung-Russell, ou « diagramme HR »), les points ne se répartissent pas au hasard mais selon des séquences définies (comme la séquence principale). L'analyse de la structure de ce diagramme a permis d'établir la théorie de l'évolution d'une étoile. Les points ayant même abscisse dans ce diagramme ont même température effective et par conséquent, par définition de cette température effective, même énergie lumineuse émise par unité de surface. Il en résulte que l'ordonnée du point représentatif d'une étoile sur une même verticale est une simple mesure du carré du rayon de l'astre. Les étoiles les plus grosses (géantes ou supergéantes) sont situées dans le haut du diagramme et les plus petites (naines), dans la partie basse.

## Processus nucléaires

Dans les étoiles se succèdent des phases de fusion (qui peuvent avoir lieu dans le cœur même de l'étoile, ou dans les couches adjacentes à celui-ci) et de contraction. Ces fusions successives vont produire tous les éléments que l'on connait, jusqu'au fer.
Au début de la vie d'une étoile se trouve un nuage de gaz (principalement de l'hydrogène), qui, pour une raison ou une autre, va commencer à se contracter. Avec la contraction, la température augmente progressivement. Quand elle dépasse 10 millions de kelvins la température est suffisante pour que les noyaux d'hydrogène aient assez d'énergie pour vaincre la barrière coulombienne et ainsi fusionner. L'énergie libérée va contrebalancer l'effet de la gravité, et l'étoile atteint un premier équilibre.
L'hydrogène disponible à ce moment-là a été synthétisé lors du Big Bang.

### Fusion de l'hydrogène

L'hydrogène est le carburant majeur des étoiles et sa fusion est le premier maillon de la chaîne de nucléosynthèse. Il y a deux manières de transformer l'hydrogène en hélium :
* Le cycle proton-proton permet de transformer de l'hydrogène en $^4$He. Ce cycle est divisé en plusieurs réactions. La première réaction de ce cycle permet la transformation d'un proton et d'un neutron en deutérium. Ensuite, selon ce qui est disponible au cœur de l'étoile (donc selon le moment de son évolution auquel se produit cette réaction), par l'intermédiaire du tritium ou de l'hélium 3 avec lesquels les atomes de deutérium vont fusionner, l'étoile va fabriquer de l'hélium 4 ;
* Le cycle CNO lui se produit à plus haute température, car les atomes qui vont fusionner lors de ce cycle ont des masses atomiques beaucoup plus élevées que l'hydrogène. Une autre condition pour que ce cycle puisse participer à la synthèse d'éléments, c'est que les éléments carbone, oxygène et azote soient présents dans le cœur de l'étoile. La figure détails du cycle CNO en précise le cycle. Comme on peut le voir ce cycle se décompose en trois sous-cycles (I, II et III), qui vont chacun produire différents isotopes qui serviront de base au sous-cycle suivant.
Globalement, la réaction de fusion de l'hydrogène peut s'écrire de la manière suivante :
    $4p\leftrightarrow^{4}{\rm {He}}+2e^{+}+2\nu +E$.
Cette réaction de fusion de l'hydrogène est la plus exothermique de toutes les réactions qui vont se produire au cœur des étoiles. Comme les étoiles sont composées majoritairement d'hydrogène, elles disposent à ce moment de leur vie d'une grande quantité de combustible, qui leur fournit donc une grande quantité d'énergie. Cela explique pourquoi les étoiles passent la plus grande partie de leur existence dans cette phase de combustion de l'hydrogène.
Au bout d'un temps dépendant de la masse de l'étoile, la quantité d'hydrogène dans le cœur de l'étoile finit par ne plus être suffisante pour entretenir un taux de réaction assez élevé pour contrebalancer l'effet de la gravitation. Le cœur de l'étoile va alors se contracter. Lors de cette contraction, sa température va augmenter (comme n'importe quel gaz que l'on comprime) de même que sa densité. Quand la température atteint quelque 10$^8$ Kelvins, la fusion de l'hélium peut s'amorcer.

### Fusion de l'hélium

C'est lorsque cette réaction démarre que la structure en couches d'une étoile se forme. En effet, la température au centre de l'étoile est telle que les couches externes au noyau sont suffisamment chaudes pour que puissent s'amorcer des réactions de fusion.
Deux réactions permettent la transformation d'hélium en éléments plus lourds :
    $_{2}^{4}{\rm {He}}+_{2}^{4}{\rm {He}}\leftrightarrow_{4}^{8}{\rm {Be}}+\gamma$ (1)
    $_{4}^{8}{\rm {Be}}+_{2}^{4}{\rm {He}}\leftrightarrow_{6}^{12}{\rm {C}}+\gamma$ (2)
La réaction (2) a posé un problème à première vue insoluble. En effet on trouve une grande quantité de carbone $^{12}$C dans l'Univers. Mais des calculs théoriques ont montré que le béryllium $^8$Be avait une demi-vie de l'ordre de 10$^{-16}$ secondes. Cela signifie qu'il est hautement improbable que l'on puisse fabriquer du $^{12}$C en quantités importantes à partir du $^8$Be. Pour tenter de résoudre ce problème, Fred Hoyle avait suggéré que la réaction entre le $^4$He et le $^8$Be devait être en résonance avec un niveau d'énergie inconnu du $^{12}$C. Si ce niveau résonnant existe, alors, la section efficace de la réaction numéro 2 sera considérablement augmentée, la rendant par là même possible. Quelques années après cette prédiction, des mesures en laboratoire ont montré qu'effectivement, cet état excité existait. Ce carbone disponible va pouvoir réagir lui aussi avec les atomes d'hélium présents selon la réaction suivante:
    $_{6}^{12}{\rm {C}}+_{2}^{4}{\rm {He}}\leftrightarrow_{8}^{16}{\rm {O}}+\gamma$
À la fin de la phase de fusion de l'hélium, le cœur de l'étoile est donc composé majoritairement de carbone $^{12}$C et d'oxygène $^{16}$O. Comme pour la fin de la fusion de l'hydrogène, ce qui se passe ensuite ne dépend que de la masse de l'étoile. Si celle-ci est suffisante, alors le cœur de l'étoile peut se contracter de nouveau. La température et la densité peuvent elles aussi augmenter, et atteindre des valeurs qui permettent successivement les fusions du carbone et du néon.

### Fusions du carbone et du néon

La fusion du carbone s'amorce quand la température au cœur de l'étoile dépasse le milliard de kelvins. Trois réactions principales se produisent, pouvant former du sodium, du néon ou du magnésium 23 :
* $_{6}^{12}{\rm {C}}+_{6}^{12}{\rm {C}}\leftrightarrow_{11}^{23}{\rm {Na}}+p$
* $_{6}^{12}{\rm {C}}+_{6}^{12}{\rm {C}}\leftrightarrow_{10}^{20}{\rm {Ne}}+\alpha$
* $_{6}^{12}{\rm {C}}+_{6}^{12}{\rm {C}}\leftrightarrow_{12}^{23}{\rm {Mg}}+n$
Si la température du cœur de l'étoile reste proche du milliard de degré, les deux premières réactions sont favorisées. Si au contraire, elle s'élève au-dessus de 1,1 × 10$^9$ Kelvins, alors, c'est la 3e réaction qui est prédominante. Et de même, quand le carbone vient à manquer, le cœur se contracte, la température s'élève.
Lorsque celle-ci dépasse 1,2 × 10$^9$, les atomes de néon ont suffisamment d'énergie pour que leur fusion ait lieu. Deux réactions principales ont lieu qui peuvent produire de l'oxygène ou du magnésium 24 :
* $_{10}^{20}{\rm {Ne}}+\gamma \leftrightarrow_{8}^{16}{\rm {O}}+\alpha$
* $_{10}^{20}{\rm {Ne}}+\alpha \leftrightarrow_{12}^{24}{\rm {Mg}}+\gamma$
La combustion du carbone et du néon produit aussi des neutrons par les réactions suivantes :
* $_{6}^{13}{\rm {C}}+_{2}^{4}{\rm {He}}\leftrightarrow_{8}^{16}{\rm {O}}+n$
* $_{10}^{22}{\rm {Ne}}+_{2}^{4}{\rm {He}}\leftrightarrow_{12}^{25}{\rm {Mg}}+n$
Lors du processus s les neutrons produits de ces réactions sont absorbés par des noyaux qui ainsi augmentent leur nombre de nucléons. Dans ces noyaux les neutrons excédentaires se transforment en protons en émettant un électron et un neutrino. Le processus s peut produire les noyaux jusqu'au bismuth 209.
La phase suivante fait intervenir cette fois les atomes d'oxygène.

### Fusion de l'oxygène

La dernière contraction du cœur de l'étoile a conduit celui-ci à atteindre une température de plus de 2 milliards de kelvins. À cette température, les atomes d'oxygène fusionnent, et cinq principales réactions ont lieu :
* $_{8}^{16}{\rm {O}}+_{8}^{16}{\rm {O}}\leftrightarrow_{14}^{28}{\rm {Si}}+\alpha +9.594MeV$
* $_{8}^{16}{\rm {O}}+_{8}^{16}{\rm {O}}\leftrightarrow_{15}^{31}{\rm {P}}+p+7.678MeV$
* $_{8}^{16}{\rm {O}}+_{8}^{16}{\rm {O}}\leftrightarrow_{16}^{31}{\rm {S}}+n+1.500MeV$
* $_{8}^{16}{\rm {O}}+_{8}^{16}{\rm {O}}\leftrightarrow_{14}^{30}{\rm {Si}}+2p+0.381MeV$
* $_{8}^{16}{\rm {O}}+_{8}^{16}{\rm {O}}\leftrightarrow_{15}^{30}{\rm {P}}+D+2.409MeV$
Alternativement :
* $_{8}^{16}{\rm {O}}+_{8}^{16}{\rm {O}}\leftrightarrow_{16}^{32}{\rm {S}}+\gamma +xMeV$
* $_{8}^{16}{\rm {O}}+_{8}^{16}{\rm {O}}\leftrightarrow_{12}^{24}{\rm {Mg}}+2\alpha +xMeV$
Les particules alpha, les neutrons et les protons libérés lors de ces réactions ont suffisamment d'énergie pour participer à la synthèse d'autres éléments. On va donc voir apparaitre lors de cette phase un grand nombre d'éléments, tels que le chlore, l'argon, le potassium, le calcium, le titane, etc.
Une fois l'oxygène épuisé, se termine la dernière phase de fusion d'éléments au cœur de l'étoile : la fusion du silicium.

### Fusion du silicium

Nous sommes là dans les tout derniers moments de la vie d'une étoile. Comme on peut le voir au tableau Temps de fusion, l'étoile n'a plus que quelques heures à vivre. Le cœur s'est à nouveau contracté, et cette fois-ci, la température atteint près de 3 milliards de kelvins. Les atomes de silicium sont brisés par les photons gamma présents, et libèrent des neutrons, des protons et des particules $\alpha$ qui vont interagir avec les atomes de $^{28}$Si présents pour former tous les éléments jusqu'au fer.
Lorsque la quantité de fer est trop importante dans le cœur, et la quantité de silicium trop faible pour soutenir la pression de rayonnement produite par la fusion, l'étoile est à la toute fin de sa vie.