eagle_6.yml 14.1 KB
Newer Older
1 2 3
MetaData:
  run_name:   EAGLE-L0006N0094-Ref

4 5
# Define the system of units to use internally. 
InternalUnitSystem:
6 7
  UnitMass_in_cgs:     1.98848e43    # 10^10 M_sun in grams
  UnitLength_in_cgs:   3.08567758e24 # Mpc in centimeters
8 9 10 11
  UnitVelocity_in_cgs: 1e5           # km/s in centimeters per second
  UnitCurrent_in_cgs:  1             # Amperes
  UnitTemp_in_cgs:     1             # Kelvin

12
# Structure finding options
13
StructureFinding:
14 15 16
  config_file_name:     stf_input.cfg # Name of the STF config file.
  basename:             ./stf         # Common part of the name of output files.
  scale_factor_first:   0.92          # Scale-factor of the first snaphot (cosmological run)
17
  time_first:           0.01        # Time of the first structure finding output (in internal units).
18
  delta_time:           1.10          # Time difference between consecutive structure finding outputs (in internal units) in simulation time intervals.
19

20 21 22 23 24 25 26
# Cosmological parameters
Cosmology:
  h:              0.6777        # Reduced Hubble constant
  a_begin:        0.9090909     # Initial scale-factor of the simulation
  a_end:          1.0           # Final scale factor of the simulation
  Omega_m:        0.307         # Matter density parameter
  Omega_lambda:   0.693         # Dark-energy density parameter
27
  Omega_b:        0.0482519     # Baryon density parameter
28
  
29
Scheduler:
30
  max_top_level_cells:    8
31
  links_per_tasks: 500
32
  
33 34 35 36 37
# Parameters governing the time integration
TimeIntegration:
  time_begin: 0.    # The starting time of the simulation (in internal units).
  time_end:   1e-2  # The end time of the simulation (in internal units).
  dt_min:     1e-10 # The minimal time-step size of the simulation (in internal units).
38
  dt_max:     1e-3  # The maximal time-step size of the simulation (in internal units).
39
  
40 41 42
# Parameters governing the snapshots
Snapshots:
  basename:            eagle # Common part of the name of output files
43
  scale_factor_first:  0.91  # Scale-factor of the first snaphot (cosmological run)
44
  time_first:          0.01  # Time of the first output (non-cosmological run) (in internal units)
45
  delta_time:          1.01  # Time difference between consecutive outputs (in internal units)
46
  compression:         1
47 48 49

# Parameters governing the conserved quantities statistics
Statistics:
50 51 52
  scale_factor_first:  0.92 # Scale-factor of the first stat dump (cosmological run)
  time_first:          0.01 # Time of the first stat dump (non-cosmological run) (in internal units)
  delta_time:          1.05 # Time between statistics output
53 54 55

# Parameters for the self-gravity scheme
Gravity:
56
  eta:                    0.025    # Constant dimensionless multiplier for time integration.
57 58
  theta:                  0.7      # Opening angle (Multipole acceptance criterion)
  mesh_side_length:       16
59
  dithering:              0
60 61 62 63
  comoving_DM_softening:         0.0026994 # Comoving DM softening length (in internal units).
  max_physical_DM_softening:     0.0007    # Max physical DM softening length (in internal units).
  comoving_baryon_softening:     0.0026994 # Comoving DM softening length (in internal units).
  max_physical_baryon_softening: 0.0007    # Max physical DM softening length (in internal units).
64

65 66
# Parameters for the hydrodynamics scheme
SPH:
67 68 69 70 71 72 73
  resolution_eta:        1.2348           # Target smoothing length in units of the mean inter-particle separation (1.2348 == 48Ngbs with the cubic spline kernel).
  h_min_ratio:           0.1              # Minimal smoothing in units of softening.
  h_max:                 0.5              # Maximal softening in co-moving internal units.
  CFL_condition:         0.1              # Courant-Friedrich-Levy condition for time integration.
  minimal_temperature:   100              # (internal units)
  particle_splitting:    1
  particle_splitting_mass_threshold: 7e-4 # Internal units (i.e. 7e6 Msun ~ 4 times the initial gas mass)
74

75 76 77 78
# Parameters of the stars neighbour search
Stars:
  resolution_eta:        1.1642   # Target smoothing length in units of the mean inter-particle separation
  h_tolerance:           7e-3
79 80
  overwrite_birth_time:  1
  birth_time:            0.33333  # Pretend all the stars were born at z = 2
81

82 83
# Parameters for the Friends-Of-Friends algorithm
FOF:
84 85 86 87
  basename:                        fof_output  # Filename for the FOF outputs.
  min_group_size:                  256         # The minimum no. of particles required for a group.
  linking_length_ratio:            0.2         # Linking length in units of the main inter-particle separation.
  black_hole_seed_halo_mass_Msun:  1.5e10      # Minimal halo mass in which to seed a black hole (in solar masses).
88 89
  scale_factor_first:              0.91        # Scale-factor of first FoF black hole seeding calls.
  delta_time:                      1.005       # Scale-factor ratio between consecutive FoF black hole seeding calls.
