Code owners
Assign users and groups as approvers for specific file changes. Learn more.
binary_c_parameters.rst 161.11 KiB
Binary_c parameters
The following chapter contains all the parameters that the current version of binary_c can handle, along with their descriptions and other properties.
- This information was obtained by the following binary_c build:
- binary_c git branch: branch_david binary_c git revision: 5834:20211216:a3146490e Built on: Dec 20 2021 20:33:01
Section: stars
Parameter: metallicity
Description: This sets the metallicity of the stars, i.e. the amount (by mass) of matter which is not hydrogen or helium. If you are using the BSE algorithm, this must be 1e-4 <= metallicity <= 0.03. See also nucsyn_metallicity and effective_metallicity.
Parameter input type: Float
Default value: 0.02
Parameter: effective_metallicity
Description: This sets effective metallicity of stars as used in routines like the Schneider wind loss. If not set, or set to DEFAULT_TO_METALLICITY (==-1, the default), this is just the same as metallicity. The main difference between effective_metallicity and metallicity is the range of validity: 0 <= effective_metallicity <= 1, while metallicity's range of validity is limited by the stellar evolution algorithm (so, for BSE, is 1e-4 <= metallicity <= 0.03).
Parameter input type: Float
Default value: -1
Macros: ['DEFAULT_TO_METALLICITY = -1']
Parameter: M_1
Description: The initial mass of star one (in solar units, internally this is star index 0).
Parameter input type: Float
Default value: 0
Parameter: M_2
Description: The initial mass of star two (in solar units, internally this is star index 1).
Parameter input type: Float
Default value: 0
Parameter: M_3
Description: The initial mass of star three (in solar units, internally this is star index 2).
Parameter input type: Float
Default value: 0
Parameter: M_4
Description: The initial mass of star four (in solar units, internally this is star index 3).
Parameter input type: Float
Default value: 0
Parameter: vrot1
Description: The initial equatorial rotational velocity of star one (in km/s, internally this is star index 0). If 0.0, the Hurley et al 2000/2002 prescription is used to set the main-sequence velocity, so for a truly non-rotating star, set to something small (e.g. 0.001). See also vrot2,3,4.
Parameter input type: Float
Default value: 0
Macros: ['VROT_BSE = 0', 'VROT_BREAKUP = -1', 'VROT_SYNC = -2', 'VROT_NON_ROTATING = -3']
Parameter: vrot2
Description: The initial equatorial rotational velocity of star two (in km/s, internally this is star index 1). If 0.0, the Hurley et al 2000/2002 prescription is used to set the main-sequence velocity, so for a truly non-rotating star, set to something small (e.g. 0.001). See also vrot1,3,4.
Parameter input type: Float
Default value: 0
Macros: ['VROT_BSE = 0', 'VROT_BREAKUP = -1', 'VROT_SYNC = -2', 'VROT_NON_ROTATING = -3']
Parameter: vrot3
Description: The initial equatorial rotational velocity of star three (in km/s, internally this is star index 2). If 0.0, the Hurley et al 2000/2002 prescription is used to set the main-sequence velocity, so for a truly non-rotating star, set to something small (e.g. 0.001). See also vrot1,2,4.
Parameter input type: Float
Default value: 0
Macros: ['VROT_BSE = 0', 'VROT_BREAKUP = -1', 'VROT_SYNC = -2', 'VROT_NON_ROTATING = -3']
Parameter: vrot4
Description: The initial equatorial rotational velocity of star four (in km/s, internally this is star index 3). If 0.0, the Hurley et al 2000/2002 prescription is used to set the main-sequence velocity, so for a truly non-rotating star, set to something small (e.g. 0.001). See also vrot1,2,3.
Parameter input type: Float
Default value: 0
Macros: ['VROT_BSE = 0', 'VROT_BREAKUP = -1', 'VROT_SYNC = -2', 'VROT_NON_ROTATING = -3']
Parameter: Prot1
Description: The initial equatorial rotational velocity of star one (in km/s, internally this is star index 0). See also Prot2,3,4.
Parameter input type: Float
Default value: 0
Parameter: Prot2
Description: The initial equatorial rotational velocity of star two (in km/s, internally this is star index 1). See also Prot1,3,4.
Parameter input type: Float
Default value: 0
Parameter: Prot3
Description: The initial equatorial rotational period of star three (in days, internally this is star index 2). See also Prot1,2,4.
Parameter input type: Float
Default value: 0
Parameter: Prot4
Description: The initial equatorial rotational period of star four (in days, internally this is star index 3). See also Prot1,2,3.
Parameter input type: Float
Default value: 0
Parameter: inclination1
Description: The initial inclination of star one (in degrees).
Parameter input type: Float
Default value: 0
Parameter: inclination2
Description: The initial inclination of star two (in degrees).
Parameter input type: Float
Default value: 0
Parameter: inclination3
Description: The initial inclination of star three (in degrees).
Parameter input type: Float
Default value: 0
Parameter: inclination4
Description: The initial inclination of star four (in degrees).
Parameter input type: Float
Default value: 0
Parameter: B_1
Description: The initial magnetic field of star one (in Gauss, internally this is star index 0).
Parameter input type: Float
Default value: 0
Parameter: B_2
Description: The initial magnetic field of star two (in Gauss, internally this is star index 1).
Parameter input type: Float
Default value: 0
Parameter: B_3
Description: The initial magnetic field of star three (in Gauss, internally this is star index 2).
Parameter input type: Float
Default value: 0
Parameter: B_4
Description: The initial magnetic field of star four (in Gauss, internally this is star index 3).
Parameter input type: Float
Default value: 0
Parameter: B_inclination1
Description: The initial inclination of the magnetic field of star one (in degrees).
Parameter input type: Float
Default value: 0
Parameter: B_inclination2
Description: The initial inclination of the magnetic field of star two (in degrees).
Parameter input type: Float
Default value: 0
Parameter: B_inclination3
Description: The initial inclination of the magnetic field of star three (in degrees).
Parameter input type: Float
Default value: 0
Parameter: B_inclination4
Description: The initial inclination of the magnetic field of star four (in degrees).
Parameter input type: Float
Default value: 0
Parameter: stellar_type_1
Description: Set the stellar type of star 1 (internal index 0), usually MAIN_SEQUENCE (main sequence). Note that setting the stellar type only works for stars with both age=0 and core_mass=0, i.e. main sequence (hydrogen or helium), white dwarfs, black holes and neutrn stars.
Parameter input type: Integer
Default value: 0
Macros: ['LOW_MASS_MS = 0', 'MS = 1', 'HG = 2', 'GIANT_BRANCH = 3', 'CHeB = 4', 'EAGB = 5', 'TPAGB = 6', 'HeMS = 7', 'HeHG = 8', 'HeGB = 9', 'HeWD = 10', 'COWD = 11', 'ONeWD = 12', 'NS = 13', 'BH = 14', 'MASSLESS_REMNANT = 15', 'LOW_MASS_MAIN_SEQUENCE = 0', 'MAIN_SEQUENCE = 1', 'HERTZSPRUNG_GAP = 2', 'FIRST_GIANT_BRANCH = 3', 'CORE_HELIUM_BURNING = 4', 'EARLY_ASYMPTOTIC_GIANT_BRANCH = 5', 'THERMALLY_PULSING_ASYMPTOTIC_GIANT_BRANCH = 6', 'NAKED_MAIN_SEQUENCE_HELIUM_STAR = 7', 'NAKED_HELIUM_STAR_HERTZSPRUNG_GAP = 8', 'NAKED_HELIUM_STAR_GIANT_BRANCH = 9', 'HELIUM_WHITE_DWARF = 10', 'CARBON_OXYGEN_WHITE_DWARF = 11', 'OXYGEN_NEON_WHITE_DWARF = 12', 'NEUTRON_STAR = 13', 'BLACK_HOLE = 14', 'STAR_WITH_NO_MASS = 15']
Parameter: stellar_type_2
Description: Set the stellar type of star 2 (internal index 1), usually MAIN_SEQUENCE (main sequence). Note that setting the stellar type only works for stars with both age=0 and core_mass=0, i.e. main sequence (hydrogen or helium), white dwarfs, black holes and neutrn stars.
Parameter input type: Integer
Default value: 0
Macros: ['LOW_MASS_MS = 0', 'MS = 1', 'HG = 2', 'GIANT_BRANCH = 3', 'CHeB = 4', 'EAGB = 5', 'TPAGB = 6', 'HeMS = 7', 'HeHG = 8', 'HeGB = 9', 'HeWD = 10', 'COWD = 11', 'ONeWD = 12', 'NS = 13', 'BH = 14', 'MASSLESS_REMNANT = 15', 'LOW_MASS_MAIN_SEQUENCE = 0', 'MAIN_SEQUENCE = 1', 'HERTZSPRUNG_GAP = 2', 'FIRST_GIANT_BRANCH = 3', 'CORE_HELIUM_BURNING = 4', 'EARLY_ASYMPTOTIC_GIANT_BRANCH = 5', 'THERMALLY_PULSING_ASYMPTOTIC_GIANT_BRANCH = 6', 'NAKED_MAIN_SEQUENCE_HELIUM_STAR = 7', 'NAKED_HELIUM_STAR_HERTZSPRUNG_GAP = 8', 'NAKED_HELIUM_STAR_GIANT_BRANCH = 9', 'HELIUM_WHITE_DWARF = 10', 'CARBON_OXYGEN_WHITE_DWARF = 11', 'OXYGEN_NEON_WHITE_DWARF = 12', 'NEUTRON_STAR = 13', 'BLACK_HOLE = 14', 'STAR_WITH_NO_MASS = 15']
Parameter: stellar_type_3
Description: Set the stellar type of star 3 (internal index 2), usually MAIN_SEQUENCE (main sequence). Note that setting the stellar type only works for stars with both age=0 and core_mass=0, i.e. main sequence (hydrogen or helium), white dwarfs, black holes and neutrn stars.
Parameter input type: Integer
Default value: 0
Macros: ['LOW_MASS_MS = 0', 'MS = 1', 'HG = 2', 'GIANT_BRANCH = 3', 'CHeB = 4', 'EAGB = 5', 'TPAGB = 6', 'HeMS = 7', 'HeHG = 8', 'HeGB = 9', 'HeWD = 10', 'COWD = 11', 'ONeWD = 12', 'NS = 13', 'BH = 14', 'MASSLESS_REMNANT = 15', 'LOW_MASS_MAIN_SEQUENCE = 0', 'MAIN_SEQUENCE = 1', 'HERTZSPRUNG_GAP = 2', 'FIRST_GIANT_BRANCH = 3', 'CORE_HELIUM_BURNING = 4', 'EARLY_ASYMPTOTIC_GIANT_BRANCH = 5', 'THERMALLY_PULSING_ASYMPTOTIC_GIANT_BRANCH = 6', 'NAKED_MAIN_SEQUENCE_HELIUM_STAR = 7', 'NAKED_HELIUM_STAR_HERTZSPRUNG_GAP = 8', 'NAKED_HELIUM_STAR_GIANT_BRANCH = 9', 'HELIUM_WHITE_DWARF = 10', 'CARBON_OXYGEN_WHITE_DWARF = 11', 'OXYGEN_NEON_WHITE_DWARF = 12', 'NEUTRON_STAR = 13', 'BLACK_HOLE = 14', 'STAR_WITH_NO_MASS = 15']
Parameter: stellar_type_4
Description: Set the stellar type of star 4 (internal index 3), usually MAIN_SEQUENCE (main sequence). Note that setting the stellar type only works for stars with both age=0 and core_mass=0, i.e. main sequence (hydrogen or helium), white dwarfs, black holes and neutrn stars.
Parameter input type: Integer
Default value: 0
Macros: ['LOW_MASS_MS = 0', 'MS = 1', 'HG = 2', 'GIANT_BRANCH = 3', 'CHeB = 4', 'EAGB = 5', 'TPAGB = 6', 'HeMS = 7', 'HeHG = 8', 'HeGB = 9', 'HeWD = 10', 'COWD = 11', 'ONeWD = 12', 'NS = 13', 'BH = 14', 'MASSLESS_REMNANT = 15', 'LOW_MASS_MAIN_SEQUENCE = 0', 'MAIN_SEQUENCE = 1', 'HERTZSPRUNG_GAP = 2', 'FIRST_GIANT_BRANCH = 3', 'CORE_HELIUM_BURNING = 4', 'EARLY_ASYMPTOTIC_GIANT_BRANCH = 5', 'THERMALLY_PULSING_ASYMPTOTIC_GIANT_BRANCH = 6', 'NAKED_MAIN_SEQUENCE_HELIUM_STAR = 7', 'NAKED_HELIUM_STAR_HERTZSPRUNG_GAP = 8', 'NAKED_HELIUM_STAR_GIANT_BRANCH = 9', 'HELIUM_WHITE_DWARF = 10', 'CARBON_OXYGEN_WHITE_DWARF = 11', 'OXYGEN_NEON_WHITE_DWARF = 12', 'NEUTRON_STAR = 13', 'BLACK_HOLE = 14', 'STAR_WITH_NO_MASS = 15']
Parameter: max_stellar_type_1
Description: The maximum stellar type of star 1 (internal index 0). Evolution is stopped when the star reaches this stellar type. If this is negative, massless remnants are allowed, and the maximum stellar type is the absolute value.
Parameter input type: Integer
Default value: 16
Macros: ['LOW_MASS_MS = 0', 'MS = 1', 'HG = 2', 'GIANT_BRANCH = 3', 'CHeB = 4', 'EAGB = 5', 'TPAGB = 6', 'HeMS = 7', 'HeHG = 8', 'HeGB = 9', 'HeWD = 10', 'COWD = 11', 'ONeWD = 12', 'NS = 13', 'BH = 14', 'MASSLESS_REMNANT = 15', 'LOW_MASS_MAIN_SEQUENCE = 0', 'MAIN_SEQUENCE = 1', 'HERTZSPRUNG_GAP = 2', 'FIRST_GIANT_BRANCH = 3', 'CORE_HELIUM_BURNING = 4', 'EARLY_ASYMPTOTIC_GIANT_BRANCH = 5', 'THERMALLY_PULSING_ASYMPTOTIC_GIANT_BRANCH = 6', 'NAKED_MAIN_SEQUENCE_HELIUM_STAR = 7', 'NAKED_HELIUM_STAR_HERTZSPRUNG_GAP = 8', 'NAKED_HELIUM_STAR_GIANT_BRANCH = 9', 'HELIUM_WHITE_DWARF = 10', 'CARBON_OXYGEN_WHITE_DWARF = 11', 'OXYGEN_NEON_WHITE_DWARF = 12', 'NEUTRON_STAR = 13', 'BLACK_HOLE = 14', 'STAR_WITH_NO_MASS = 15']
Parameter: max_stellar_type_2
Description: The maximum stellar type of star 2 (internal index 1). Evolution is stopped when the star reaches this stellar type. If this is negative, massless remnants are allowed, and the maximum stellar type is the absolute value.
Parameter input type: Integer
Default value: 16
Macros: ['LOW_MASS_MS = 0', 'MS = 1', 'HG = 2', 'GIANT_BRANCH = 3', 'CHeB = 4', 'EAGB = 5', 'TPAGB = 6', 'HeMS = 7', 'HeHG = 8', 'HeGB = 9', 'HeWD = 10', 'COWD = 11', 'ONeWD = 12', 'NS = 13', 'BH = 14', 'MASSLESS_REMNANT = 15', 'LOW_MASS_MAIN_SEQUENCE = 0', 'MAIN_SEQUENCE = 1', 'HERTZSPRUNG_GAP = 2', 'FIRST_GIANT_BRANCH = 3', 'CORE_HELIUM_BURNING = 4', 'EARLY_ASYMPTOTIC_GIANT_BRANCH = 5', 'THERMALLY_PULSING_ASYMPTOTIC_GIANT_BRANCH = 6', 'NAKED_MAIN_SEQUENCE_HELIUM_STAR = 7', 'NAKED_HELIUM_STAR_HERTZSPRUNG_GAP = 8', 'NAKED_HELIUM_STAR_GIANT_BRANCH = 9', 'HELIUM_WHITE_DWARF = 10', 'CARBON_OXYGEN_WHITE_DWARF = 11', 'OXYGEN_NEON_WHITE_DWARF = 12', 'NEUTRON_STAR = 13', 'BLACK_HOLE = 14', 'STAR_WITH_NO_MASS = 15']
Parameter: max_stellar_type_3
Description: The maximum stellar type of star 3 (internal index 2). Evolution is stopped when the star reaches this stellar type. If this is negative, massless remnants are allowed, and the maximum stellar type is the absolute value.
Parameter input type: Integer
Default value: 16
Macros: ['LOW_MASS_MS = 0', 'MS = 1', 'HG = 2', 'GIANT_BRANCH = 3', 'CHeB = 4', 'EAGB = 5', 'TPAGB = 6', 'HeMS = 7', 'HeHG = 8', 'HeGB = 9', 'HeWD = 10', 'COWD = 11', 'ONeWD = 12', 'NS = 13', 'BH = 14', 'MASSLESS_REMNANT = 15', 'LOW_MASS_MAIN_SEQUENCE = 0', 'MAIN_SEQUENCE = 1', 'HERTZSPRUNG_GAP = 2', 'FIRST_GIANT_BRANCH = 3', 'CORE_HELIUM_BURNING = 4', 'EARLY_ASYMPTOTIC_GIANT_BRANCH = 5', 'THERMALLY_PULSING_ASYMPTOTIC_GIANT_BRANCH = 6', 'NAKED_MAIN_SEQUENCE_HELIUM_STAR = 7', 'NAKED_HELIUM_STAR_HERTZSPRUNG_GAP = 8', 'NAKED_HELIUM_STAR_GIANT_BRANCH = 9', 'HELIUM_WHITE_DWARF = 10', 'CARBON_OXYGEN_WHITE_DWARF = 11', 'OXYGEN_NEON_WHITE_DWARF = 12', 'NEUTRON_STAR = 13', 'BLACK_HOLE = 14', 'STAR_WITH_NO_MASS = 15']
Parameter: max_stellar_type_4
Description: The maximum stellar type of star 4 (internal index 3). Evolution is stopped when the star reaches this stellar type. If this is negative, massless remnants are allowed, and the maximum stellar type is the absolute value.
Parameter input type: Integer
Default value: 16
Macros: ['LOW_MASS_MS = 0', 'MS = 1', 'HG = 2', 'GIANT_BRANCH = 3', 'CHeB = 4', 'EAGB = 5', 'TPAGB = 6', 'HeMS = 7', 'HeHG = 8', 'HeGB = 9', 'HeWD = 10', 'COWD = 11', 'ONeWD = 12', 'NS = 13', 'BH = 14', 'MASSLESS_REMNANT = 15', 'LOW_MASS_MAIN_SEQUENCE = 0', 'MAIN_SEQUENCE = 1', 'HERTZSPRUNG_GAP = 2', 'FIRST_GIANT_BRANCH = 3', 'CORE_HELIUM_BURNING = 4', 'EARLY_ASYMPTOTIC_GIANT_BRANCH = 5', 'THERMALLY_PULSING_ASYMPTOTIC_GIANT_BRANCH = 6', 'NAKED_MAIN_SEQUENCE_HELIUM_STAR = 7', 'NAKED_HELIUM_STAR_HERTZSPRUNG_GAP = 8', 'NAKED_HELIUM_STAR_GIANT_BRANCH = 9', 'HELIUM_WHITE_DWARF = 10', 'CARBON_OXYGEN_WHITE_DWARF = 11', 'OXYGEN_NEON_WHITE_DWARF = 12', 'NEUTRON_STAR = 13', 'BLACK_HOLE = 14', 'STAR_WITH_NO_MASS = 15']
Parameter: probability
Description: The probability is a weighting applied to the star based on, say, the initial mass function. When running a grid of stars to simulate all stars, the summed probability of all the stars should be 1.0.
Parameter input type: Float
Default value: 1
Parameter: phasevol
Description: The system's phase volume, used by binary_grid.
Parameter input type: Float
Default value: NULL
Parameter: stellar_structure_algorithm
Description: Set the stellar structure algorithm. 0=modified BSE (default), 1=none, 2=external function (must be defined by the calling code), 3=binary_c (not yet implemented).
Parameter input type: Integer
Default value: 0
Macros: ['STELLAR_STRUCTURE_ALGORITHM_MODIFIED_BSE = 0', 'STELLAR_STRUCTURE_ALGORITHM_NONE = 1', 'STELLAR_STRUCTURE_ALGORITHM_EXTERNAL_FUNCTION = 2', 'STELLAR_STRUCTURE_ALGORITHM_MINT = 3']
Parameter: solver
Description: The type of solver. Default is the Forward-Euler (0), but could be RK2 (1), RK4 (2) or a predictor-corretor (3).
Parameter input type: Integer
Default value: 0
Macros: ['SOLVER_FORWARD_EULER = 0', 'SOLVER_RK2 = 1', 'SOLVER_RK4 = 2', 'SOLVER_PREDICTOR_CORRECTOR = 3']
Parameter: max_evolution_time
Description: Set the maximum age for the stars (Myr).
Parameter input type: Float
Default value: 15000
Parameter: max_model_number
Description: Set the maximum number of models, ignored if 0 (default is 0).
Parameter input type: Integer
Default value: 0
Parameter: monte_carlo_kicks
Description: Turn on Monte-Carlo SN kicks. On (True) by default, and indeed other algorithms are probably broken.
Parameter input type: True|False
Default value: True
Parameter: disable_debug
Description: Disables debug output. Only has an effect when DEBUG is 1, which probably requires a rebuild. Default FALSE.
Parameter input type: True|False
Default value: False
Parameter: timestep_logging
Description: Turn on timestep logging (default is False).
Parameter input type: True|False
Default value: False
Parameter: rejects_in_log
Description: Show timestep rejections in the main log (default is False).
Parameter input type: True|False
Default value: False
Parameter: vandenHeuvel_logging
Description: Turn on van den Heuvel logging (default is False).
Parameter input type: True|False
Default value: False
Parameter: evolution_splitting
Description: If True, turn on splitting of an evolutionary run if splitpoint (e.g. supernova) occurs.
Parameter input type: True|False
Default value: False
Parameter: disable_events
Description: Whether to disable the new events code (defaults to False, so we use events by default)
Parameter input type: True|False
Default value: False
Parameter: evolution_splitting_sn_eccentricity_threshold
Description: Threshold eccentricity above which evolution splitting happens in a system with no SN kick. (0.01)
Parameter input type: Float
Default value: 0.01
Parameter: evolution_splitting_sn_n
Description: Number of runs to split into when a SN occurs.
Parameter input type: Integer
Default value: 10
Parameter: evolution_splitting_maxdepth
Description: Max number of splits in an evolutionary run.
Parameter input type: Integer
Default value: 2
Parameter: equation_of_state_algorithm
Description: Set the equation of state algorithm. 0 = Paczynski.
Parameter input type: Integer
Default value: NULL
Macros: ['EQUATION_OF_STATE_PACZYNSKI = 0']
Parameter: opacity_algorithm
Description: Set the opacity algorithm. 0 = Paczynski, 1 = Ferguson/Opal.
Parameter input type: Integer
Default value: NULL
Macros: ['OPACITY_ALGORITHM_PACZYNSKI = 0', 'OPACITY_ALGORITHM_FERGUSON_OPAL = 1', 'OPACITY_ALGORITHM_STARS = 2']
Parameter: wind_mass_loss
Description: Defines the algorithm used for stellar winds. 0 = none, 1 = Hurley et al. (2002), 2 = Schneider (2018), 3 = Schneider + Sander and Wink (2020) for helium stars, 4 = Schneider + Hurley for LBV mass loss
Parameter input type: Unsigned integer
Default value: 3
Macros: ['WIND_ALGORITHM_NONE = 0', 'WIND_ALGORITHM_HURLEY2002 = 1', 'WIND_ALGORITHM_SCHNEIDER2018 = 2', 'WIND_ALGORITHM_BINARY_C_2020 = 3', 'WIND_ALGORITHM_HENDRIKS_2022 = 4']
Extra: 0
Parameter: gbwind
Description: Wind prescription for first red giant branch stars. 0=Reimers (Hurley et al 2000/2002; choose gb_reimers_eta=0.5 for their mass loss rate), 1=Schroeder+Cuntz 2005 (set gb_reimers_eta=1.0 for their mass loss rate).
Parameter input type: Integer
Default value: 0
Macros: ['GB_WIND_REIMERS = 0', 'GB_WIND_SCHROEDER_CUNTZ_2005 = 1', 'GB_WIND_GOLDMAN_ETAL_2017 = 2', 'GB_WIND_BEASOR_ETAL_2020 = 3']
Parameter: postagbwind
Description: Apply special post-(A)GB prescription. Default is POSTAGB_WIND_USE_GIANT which means we just use whatever is prescribed on the giant branch. Other options include: POSTAGB_WIND_NONE = 1 (no wind on the post (A)GB), POSTAGB_WIND_KRTICKA2020 = 2 which uses Krticka, Kubát and Krticková (2020, A&A 635, A173).
Parameter input type: Integer
Default value: 0
Macros: ['POSTAGB_WIND_GIANT = 0', 'POSTAGB_WIND_NONE = 1', 'POSTAGB_WIND_KRTICKA2020 = 2']
Parameter: Teff_postAGB_min
Description: The minimum temperature for which we apply post-(A)GB winds. See also Teff_postAGB_max. (6000 K)
Parameter input type: Float
Default value: 6000
Parameter: Teff_postAGB_max
Description: The maximum temperature for which we apply post-(A)GB winds. See also Teff_postAGB_min. (120000 K)
Parameter input type: Float
Default value: 120000
Parameter: mattsson_Orich_tpagbwind
Description: Experimental: turns on Mattsson's TPAGB wind when the star is oxygen rich. Requires MATTSSON_MASS_LOSS.
Parameter input type: Integer
Default value: NULL
Parameter: magnetic_braking_factor
Description: Multiplier for the magnetic braking angular momentum loss rate.
Parameter input type: Float
Default value: 1
Parameter: magnetic_braking_gamma
Description: gamma factor in Rappaport style magnetic braking expression.
Parameter input type: Float
Default value: 3
Parameter: magnetic_braking_algorithm
Description: Algorithm for the magnetic braking angular momentum loss rate. 0 = Hurley et al. 2002, 1 = Andronov, Pinnsonneault and Sills 2003, 2 = Barnes and Kim 2010
Parameter input type: Integer
Default value: 0
Macros: ['MAGNETIC_BRAKING_ALGORITHM_HURLEY_2002 = 0', 'MAGNETIC_BRAKING_ALGORITHM_ANDRONOV_2003 = 1', 'MAGNETIC_BRAKING_ALGORITHM_BARNES_2010 = 2', 'MAGNETIC_BRAKING_ALGORITHM_RAPPAPORT_1983 = 3']
Parameter: helium_flash_mass_loss
Description: Mass to be lost at the helium flash.
Parameter input type: Float
Default value: 0
Parameter: gb_reimers_eta
Description: First red giant branch wind multiplication factor, cf. eta in Reimers' mass loss formula. (This multiplies the 4e-13 in Reimers' formula, or the 8e-14 in Schroeder and Cuntz.)
Parameter input type: Float
Default value: 0.4
Parameter: gbwindfac
Description: Multiplier for the giant branch wind mass loss rate
Parameter input type: Float
Default value: 1
Parameter: tpagbwindfac
Description: Multiplier for the TPAGB wind mass loss rate
Parameter input type: Float
Default value: 1
Parameter: eagbwindfac
Description: Multiplier for the EAGB wind mass loss rate
Parameter input type: Float
Default value: 1
Parameter: nieuwenhuijzen_windfac
Description: Multiplier for the Nieuwenhuijzen & de Jager wind mass loss rate
Parameter input type: Float
Default value: 1
Parameter: tpagbwind
Description: Wind prescription during the TPAGB. 0=Karakas 2002 (a modified Vassiliadis and Wood 1993), 1=Hurley et al 2000/2002 (Vassiliadis and Wood 1993), 2=Reimers, 3=Bloecker, 4=Van Loon, 5=Rob's C-wind (broken?), 6,7=Vassiliadis and Wood 1993 (Karakas,Hurley variants respectively) when C/O>1, 8=Mattsson, 9 = Goldman et al. (2017), 10 = Beasor et al. (2020).
Parameter input type: Integer
Default value: 0
Macros: ['TPAGB_WIND_VW93_KARAKAS = 0', 'TPAGB_WIND_VW93_ORIG = 1', 'TPAGB_WIND_REIMERS = 2', 'TPAGB_WIND_BLOECKER = 3', 'TPAGB_WIND_VAN_LOON = 4', 'TPAGB_WIND_ROB_CWIND = 5', 'TPAGB_WIND_VW93_KARAKAS_CARBON_STARS = 6', 'TPAGB_WIND_VW93_ORIG_CARBON_STARS = 7', 'TPAGB_WIND_MATTSSON = 8', 'TPAGB_WIND_GOLDMAN_ETAL_2017 = 9', 'TPAGB_WIND_BEASOR_ETAL_2020 = 10']
Parameter: eagbwind
Description: Wind prescription during the EAGB. 0=BSE (Hurley+2002, based on VW93), 1 = Goldman et al. (2017), 2 = Beasor et al. (2020).
Parameter input type: Integer
Default value: 0
Macros: ['EAGB_WIND_BSE = 0', 'EAGB_WIND_GOLDMAN_ETAL_2017 = 1', 'EAGB_WIND_BEASOR_ETAL_2020 = 2']
Parameter: wind_gas_to_dust_ratio
Description: Gas to dust ratio used in wind calculations (e.g. Goldman et al. 2017). Typically 200 (Milky Way)-500 (Magellanic Clouds). Default is 200, approximately as in MW stars.
Parameter input type: Float
Default value: 200
Parameter: vwind_multiplier
Description: Multiplier for the stellar wind velocity.
Parameter input type: Float
Default value: 1
Parameter: vwind_beta
Description: Beta for stellar wind speed calculations, where vwind=sqrt(beta) * escape velocity. Default 0.125 (from BSE, Hurley et al. 2002).
Parameter input type: Float
Default value: 0.125
Parameter: superwind_mira_switchon
Description: In the Vassiliadis and Wood (1993) AGB wind prescription, the superwind is turned on at a given Mira period, usually 500 days. You can vary when this switch-on happens with this parameter.
Parameter input type: Float
Default value: 500
Parameter: vw93_mira_shift
Description: In the Vassiliadis and Wood (1993) AGB wind prescription, the wind loss rate depends on the Mira period plus this offset. Requires VW93_MIRA_SHIFT.
Parameter input type: Float
Default value: 0
Parameter: vw93_multiplier
Description: In the Vassiliadis and Wood (1993) AGB wind prescription, the wind loss rate is multiplied by this factor. Requires VW93_MULTIPLIER.
Parameter input type: Float
Default value: 1
Parameter: tpagb_reimers_eta
Description: TPAGB Reimers wind multiplication factor, cf. eta in Reimers' mass loss formula. (This multiplies the 4e-13 in Reimers' formula, or the 8e-14 in Schroeder and Cuntz.) Note that Reimers is not the default TPAGB wind prescription. See also tpagbwind.
Parameter input type: Float
Default value: 1
Parameter: Tout_Pringle_1992_multiplier
Description: Multiplier for the Tout & Pringle (1992) magnetic wind. (0.0)
Parameter input type: Float
Default value: 0
Parameter: artificial_mass_accretion_rate%d
Description: Constant mass accretion rate for star <n>.
Parameter input type: Float(scanf)
Default value: NULL
Parameter: artificial_mass_accretion_rate_by_stellar_type%d
Description: Constant mass accretion rate for stellar type <n>.
Parameter input type: Float(scanf)
Default value: NULL
Parameter: artificial_angular_momentum_accretion_rate%d
Description: Constant angular momentum accretion for star <n>.
Parameter input type: Float(scanf)
Default value: NULL
Parameter: artificial_orbital_angular_momentum_accretion_rate
Description: Constant angular momentum accretion rate on the orbit.
Parameter input type: Float
Default value: 0
Parameter: artificial_accretion_start_time
Description: Time at which artificial accretion stars. Ignored if <0 (default is -1).
Parameter input type: Float
Default value: -1
Macros: ['ARTIFICIAL_ACCRETION_IGNORE = -1']
Parameter: artificial_accretion_end_time
Description: Time at which artificial accretion ends. Ignored if <0 (default is -1).
Parameter input type: Float
Default value: -1
Macros: ['ARTIFICIAL_ACCRETION_IGNORE = -1']
Parameter: wr_wind
Description: Massive-star (WR) wind prescription. 0 = Hurley et al 2000/2002, 1=Maeder and Meynet, 2=Nugis and Lamers, 3=John Eldridge's version of Vink's early-2000s wind (See Lynnette Dray's thesis, or John Eldridge's thesis)
Parameter input type: Integer
Default value: 0
Macros: ['WR_WIND_BSE = 0', 'WR_WIND_MAEDER_MEYNET = 1', 'WR_WIND_NUGIS_LAMERS = 2', 'WR_WIND_ELDRIDGE = 3']
Parameter: wr_wind_fac
Description: Massive-star (WR) wind multiplication factor.
Parameter input type: Float
Default value: 1
Parameter: wrwindfac
Description: Massive-star (WR) wind multiplication factor. Synonymous with wr_wind_fac (which you should use instead).
Parameter input type: Float
Default value: 1
Parameter: BH_prescription
Description: Black hole mass prescrition: relates the mass of a newly formed black hole to its progenitor's (CO) core mass. BH_HURLEY2002 = 0 = Hurley et al 2000/2002, BH_BELCZYNSKI = 1 = Belczynski (early 2000s), BH_SPERA2015 = Spera+ 2015, BH_FRYER12_DELAYED = 3 = Fryer et al. (2012) delayed prescription, BH_FRYER12_RAPID = 4 = Fryer et al. (2012) rapid prescription, BH_FRYER12_STARTRACK = 5 = Fryer et al. (2012) startrack prescription.
Parameter input type: Integer
Default value: 0
Macros: ['BH_HURLEY2002 = 0', 'BH_BELCZYNSKI = 1', 'BH_SPERA2015 = 2', 'BH_FRYER12_DELAYED = 3', 'BH_FRYER12_RAPID = 4', 'BH_FRYER12_STARTRACK = 5']
Parameter: PPISN_prescription
Description: (Pulsational) Pair-Instability Supernova prescription: Relates initial helium core mass of star to whether the star undergoes PPISN or PISN. Requires PPISN flag to be True (see binary_c_parameters.h). 0=no ppisn, 1=Farmer et al 2019. 2=Hendriks et al 21 (in prep), which is different fit based on the same data as Farmer et al 2019
Parameter input type: Integer
Default value: 1
Macros: ['PPISN_NONE = 0', 'PPISN_FARMER19 = 1', 'PPISN_NEW_FIT_21 = 2']
Extra: Ignore
Parameter: PPISN_additional_massloss
Description: (Pulsational) Pair-Instability Supernova additional massloss (on top of the result of the prescription). Value currently is not required to be positive
Parameter input type: Float
Default value: 0
Extra: Ignore
Parameter: PPISN_core_mass_range_shift
Description: (Pulsational) Pair-Instability Supernova CO core mass range shift: Value by which we shift the range of CO core mass values that undergo PPISN. This is taken into account in the fit that calculates the mass removal. Negative number is shifting the range to lower masses, positive to higher masses.
Parameter input type: Float
Default value: 0
Extra: Ignore
Parameter: use_LWN_fallback
Description: Flag whether to use the Nazehdin-Lovegrove-Woosley hydrogen envelope ejection at failed supernova.
Parameter input type: Integer
Default value: 1
Extra: Ignore
Parameter: sn_kick_distribution_II
Description: Set the distribution of speeds applied to kick type II core collapse supernova systems. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Integer
Default value: 1
Macros: ['KICK_VELOCITY_FIXED = 0', 'KICK_VELOCITY_MAXWELLIAN = 1', 'KICK_VELOCITY_CUSTOM = 2']
Parameter: sn_kick_distribution_ECAP
Description: Set the distribution of speeds applied to the remnants of electron-capture supernovae. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Integer
Default value: 0
Macros: ['KICK_VELOCITY_FIXED = 0', 'KICK_VELOCITY_MAXWELLIAN = 1', 'KICK_VELOCITY_CUSTOM = 2']
Parameter: sn_kick_distribution_NS_NS
Description: Set the distribution of speeds applied to kick neutron stars and black holes that survive a NS-NS merger. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Integer
Default value: 0
Macros: ['KICK_VELOCITY_FIXED = 0', 'KICK_VELOCITY_MAXWELLIAN = 1', 'KICK_VELOCITY_CUSTOM = 2']
Parameter: sn_kick_distribution_IBC
Description: Set the distribution of speeds applied to kick newly-born neutron stars and black holes after a type Ib/c core-collapse supernova. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Integer
Default value: 1
Macros: ['KICK_VELOCITY_FIXED = 0', 'KICK_VELOCITY_MAXWELLIAN = 1', 'KICK_VELOCITY_CUSTOM = 2']
Parameter: sn_kick_distribution_GRB_COLLAPSAR
Description: Set the distribution of speeds applied to kick newly-born neutron stars and black holes after a type Ib/c core-collapse supernova which is also a collapsar. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Integer
Default value: 0
Macros: ['KICK_VELOCITY_FIXED = 0', 'KICK_VELOCITY_MAXWELLIAN = 1', 'KICK_VELOCITY_CUSTOM = 2']
Parameter: sn_kick_distribution_TZ
Description: Set the distribution of speeds applied to kick newly-born neutron stars and black holes at the death of a Thorne-Zytkow object. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Integer
Default value: 0
Macros: ['KICK_VELOCITY_FIXED = 0', 'KICK_VELOCITY_MAXWELLIAN = 1', 'KICK_VELOCITY_CUSTOM = 2']
Parameter: sn_kick_distribution_AIC_BH
Description: Set the distribution of speeds applied to kick newly-born neutron stars black holes after accretion induced collapse of a neutron star. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Integer
Default value: 0
Macros: ['KICK_VELOCITY_FIXED = 0', 'KICK_VELOCITY_MAXWELLIAN = 1', 'KICK_VELOCITY_CUSTOM = 2']
Parameter: sn_kick_distribution_BH_BH
Description: Set the distribution of speeds applied to black holes formed by the merger of two black holes. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Integer
Default value: 0
Macros: ['KICK_VELOCITY_FIXED = 0', 'KICK_VELOCITY_MAXWELLIAN = 1', 'KICK_VELOCITY_CUSTOM = 2']
Parameter: sn_kick_distribution_BH_NS
Description: Set the distribution of speeds applied to black holes formed by the merger of a neutron star and a black holes. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Integer
Default value: 0
Macros: ['KICK_VELOCITY_FIXED = 0', 'KICK_VELOCITY_MAXWELLIAN = 1', 'KICK_VELOCITY_CUSTOM = 2']
Parameter: sn_kick_distribution_IA_Hybrid_HeCOWD
Description: Set the distribution of speeds applied to any survivor of a hybrid He-COWD SNIa explosion. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Integer
Default value: 0
Macros: ['KICK_VELOCITY_FIXED = 0', 'KICK_VELOCITY_MAXWELLIAN = 1', 'KICK_VELOCITY_CUSTOM = 2']
Parameter: sn_kick_distribution_IA_Hybrid_HeCOWD_subluminous
Description: Set the distribution of speeds applied to any survivor of a subluminous hybrid He-COWD SNIa explosion. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Integer
Default value: 0
Macros: ['KICK_VELOCITY_FIXED = 0', 'KICK_VELOCITY_MAXWELLIAN = 1', 'KICK_VELOCITY_CUSTOM = 2']
Parameter: sn_kick_distribution_PPISN
Description: Set the distribution of speeds applied to PPISN supernovae. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Integer
Default value: 0
Macros: ['KICK_VELOCITY_FIXED = 0', 'KICK_VELOCITY_MAXWELLIAN = 1', 'KICK_VELOCITY_CUSTOM = 2']
Parameter: sn_kick_distribution_PISN
Description: Set the distribution of speeds applied to PISN supernovae. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Integer
Default value: 0
Macros: ['KICK_VELOCITY_FIXED = 0', 'KICK_VELOCITY_MAXWELLIAN = 1', 'KICK_VELOCITY_CUSTOM = 2']
Parameter: sn_kick_distribution_PHDIS
Description: Set the distribution of speeds applied to PHDIS supernovae. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Integer
Default value: 0
Macros: ['KICK_VELOCITY_FIXED = 0', 'KICK_VELOCITY_MAXWELLIAN = 1', 'KICK_VELOCITY_CUSTOM = 2']
Parameter: sn_kick_dispersion_II
Description: Set the dispersion of speeds applied to kick type II core collapse supernova systems. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Float
Default value: 190
Parameter: sn_kick_dispersion_ECAP
Description: Set the dispersion of speeds applied to the remnants of electron-capture supernovae. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Float
Default value: 0
Parameter: sn_kick_dispersion_NS_NS
Description: Set the dispersion of speeds applied to kick neutron stars and black holes that survive a NS-NS merger. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Float
Default value: 0
Parameter: sn_kick_dispersion_IBC
Description: Set the dispersion of speeds applied to kick newly-born neutron stars and black holes after a type Ib/c core-collapse supernova. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Float
Default value: 190
Parameter: sn_kick_dispersion_GRB_COLLAPSAR
Description: Set the dispersion of speeds applied to kick newly-born neutron stars and black holes after a type Ib/c core-collapse supernova which is also a collapsar. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Float
Default value: 0
Parameter: sn_kick_dispersion_TZ
Description: Set the dispersion of speeds applied to kick newly-born neutron stars and black holes at the death of a Thorne-Zytkow object. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Float
Default value: 0
Parameter: sn_kick_dispersion_AIC_BH
Description: Set the dispersion of speeds applied to kick newly-born neutron stars black holes after accretion induced collapse of a neutron star. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Float
Default value: 0
Parameter: sn_kick_dispersion_BH_BH
Description: Set the dispersion of speeds applied to black holes formed by the merger of two black holes. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Float
Default value: 0
Parameter: sn_kick_dispersion_BH_NS
Description: Set the dispersion of speeds applied to black holes formed by the merger of a neutron star and a black holes. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Float
Default value: 0
Parameter: sn_kick_dispersion_IA_Hybrid_HeCOWD
Description: Set the dispersion of speeds applied to the survivor of a SNIa explosion of a hybrid He-COWD. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Float
Default value: 0
Parameter: sn_kick_dispersion_IA_Hybrid_HeCOWD_subluminous
Description: Set the dispersion of speeds applied to the survivor of a subluminous SNIa explosion of a hybrid He-COWD. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Float
Default value: 0
Parameter: sn_kick_dispersion_PPISN
Description: Set the dispersion of speeds applied to the survivor of a PPISN supernova. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Float
Default value: 0
Parameter: sn_kick_dispersion_PISN
Description: Set the dispersion of speeds applied to the survivor of a PISN supernova. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Float
Default value: 0
Parameter: sn_kick_dispersion_PHDIS
Description: Set the dispersion of speeds applied to the survivor of a PHDIS supernova. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).
Parameter input type: Float
Default value: 0
Parameter: sn_kick_companion_IA_He
Description: Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when a Ia He supernova occurs. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.
Parameter input type: Float
Default value: 0
Macros: ['SN_IMPULSE_NONE = 0', 'SN_IMPULSE_LIU2015 = 1', 'SN_IMPULSE_WHEELER1975 = 2']
Parameter: sn_kick_companion_IA_ELD
Description: Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when a Ia ELD (sub-Mch) supernova occurs. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.
Parameter input type: Float
Default value: 0
Macros: ['SN_IMPULSE_NONE = 0', 'SN_IMPULSE_LIU2015 = 1', 'SN_IMPULSE_WHEELER1975 = 2']
Parameter: sn_kick_companion_IA_CHAND
Description: Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when a Ia Mch supernova occurs. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.
Parameter input type: Float
Default value: 0
Macros: ['SN_IMPULSE_NONE = 0', 'SN_IMPULSE_LIU2015 = 1', 'SN_IMPULSE_WHEELER1975 = 2']
Parameter: sn_kick_companion_AIC
Description: Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when an accretion induced collapse (supernova) occurs. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.
Parameter input type: Float
Default value: 0
Macros: ['SN_IMPULSE_NONE = 0', 'SN_IMPULSE_LIU2015 = 1', 'SN_IMPULSE_WHEELER1975 = 2']
Parameter: sn_kick_companion_ECAP
Description: Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when an electron capture supernova occurs. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.
Parameter input type: Float
Default value: 0
Macros: ['SN_IMPULSE_NONE = 0', 'SN_IMPULSE_LIU2015 = 1', 'SN_IMPULSE_WHEELER1975 = 2']
Parameter: sn_kick_companion_IA_He_Coal
Description: Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when a Ia helium merger supernova occurs. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.
Parameter input type: Float
Default value: 0
Macros: ['SN_IMPULSE_NONE = 0', 'SN_IMPULSE_LIU2015 = 1', 'SN_IMPULSE_WHEELER1975 = 2']