TYCHO Users Manual

TYCHO Users Manual

Macroscopic Physics basis

Microscopic Physics basis

Abundances are user defined. Tables for Anders-Grevesse, Lodders, Asplund, Big-Bang and interpolated metal poor abundances (and alpha-rich) are standard.

TYCHO Data

TYCHO produces data is several ways:
  • An online log in a terminal window
  • Online graphics in 3 PGPLOT windows
  • Regular "dumps" of snapshots (which are re-startable)
  • Time sequences for position in the HR diagram (hr.??) and for convective status (cv.??)
    • Here "??" represents a two-character symbol identifying the computational sequence, and is user defined in the control file "params.d".
    Control file: params.d

    Programs for Model modification

    OMIT: This is an auxillary program to "slice" layers off the grid (top down, or bottom up) to prepare a more efficient computational domain for explicit hydrodynamics. Generally this will be the removal of a core region, to examine hydrodynamics of the mantle, or visa versa.

    Control file: omit.in

    GENEX: This auxillary program takes a model (in file "old.model") and maps it into "new.model", with optional changes:

  • mass and radius scaling
  • several frequently used abundance redefinitions
    • Control file: genex.in

    Mapping hydrostatic --> explicit hydrodynamic models

    Analysis Suite

    In additon to the basic hydrodynamics/evolution code TYCHO, a suite of analysis programs are included. These use PGPLOT graphics and allow rapid investigation of both snapshots and time sequences.

    snapshots

  • GENPLOT (workhorse graphics of snapshots)
  • GENNUC (rund of abundances in snapshots)
  • GENRATE (details of reactions in individual zones in snapshots)

    time sequences

  • HRPLOT (HR diagram)
  • CVPLOT (evolution of convective structure)
  • RATIOS (evolution of abundance pairs versus Teff, time, or other variable)

    The control files for all these programs are named according to the program, plus a postfix: ".in"

    A 3D linear stability tool is included, and a nonlinear stability tool is being developed.

    MODIFIED UP TO HERE-------------------------------------

    THIS PAGE IS WORK IN PROGRESS (4/10/2008)
    DO NOT USE ITEMS BELOW HERE

    Getting started: the first run

    Get the OPAL2 opacity tarfile here

    Un-tar this file and move it to your tycho-6.0 directory. Its use is explained in opal2.dir/type2tab/README. It allows you to replace the default opal tables codataa, codatab, codatac, codatad, and codatae with tables for different metallicity.

    Get the genex tarfile here

    Un-tar this file and move it to your tycho-6.0 directory. Type "make genex" in its genex.dir/src subdirectory.

    The use of genex is explained in genex.dir/README. It allows you to modify an imodel-type file by scaling the mass and radius, or changing the composition, and provides the result in a file called new.model.

    To scale the metallicity to half of solar, set ifm = 2 and zpop = 0.5 in genex.dir/genex.in, for example.

    To increase the mass and radius by 10 percent (at roughly constant temperature), set ifm = 1, fm=1.1 and fr=1.1.

    Satisfactory results may required repeated small modifications, followed by a few timesteps with tycho6.

    Get the burn.dir tarfile here

    Un-tar this file and move it to your tycho-6.0 directory. Type "make Burnf" in its burn.dir/src subdirectory. The use of Burnf will be discussed in class. There will now be a Burnf executeable in /burn.dir. Its parameter file in called nuc.in. Type Burnf and see how many burning stages you can identify. Hint: look for the temperatures and compare with my book. Expect several in a few seconds!

    A more detailed discussion of running Burnf is here.

    Get asci file nuc.in.bb2 here

    Move this file to your tycho-6.0/burn.dir directory. If you have already made Burnf, you can do a cosmological nucleosynthesis run: 

    cp nuc.in.bb2 nuc.in
Burnf
    After about 13 seconds (2.77e5 seconds for the cosmology), you should have an output file "nuc.out.d" which gives details. You may want to modify fnuclear.f in /burn.dir/src to give output matched to your graphics package. See Chapter 5 of Supernovae and Nucleosynthesis, for a discussion of what you have just calculated. You might try modifying the initial density in line 21 in nuc.in 
    Intial density (overwritten for ievol=2)..........  rho      1.5000E-01
    to generate a version of Figure 5.4.

    Get Audra's pgburn.dir tarfile here

    In addition to pgplot graphics, the /pgburn.dir/src/fnuclear.f file has a pause inserted at line 1096. Burnf will proceed for nwrit steps (see nuc.in), bring up three graphics windows, and pause. If you type "go", it will continue. Otherwise it stops. This allows the windows to appear for enough time to study them. You can step through all the nucleosynthesis stages this way. You can comment out the pause to avoid this feature entirely.

    Read Audra's manual here

    Get tycho-6.11 tarfile here

    Version tycho-6.11 has an 80 nuclei network, Alexander surface opacities, and several fixes. This directory is set up for population III (z=0 metallicity) stars; 4 and 23 solar mass sequences have been run and are included.

    Running TYCHO

    Analysis

    Return to main TYCHO page.

    E-mail:darnett@as.arizona.edu

    TYCHO is not public domain software, but is freely available for non-commercial purposes. Absolutely no guarantee is given; use at your own risk...

    Last updated 2-20-02