Library of EZ Star Plots

ZAMS for Z=0.02 to Z=0.0001

Plots for Star Masses from 0.6 Msun to 100 Msun
Tracks for Different Mass Ranges at Z=0.02

Comparison of Tracks at Various Z's

Some Stars at Z's from 0.0001 to 0.03

Details for Lots of Individual Stars at Z=0.02
Helium Flash Delayed by Mass Loss at Tip of Giant Branch

     The plots here primarily come from the demo programs I created to give EZ a test drive.  The programs, which can be downloaded as part of the EZ software distribution (available at my home page), illustrate the use of EZ and can provide templates for your own experiments.   They use the standard interfaces and were not (too) difficult to write.   I hope they'll show that the system gives access to some interesting phenomena.

    The links below let you download pdf files, with each file containing several figures on separate pages.   Rather than have a lot of jpegs floating around, I decided just to show some example thumbnails here and provide the pdf files so you can view or print as you desire (the Acrobat pdf reader is available from Adobe at no charge).



    Demo 1 computes the zero age main sequences for Z=0.02, 0.001, 0.004, and 0.0001.   The variation with mass along the ZAMS is shown for luminosity, surface temperature, central temperature, central density, radius, P-P versus CNO nuclear reactions, opacity, nuclear burning time scale, and convection zones.
ZAMS for 0.1 Msun to 100 Msun, Z=0.02 and Z=0.0001
Click to download pdf file.







    Demo 2 begins with plots of the Hertzsprung-Russell diagram tracks and central temperature and density tracks for a sample of Z=0.02 stars over a wide range of initial mass.   The tracks are marked to show the location of break-even for the net power from nuclear reactions beyond hydrogen burning minus the total power lost to neutrino cooling (i.e., the onset of significant helium burning).   The temperature-density plots include lines showing the degeneracy contours for Psi equal 0, 5, and 100.   NOTE: For the massive stars, late stages of evolution depend critically on the details of wind-driven mass loss.   I've made no attempt to address this issue and simply use a Reimers' wind with eta=1 as a zeroth order approximation.

H-R Tracks and Central T-Rho Tracks for 0.3 Msun to 100 Msun
Click to download pdf file.




    Demo 2 also also looks in detail at a number of stars with masses from 0.3 Msun up to 100 Msun.   For each one, it traces the evolution of convection zones, burning zones, central abundances, power sources, neutrino cooling, and shows profiles taken at key moments along the way.   These are all simulated with convective overshooting and a Reimers' wind (eta=1).   The tests start at the Z=0.02 ZAMS and are terminated when either the core becomes highly degenerate or something such as a helium flash requires a timestep smaller than the dynamic timescale.

    In addition to the usual HR diagram and central temperature-density tracks, the figures include histories of the later values for radius, neutino losses, power from triple-alpha, power from alpha-capture, power from carbon burning, center degeneracy, and total metal fraction.   There are plots that trace the evolution of convection zones, burning zones, central abundances (of helium, carbon, oxygen, nitrogen, and neon), nuclear power sources, and neutrino cooling.   The burning zones are subdivided into regions producing at least 1 erg/g/sec.   The power is subdivided into sources: PP+CNO, triple-alpha, alpha-capture, other nuclear burning, and total neutrino cooling. The neutrino cooling is subdivided into its various sources: plasmon decay, bremsstrahlung, photoneutrinos, and pair production.   Finally, there are figures showing sets of profiles by mass coordinate taken at key moments along the entire evolution.   These have profiles for abundances, power sources, temperature, and density.

    Plots combining various Z's ­ It can also be instructive to look at a single plot combining the tracks for several metallicities at a given initial masss.   The following files show such combined H-R and T-Rho tracks for a selection of masses.
Click to download pdf file.
  1  
  2  
  3  
  4  
  6  
  12  
  25  

     The results are illustrated here by the H-R and T-Rho tracks for stars of 1 Msun and 6 Msun.




    Different Z's   ­   The initial ZAMS metallicity has a big effect on the evolution.   Here is a selection of several cases for comparison.
Click to download pdf file.
  Z=0.03  
  1  
  2  
  3  
  4  
  6  
  12  
  25  
  Z=0.02  
  1  
  2  
  3  
  4  
  6  
  12  
  25  
  Z=0.01  
  1  
  2  
  3  
  4  
  6  
  12  
  25  
  Z=0.004  
  1  
  2  
  3  
  4  
  6  
  12  
  25  
  Z=0.001  
  1  
  2  
  3  
  4  
  6  
  12  
  25  
  Z=0.0001  
  1  
  2  
  3  
  4  
  6  
  12  
  25  
     For example, the following images show the H-R and T-Rho tracks for a 6 Msun star with Z's of 0.0001 and 0.03. The files have the complete sets of plots.




    Individual Stars ­  The wide range of initial masses provide illustrations of white dwarf formation (up to 0.7 Msun), helium core flash (up to 2 Msun), and helium core burning occurring under a hydrogen burning shell.   For the more massive stars (over 2 Msun), the period of helium burning in the core is followed by a stage of two burning shells, one consuming helium and the other hydrogen, over a growing carbon/oxygen core.   NOTE: As mentioned above, for the massive stars, late stages of evolution depend critically on the details of wind-driven mass loss.   I've made no attempt to address this issue and simply use a Reimers' wind with eta=1 as a zeroth order approximation.
Click to download pdf file for the indicated initial mass.
  0.6  
  0.7  
  0.8  
  0.9  
  1.0  
  1.1  
  1.2  
  1.3  
  1.4  
  1.5  
  1.6  
  1.7  
  1.8  
  1.9  
  2.0  
  2.1  
  2.2  
  2.3  
  2.4  
  2.5  
  2.6  
  2.7  
  2.8  
  2.9  
  3.0  
  3.1  
  3.2  
  3.3  
  3.4  
  3.5  
  3.6  
  3.7  
  3.8  
  3.9  
  4  
  5  
  6  
  7  
  8  
  9  
  10  
  12  
  16  
  25  
  40  
  100  


     The 0.8 to 2.0 Msun cases include information for the evolution after the helium flash.   Since the code doesn't actually follow the star through the details of the flash, this is only meant to be suggestive.   The post-flash evolution is done with a star that matches the pre-flash star in total mass, helium core mass, and metallicity.   It has recovered from the trauma of the flash and is happily burning helium in its core.



    Late Helium Flash ­ This example considers a 1 Msun star with Z=0.02.   It takes the star from the main sequence and saves state near the tip of the giant branch.   Then it continues the evolution in three different ways.   First, it evolves with no mass loss up to the helium flash.   Then, it restores the saved state and evolves the star again with the envelope ejected quickly.   This gives a smaller core and produces a white dwarf.   Finally, it restores state and tries again with a reduced rate of mass loss, which gives a slightly larger core and a late helium flash.
1.0 Msun with Late Helium Flash after Envelope Ejection
Click to download pdf file.