MetaData: run_name: StromgrenSphere-3D-singlebin # Define the system of units to use internally. InternalUnitSystem: UnitMass_in_cgs: 1.98841586e+33 # 1 M_Sol UnitLength_in_cgs: 3.08567758e21 # kpc in cm UnitVelocity_in_cgs: 1.e5 # km/s UnitCurrent_in_cgs: 1. UnitTemp_in_cgs: 1. # K # Parameters governing the time integration TimeIntegration: max_nr_rt_subcycles: 1 time_begin: 0. # The starting time of the simulation (in internal units). time_end: 0.512 # end time dt_min: 1.e-12 # The minimal time-step size of the simulation (in internal units). dt_max: 1.e-03 # The maximal time-step size of the simulation (in internal units). # Parameters governing the snapshots Snapshots: basename: output_singlebin # Common part of the name of output files time_first: 0. # Time of the first output (in internal units) delta_time: 0.0102 # Time between snapshots # Parameters governing the conserved quantities statistics Statistics: time_first: 0. delta_time: 1e-3 # Time between statistics output # 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). CFL_condition: 0.8 # Courant-Friedrich-Levy condition for time integration. minimal_temperature: 10 # Kelvin # Parameters related to the initial conditions InitialConditions: file_name: ./stromgrenSphere-3D.hdf5 # The file to read periodic: 1 # periodic ICs. Stars: resolution_eta: 2.2348 # (Optional) Target smoothing length in units of the mean inter-particle separation (1.2348 == 48Ngbs with the cubic spline kernel). Defaults to the SPH value. GEARRT: f_reduce_c: 0.01 # reduce the speed of light for the RT solver by multiplying c with this factor CFL_condition: 0.9 photon_groups_Hz: [3.288e15] # Lower photon frequency group bin edges in Hz. Needs to have exactly N elements, where N is the configured number of bins --with-RT=GEAR_N stellar_luminosity_model: const # Which model to use to determine the stellar luminosities. const_stellar_luminosities_LSol: [6.198024e+04] # # (Conditional) constant star luminosities for each photon frequency group to use if stellar_luminosity_model:const is set, in units of Solar Luminosity. hydrogen_mass_fraction: 1.00 # total hydrogen (H + H+) mass fraction in the metal-free portion of the gas set_equilibrium_initial_ionization_mass_fractions: 0 # (Optional) set the initial ionization fractions depending on gas temperature assuming ionization equilibrium. set_initial_ionization_mass_fractions: 1 # (Optional) manually overwrite initial mass fraction of each species (using the values you set below) mass_fraction_HI: 0.999999 # (Conditional) If set_initial_ionization_fractions=1, needed to set initial HI mass fractions to this value mass_fraction_HII: 1.e-6 # (Conditional) If set_initial_ionization_fractions=1, needed to set initial HII mass fractions to this value mass_fraction_HeI: 0. # (Conditional) If set_initial_ionization_fractions=1, needed to set initial HeI mass fractions to this value mass_fraction_HeII: 0. # (Conditional) If set_initial_ionization_fractions=1, needed to set initial HeII mass fractions to this value mass_fraction_HeIII: 0. # (Conditional) If set_initial_ionization_fractions=1, needed to set initial HeIII mass fractions to this value stellar_spectrum_type: 1 # Which radiation spectrum to use. 0: constant from 0 until some max frequency set by stellar_spectrum_const_max_frequency_Hz. 1: blackbody spectrum. stellar_spectrum_blackbody_temperature_K: 1.e5 # (Conditional) if stellar_spectrum_type=1, use this temperature (in K) for the blackbody spectrum. stars_max_timestep: 6.250000e-05 # (Optional) restrict the maximal timestep of stars to this value (in internal units) # commented parameters are for isothermal single-frequency bin Stromgren Sphere SPHM1RT: cred: 2.99792458e3 # value of reduced speed of light for the RT solver in code unit CFL_condition: 0.1 # CFL condition for RT, independent of hydro chi: [0, 0, 0, 0] # (Optional) initial opacity in code unit for all gas particles photon_groups_Hz: [3.288e15, 5.945e19, 5.946e19] # Photon frequency group bin edges in Hz. Needs to be 1 less than the number of groups (N) requested during the configuration (--with-RT=SPHM1RT_N). Outer edges of zero and infinity are assumed. use_const_emission_rates: 1 # (Optional) use constant emission rates for stars as defined with star_emission_rates parameter star_emission_rates: [1e-32, 2.34e11, 1e-32, 1e-32] # (Optional) constant star emission rates (internal unit: energy/time) for each photon frequency group to use if use_constant_emission_rates is set. # = 5e48 photons/sec stellar_spectrum_type: 1 # Which radiation spectrum to use. 0: constant from 0 until some max frequency set by stellar_spectrum_const_max_frequency_Hz. 1: blackbody spectrum. stars_max_timestep: 7.8125e-6 # (Optional) restrict the maximal timestep of stars to this value (in internal units). Set to negative to turn off. stellar_spectrum_blackbody_temperature_K: 1.e5 # (Conditional) if stellar_spectrum_type=1, use this temperature (in K) for the blackbody spectrum. skip_thermochemistry: 0 # (Optional) skip the thermochemistry. This is intended only for debugging and testing the radiation transport, as it breaks the purpose of RT. init_mass_fraction_metal: 0. # (Optional) Inital mass fraction of particle mass in *all* metals (if it is set, the initial fraction will be over-written.) init_mass_fraction_Hydrogen: 1.0 # (Conditional) (if init_mass_fraction_metal != -1.0f) Inital mass fraction of particle mass in Hydrogen init_mass_fraction_Helium: 0.0 # (Conditional) (if init_mass_fraction_metal != -1.0f) Inital mass fraction of particle mass in Helium coolingon: 0 # (Optional) switch for cooling (and photoheating), but photo-ionization will be ongoing even if coolingon==0 useabundances: 1 # (Optional) use the species abundances below, instead of reading from initial condition init_species_abundance_e: 1e-5 # (Conditional) (if useabundances==1) free electron abundances (in unit hydrogen number density:nH) init_species_abundance_HI: 0.99999 # (Conditional) (if useabundances==1) HI abundances (in unit hydrogen number density:nH) init_species_abundance_HII: 1e-5 # (Conditional) (if useabundances==1) HII abundances (in unit hydrogen number density:nH) init_species_abundance_HeI: 0.0 # (Conditional) (if useabundances==1) HeI abundances (in unit hydrogen number density:nH) init_species_abundance_HeII: 0.0 # (Conditional) (if useabundances==1) HeII abundances (in unit hydrogen number density:nH) init_species_abundance_HeIII: 0.0 # (Conditional) (if useabundances==1) HeIII abundances (in unit hydrogen number density:nH) relativeTolerance: 1e-3 # (Optional) Relative tolerance for SPHM1RT thermo-chemistry intergration absoluteTolerance: 1e-10 # (Optional) Absolute tolerance for SPHM1RT thermo-chemistry integration explicitTolerance: 0.1 # (Optional) Tolerance below which we use the explicit solution in SPHM1RT thermo-chemistry ionizing_photon_energy_erg: [3.0208e-11, 5.61973e-11, 1.05154e-10] # (Optional) ionizing photon energy in erg averaged within frequency bins useparams: 1 # (Optional) switch to use thermo-chemistry parameters from the parameter file sigma_cross: [8.13e-18, 1e-32, 1e-32] # (Conditional) (if useparams=1) The cross section of ionizing photons for hydrogen (cm^2) alphaB: 2.59e-13 # (Conditional) (if useparams=1) The case B recombination coefficient for hydrogen (cgs) beta: 3.1e-16 # (Conditional) (if useparams=1) The collisional ionization coefficient for hydrogen (cgs) reinject: 1 # (Optional) gather energy around injection radius and re-inject the energy