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

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

17
18
19
20
21
22
23
24
25
# 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
  Omega_b:        0.0455        # Baryon density parameter
  
26
27
28
# Parameters governing the time integration
TimeIntegration:
  time_begin: 0.    # The starting time of the simulation (in internal units).
29
30
  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).
31
  dt_max:     1e-2  # The maximal time-step size of the simulation (in internal units).
32

33
34
35
# Parameters governing the snapshots
Snapshots:
  basename:            eagle # Common part of the name of output files
36
  scale_factor_first:  0.91  # Scale-factor of the first snaphot (cosmological run)
37
  time_first:          0.01  # Time of the first output (non-cosmological run) (in internal units)
38
  delta_time:          1.01  # Time difference between consecutive outputs (in internal units)
39
40
41

# Parameters governing the conserved quantities statistics
Statistics:
42
43
44
  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
45

46
47
# Parameters for the self-gravity scheme
Gravity:
48
  eta:                    0.025    # Constant dimensionless multiplier for time integration. 
Matthieu Schaller's avatar
Matthieu Schaller committed
49
  theta:                  0.7     # Opening angle (Multipole acceptance criterion)
50
51
  comoving_softening:     0.0026994 # Comoving softening length (in internal units).
  max_physical_softening: 0.0007    # Physical softening length (in internal units).
52
  mesh_side_length:       256
53

54
55
56
# Parameters for the hydrodynamics scheme
SPH:
  resolution_eta:        1.2348   # Target smoothing length in units of the mean inter-particle separation (1.2348 == 48Ngbs with the cubic spline kernel).
57
  h_min_ratio:           0.1      # Minimal smoothing in units of softening.
58
  CFL_condition:         0.1      # Courant-Friedrich-Levy condition for time integration.
59
  minimal_temperature:   100      # (internal units)
60
61
62

# Parameters related to the initial conditions
InitialConditions:
63
  file_name:  ./EAGLE_ICs_25.hdf5    # The file to read
64
  periodic:   1
65
  cleanup_h_factors: 1               # Remove the h-factors inherited from Gadget
66
67
  cleanup_velocity_factors: 1        # Remove the sqrt(a) factor in the velocities inherited from Gadget

68
69
70
71
72
73
74
75
76
77
78
79
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
  
80
EAGLECooling:
81
  dir_name:                ./coolingtables/
82
  H_reion_z:               11.5 
83
  H_reion_eV_p_H:          2.0
84
85
  He_reion_z_centre:       3.5
  He_reion_z_sigma:        0.5
86
  He_reion_eV_p_H:         2.0
87

88
89
90
Scheduler:
  links_per_tasks:    25

91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
# 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.
  KS_min_over_density:               57.7      # The over-density above which star-formation is allowed.
  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.
  KS_temperature_margin_dex:         0.5       # Logarithm base 10 of the maximal temperature difference above the EOS allowed to form stars.
  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
117
118
119

# EAGLE feedback model
EAGLEFeedback:
120
121
  use_SNII_feedback:                1               # Global switch for SNII thermal (stochastic) feedback.
  use_SNIa_feedback:                1               # Global switch for SNIa thermal (continuous) feedback.
122
123
124
125
126
127
128
129
130
131
  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.
  SNII_min_mass_Msun:               6.0             # Minimal mass considered for SNII feedback (not SNII enrichment!) in solar masses.
  SNII_max_mass_Msun:             100.0             # Maximal mass considered for SNII feedback (not SNII enrichment!) in solar masses.
  SNII_wind_delay_Gyr:              0.03            # Time in Gyr between a star's birth and the SNII thermal feedback event.
  SNII_delta_T_K:                   3.16228e7       # Change in temperature to apply to the gas particle in a SNII thermal feedback event in Kelvin.
132
133
134
135
136
137
138
  SNII_energy_erg:                  1.0e51          # Energy of one SNII explosion in ergs.
  SNII_energy_fraction_min:         3.0             # Maximal fraction of energy applied in a SNII feedback event.
  SNII_energy_fraction_max:         0.3             # Minimal fraction of energy applied in a SNII feedback event.
  SNII_energy_fraction_Z_0:         0.0012663729    # Pivot point for the metallicity dependance of the SNII energy fraction (metal mass fraction).
  SNII_energy_fraction_n_0_H_p_cm3: 0.67            # Pivot point for the birth density dependance of the SNII energy fraction in cm^-3.
  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.
139
140
141
  SNIa_max_mass_Msun:              8.0              # Maximal mass considered for SNIa feedback and enrichment in solar masses.
  SNIa_timescale_Gyr:              2.0              # Time-scale of the exponential decay of the SNIa rates in Gyr.
  SNIa_efficiency_p_Msun:          0.002            # Normalisation of the SNIa rates in inverse solar masses.
142
  SNIa_energy_erg:                 1.0e51           # Energy of one SNIa explosion in ergs.
143
144
145
146
147
148
149
150
151
152
  AGB_ejecta_velocity_km_p_s:      10.0             # Velocity of the AGB ejectas in km/s.
  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.
153
154
155
156
157
158
159
160

# EAGLE AGN model
EAGLEAGN:
  max_eddington_fraction:           1.         # Maximal allowed accretion rate in units of the Eddington rate.
  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.