90

91 92 93
# Parameters related to the initial conditions
InitialConditions:
  file_name:  ./EAGLE_ICs_6.hdf5     # The file to read
94
  periodic:   1
95
  cleanup_h_factors: 1               # Remove the h-factors inherited from Gadget
96 97
  cleanup_velocity_factors: 1        # Remove the sqrt(a) factor in the velocities inherited from Gadget

98 99 100 101 102 103 104 105 106 107 108
EAGLEChemistry: 	     # Solar abundances
  init_abundance_metal:      0.014
  init_abundance_Hydrogen:   0.70649785
  init_abundance_Helium:     0.28055534
  init_abundance_Carbon:     2.0665436e-3
  init_abundance_Nitrogen:   8.3562563e-4
  init_abundance_Oxygen:     5.4926244e-3
  init_abundance_Neon:       1.4144605e-3
  init_abundance_Magnesium:  5.907064e-4
  init_abundance_Silicon:    6.825874e-4
  init_abundance_Iron:       1.1032152e-3
109

110
EAGLECooling:
111
  dir_name:                ./coolingtables/
112
  H_reion_z:               11.5 
113
  H_reion_eV_p_H:          2.0
114 115
  He_reion_z_centre:       3.5
  He_reion_z_sigma:        0.5
116
  He_reion_eV_p_H:         2.0
117

118 119 120 121 122 123 124
# EAGLE star formation parameters
EAGLEStarFormation:
  EOS_density_norm_H_p_cm3:          0.1       # Physical density used for the normalisation of the EOS assumed for the star-forming gas in Hydrogen atoms per cm^3.
  EOS_temperature_norm_K:            8000      # Temperature om the polytropic EOS assumed for star-forming gas at the density normalisation in Kelvin.
  EOS_gamma_effective:               1.3333333 # Slope the of the polytropic EOS assumed for the star-forming gas.
  KS_normalisation:                  1.515e-4  # The normalization of the Kennicutt-Schmidt law in Msun / kpc^2 / yr.
  KS_exponent:                       1.4       # The exponent of the Kennicutt-Schmidt law.
125
  min_over_density:                  57.7      # The over-density above which star-formation is allowed.
126 127
  KS_high_density_threshold_H_p_cm3: 1e3       # Hydrogen number density above which the Kennicut-Schmidt law changes slope in Hydrogen atoms per cm^3.
  KS_high_density_exponent:          2.0       # Slope of the Kennicut-Schmidt law above the high-density threshold.
128
  EOS_entropy_margin_dex:            0.5       # Logarithm base 10 of the maximal entropy above the EOS at which stars can form.
129 130 131 132 133 134 135 136 137 138 139 140 141 142 143
  threshold_norm_H_p_cm3:            0.1       # Normalisation of the metal-dependant density threshold for star formation in Hydrogen atoms per cm^3.
  threshold_Z0:                      0.002     # Reference metallicity (metal mass fraction) for the metal-dependant threshold for star formation.
  threshold_slope:                   -0.64     # Slope of the metal-dependant star formation threshold
  threshold_max_density_H_p_cm3:     10.0      # Maximal density of the metal-dependant density threshold for star formation in Hydrogen atoms per cm^3.
  
# Parameters for the EAGLE "equation of state"
EAGLEEntropyFloor:
  Jeans_density_threshold_H_p_cm3: 0.1       # Physical density above which the EAGLE Jeans limiter entropy floor kicks in expressed in Hydrogen atoms per cm^3.
  Jeans_over_density_threshold:    10.       # Overdensity above which the EAGLE Jeans limiter entropy floor can kick in.
  Jeans_temperature_norm_K:        8000      # Temperature of the EAGLE Jeans limiter entropy floor at the density threshold expressed in Kelvin.
  Jeans_gamma_effective:           1.3333333 # Slope the of the EAGLE Jeans limiter entropy floor
  Cool_density_threshold_H_p_cm3: 1e-5       # Physical density above which the EAGLE Cool limiter entropy floor kicks in expressed in Hydrogen atoms per cm^3.
  Cool_over_density_threshold:    10.        # Overdensity above which the EAGLE Cool limiter entropy floor can kick in.
  Cool_temperature_norm_K:        8000       # Temperature of the EAGLE Cool limiter entropy floor at the density threshold expressed in Kelvin.
  Cool_gamma_effective:           1.         # Slope the of the EAGLE Cool limiter entropy floor
144

145 146
# EAGLE feedback model
EAGLEFeedback:
147 148
  use_SNII_feedback:                1               # Global switch for SNII thermal (stochastic) feedback.
  use_SNIa_feedback:                1               # Global switch for SNIa thermal (continuous) feedback.
149 150 151 152 153 154
  use_AGB_enrichment:               1               # Global switch for enrichement from AGB stars.
  use_SNII_enrichment:              1               # Global switch for enrichement from SNII stars.
  use_SNIa_enrichment:              1               # Global switch for enrichement from SNIa stars.
  filename:                         ./yieldtables/  # Path to the directory containing the EAGLE yield tables.
  IMF_min_mass_Msun:                0.1             # Minimal stellar mass considered for the Chabrier IMF in solar masses.
  IMF_max_mass_Msun:              100.0             # Maximal stellar mass considered for the Chabrier IMF in solar masses.
