Sample Parameter Files from Publications¶
Here is a collection of parameter files that are equivalent or similar to ones used in publications. Each example has a description and motivation behind the choice of important parameters.
Wise et al. (2014)¶
(Publication link)
(Parameter file
)
(MUSIC parameter file
)
(Lagrangian data file 622kB)
This parameter file will set up a simulation similar to the paper’s simulation. Some parameter choices are updated to be consistent with papers written after the original simulation was run in 2012. It models Pop III star formation and the transition to galaxy formation. While the production simulation was a 2563 AMR everywhere simulation, this sample file is a zoom-in calculation (the user needs to generate the ICs with MUSIC) focused on a single 108 solar mass halo at z = 10.
Hierarchy setup¶
1 Mpc comoving box with a 643 root grid resolution and 2 nested grids
MaximumRefinementLevel = 14
: maximal comoving resolution of ~1 pcMaximumParticleRefinementLevel = 11
: to remove the discreteness of DM particles, smooth them over ~1 proper pc below which baryons become dominantCellFlaggingMethod = 2 4 6 8
: Refine by baryon/DM mass, Jeans length, and must-refine cosmological particlesRefineByJeansLengthSafetyFactor = 8
: Resolve the local Jeans length by at least 8 cellsMustRefineParticlesCreateParticles = 3
: Use a dynamical nested region that’s specified by the particle types given by MUSICMustRefineParticlesRefineToLevel = 2
: Refine the innermost region to at least level 2
Hydrodynamics¶
RiemannSolver = 4
: use the HLLC solver for more stability in supernova blastwavesRiemannSolverFallback = 1
: use the HLL solver when HLLC fails
Cooling and chemistry¶
MultiSpecies = 3
: 12-species primordial radiative coolingMetalCooling = 1
: tabulated metal coolinggrackle_data_file = cloudy_metals_2008_3D-lwb.h5
: metal cooling with a Lyman-Werner (LW) background from Qin+ (2020)UVbackground = 1
: Use the LW background the table aboveH2_self_shielding = 1
: Wolcott-Green et al. (2011) model
Star formation and feedback¶
Metal-free and metal-enriched star formation with radiative and supernova feedback
StarParticleCreation = 40
: Use model #3 (Pop III stars) and #5 (radiating star clusters)StarParticleFeedback = 40
: Use model #3 and #5
Pop III particles¶
PopIIIOverDensityThreshold = -1e6
: The negative indicates units in particles per cm3. This number density is consistent with Pop III formation at the resolution level of 0.1 proper pc.PopIIIH2CriticalFraction = 1e-3
: This is an appropriate H2 fraction at the above density thresholdPopIIIMetalCriticalFraction = 1.295e-6
: Use a critical transition metallicity of 10-4 of solarPopIIIUseHypernova = 0
: Use a typical 1051 erg core-collapse supernova energy for stars between 11 and 40 solar massesPopIIIInitialMassFunction = 1
: Use a Pop III IMF defined in Wise et al. (2014)PopIIIStarMass = 20
: Use a characteristic mass of 20 solar masses in the above IMF. Motivated by Hirano et al. (2015) and Susa (2019)
Star cluster particles¶
StarClusterMinDynamicalTime = 3e+06
: Motivated by Tan et al. (2006) who found that molecular clouds have dynamical times around 700 kyr and form stars over several dynamical times. There could be newer references.StarClusterIonizingLuminosity = 1.9e+46
: Use the average over 20 Myr from the metal-poor binary model of Rosdahl et al. (2018)StarClusterSNEnergy = 1e49
: Calculated from a Salpeter IMF, i.e. 1 SN per 100 solar masses of starsStarClusterFormEfficiency = 0.07
: Motivated by Krumholz et al. (2007). This converts 7% of the cold gas within the star forming cloud into stars. There could be newer references.StarClusterMinimumMass = 1000
: (in solar masses) equivalent to the smallest observed star clusters
Radiative transfer¶
Uses the adaptive ray tracing machinery that is sourced by both types of star particles
RadiativeTransferOpticallyThinH2 = 1
: Source a (1/r2) LW radiation field from the star particles. Important for star formation regulationRadiativeTransferRadiationPressure = 1
: Consider momentum transfer between the ionizing photons and absorbing gas. Plays a role at small distances in dense gas and high fluxRadiativeTransferSourceClustering = 1
: Use ray merging to conserve memory (and maybe speed it up)RadiativeTransferPhotonMergeRadius = 3.0
: Merge rays at 3 times the separation between sources that are paired in a tree