First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 1
Explosive nucleosynthesis Margarita Hernanz
Institut d’Estudis Espacials de Catalunya, IEEC-CSIC Barcelona
Themes
Scenarios of explosive nucleosynthesis• Accreting white dwarfs: classical novae, thermonuclear supernovae• Massive stars: core collapse supernovae• Accreting neutron stars: x-ray bursts• (Mergings of neutron stars..)
Synthesis of nuclides up to Fe: supernovae of all types, novaeSynthesis of heavier nuclei: r (and s), p and rp processesObservational clues:
• Chemical evolution of the Galaxy • Individual stars• Galactic radioactivity• Primitive meteorites
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 2
Solar System abundances (by mass number)
α-nuclei
Fe-peak
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 3
Sites of explosive nucleosynthesis
• SUPERNOVAE:* Thermonuclear supernovae (SN Ia): exploding white dwarfs in binary systems (no remnant)* Core collapse supernovae (SN II, SN Ib/c): exploding massive stars (M ≥ 10 M ) (neutron star or black hole remnant)
v ∼ 104 km/s, E ∼ 1051 erg, Mej ∼ M
• CLASSICAL NOVAE:* Explosion of the external H-rich accreted shells of a white dwarf in a binary system
v ∼ 102 - 103 km/s, E ∼ 1045 erg, Mej ∼ 10-4 - 10-5 M
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 4
Effect of accretion: depends on L and M of the WD, accretion rate and chemical comp. of accreted matter properties of the binary system: M1, M2, Porb –A, (dM/dt)accr
explosive H-nuclear burning on top of the WD envelope ejection, Nova explosion: E~1045 erg, Mejected ~10-4-5 M , vejec ~ 102-3 km/s, (H, He, CNO, Ne, ...), ~30/yr in the Galaxy
“central” explosive C burning total disruption of the star, thermonuclear Supernova (Ia): E~1051 erg, Mejected ~MWD (1033 M ) , vejec ~ 104 km/s (Fe, Si, Ca...), ~1/100 yr in the Galaxy
(Can novae accumulate enough mass to finally explode as SNIa?)
White dwarfs in close binary systems
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 5
• X-ray bursts:* Thermonuclear flashes on accreting neutron stars: unstable burning of H and He. Very short rise times and durations (~10 s); recurrence times ~ hrs - days
E ∼ 1039 erg, Mej ∼ 0
* Superbursts: X-ray bursts lasting 1000 times longer and releasing 1000 times more energy (1042 erg); recurrence times uncertain (up to ~5 years). Probably related to C burning
* Some ultracompact LMXBs with high Ne/O ratios (observed in X-rays with ASCA, Chandra...) exhibit X-ray bursts, in apparent contradiction with H and He depletion in donor stars (Juett et al. 2003): spallation of accreted elements (Bildsten et al. 1992)?
Other sites of explosive nucleosynthesis
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 7
- Defining characteristic of SN Ia : lack of H and presence of Si II (λ6355) inspectrum
- Other observational properties:• homogeneity: ∼ 90 % of all SN Ia have similar spectra, light curves and peak
absolute magnitudes• SN Ia appear in both elliptical and spiral galaxies
thermonuclear disruption of mass accreting CO white dwarfs
-Scenario?
Explosive nucleosynthesis in thermonuclear supernovae
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 8
• Single degenerate scenario
WD + Normal companion
(H or He accretion)
• Double degenerate scenario
WD + WD merging
(He or C-O accretion)
Scenario of thermonuclear supernovae
What’s the binary system progenitor of type Ia supernovae?
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 9
- Defining characteristic of SN Ia : lack of H and presence of Si II (λ6355) inspectrum
- Other observational properties:• homogeneity: ∼ 90 % of all SN Ia have similar spectra, light curves and peak
absolute magnitudes• SN Ia appear in both elliptical and spiral galaxies
thermonuclear disruption of mass accreting CO white dwarfs
-Scenario?
- Mass of the white dwarf: Chandrasekhar or sub-Chandrasekhar ?relation with the explosion mechanism
Explosive nucleosynthesis in thermonuclear supernovae
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 10
Thermonuclear supernovae: explosion mechanisms
Detonation: supersonic flameIf ρ > 107 g/cm3 ⇒ C,O → NiIf ρ ≤107 g/cm3 ⇒ C,O → Si,Ca, S, ...
Central carbon ignition : Chandrasekhar mass model
Deflagration: subsonic velocitylaminar flame: v ~ 0.01 cs
The laminar flame becomes turbulent (Rayleigh-Taylor instability)
Flame surface increaseseffective velocity increasesTurbulent flame: v ~ 0.1 - 0.3 cs
Deflagration + detonation:delayed detonation
Courtesy Jordi Isern
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 11
Thermonuclear supernovae: explosion mechanisms
MHe ~ 0.2 - 0.3 M
He-detonation
Pressure wave
MCO ~ 0.5 - 1.1 M
C-detonation
Off-center ignition: sub Chandrasekhar mass model
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 12
Explosive nucleosynthesis in thermonuclear supernovae
Nucleosynthesis in multi-dimensional SNIa explosions
Travaglio et al. 2004, A&A
coupling of a tracer particle method to 2D and 3D Eulerian hydrodynamic SNIa calculations
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 13
Explosive nucleosynthesis in thermonuclear supernovae
Travaglio et al. 2004
Chandrasekhar mass model pure deflagration model in 3D
Main differences w.r.t.
W7:
• 56Fe smaller
• unburnt 12C, 16O, 22Ne
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 14
Important diagnostic of explosion mechanism: observation of gamma-rays
Isotope Decay chain Lifetime Line energy (keV)
56Ni 56Ni 56Co 8.8 d 158, 812, 750, 480
56Co 56Co 56Fe 111 d 847, 1238
57Ni 57Ni 57Co 57Fe (52 h) 390 d 122
44Ti 44Ti 44Sc 44Ca 89 y (5.4 h) 78, 68, 1157
26Al 26Al 26Mg 1.0 x 106 y 1809
60Fe 60Fe 60Co 60Ni 2.0 x 106 y (7.6 y) 1173, 1332
7Be 7Be 7Li 77d 478
22Na 22Na 22Ne 3.8 y 1275
e- capture β+ β-
Sup
erno
v ae
mai
nly
Nov
ae
mai
nly
Radioactive isotopes relevant for γ-ray line astronomy
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 15
20d5 Mpc
SNe Ia: γ-ray spectral evolution Deflagration Delayed detonation Detonation Sub-Chandrasekhar
• DEF only shows the continuum• DEL,DET,SUB display strong lines• 56Ni still present• 56Ni & 56 Co prominent in SUB/DET
Gómez-Gomar, Isern, Jean (1998)
Presence of low Z elements: continuum extends to lower energies in DEF & DEL
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 16
SNe Ia: γ-ray spectral evolution Deflagration Delayed detonation Detonation Sub-Chandrasekhar
60d5Mpc
• 56Ni lines have disappeared• 122 - 136 keV 57Co lines appear• the energy cut-off of DEF is still low
Gómez-Gomar, Isern, Jean (1998)
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 17
Gamma-rays provide diagnostics of the explosion mechanism
( )( )
56
200 max 56(847) / (158)
tot
d
surface
NiF F
Ni∝
F(847)200d/F(158)max
DEF 8 DEL 2.2
SUB 1.3 DET 0.7
Gómez-Gomar, Isern, Jean (1998)
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 18
INTEGRAL/SPI (847 keV line):
• 5 Mpc (DEFa, ΔE=20 keV, 100days)
• 8 Mpc (DET, ΔE=40 keV, 70d)
d=1 Mpc
WARNING: lines are very broad: detectability very difficult
Isern, Bravo, Hirschmann, García-Senz, 2004
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 20
Synthesis of radioactive isotopes in core collapse supernovae
Models* Thielemann, Nomoto and Hashimoto (1996): M = 13 - 25 M , Z
(Nakamura et al. 1999: effect of low Z)- Explosion induced by depositing thermal E in the Fe core “thermal bomb”
* Woosley and Weaver (1995): M = 11 - 40 M , Z = 0, 10-4 - Z- Explosion induced by piston located at the outer edge of the Fe core
* Limongi and Chieffi (2003), Rauscher et al. (2002) ...Some results: “3 coarse classes of nuclei” (Rauscher 2004)• determined by stellar evolution only (mainly hydrostatic); Y=f(M)
sensitive to unceratinties in the reaction rates and to mixing ...He, C, O, Ne, Mg, 26Al, 60Fe (very dependent on 59Fe(n,γ)60Fe)
• depend on stellar evolution and explosion energy; weakly dependent on MSi, S, Ar, Ca
• probe the explosion mechanism: depend on size of pre-SN Fe core, Mcut, Eexplosion, electron abundance Ye (neutronization)44Ti, Fe group (56,57Ni), r-process elements
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 21
56,57Ni and 44Ti production in core collapse supernovae
56,57Ni: explosive Si burning
44Ti: produced during α-rich freezout from NSE in the hottest and deepest layers ejected during the explosion (also in sub-Ch SNIa)
large sensitivity to mass cut location and amount of mass fallback
M(44Ti) < 10-4 M (M ≥ 30 M ==> no 44Ti ejection)44Ti ejection accompanied by 56,57Ni
Diagnostic of models: observation of individual objects• γ-ray lines: 56,57Ni and 44Ti yields (better in SN Ia; detected in SN 1987A)• bolometric light curves
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 22
Diehl & Timmes, 1998
56,57Ni and 44Tiproduction sites
Explosive burning inmassive stars (core collapse supernovae56,57Ni and 44Ti areproduced in the same internal zones
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 23
56,57Ni and 44Ti production in core collapse supernovae
56,57Ni: explosive Si burning
44Ti: produced during α-rich freezout from NSE in the hottest and deepest layers ejected during the explosion
large sensitivity to mass cut location and amount of mass fallback
M(44Ti) < 10-4 M (M ≥ 30 M ==> no 44Ti ejection)44Ti ejection accompanied by 56,57Ni
Diagnostic of models: observation of individual objects• γ-ray lines: 56,57Ni (better in SN Ia; but detected in SN 1987A), and 44Ti yields • bolometric light curves: idem (τ: 111d, 390d, 89 yrs)
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 25
Consequences of nuclear electron captures in core collapse supernovae
Realistic treatment of e-capture on heavy nuclei (A>60) implies significant changes in the hydrodynamics of core collapse and bounce, and also on postbounce evolution
see for instance Hix et al. 2003, Phys. Rev Let. (ask Martinez Pinedo)
Still fail to produce an explosion in the spherically symmetric case, but mass behind the shock when it is launched is reduced by 20%, neutrino L is boosted by 15% ...
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 27
Mass transfer from thecompanion star onto
the white dwarf (cataclysmic variable)
Hydrogen burning indegenerate conditions on
top of the white dwarf
Thermonuclear runaway
Explosive H-burning
Decay of short-lived radioactive nuclei in the outer envelope(transported by convection)
Envelope expansion, L increaseand mass ejection
Scenario of classical novae
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 28
Nova Models: Thermonuclear Burning of Hydrogen. CNO cycle
13N
15O 16O
14N
12C 13C
15N
17O
17F
AX(p,γ)
(p,α)(β+)
14O
18Ne
18F 19F
18O102s 176s
93s
862s
158min
• Start: τβ+ < τ(p,γ)
CNO cycle operates in equ.• T ~108 K: τβ+ > τ(p,γ)
CNO cycle β+-limited(bottle neck)
• Convection: fresh fuel brought to the burning
shellτconv < τβ+: β+-unstable nuclei to
external cooler regions where they are preserved from destructionLater decay on the surface leads to expansion and luminosity increase
~
Z
N
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 29Gehrz et al 1998, PASP
V1370 Aql 1982 Z=0.86=45 Z ; Ne=0.56=296 Ne
QU Vul 1984 Z=0.44=23 Z ; Ne=0.18=100 Ne
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 30
Nova Models: need of core-envelope mixing
• Z observed >> solar mixing CO or ONe core – solar envelope accreted
• Explosion itself (fast nova) initial overabundance of CNO mixing
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 31
Many classical nova ejecta are enriched in CNO and Ne. Rosner and coworkers recently suggested that the enrichment might originate in the resonant interaction between large-scale shear flows in the accreted H/He envelope and gravity waves at the interface between the envelope and the underlying C/O white dwarf . The shear flow amplifies the waves, which eventually form cusps and break. This wave breaking injects a spray of C/O into the superincumbent H/He.
In the absence of enrichment prior to ignition, the base of the convective zone, does not reach the C/O interface. As a result, there is no additional mixing, and the runaway is slow. In contrast, the formation of a mixed layer during the accretion of H/He, prior to ignition, causes a more violent runaway. The envelope can be enriched by ≤ 25% of C/O by mass (consistent with that observed in some ejecta) for shear velocities, over the surface, with Mach numbers ≤ 0.4.
ApJ, 602 (2004)
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 34
The underlying white dwarf
White dwarfs are the endpoints of the stellar evolution of stars with masses below 11-12 M .
M≤ 8-10 M CO white dwarfs (He burning)
8-10 M ≤M ≤ 12 M ONe white dwarfs (C burning)
10 M 1.2 M ONe core
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 35
The underlying white dwarf
Ritossa, García-Berro & Iben, 1996, ApJ
see also Domínguez, Tornambè & Isern 1993
10M mass Population I star evolved from the H-burning main sequence through carbon burning
1.2M ONe core
≠ONeMg core predicted by hydrostatic C-burning (Arnett & Truran, 1969)
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 36
The underlying white dwarf
Size of the CO core at the beginning of C burning, for single and binary evolution
Mass point at which C is ignited
Minimum mass required for C-ignition to take place (*): 8.1 M (single) and 8.7 M (binary)
Off-center C-ignition
Central C ignition:
11 M for single evolution
12 M for binary evolution
Be careful when adopting M>8Mas “massive stars”
**
Gil Pons, García-Berro, José, Hernanz & Truran, 2003, A&A
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 37
The underlying white dwarf
Gil-Pons, García-Berro, José, Hernanz, Truran, 2003, A&A
Size of the final core for single and binary evolution: relevance of new Minitial-Mfinalmass relation for the fraction of novae hosting ONe white dwarfs: smaller number but still around 30%
ONe core mass with a “CO buffer” (binary evolution)
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 38
The underlying White Dwarf
MZMAS=9M MZMAS=10M
“CO buffer” on top of an ONe core
(Gil-Pons et al., 2003, A&A)
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 39
The underlying White Dwarf
José, Hernanz, García-Berro, Gil-Pons, 2003, ApJL
A
CO buffer on top of ONe core: weird nuclesoynthesis potentially leading to missclassification of novae
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 40
The underlying White Dwarf
Relevance of CO buffer on top of ONe WD for nova nucleosynthesis: lack of Ne in the ejecta: misclassification of novae (non-Ne nova ≠ CO nova)
José, Hernanz, García-Berro, Gil-Pons, 2003, ApJL
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 41
Relevance of nucleosynthesis in classical novae
interpretation of elemental abundances observed in individual objectschemical evolution of the Galaxypresolar meteoritic grains dust formation in novae (IR obs.)gamma-ray emission
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 43
Nova nucleosynthesis and chemical evolution of the Galaxy
Mejec(theor.) ~ 2x10-5 M /nova
R(novae) ~ 35 novae/yr
Age of the Galaxy ~ 1010 yrs
Mejec,total(novae) ~ 7x106 M = (7x10-4 M /yr) ≈ 1/3000 Mgal(gas+dust)
Novae can account for the galactic abundances of the isotopes they overproduce (w.r.t. sun) by factors ≥ 3000
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 44
Novae nucleosynthesis: overproductions w.r.t. solar
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 45
Novae nucleosynthesis: overproductions w.r.t. solar
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 46
ONe
1.35M
Novae nucleosynthesis: overproductions w.r.t. solar
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 47
MNRAS, 2004
7Li CNO: 13C, 15N, 17O
See as well Alibés, Labay & Canal, 2001, A&A
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 48
IR observations indicate that dust grains are formed in many novae
Dust in novae
Nova Cyg 1978
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 49
Novae and presolar meteoritic grainsPrimitive meteorites contain presolar grains, which condensed
in stellar atmospheres or in supernova or nova ejecta, and survived their “interstellar trip” and solar system formation
Isotopic abundances measurements in lab
allow to ascertain their origin
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 50
Five SiC and two graphite grains from the Murchison and Acfer 094 meteorites show isotopic compositions indicating a nova origin: Amari, Gao, Nittler, Zinner, José, Hernanz & Lewis (2001); José, Hernanz, Amari, Lodders & Zinner (2004)
Novae and presolar meteoritic grains
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 51
Novae and presolar meteoritic grains
Amari et al. 2001
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 52
Main radioactive isotopes synthesized in classical novae
Nucleus τ Type of emission Nova type
13N 862 s511 keV linecontinuum (E<511 keV)
CO and ONe
18F 158 min511 keV linecontinuum (E<511 keV)
CO and ONe
7Be 77 days 478 keV line CO mainly
22Na 3.75 yr 1275 keV line ONe
26Al 1.0X106 yr 1809 keV line ONe
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 53
Spectra of CO and ONe novaeMWD = 1.15 M
d=1 kpc
Gómez-Gomar, Hernanz, José, Isern,1998, MNRAS
d=1 kpc
coONe
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 54
Light curves: 1275 keV (22Na) line
Rise phase
Exponentialdecline
tmax: 20 days (1.15M ),12 days (1.25 M ), line width ∼ 20 keV; duration: months
Flux (max) ∼ 2x10-5 ph/cm2/s; Mejected(22Na) ~ (6-7)x10-9M
Only in ONe
novae
d=1 kpc
predicted theoretically by Clayton & Hoyle, 1974
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 55
Observations : 1275 keV line (22Na)
CGRO/COMPTEL searched for 1275 keV emission in many novae: no detection, upper limits
CGRO/COMPTEL most constraining upper limit (Nova Cyg 1992, d=2.3 kpc) in agreement with current theoretical predictions:
F<2.3x10-5phot/cm2/s Mej(22Na)< 3.0x10-8M
Iyudin et al. 1995, A&A
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 56
Nuclear uncertainties: 1275 keV (22Na) line
Uncertainties in the rates 21Na(p,γ)22Mg and 22Na(p,γ)23Mgtranslated into uncertainties by factors around 3 in the 22Na yields (José, Coc & Hernanz, 1999) – and therefore in the flux of the 1275 keV line.
These uncertainties have been reduced by recent measurements:
Jenkins et al 2004: 22Na(p,γ)23Mg
Bishop et al. 2003, D’Auria et al. 2004: 21Na(p,γ)22Mg
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 57
Galactic distribution of γ-ray emission from novae
Observations (upper limits): Leising et al. 1988, Harris et al 1991, 1996
Theoretical predictions:Jean, Hernanz, Gómez-Gomar, José, 2000, MNRAS
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 58
Light curves: 478 keV (7Be) line
Only in CO novae
tmax: 13 days (0.8M )
5 days (1.15 M )
duration: some weeks
Flux ∼ (1-2)x10-6 ph/cm2/s
Mejected(7Be) ~ (0.7-1.1)x10-10M
Line width: 3-7 keV
d=1 kpc
predicted theoretically by Clayton 1981
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 59
Light curves: 511 keV line
• 511 keV line in ONe novae remains after 2 days until ∼1 week because of e+ from22Na• Intense (but short duration)• Very early appearence, before visual maximum (i.e, before discovery)
Model tmax* (h) Fmax (ph/cm2/s)**
CO, 0.8 M - - - 2.6 x 10-5
CO, 1.15 M 6.5 5.3 x 10-4
ONe, 1.15 M 6 1.0 x 10-3
ONe, 1.25 M 5 1.9 x 10-3
In CO and ONe novae
d=1 kpc
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 60
Nuclear uncertainties related with 18F synthesis (511 keV & continuum emission)
Coc, Hernanz, José, Thibaud, 2000, A&A
Rates obtained including the latest experimental data up to the end of 1999
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 61
18F(p,α)15O• de Séréville, Coc, Angulo et al. 2003,Phys. Rev. C: reduction of the uncertainty and nominal rate similar to Coc et al. 2000• Other measurements: Bardayan et al.; Kozub et al.: rates 3-5 times lower than expected
17O(p,γ)18F - 17O(p,α)14N• Fox et al. 2004: uncertainties reduced• Orsay results: less 17O, less 18F
Nuclear uncertainties related with 18F synthesis (511 keV & continuum emission)
Uncertainties in the 18F yields reduced by a factor of ~5 (18F+p) and by
a factor of ~3 (17O+p) in the nova T range
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 62
End point of nova nucleosynthesis: 40Ca
First CARINA Workshop, 8-10 June 2005, Aiguablava M. Hernanz 63
T=6x108 K
T=3x108 K T=109 K
T=1.7x109 K
José & Moreno (2002)
rp process on X-ray bursts