155 156 157 158
  SNII_min_mass_Msun:               6.0             # Minimal mass considered for SNII stars in solar masses.
  SNII_max_mass_Msun:             100.0             # Maximal mass considered for SNII stars in solar masses.
  SNII_sampled_delay:               0               # Sample the SNII lifetimes to do feedback.
  SNII_wind_delay_Gyr:              0.03            # Time in Gyr between a star's birth and the SNII thermal feedback event when not sampling.
159
  SNII_delta_T_K:                   3.16228e7       # Change in temperature to apply to the gas particle in a SNII thermal feedback event in Kelvin.
160
  SNII_energy_erg:                  1.0e51          # Energy of one SNII explosion in ergs.
161 162
  SNII_energy_fraction_min:         0.3             # Minimal fraction of energy applied in a SNII feedback event.
  SNII_energy_fraction_max:         3.0             # Maximal fraction of energy applied in a SNII feedback event.
163
  SNII_energy_fraction_Z_0:         0.0012663729    # Pivot point for the metallicity dependance of the SNII energy fraction (metal mass fraction).
164
  SNII_energy_fraction_n_0_H_p_cm3: 1.4588          # Pivot point for the birth density dependance of the SNII energy fraction in cm^-3.
165 166
  SNII_energy_fraction_n_Z:         0.8686          # Power-law for the metallicity dependance of the SNII energy fraction.
  SNII_energy_fraction_n_n:         0.8686          # Power-law for the birth density dependance of the SNII energy fraction.
167 168 169
  SNIa_DTD:                         PowerLaw        # Functional form of the SNIa delay time distribution.
  SNIa_DTD_delay_Gyr:               0.04            # Stellar age after which SNIa start in Gyr (40 Myr corresponds to stars ~ 8 Msun).
  SNIa_DTD_power_law_norm_p_Msun:   0.0012          # Normalization of the SNIa delay time distribution (in Msun^-1).
170
  SNIa_energy_erg:                 1.0e51           # Energy of one SNIa explosion in ergs.
171
  AGB_ejecta_velocity_km_p_s:      10.0             # Velocity of the AGB ejectas in km/s.
172 173
  stellar_evolution_age_cut_Gyr:    0.1             # Stellar age in Gyr above which the enrichment is down-sampled.
  stellar_evolution_sampling_rate:   10             # Number of time-steps in-between two enrichment events for a star above the age threshold.
174 175 176 177 178 179 180 181 182
  SNII_yield_factor_Hydrogen:       1.0             # (Optional) Correction factor to apply to the Hydrogen yield from the SNII channel.
  SNII_yield_factor_Helium:         1.0             # (Optional) Correction factor to apply to the Helium yield from the SNII channel.
  SNII_yield_factor_Carbon:         0.5             # (Optional) Correction factor to apply to the Carbon yield from the SNII channel.
  SNII_yield_factor_Nitrogen:       1.0             # (Optional) Correction factor to apply to the Nitrogen yield from the SNII channel.
  SNII_yield_factor_Oxygen:         1.0             # (Optional) Correction factor to apply to the Oxygen yield from the SNII channel.
  SNII_yield_factor_Neon:           1.0             # (Optional) Correction factor to apply to the Neon yield from the SNII channel.
  SNII_yield_factor_Magnesium:      2.0             # (Optional) Correction factor to apply to the Magnesium yield from the SNII channel.
  SNII_yield_factor_Silicon:        1.0             # (Optional) Correction factor to apply to the Silicon yield from the SNII channel.
  SNII_yield_factor_Iron:           0.5             # (Optional) Correction factor to apply to the Iron yield from the SNII channel.
183

184 185
# EAGLE AGN model
EAGLEAGN:
186
  subgrid_seed_mass_Msun:           1.5e5      # Black hole subgrid mass at creation time in solar masses.
187
  max_eddington_fraction:           1.5        # Maximal allowed accretion rate in units of the Eddington rate.
188
  eddington_fraction_for_recording: 0.1        # Record the last time BHs reached an Eddington ratio above this threshold.
189
  viscous_alpha:                    1e6        # Normalisation constant of the Bondi viscuous time-scale accretion reduction term
190 191 192 193
  radiative_efficiency:             0.1        # Fraction of the accreted mass that gets radiated.
  coupling_efficiency:              0.15       # Fraction of the radiated energy that couples to the gas in feedback events.
  AGN_delta_T_K:                    3.16228e8  # Change in temperature to apply to the gas particle in an AGN feedback event in Kelvin.
  AGN_num_ngb_to_heat:              1.         # Target number of gas neighbours to heat in an AGN feedback event.
194 195 196
  max_reposition_mass:              2e8        # Maximal BH mass considered for BH repositioning in solar masses.
  threshold_major_merger:           0.333      # Mass ratio threshold to consider a BH merger as 'major'
  threshold_minor_merger:           0.1        # Mass ratio threshold to consider a BH merger as 'minor'