N-Body Simulation Data Operators¶

Class providing operators that act on data from N-body simulations—post-processing transformations applied sequentially to the nBodyData objects imported by nbodyImporterClass. Typical operations include computing halo masses, spin parameters, concentrations, merger histories, or galaxy occupation fractions from the raw particle data. Multiple operators can be chained to build a complete analysis pipeline, with the null operator leaving the data unchanged for pass-through use cases.

Default implementation: nbodyOperatorNull

Methods¶

operate → void

Apply this operator to the provided array of N-body simulation datasets, transforming or augmenting particle and halo properties in place.

type(nBodyData), intent(inout), dimension(:) :: simulations

`nbodyOperatorAddAttributes`¶

An N-body data operator which attaches named scalar attributes to the simulation dataset, enabling metadata such as simulation parameters or analysis flags to be stored alongside particle data. Attribute names and corresponding values are specified as arrays via the [names] and [values] parameters.

Parameters

[names] — A list of attribute names.
[values] — A list of attribute values.

`nbodyOperatorAngularMomentum`¶

An N-body data operator which determines the mean angular momentum of particles. Also finds the angular velocity vector for the rotating frame in which the angular momentum would be zero.

Parameters

[selfBoundParticlesOnly] — If true, the mean angular momentum is computed only for self-bound particles
[bootstrapSampleCount] (default 30_c_size_t) — The number of bootstrap resamples of the particles that should be used.

`nbodyOperatorConcentrationDistributionFunction`¶

An N-body data operator which computes the halo concentration distribution function by binning halos as a function of concentration within a specified mass and concentration range. Mass limits and binning are set by [massMinimum] , [massMaximum], and [massCountPerDecade], concentration limits and binning by [concentrationMinimum], [concentrationMaximum], and [concentrationCountPerDecade].

Parameters

[massMinimum] — The minimum halo mass (in \(\mathrm{M}_\odot\)) below which halos are excluded from the concentration distribution function.
[massMaximum] — The maximum halo mass (in \(\mathrm{M}_\odot\)) above which halos are excluded from the concentration distribution function.
[massCountPerDecade] — The number of logarithmic bins per decade of mass used when constructing the concentration distribution function.
[concentrationMinimum] — The minimum halo concentration parameter below which halos are excluded from the distribution function histogram.
[concentrationMaximum] — The maximum halo concentration parameter above which halos are excluded from the distribution function histogram.
[concentrationCountPerDecade] — The number of logarithmic bins per decade of concentration parameter used when constructing the concentration distribution function.
[description] — A human-readable description of this concentration distribution function dataset, stored as metadata in the output file.
[simulationReference] — A bibliographic reference for the N-body simulation from which this concentration distribution is derived, stored as output metadata.
[simulationURL] — A URL pointing to the publicly accessible dataset or documentation for the N-body simulation, stored as output metadata.

`nbodyOperatorConvexHull`¶

An N-body data operator which constructs the convex hull of the particles.

Methods

tabulate — Tabulate the virial density contrast as a function of mass and time.
restoreTable — Restore a tabulated solution from file.
storeTable — Store a tabulated solution to file.

Parameters

[velocityCharacteristic] (default 250.0d0) — The velocity scale at which the SNe-driven outflow rate equals the star formation rate in disks.
[exponent] (default 3.5d0) — The velocity scaling of the SNe-driven outflow rate in disks.
[fraction] (default 0.01d0) — The normalization \(f\) of the outflow rate relative to the star formation rate at a reference halo velocity of 200 km/s and expansion factor of 1, setting the overall mass-loading amplitude of the halo-scaling feedback model.
[exponentVelocity] (default -2.0d0) — The exponent of virial velocity in the outflow rate in disks.
[exponentRedshift] (default 0.0d0) — The power-law exponent of the cosmological expansion factor \((1+z)\) in the halo-scaling outflow rate, allowing the mass-loading factor to evolve with redshift; a value of zero gives no redshift evolution.
[toleranceRelativeVelocityDispersion] (default 1.0d-6) — The relative tolerance to use in numerical solutions for the velocity dispersion in dark-matter-only density profiles.
[toleranceRelativeVelocityDispersionMaximum] (default 1.0d-3) — The maximum relative tolerance to use in numerical solutions for the velocity dispersion in dark-matter-only density profiles.
[radiusNormalization] (default 3.3d-6) — The initial value appearing in the radius-mass relation
[toleranceAbsoluteMass] (default 1.0d-6) — The mass tolerance used to judge whether the nuclear star cluster is physically plausible.
[toleranceRelativeMetallicity] (default 1.0d-4) — The metallicity tolerance for ODE solution.
[inactiveLuminositiesStellar] (default .false.) — Specifies whether or not nuclear star cluster stellar luminosities are inactive properties (i.e. do not appear in any ODE being solved).
[scaleRelativeMass] (default 1.0d-2) — The mass scale, relative to the total mass of the node, below which calculations in the delayed very simple hot halo component are allowed to become inaccurate.
[starveSatellites] (default .false.) — Specifies whether or not the hot halo should be removed (“starved”) when a node becomes a satellite.
[starveSatellitesOutflowed] (default .false.) — Specifies whether or not the outflowed hot halo should be removed (“starved”) when a node becomes a satellite.
[outflowReturnOnFormation] (default .false.) — Specifies whether or not outflowed gas should be returned to the hot reservoir on halo formation events.
[angularMomentumAlwaysGrows] (default .false.) — Specifies whether or not negative rates of accretion of angular momentum into the hot halo will be treated as positive for the purposes of computing the hot halo angular momentum.
[fractionBaryonLimitInNodeMerger] (default .false.) — Controls whether the hot gas content of nodes should be limited to not exceed the universal baryon fraction at node merger events. If set to true, hot gas (and angular momentum, abundances, and chemicals proportionally) will be removed from the merged halo to the unaccreted gas reservoir to limit the baryonic mass to the universal baryon fraction where possible.
[scaleAbsoluteMass] (default 100.0d0) — The absolute mass scale below which calculations in the very simple disk component are allowed to become inaccurate.
[toleranceAbsoluteMass] (default 1.0d-6) — The mass tolerance used to judge whether the disk is physically plausible.
[toleranceAbsoluteMass] (default 1.0d-6) — The mass tolerance used to judge whether the disk is physically plausible.
[toleranceRelativeMetallicity] (default 1.0d-4) — The metallicity tolerance for ODE solution.
[radiusStructureSolver] (default 1.0d0) — The radius (in units of the standard scale length) to use in solving for the size of the disk.
[structureSolverUseCole2000Method] (default .false.) — If true, use the method described in Cole et al. (2000) to correct for difference between thin disk and spherical mass distributions when solving for disk radii.
[diskNegativeAngularMomentumAllowed] (default .true.) — Specifies whether or not negative angular momentum is allowed for the disk.
[inactiveLuminositiesStellar] (default .false.) — Specifies whether or not disk stellar luminosities are inactive properties (i.e. do not appear in any ODE being solved).
[postStepZeroNegativeMasses] (default .true.) — If true, negative masses will be zeroed after each ODE step. Note that this can lead to non-conservation of mass.
[ratioAngularMomentumSolverRadius] (default ratioAngularMomentumSolverRadiusDefault) — The assumed ratio of the specific angular momentum at the structure solver radius to the mean specific angular momentum of the standard disk component.
[scaleAbsoluteMass] (default 100.0d0) — The absolute mass scale below which calculations in the very simple spheroid component are allowed to become inaccurate.
[toleranceAbsoluteMass] (default 1.0d-6) — The mass tolerance used to judge whether the spheroid is physically plausible.
[efficiencyEnergeticOutflow] (default 1.0d-2) — The proportionality factor relating mass outflow rate from the spheroid to the energy input rate divided by \(V_\mathrm{spheroid}^2\).
[toleranceRelativeMetallicity] (default 1.0d-4) — The metallicity tolerance for ODE solution.
[toleranceAbsoluteMass] (default 1.0d-6) — The mass tolerance used to judge whether the spheroid is physically plausible.
[inactiveLuminositiesStellar] (default .false.) — Specifies whether or not spheroid stellar luminosities are inactive properties (i.e. do not appear in any ODE being solved).
[postStepZeroNegativeMasses] (default .true.) — If true, negative masses will be zeroed after each ODE step. Note that this can lead to non-conservation of mass.
[ratioAngularMomentumScaleRadius] (default ratioAngularMomentumScaleRadiusDefault) — The assumed ratio of the specific angular momentum at the scale radius to the mean specific angular momentum of the standard spheroid component.
[outputMergers] (default .false.) — Determines whether or not properties of black hole mergers will be output.
[fileNames] — The name of the file(s) from which merger tree data should be read when using the [mergerTreeConstruct]\(=\)read tree construction method.
[forestSizeMaximum] (default 0_c_size_t) — The maximum number of nodes allowed in a forest before it will be broken up into trees and processed individually. A value of 0 implies that forests should never be split.
[presetMergerTimes] (default .true.) — Specifies whether merging times for subhalos should be preset when reading merger trees from a file.
[presetMergerNodes] (default .true.) — Specifies whether the target nodes for mergers should be preset (i.e. determined from descendant nodes). If they are not, merging will be with each satellite’s host node.
[presetSubhaloMasses] (default .true.) — Specifies whether subhalo mass should be preset when reading merger trees from a file.
[subhaloAngularMomentaMethod] (default var_str('summation')) — Specifies how to account for subhalo angular momentum when adding subhalo mass to host halo mass.
[presetSubhaloIndices] (default .true.) — Specifies whether subhalo indices should be preset when reading merger trees from a file.
[presetPositions] (default .true.) — Specifies whether node positions should be preset when reading merger trees from a file.
[presetScaleRadii] (default .true.) — Specifies whether node scale radii should be preset when reading merger trees from a file.
[scaleRadiiFailureIsFatal] (default .true.) — Specifies whether failure to set a node scale radii should be regarded as a fatal error. (If not, a fallback method to set scale radius is used in such cases.)
[presetScaleRadiiConcentrationMinimum] (default 3.0d0) — The lowest concentration (\(c=r_\mathrm{vir}/r_\mathrm{s}\)) allowed when setting scale radii, \(r_\mathrm{s}\).
[presetScaleRadiiConcentrationMaximum] (default 60.0d0) — The largest concentration (\(c=r_\mathrm{vir}/r_\mathrm{s}\)) allowed when setting scale radii, \(r_\mathrm{s}\).
[presetScaleRadiiMinimumMass] (default 0.0d0) — The minimum halo mass for which scale radii should be preset (if [presetScaleRadii]\(=\)true).
[presetUnphysicalAngularMomenta] (default .false.) — When reading merger trees from file and presetting halo angular momenta, detect unphysical (<=0) angular momenta and preset them using the selected halo spin method.
[presetAngularMomenta] (default .true.) — Specifies whether node angular momenta should be preset when reading merger trees from a file.
[presetAngularMomenta3D] (default .false.) — Specifies whether node 3-D angular momenta vectors should be preset when reading merger trees from a file.
[presetOrbits] (default .true.) — Specifies whether node orbits should be preset when reading merger trees from a file.
[presetOrbitsSetAll] (default .true.) — Forces all orbits to be set. If the computed orbit does not cross the virial radius, then select one at random instead.
[presetOrbitsAssertAllSet] (default .true.) — Asserts that all virial orbits must be preset. If any can not be set, Galacticus will stop.
[presetOrbitsBoundOnly] (default .true.) — Specifies whether only bound node orbits should be set.
[beginAt] (default -1_kind_int8) — Specifies the index of the tree to begin at. (Use -1 to always begin with the first tree.)
[outputTimeSnapTolerance] (default 0.0d0) — The relative tolerance required to “snap” a node time to the closest output time.
[missingHostsAreFatal] (default .true.) — Specifies whether nodes with missing host nodes should be considered to be fatal—see the discussion of missing host nodes in the class description above.
[treeIndexToRootNodeIndex] (default .false.) — Specifies whether tree indices should always be set to the index of their root node.
[allowBranchJumps] (default .true.) — Specifies whether nodes are allowed to jump between branches.
[allowSubhaloPromotions] (default .true.) — Specifies whether subhalos are permitted to be promoted to being isolated halos.
[alwaysPromoteMostMassive] (default .false.) — If true, the most massive progenitor is always promoted to be the primary progenitor even if it is a subhalo. Otherwise, isolated progenitors are given priority over subhalo progenitors, even if they are less massive.
[presetNamedReals] — Names of real datasets to be additionally read and stored in the nodes of the merger tree when using the [mergerTreeConstruct]\(=\)read tree construction method.
[presetNamedIntegers] — Names of integer datasets to be additionally read and stored in the nodes of the merger tree when using the [mergerTreeConstruct]\(=\)read tree construction method.
[fatalMismatches] (default .true.) — Specifies whether mismatches in cosmological parameter values between Galacticus and “Sussing Merger Trees” format (Srisawat et al., 2013) merger tree files should be considered fatal.
[fatalNonTreeNode] (default .true.) — Specifies whether nodes in snapshot files but not in the merger tree file should be considered fatal when importing from the “Sussing Merger Trees” format (Srisawat et al., 2013).
[subvolumeCount] (default 1) — Specifies the number of subvolumes along each axis into which a “Sussing Merger Trees” format (Srisawat et al., 2013) merger tree files should be split for processing through Galacticus.
[subvolumeBuffer] (default 0.0d0) — Specifies the buffer region (in units of Mpc\(/h\) to follow the format convention) around subvolumes of a “Sussing Merger Trees” format (Srisawat et al., 2013) merger tree file which should be read in to ensure that no halos are missed from trees.
[subvolumeIndex] (default [0,0,0]) — Specifies the index (in each dimension) of the subvolume of a “Sussing Merger Trees” format (Srisawat et al., 2013) merger tree file to process. Indices range from 0 to [subvolumeCount]\(-1\).
[badValue] (default -0.5d0) — Use for bad value detection in “Sussing” merger trees. Values for scale radius and halo spin which exceed this threshold are assumed to be bad.
[badValueTest] (default var_str('lessThan')) — Use for bad value detection in “Sussing” merger trees. Values which exceed the threshold in ths specified direction are assumed to be bad.
[treeSampleRate] (default 1.0d0) — Specify the probability that any given tree should processed (to permit subsampling).
[massOptions] (default var_str('default')) — Mass option for Sussing merger trees.
[mergeProbability] (default 0.1d0) — The largest probability of branching allowed in a timestep in merger trees built by the Cole et al. (2000) method.
[accretionLimit] (default 0.1d0) — The largest fractional mass change due to subresolution accretion allowed in a timestep in merger trees built by the Cole et al. (2000) method.
[redshiftMaximum] (default 1.0d5) — The highest redshift to which merger trees will be built in the Cole et al. (2000) method.
[toleranceTimeEarliest] (default 2.0d-6) — The fractional tolerance used to judge if a branch is at the earliest allowed time in the tree.
[branchIntervalStep] (default .true.) — If false use the original Cole et al. (2000) method to determine whether branching occurs in a timestep. If true draw branching intervals from a negative exponential distribution.
[toleranceResolutionSelf] (default 1.0d-6) — The fractional tolerance in node mass at the resolution limit below which branch mis-orderings will be ignored.
[toleranceResolutionParent] (default 1.0d-3) — The fractional tolerance in parent node mass at the resolution limit below which branch mis-orderings will be ignored.
[ignoreNoProgress] (default .false.) — If true, failure to make progress on a branch will be ignored (and the branch terminated).
[ignoreWellOrdering] (default .false.) — If true, non-well-ordered tree branches are pruned away instead of causing errors..
[redshiftBase] (default 0.0d0) — The redshift at which to plant the base node when building merger trees.
[timeSnapTolerance] (default 1.0d-6) — The fractional tolerance within which the tree base time will be snapped to a nearby output time.
[treeBeginAt] (default 0) — The index (in order of increasing base halo mass) of the tree at which to begin when building merger trees. A value of “0” means to begin with tree number 1 (if processing trees in ascending order), or equal to the number of trees (otherwise).
[processDescending] (default .true.) — If true, causes merger trees to be processed in order of decreasing mass.
[splitTrees] (default .false.) — If true, prune away any nodes of the tree that are not needed to determine evolution up to the latest time at which a node is present inside the lightcone. This typically leads to a tree splitting into a forest of trees.
[label] — A label for the mass function.
[comment] — A descriptive comment for the mass function.
[starFormationRates] — The star formation rates corresponding to bin centers.
[covarianceBinomialBinsPerDecade] (default 10) — The number of bins per decade of star formation rate to use when constructing star formation rate function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMinimum] (default 1.0d10) — The star formation rate to consider when constructing star formation rate function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMaximum] (default 1.0d12) — The maximum star formation rate to consider when constructing star formation rate function covariance matrices for main branch galaxies.
[targetLabel] — Label for the target dataset.
[functionValueTarget] — The target function for likelihood calculations.
[functionCovarianceTarget] — The target function covariance for likelihood calculations.
[likelihoodBins] — Controls which bins in the stellar mass–halo mass relation will be used in computing the likelihood: * not present: all bins are included in the likelihood calculation; * list of integers: use only the mass bin(s) given in this list in the likelihood calculation; * auto: use only bins which have a non-zero number of halos contributing to them in the likelihood calculation.
[fileNameTarget] — The name of the file containing the target data.
[redshiftInterval] (default 1) — The redshift interval to use.
[likelihoodNormalize] (default .false.) — If true, then normalize the likelihood to make it a probability density.
[computeScatter] (default .false.) — If true, the scatter in log10(stellar mass) is computed. Otherwise, the mean is computed.
[systematicErrorPolynomialCoefficient] (default [0.0d0]) — The coefficients of the systematic error polynomial for stellar mass in the stellar vs halo mass relation.
[systematicErrorMassHaloPolynomialCoefficient] (default [0.0d0]) — The coefficients of the systematic error polynomial for halo mass in the stellar vs halo mass relation.
[errorTolerant] (default .false.) — Error tolerance for the N-body spin distribution operator.
[logNormalRange] (default 100.0d0) — The multiplicative range of the log-normal distribution used to model the distribution of the mass and energy terms in the spin parameter. Specifically, the lognormal distribution is truncated outside the range \((\lambda_\mathrm{m}/R,\lambda_\mathrm{m} R\), where \(\lambda_\mathrm{m}\) is the measured spin, and \(R=\)[logNormalRange]
[fileName] — The name of the file from which to read spin distribution function parameters.
[comment] — A comment describing this analysis.
[label] — A label for this analysis.
[label] — A label for the spin distribution function.
[comment] — A descriptive comment for the spin distribution function.
[redshift] — The redshift at which to compute the spin distribution function.
[massMinimum] — Minimum halo mass for the spin distribution function.
[massMaximum] — Maximum halo mass for the spin distribution function.
[spinMinimum] — Minimum spin for the spin distribution function.
[spinMaximum] — Maximum spin for the spin distribution function.
[countSpinsPerDecade] — Number of spins per decade at which to compute the spin distribution function.
[timeRecent] — Halos which experienced a major node merger within a time \(\Delta t=\)[timeRecent] of the analysis time will be excluded from the analysis.
[particleCountMinimum] — The minimum particle count to assume when computing N-body errors on spins.
[massParticle] — The mass of the particle used in the N-body simulation from which spins were measured.
[energyEstimateParticleCountMaximum] — The maximum number of particles used in estimating halo energies when measuring spins from the N-body simulation.
[targetLabel] — Label for the target dataset.
[functionValueTarget] — The target function for likelihood calculations.
[functionCovarianceTarget] — The target function covariance for likelihood calculations.
[label] — A label for the mass function.
[comment] — A descriptive comment for the mass function.
[masses] — The masses corresponding to bin centers.
[covarianceBinomialBinsPerDecade] (default 10) — The number of bins per decade of halo mass to use when constructing HI mass function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMinimum] (default 1.0d8) — The minimum halo mass to consider when constructing HI mass function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMaximum] (default 1.0d16) — The maximum halo mass to consider when constructing HI mass function covariance matrices for main branch galaxies.
[targetLabel] — Label for the target dataset.
[functionValueTarget] — The target function for likelihood calculations.
[functionCovarianceTarget] — The target function covariance for likelihood calculations.
[label] — A label for the luminosity function.
[comment] — A descriptive comment for the luminosity function.
[magnitudesAbsolute] — The absolute magnitudes corresponding to bin centers.
[covarianceBinomialBinsPerDecade] (default 10) — The number of bins per decade of halo mass to use when constructing luminosity function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMinimum] (default 1.0d8) — The minimum halo mass to consider when constructing luminosity function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMaximum] (default 1.0d16) — The maximum halo mass to consider when constructing luminosity function covariance matrices for main branch galaxies.
[targetLabel] — Label for the target dataset.
[functionValueTarget] — The target function for likelihood calculations.
[functionCovarianceTarget] — The target function covariance for likelihood calculations.
[label] — A label for the luminosity function.
[comment] — A descriptive comment for the luminosity function.
[luminosities] — The luminosities corresponding to bin centers.
[covarianceBinomialBinsPerDecade] (default 10) — The number of bins per decade of halo mass to use when constructing luminosity function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMinimum] (default 1.0d8) — The minimum halo mass to consider when constructing luminosity function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMaximum] (default 1.0d16) — The maximum halo mass to consider when constructing luminosity function covariance matrices for main branch galaxies.
[includeNitrogenII] (default .false.) — If true, include contamination by the [NII] (6548AA \(+\) 6584AA) doublet.
[depthOpticalISMCoefficient] (default 1.0d0) — Multiplicative coefficient for optical depth in the ISM.
[targetLabel] — Label for the target dataset.
[functionValueTarget] — The target function for likelihood calculations.
[functionCovarianceTarget] — The target function covariance for likelihood calculations.
[label] — A label for the mass function.
[comment] — A descriptive comment for the mass function.
[masses] — The masses corresponding to bin centers.
[covarianceBinomialBinsPerDecade] (default 10) — The number of bins per decade of halo mass to use when constructing stellar mass function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMinimum] (default 1.0d8) — The minimum halo mass to consider when constructing stellar mass function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMaximum] (default 1.0d16) — The maximum halo mass to consider when constructing stellar mass function covariance matrices for main branch galaxies.
[targetLabel] — Label for the target dataset.
[functionValueTarget] — The target function for likelihood calculations.
[functionCovarianceTarget] — The target function covariance for likelihood calculations.
[rootVarianceFractionalMinimum] (default 0.0d0) — The minimum fractional root variance (relative to the target dataset).
[fileName] — The name of the file from which to read concentration distribution function parameters.
[comment] — A comment describing this analysis.
[label] — A label for this analysis.
[label] — A label for the concentration distribution function.
[comment] — A descriptive comment for the concentration distribution function.
[redshift] — The redshift at which to compute the concentration distribution function.
[massMinimum] — Minimum halo mass for the concentration distribution function.
[massMaximum] — Maximum halo mass for the concentration distribution function.
[concentrationMinimum] — Minimum concentration for the concentration distribution function.
[concentrationMaximum] — Maximum concentration for the concentration distribution function.
[countConcentrationsPerDecade] — Number of concentrations per decade at which to compute the concentration distribution function.
[timeRecent] — Halos which experienced a major node merger within a time \(\Delta t=\)[timeRecent] of the analysis time will be excluded from the analysis.
[massParticle] — The particle mass in the source N-body simulation.
[targetLabel] — Label for the target dataset.
[functionValueTarget] — The target function for likelihood calculations.
[functionCovarianceTarget] — The target function covariance for likelihood calculations.
[fileName] — The name of the file from which to read star forming main sequence function parameters.
[comment] — A comment describing this analysis.
[label] — A label for this analysis.
[label] — A label for the star forming main sequence function.
[comment] — A descriptive comment for the star forming main sequence function.
[massMinimum] — Minimum stellar mass for the star forming main sequence function.
[massMaximum] — Maximum stellar mass for the star forming main sequence function.
[countMassesPerDecade] — Number of masses per decade at which to compute the star forming main sequence function.
[targetLabel] — Label for the target dataset.
[meanValueTarget] — The target function for likelihood calculations.
[meanCovarianceTarget] — The target function covariance for likelihood calculations.
[label] — A label for the mass function.
[comment] — A descriptive comment for the mass function.
[separations] — The separations corresponding to bin centers.
[massMinima] — The minimum mass of each mass sample.
[massMaxima] — The maximum mass of each mass sample.
[massHaloBinsPerDecade] (default 10) — The number of bins per decade of halo mass to use when constructing the mass function covariance matrix for main branch galaxies.
[massHaloMinimum] (default 1.0d8) — The minimum halo mass to consider when constructing the mass function covariance matrix for main branch galaxies.
[massHaloMaximum] (default 1.0d16) — The maximum halo mass to consider when constructing the mass function covariance matrix for main branch galaxies.
[wavenumberCount] (default 60_c_size_t) — The number of bins in wavenumber to use in computing the correlation function.
[wavenumberMinimum] (default 1.0d-3) — The minimum wavenumber to use when computing the correlation function.
[wavenumberMaximum] (default 1.0d4) — The maximum wavenumber to use when computing the correlation function.
[integralConstraint] — The integral constraint for these correlation functions.
[depthLineOfSight] — The line-of-sight depth over which the correlation function was projected.
[halfIntegral] — Set to true if the projection integrand should be over line-of-sight depths greater than zero.
[binnedProjectedCorrelationTarget] — The target function for likelihood calculations.
[binnedProjectedCorrelationCovarianceTarget] — The target function covariance for likelihood calculations.
[targetLabel] (default var_str('')) — A label for the target dataset in a plot of this analysis.
[starFormationRateSpecificQuiescentLogarithmic] — The base-10 logarithm specific star formation rate (in units of Gyr\(^{-1}\)) separating quiescent and star-forming galaxies.
[starFormationRateSpecificLogarithmicError] — The observational fractional error in specific star formation rate (in units of dex) of galaxies.
[fileName] — The name of the file from which to read quiescent fraction function parameters.
[comment] — A comment describing this analysis.
[label] — A label for this analysis.
[label] — A label for the star forming main sequence function.
[comment] — A descriptive comment for the star forming main sequence function.
[massMinimum] — Minimum stellar mass for the star forming main sequence function.
[massMaximum] — Maximum stellar mass for the star forming main sequence function.
[countMassesPerDecade] — Number of masses per decade at which to compute the star forming main sequence function.
[targetLabel] — Label for the target dataset.
[meanValueTarget] — The target function for likelihood calculations.
[meanCovarianceTarget] — The target function covariance for likelihood calculations.
[nonAnalyticSolver] (default var_str('fallThrough')) — Selects how solutions are computed when no analytic solution is available. If set to “fallThrough” then the solution ignoring heating is used, while if set to “numerical” then numerical solvers are used to find solutions.
[radiusFractionalTruncateMinimum] (default 2.0d0) — The minimum radius (in units of the virial radius) to begin truncating the density profile.
[radiusFractionalTruncateMaximum] (default 4.0d0) — The maximum radius (in units of the virial radius) to finish truncating the density profile.
[nonAnalyticSolver] (default var_str('fallThrough')) — Selects how solutions are computed when no analytic solution is available. If set to “fallThrough” then the solution ignoring heating is used, while if set to “numerical” then numerical solvers are used to find solutions.
[velocityDispersionApproximate] (default .true.) — If true, radial velocity dispersion is computed using an approximate method in which we assume that \(\sigma_\mathrm{r}^2(r) \rightarrow \sigma_\mathrm{r}^2(r) - (2/3) \epsilon(r)\), where \(\epsilon(r)\) is the specific heating energy. If false then radial velocity dispersion is computed by numerically solving the Jeans equation.
[tolerateEnclosedMassIntegrationFailure] (default .false.) — If true, tolerate failures to find the mass enclosed as a function of radius.
[tolerateVelocityDispersionFailure] (default .false.) — If true, tolerate failures to compute the velocity dispersion.
[tolerateVelocityMaximumFailure] (default .false.) — If true, tolerate failures to find the radius of the maximum circular velocity.
[toleratePotentialIntegrationFailure] (default .false.) — If true, tolerate numerical failures when computing the gravitational potential of a heated dark matter profile, allowing the calculation to continue with a fallback result rather than aborting.
[toleranceRelativeVelocityDispersion] (default 1.0d-6) — The relative tolerance to use in numerical solutions for the velocity dispersion in dark-matter-only density profiles.
[toleranceRelativeVelocityDispersionMaximum] (default 1.0d-3) — The maximum relative tolerance to use in numerical solutions for the velocity dispersion in dark-matter-only density profiles.
[fractionRadiusFinalSmall] (default 1.0d-3) — The initial radius is limited to be no smaller than this fraction of the final radius. This can help avoid problems in profiles that are extremely close to being disrupted.
[toleranceRelativePotential] (default 1.0d-3) — The maximum allowed relative tolerance to use in numerical solutions for the gravitational potential in dark-matter-only density profiles before aborting.
[tolerateVelocityMaximumFailure] (default .true.) — If true, tolerate failures to find the radius of the peak in the rotation curve.
[lengthResolution] — The gravitational softening length \(\Delta x\) (in Mpc) of the N-body simulation, which sets the minimum spatial scale below which the dark matter profile is smoothed to avoid artificial two-body effects.
[massResolution] — The mass resolution \(\Delta M\) (in \(\mathrm{M}_\odot\)) of the N-body simulation, representing the minimum halo mass that can be resolved; profiles of halos near this limit are softened to account for particle discreteness effects.
[resolutionIsComoving] — If true, the resolution length is assumed to be fixed in comoving coordinates, otherwise in physical coordinates.
[nonAnalyticSolver] (default var_str('fallThrough')) — Selects how solutions are computed when no analytic solution is available. If set to “fallThrough” then the solution ignoring heating is used, while if set to “numerical” then numerical solvers are used to find solutions.
[C] (default 400.0d0) — The parameter \(C\) appearing in the halo concentration algorithm of Ludlow et al. (2016).
[f] (default 0.02d0) — The parameter \(f\) appearing in the halo concentration algorithm of Ludlow et al. (2016).
[timeFormationSeekDelta] (default 0.0d0) — The parameter \(\Delta \log t\) by which the logarithm of the trial formation time is incremented when stepping through the formation history of a node to find the formation time. If set to zero (or a negative value) the cumulative mass histories of nodes are assumed to be monotonic functions of time, and the formation time is instead found by a root finding algorithm,
[massBoundIsInactive] (default .false.) — Specifies whether or not the bound mass of the satellite component is inactive (i.e. does not appear in any ODE being solved).
[useLastIsolatedTime] (default .false.) — If true, evaluate the halo virial radius using a the virial density definition at the last isolated time of the halo.
[filterName] — The filter to select.
[filterType] — The filter type (rest or observed) to select.
[redshiftBand] — The redshift of the band (if not the output redshift).
[postprocessChain] — The postprocessing chain to use.
[cloudyTableFileName] (default var_str('%DATASTATICPATH%/hiiRegions/emissionLineLuminosities_BC2003_highResolution_imfChabrier.hdf5')) — The file of emission line luminosities to use.
[lineNames] — The emission lines to extract.
[component] — The component from which to extract star formation rate.
[toleranceRelative] (default 1.0d-3) — The relative tolerance used in integration over stellar population spectra.
[component] — The component from which to extract star formation rate.
[radiusCore] — The soliton core radius (in Mpc) characterizing the size of the quantum pressure-supported central core of the fuzzy dark matter halo; the density profile flattens inside this scale.
[densitySolitonCentral] — The central density (in \(\mathrm{M}_\odot\)/Mpc\(^3\)) of the solitonic core at \(r=0\), which sets the overall normalization of the density profile \(\rho(r) = \rho_\mathrm{c} [1+(r/r_c)^2]^{-8}\).
[toleranceRelativePotential] (default 1.0d-3) — The relative tolerance used in numerical ODE solutions for the gravitational potential of the solitonic core profile.
[dimensionless] (default .true.) — If true the soliton profile is treated as dimensionless (scale-free), allowing its radial and density quantities to be specified in arbitrary units.
[componentType] (default var_str('unknown')) — The galactic structure component type (e.g.dark matter halo, disk, spheroid) represented by this mass distribution, used for component-specific queries.
[massType] (default var_str('unknown')) — The mass type (e.g.dark matter, baryonic, total) represented by this mass distribution, used for mass-type-specific queries.
[radiusTransition] — The transition radius (in Mpc) at which the density profile smoothly switches from the halo profile to the accretion flow, controlled by the fourth-order transition function \(f_\mathrm{trans}(r)\).
[nonAnalyticSolver] (default var_str('fallThrough')) — Selects how solutions are computed when no analytic solution is available. If set to “fallThrough” then the solution ignoring heating is used, while if set to “numerical” then numerical solvers are used to find solutions.
[componentType] (default var_str('unknown')) — The component type that this mass distribution represents.
[massType] (default var_str('unknown')) — The mass type that this mass distribution represents.
[timeAge] — The age of the halo (in Gyr) since its formation, determining the total time available for SIDM self-interactions to thermalize the inner halo and produce an isothermal core.
[velocityRelativeMean] — Mean relative velocity to calculate self interaction cross section.
[nonAnalyticSolver] (default var_str('fallThrough')) — Selects how solutions are computed when no analytic solution is available. If set to “fallThrough” then the solution ignoring heating is used, while if set to “numerical” then numerical solvers are used to find solutions.
[componentType] (default var_str('unknown')) — The component type that this mass distribution represents.
[massType] (default var_str('unknown')) — The mass type that this mass distribution represents.
[nonAnalyticSolver] (default var_str('fallThrough')) — Selects how solutions are computed when no analytic solution is available. If set to “fallThrough” then the solution ignoring heating is used, while if set to “numerical” then numerical solvers are used to find solutions.
[componentType] (default var_str('unknown')) — The component type that this mass distribution represents.
[massType] (default var_str('unknown')) — The mass type that this mass distribution represents.
[tolerateVelocityMaximumFailure] (default .false.) — If true, tolerate failures to find the radius of the peak in the rotation curve.
[tolerateEnclosedMassIntegrationFailure] (default .false.) — If true, tolerate failures to find the mass enclosed as a function of radius.
[toleratePotentialIntegrationFailure] (default .false.) — If true, tolerate failures to compute the potential.
[fractionRadiusFinalSmall] (default 1.0d-3) — The initial radius is limited to be no smaller than this fraction of the final radius. This can help avoid problems in profiles that are extremely close to being disrupted.
[toleranceRelativePotential] (default 1.0d-3) — The maximum allowed relative tolerance to use in numerical solutions for the gravitational potential in dark-matter-only density profiles before aborting.
[lengthResolution] — The spatial resolution length scale (in Mpc) below which the underlying density profile is softened to a flat core, mimicking the finite force resolution of an N-body simulation.
[nonAnalyticSolver] (default var_str('fallThrough')) — Selects how solutions are computed when no analytic solution is available. If set to “fallThrough” then the solution ignoring heating is used, while if set to “numerical” then numerical solvers are used to find solutions.
[componentType] (default var_str('unknown')) — The component type that this mass distribution represents.
[massType] (default var_str('unknown')) — The mass type that this mass distribution represents.
[massMinimum] — The minimum halo mass (in \(\mathrm{M}_\odot\)) below which halos are excluded from the mass function histogram.
[massMaximum] — The maximum halo mass (in \(\mathrm{M}_\odot\)) above which halos are excluded from the mass function histogram.
[massCountPerDecade] — The number of logarithmic bins per decade of halo mass used when constructing the halo mass function.
[description] — A human-readable description of this mass function dataset, stored as metadata in the output file.
[simulationReference] — A bibliographic reference for the N-body simulation from which this mass function is derived, stored as metadata.
[simulationURL] — A URL pointing to the publicly accessible dataset or documentation for the N-body simulation, stored as metadata.
[bootstrapSampleCount] (default 30_c_size_t) — The number of bootstrap resamples of the particles that should be used.
[representativeMinimumCount] (default 10_c_size_t) — Minimum number of representative particles used to compute the center of a halo.
[tolerance] (default 1.0d-2) — The tolerance in the summed weight of bound particles which must be attained to declare convergence.
[bootstrapSampleRate] (default 1.0d0) — The sampling rate for particles.
[representativeFraction] (default 0.05d0) — Fraction of bound particles used to compute the center of a halo.
[analyzeAllParticles] (default .true.) — If true, all particles are assumed to be self-bound at the beginning of the analysis. Unbound particles at previous times are allowed to become bound in the current snapshot. If false and the self-bound information from the previous snapshot is available, only the particles that are self-bound at the previous snapshot are assumed to be bound at the beginning of the analysis.
[useVelocityMostBound] (default .false.) — If true, the velocity of the most bound particle in velocity space is used as the representative velocity of the satellite. If false, use the mass weighted mean velocity (center-of-mass velocity) of self-bound particles instead.
[orderRotation] (default var_str('none')) — The order in which evaluation of likelihoods should be rotated as a function of process number.
[logLikelihoodAccept] (default huge(0.0d0)) — The log-likelihood which should be “accepted”—once the log-likelihood reaches this value (or larger) no further updates to the chain will be made.
[report] (default .false.) — If true, report on the log-likelihood obtained.
[means] — The mean of the multivariate normal distribution.
[covariance] — The covariance matrix for the of the multivariate normal distribution.
[countForestsMaximum] (default -1_c_size_t) — If set to a positive number, this is the maximum number of forests that will be evolved.
[walltimeMaximum] (default -1_kind_int8) — If set to a positive number, this is the maximum wall time for which forest evolution is allowed to proceed before the task gives up.
[tolerateFailures] (default .false.) — If true then failures to evolve a forest are tolerated. The forest is evolved no further, but evolution of other forests continues.
[evolveForestsInParallel] (default .true.) — If true then each forest is evolved by a separate OpenMP thread. Otherwise, a single thread evolves all forests.
[suspendToRAM] (default .true.) — Specifies whether trees should be suspended to RAM (otherwise they are suspend to file).
[suspendPath] — The path to which tree suspension files will be stored.
[timeIntervalCheckpoint] (default -1_kind_int8) — If positive, gives the time in seconds between storing of checkpoint files. If zero or negative, no checkpointing is performed..
[fileNameCheckpoint] — The path to which checkpoint data will be stored.
[logM0] (default 10.0d0) — The parameter \(\log_{10} M_0\) (with \(M_0\) in units of \(\mathrm{M}_\odot\)) appearing in the star formation rate threshold expression for the star formation rate galactic filter class.
[logSFR0] (default 9.0d0) — The parameter \(\alpha_0\) appearing in the star formation rate threshold expression for the star formation rate galactic filter class.
[logSFR1] (default 0.0d0) — The parameter \(\alpha_1\) appearing in the star formation rate threshold expression for the star formation rate galactic filter class.
[cW] (default 3.78062835d0) — The parameter \(c_\mathrm{W}\) in the Bohr et al. (2021) power spectrum window function.
[beta] (default 3.4638743d0) — The parameter \(\beta\) in the Bohr et al. (2021) power spectrum window function.
[transferFunctionType] (default var_str('darkMatter')) — Specifies whether to use the darkMatter or total transfer function.
[fileName] — The name of the file from which to read a tabulated transfer function.
[redshift] (default 0.0d0) — The redshift of the transfer function to read.
[factorWavenumberSmoothExtrapolation] (default 0.0d0) — If positive, and extrapolation is used at high wavenumbers, the slope for extrapolation will be set by averaging over wavenumbers from \(k_\mathrm{max}/f\) to \(k_\mathrm{max}\), where \(f=\)[factorWavenumberSmoothExtrapolation] and \(k_\mathrm{max}\) is the highest wavenumber tabulated. This avoids spurious extrapolation for highly oscillatory transfer functions.
[acceptNegativeValues] (default .false.) — If true, negative values in the transfer function are allowed (and the absolute value is taken prior to interpolation). Otherwise, negative values result in an error.
[fractionalTimeStep] (default 0.01d0) — The fractional time step used when computing barrier crossing rates (i.e. the step used in finite difference calculations).
[fileName] (default var_str('none')) — The name of the file to/from which tabulations of barrier first crossing probabilities should be written/read. If set to “none” tables will not be stored.
[fractionalTimeStep] (default 0.01d0) — The fractional time step used when computing barrier crossing rates (i.e. the step used in finite difference calculations).
[varianceNumberPerUnitProbability] (default 1000) — The number of points to tabulate per unit variance for first crossing probabilities.
[varianceNumberPerUnit] (default 40) — The number of tabulation points per unit of \(\sigma^2\) used when building the rate look-up table for the Farahi excursion-set first-crossing distribution; higher values improve interpolation accuracy at the cost of memory and initialization time.
[varianceNumberPerDecade] (default 400) — The number of points to tabulate per decade of progenitor variance for first crossing rates.
[varianceNumberPerDecadeNonCrossing] (default 40) — The number of points to tabulate per decade of progenitor variance for non-crossing rates.
[timeNumberPerDecade] (default 10) — The number of tabulation points per decade of cosmic time used when building the first-crossing rate look-up table as a function of time; higher values improve temporal interpolation accuracy for rapidly evolving cosmologies.
[varianceIsUnlimited] (default .false.) — If true, the variance is assumed to have no upper limit (e.g. as in the case of CDM). This allows the tabulated solutions to be extended arbitrarily. Otherwise, tables are extended to encompass just the range of variance requested.
[linkingLength] (default 0.2d0) — The friends-of-friends linking length to use in computing virial density contrasts with the percolation analysis of More et al. (2011).

`nbodyOperatorConvexHullOverdensity`¶

An N-body data operator which computes the overdensity within the convex hull of the particles.

`nbodyOperatorConvexHullVolume`¶

An N-body data operator which computes the convex hull volume of the particles.

`nbodyOperatorDeleteProperties`¶

An N-body data operator which removes one or more named properties from the N-body simulation dataset, reducing memory usage or eliminating unwanted quantities before further processing. The list of properties to be deleted is specified as an array of names via the [names] parameter.

Parameters

[propertyNames] — A list of named properties to be deleted from the simulation.

`nbodyOperatorDistanceFromPoint`¶

An N-body data operator which computes and stores the Euclidean distance of each particle from a fixed reference point in three-dimensional Cartesian space, adding this as a new particle property. The reference point coordinates are specified as a three-element array via the [point] parameter.

Parameters

[point] — The Cartesian coordinates of the point from which to compute the distance.

`nbodyOperatorEnergyTensors`¶

An N-body data operator which determines the kinetic and Chandrasekhar potential energy tensors of a particle distribution, characterizing the dynamical state and shape of N-body halos. Self-bound particle restriction is controlled by [selfBoundOnly], and statistical uncertainties can be estimated using bootstrap resampling via [bootstrapSampleCount].

Parameters

[selfBoundParticlesOnly] — If true, the mean position and velocity are computed only for self-bound particles
[bootstrapSampleCount] (default 30_c_size_t) — The number of bootstrap resamples of the particles that should be used.

`nbodyOperatorEnvironmentalOverdensity`¶

An N-body data operator which measures the local environmental overdensity around each particle by counting neighbors within a spherical aperture and comparing to the mean simulation density. The aperture radius is set by [radiusEnvironment], periodic boundary conditions by [usePeriodic], and sampling efficiency by [sampleRate].

Parameters

[radiusSphere] — The radius of the sphere within which to measure environmental overdensity.
[densityParticleMean] — The mean density of particles in the simulation.
[sampleRate] (default 1_c_size_t) — One in [sampleRate] particles will be sampled when computed environmental overdensities.
[lengthBox] (default 0.0d0) — The length of the periodic box.
[periodic] (default .false.) — If true, periodic boundary conditions will be used.

`nbodyOperatorExportIRATE`¶

An N-body data operator which exports N-body simulation data to the IRATE (IRvine Astrophysical simulaTion structurE) HDF5 format for interoperability with other analysis tools. The output file path is set by [fileName], with the snapshot index and redshift specified by [snapshot] and [redshift].

Methods

readHalos — Read a snapshot from the IRATE format file.
writeHalos — Write a snapshot to a IRATE format file.
readSimulation — Read the requested properties of the simulation from an IRATE format file.
writeSimulation — Write the requested properties of the simulation from an IRATE format file.
copySimulation — Copy “SimulationProperties” group from one IRATE file to another.
copyCosmology — Copy “Cosmology” group from one IRATE file to another.

Parameters

[fileName] — The name of the file to which data should be exported.
[snapshot] — The snapshot index of the data.
[redshift] — The redshift of the data.
[fileName] — The name of the file to read.
[snapshot] — The snapshot number to read.
[label] (default var_str('primary')) — A label for the simulation.
[properties] — Properties to read from the simulation.

`nbodyOperatorFilterBox`¶

An N-body data operator which retains only particles whose positions fall within a specified axis-aligned cuboid region, enabling spatially localized analysis of simulation volumes. The lower and upper boundaries of the cuboid are specified as three-element coordinate arrays via [cornerLower] and [cornerUpper].

Parameters

[boundLow] — The lower boundary of the cuboid region.
[boundHigh] — The upper boundary of the cuboid region.

`nbodyOperatorFilterConvexHull`¶

An N-body data operator which retains only particles whose positions fall inside the convex hull of another simulation dataset, enabling filtering to an irregular survey geometry or simulation subvolume boundary. The convex hull is taken from the simulation identified by the [simulationIndex] parameter.

Parameters

[hullFromSimulation] — The index of the simulation from which to take the convex hull.

`nbodyOperatorFilterID`¶

An N-body data operator which retains only the subset of particles whose unique IDs match a specified list, enabling targeted analysis of individual halos or particle subsets. The IDs to retain can be specified directly via [IDs] or loaded from a file using [fileName].

Parameters

[idSelection] — The IDs of particles to retain.
[idSelectionFileName] — The name of a file containing the IDs of particles to retain.
[idSelectionDatasetName] (default var_str('id')) — The name of the dataset containing the IDs of particles to retain.

`nbodyOperatorFilterProperties`¶

An N-body data operator which retains only particles whose named property values fall within specified minimum and maximum bounds, enabling selection by mass, concentration, spin, or any other stored property. Property names and their lower and upper limits are specified as arrays via [names], [minimum], and [maximum].

Parameters

[propertyNames] — A list of named properties on which to filter.
[rangeLow] — The lowest value of each property to pass (”-infinity” is interpreted as the lowest possible value for the property).
[rangeHigh] — The highest value of each property to pass (”+infinity” is interpreted as the highest possible value for the property).

`nbodyOperatorFilterUncontaminatedSphere`¶

An N-body data operator which filters particles to select those within a sphere such that the contamination by particles of non-preferred type is below a specified level.

Parameters

[point] — The point at which to center the filtered sphere.
[fractionContamination] — The contamination fraction allowed within the sphere.
[particleType] — The preferred particle type for filtering.
[massWeighted] — If true, contamination is weighted by particle mass, otherwise it is unweighted (i.e. depends on the number of particles only).

`nbodyOperatorFlagAlwaysIsolated`¶

An N-body data operator which identifies and flags particles that have remained isolated throughout the simulation, never having been part of a larger halo, by tracking their merger history. A particle is considered always-isolated unless its virial mass increased by more than the factor specified by [massGrowthFactor], which ignores transient accretion events.

Parameters

[massFactor] (default 2.0d0) — The factor by which virial mass must increase for previous non-isolation to be ignored.

`nbodyOperatorHaloCrossingTime`¶

An N-body data operator which computes and stores the crossing time of halos defined as \(t_\mathrm{cross} = 2 r_\mathrm{vir}/V_\mathrm{vir}\).

`nbodyOperatorHostedRootID`¶

An N-body data operator which determines the ID of the topmost root halo in the host hierarchy by iteratively following host links upward from each particle. The errorIfNotFound parameter controls whether a missing host causes a fatal error or instead assigns a sentinel value of \(-1\).

Parameters

[missingHalosAreFatal] (default .true.) — If true, if a halo is not found during the search through hosts and descendants then a fatal error occurs. Otherwise, such missing halos are ignored, and a hostedRootID value of \(-1\) is assigned to the particle.

`nbodyOperatorIdentifyFlybysMansfieldKravtsov2020`¶

An N-body data operator which identifies flyby halos using the algorithm of Mansfield and Kravtsov (2020), classifying halos that pass through a larger host without merging as flybys. The errorIfNotFound parameter controls how missing host halos are handled during the flyby search.

Parameters

[missingHostsAreFatal] — If true, missing hosts cause a fatal error. Otherwise, a missing host causes the halo to be ignored in flyby analysis of its descendants.

`nbodyOperatorInertiaTensor`¶

An N-body data operator which computes the inertia tensor eigenvalues, eigenvectors, and resulting axis ratios of an N-body halo, characterizing its three-dimensional shape. Parameters control whether only self-bound particles are used, the maximum radius for inclusion, and the number of bootstrap resamples.

Parameters

[selfBoundParticlesOnly] — If true, the maximum velocity is computed only for self-bound particles
[bootstrapSampleCount] (default 30_c_size_t) — The number of bootstrap resamples of the particles that should be used.
[radiusMaximum] — The maximum radius (relative to the position of the most-bound particle) to include when computing the inertia tensor.

`nbodyOperatorMassFunction`¶

An N-body data operator which computes the halo mass function from N-body simulation data, binning halos logarithmically in mass. Parameters specify the minimum and maximum halo mass, the number of bins per decade, and metadata such as a description, reference, and URL for the simulation.

Methods

tabulate — Tabulate the virial density contrast as a function of mass and time.
restoreTable — Restore a tabulated solution from file.
storeTable — Store a tabulated solution to file.

Parameters

[velocityCharacteristic] (default 250.0d0) — The velocity scale at which the SNe-driven outflow rate equals the star formation rate in disks.
[exponent] (default 3.5d0) — The velocity scaling of the SNe-driven outflow rate in disks.
[fraction] (default 0.01d0) — The normalization \(f\) of the outflow rate relative to the star formation rate at a reference halo velocity of 200 km/s and expansion factor of 1, setting the overall mass-loading amplitude of the halo-scaling feedback model.
[exponentVelocity] (default -2.0d0) — The exponent of virial velocity in the outflow rate in disks.
[exponentRedshift] (default 0.0d0) — The power-law exponent of the cosmological expansion factor \((1+z)\) in the halo-scaling outflow rate, allowing the mass-loading factor to evolve with redshift; a value of zero gives no redshift evolution.
[toleranceRelativeVelocityDispersion] (default 1.0d-6) — The relative tolerance to use in numerical solutions for the velocity dispersion in dark-matter-only density profiles.
[toleranceRelativeVelocityDispersionMaximum] (default 1.0d-3) — The maximum relative tolerance to use in numerical solutions for the velocity dispersion in dark-matter-only density profiles.
[radiusNormalization] (default 3.3d-6) — The initial value appearing in the radius-mass relation
[toleranceAbsoluteMass] (default 1.0d-6) — The mass tolerance used to judge whether the nuclear star cluster is physically plausible.
[toleranceRelativeMetallicity] (default 1.0d-4) — The metallicity tolerance for ODE solution.
[inactiveLuminositiesStellar] (default .false.) — Specifies whether or not nuclear star cluster stellar luminosities are inactive properties (i.e. do not appear in any ODE being solved).
[scaleRelativeMass] (default 1.0d-2) — The mass scale, relative to the total mass of the node, below which calculations in the delayed very simple hot halo component are allowed to become inaccurate.
[starveSatellites] (default .false.) — Specifies whether or not the hot halo should be removed (“starved”) when a node becomes a satellite.
[starveSatellitesOutflowed] (default .false.) — Specifies whether or not the outflowed hot halo should be removed (“starved”) when a node becomes a satellite.
[outflowReturnOnFormation] (default .false.) — Specifies whether or not outflowed gas should be returned to the hot reservoir on halo formation events.
[angularMomentumAlwaysGrows] (default .false.) — Specifies whether or not negative rates of accretion of angular momentum into the hot halo will be treated as positive for the purposes of computing the hot halo angular momentum.
[fractionBaryonLimitInNodeMerger] (default .false.) — Controls whether the hot gas content of nodes should be limited to not exceed the universal baryon fraction at node merger events. If set to true, hot gas (and angular momentum, abundances, and chemicals proportionally) will be removed from the merged halo to the unaccreted gas reservoir to limit the baryonic mass to the universal baryon fraction where possible.
[scaleAbsoluteMass] (default 100.0d0) — The absolute mass scale below which calculations in the very simple disk component are allowed to become inaccurate.
[toleranceAbsoluteMass] (default 1.0d-6) — The mass tolerance used to judge whether the disk is physically plausible.
[toleranceAbsoluteMass] (default 1.0d-6) — The mass tolerance used to judge whether the disk is physically plausible.
[toleranceRelativeMetallicity] (default 1.0d-4) — The metallicity tolerance for ODE solution.
[radiusStructureSolver] (default 1.0d0) — The radius (in units of the standard scale length) to use in solving for the size of the disk.
[structureSolverUseCole2000Method] (default .false.) — If true, use the method described in Cole et al. (2000) to correct for difference between thin disk and spherical mass distributions when solving for disk radii.
[diskNegativeAngularMomentumAllowed] (default .true.) — Specifies whether or not negative angular momentum is allowed for the disk.
[inactiveLuminositiesStellar] (default .false.) — Specifies whether or not disk stellar luminosities are inactive properties (i.e. do not appear in any ODE being solved).
[postStepZeroNegativeMasses] (default .true.) — If true, negative masses will be zeroed after each ODE step. Note that this can lead to non-conservation of mass.
[ratioAngularMomentumSolverRadius] (default ratioAngularMomentumSolverRadiusDefault) — The assumed ratio of the specific angular momentum at the structure solver radius to the mean specific angular momentum of the standard disk component.
[scaleAbsoluteMass] (default 100.0d0) — The absolute mass scale below which calculations in the very simple spheroid component are allowed to become inaccurate.
[toleranceAbsoluteMass] (default 1.0d-6) — The mass tolerance used to judge whether the spheroid is physically plausible.
[efficiencyEnergeticOutflow] (default 1.0d-2) — The proportionality factor relating mass outflow rate from the spheroid to the energy input rate divided by \(V_\mathrm{spheroid}^2\).
[toleranceRelativeMetallicity] (default 1.0d-4) — The metallicity tolerance for ODE solution.
[toleranceAbsoluteMass] (default 1.0d-6) — The mass tolerance used to judge whether the spheroid is physically plausible.
[inactiveLuminositiesStellar] (default .false.) — Specifies whether or not spheroid stellar luminosities are inactive properties (i.e. do not appear in any ODE being solved).
[postStepZeroNegativeMasses] (default .true.) — If true, negative masses will be zeroed after each ODE step. Note that this can lead to non-conservation of mass.
[ratioAngularMomentumScaleRadius] (default ratioAngularMomentumScaleRadiusDefault) — The assumed ratio of the specific angular momentum at the scale radius to the mean specific angular momentum of the standard spheroid component.
[outputMergers] (default .false.) — Determines whether or not properties of black hole mergers will be output.
[fileNames] — The name of the file(s) from which merger tree data should be read when using the [mergerTreeConstruct]\(=\)read tree construction method.
[forestSizeMaximum] (default 0_c_size_t) — The maximum number of nodes allowed in a forest before it will be broken up into trees and processed individually. A value of 0 implies that forests should never be split.
[presetMergerTimes] (default .true.) — Specifies whether merging times for subhalos should be preset when reading merger trees from a file.
[presetMergerNodes] (default .true.) — Specifies whether the target nodes for mergers should be preset (i.e. determined from descendant nodes). If they are not, merging will be with each satellite’s host node.
[presetSubhaloMasses] (default .true.) — Specifies whether subhalo mass should be preset when reading merger trees from a file.
[subhaloAngularMomentaMethod] (default var_str('summation')) — Specifies how to account for subhalo angular momentum when adding subhalo mass to host halo mass.
[presetSubhaloIndices] (default .true.) — Specifies whether subhalo indices should be preset when reading merger trees from a file.
[presetPositions] (default .true.) — Specifies whether node positions should be preset when reading merger trees from a file.
[presetScaleRadii] (default .true.) — Specifies whether node scale radii should be preset when reading merger trees from a file.
[scaleRadiiFailureIsFatal] (default .true.) — Specifies whether failure to set a node scale radii should be regarded as a fatal error. (If not, a fallback method to set scale radius is used in such cases.)
[presetScaleRadiiConcentrationMinimum] (default 3.0d0) — The lowest concentration (\(c=r_\mathrm{vir}/r_\mathrm{s}\)) allowed when setting scale radii, \(r_\mathrm{s}\).
[presetScaleRadiiConcentrationMaximum] (default 60.0d0) — The largest concentration (\(c=r_\mathrm{vir}/r_\mathrm{s}\)) allowed when setting scale radii, \(r_\mathrm{s}\).
[presetScaleRadiiMinimumMass] (default 0.0d0) — The minimum halo mass for which scale radii should be preset (if [presetScaleRadii]\(=\)true).
[presetUnphysicalAngularMomenta] (default .false.) — When reading merger trees from file and presetting halo angular momenta, detect unphysical (<=0) angular momenta and preset them using the selected halo spin method.
[presetAngularMomenta] (default .true.) — Specifies whether node angular momenta should be preset when reading merger trees from a file.
[presetAngularMomenta3D] (default .false.) — Specifies whether node 3-D angular momenta vectors should be preset when reading merger trees from a file.
[presetOrbits] (default .true.) — Specifies whether node orbits should be preset when reading merger trees from a file.
[presetOrbitsSetAll] (default .true.) — Forces all orbits to be set. If the computed orbit does not cross the virial radius, then select one at random instead.
[presetOrbitsAssertAllSet] (default .true.) — Asserts that all virial orbits must be preset. If any can not be set, Galacticus will stop.
[presetOrbitsBoundOnly] (default .true.) — Specifies whether only bound node orbits should be set.
[beginAt] (default -1_kind_int8) — Specifies the index of the tree to begin at. (Use -1 to always begin with the first tree.)
[outputTimeSnapTolerance] (default 0.0d0) — The relative tolerance required to “snap” a node time to the closest output time.
[missingHostsAreFatal] (default .true.) — Specifies whether nodes with missing host nodes should be considered to be fatal—see the discussion of missing host nodes in the class description above.
[treeIndexToRootNodeIndex] (default .false.) — Specifies whether tree indices should always be set to the index of their root node.
[allowBranchJumps] (default .true.) — Specifies whether nodes are allowed to jump between branches.
[allowSubhaloPromotions] (default .true.) — Specifies whether subhalos are permitted to be promoted to being isolated halos.
[alwaysPromoteMostMassive] (default .false.) — If true, the most massive progenitor is always promoted to be the primary progenitor even if it is a subhalo. Otherwise, isolated progenitors are given priority over subhalo progenitors, even if they are less massive.
[presetNamedReals] — Names of real datasets to be additionally read and stored in the nodes of the merger tree when using the [mergerTreeConstruct]\(=\)read tree construction method.
[presetNamedIntegers] — Names of integer datasets to be additionally read and stored in the nodes of the merger tree when using the [mergerTreeConstruct]\(=\)read tree construction method.
[fatalMismatches] (default .true.) — Specifies whether mismatches in cosmological parameter values between Galacticus and “Sussing Merger Trees” format (Srisawat et al., 2013) merger tree files should be considered fatal.
[fatalNonTreeNode] (default .true.) — Specifies whether nodes in snapshot files but not in the merger tree file should be considered fatal when importing from the “Sussing Merger Trees” format (Srisawat et al., 2013).
[subvolumeCount] (default 1) — Specifies the number of subvolumes along each axis into which a “Sussing Merger Trees” format (Srisawat et al., 2013) merger tree files should be split for processing through Galacticus.
[subvolumeBuffer] (default 0.0d0) — Specifies the buffer region (in units of Mpc\(/h\) to follow the format convention) around subvolumes of a “Sussing Merger Trees” format (Srisawat et al., 2013) merger tree file which should be read in to ensure that no halos are missed from trees.
[subvolumeIndex] (default [0,0,0]) — Specifies the index (in each dimension) of the subvolume of a “Sussing Merger Trees” format (Srisawat et al., 2013) merger tree file to process. Indices range from 0 to [subvolumeCount]\(-1\).
[badValue] (default -0.5d0) — Use for bad value detection in “Sussing” merger trees. Values for scale radius and halo spin which exceed this threshold are assumed to be bad.
[badValueTest] (default var_str('lessThan')) — Use for bad value detection in “Sussing” merger trees. Values which exceed the threshold in ths specified direction are assumed to be bad.
[treeSampleRate] (default 1.0d0) — Specify the probability that any given tree should processed (to permit subsampling).
[massOptions] (default var_str('default')) — Mass option for Sussing merger trees.
[mergeProbability] (default 0.1d0) — The largest probability of branching allowed in a timestep in merger trees built by the Cole et al. (2000) method.
[accretionLimit] (default 0.1d0) — The largest fractional mass change due to subresolution accretion allowed in a timestep in merger trees built by the Cole et al. (2000) method.
[redshiftMaximum] (default 1.0d5) — The highest redshift to which merger trees will be built in the Cole et al. (2000) method.
[toleranceTimeEarliest] (default 2.0d-6) — The fractional tolerance used to judge if a branch is at the earliest allowed time in the tree.
[branchIntervalStep] (default .true.) — If false use the original Cole et al. (2000) method to determine whether branching occurs in a timestep. If true draw branching intervals from a negative exponential distribution.
[toleranceResolutionSelf] (default 1.0d-6) — The fractional tolerance in node mass at the resolution limit below which branch mis-orderings will be ignored.
[toleranceResolutionParent] (default 1.0d-3) — The fractional tolerance in parent node mass at the resolution limit below which branch mis-orderings will be ignored.
[ignoreNoProgress] (default .false.) — If true, failure to make progress on a branch will be ignored (and the branch terminated).
[ignoreWellOrdering] (default .false.) — If true, non-well-ordered tree branches are pruned away instead of causing errors..
[redshiftBase] (default 0.0d0) — The redshift at which to plant the base node when building merger trees.
[timeSnapTolerance] (default 1.0d-6) — The fractional tolerance within which the tree base time will be snapped to a nearby output time.
[treeBeginAt] (default 0) — The index (in order of increasing base halo mass) of the tree at which to begin when building merger trees. A value of “0” means to begin with tree number 1 (if processing trees in ascending order), or equal to the number of trees (otherwise).
[processDescending] (default .true.) — If true, causes merger trees to be processed in order of decreasing mass.
[splitTrees] (default .false.) — If true, prune away any nodes of the tree that are not needed to determine evolution up to the latest time at which a node is present inside the lightcone. This typically leads to a tree splitting into a forest of trees.
[label] — A label for the mass function.
[comment] — A descriptive comment for the mass function.
[starFormationRates] — The star formation rates corresponding to bin centers.
[covarianceBinomialBinsPerDecade] (default 10) — The number of bins per decade of star formation rate to use when constructing star formation rate function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMinimum] (default 1.0d10) — The star formation rate to consider when constructing star formation rate function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMaximum] (default 1.0d12) — The maximum star formation rate to consider when constructing star formation rate function covariance matrices for main branch galaxies.
[targetLabel] — Label for the target dataset.
[functionValueTarget] — The target function for likelihood calculations.
[functionCovarianceTarget] — The target function covariance for likelihood calculations.
[likelihoodBins] — Controls which bins in the stellar mass–halo mass relation will be used in computing the likelihood: * not present: all bins are included in the likelihood calculation; * list of integers: use only the mass bin(s) given in this list in the likelihood calculation; * auto: use only bins which have a non-zero number of halos contributing to them in the likelihood calculation.
[fileNameTarget] — The name of the file containing the target data.
[redshiftInterval] (default 1) — The redshift interval to use.
[likelihoodNormalize] (default .false.) — If true, then normalize the likelihood to make it a probability density.
[computeScatter] (default .false.) — If true, the scatter in log10(stellar mass) is computed. Otherwise, the mean is computed.
[systematicErrorPolynomialCoefficient] (default [0.0d0]) — The coefficients of the systematic error polynomial for stellar mass in the stellar vs halo mass relation.
[systematicErrorMassHaloPolynomialCoefficient] (default [0.0d0]) — The coefficients of the systematic error polynomial for halo mass in the stellar vs halo mass relation.
[errorTolerant] (default .false.) — Error tolerance for the N-body spin distribution operator.
[logNormalRange] (default 100.0d0) — The multiplicative range of the log-normal distribution used to model the distribution of the mass and energy terms in the spin parameter. Specifically, the lognormal distribution is truncated outside the range \((\lambda_\mathrm{m}/R,\lambda_\mathrm{m} R\), where \(\lambda_\mathrm{m}\) is the measured spin, and \(R=\)[logNormalRange]
[fileName] — The name of the file from which to read spin distribution function parameters.
[comment] — A comment describing this analysis.
[label] — A label for this analysis.
[label] — A label for the spin distribution function.
[comment] — A descriptive comment for the spin distribution function.
[redshift] — The redshift at which to compute the spin distribution function.
[massMinimum] — Minimum halo mass for the spin distribution function.
[massMaximum] — Maximum halo mass for the spin distribution function.
[spinMinimum] — Minimum spin for the spin distribution function.
[spinMaximum] — Maximum spin for the spin distribution function.
[countSpinsPerDecade] — Number of spins per decade at which to compute the spin distribution function.
[timeRecent] — Halos which experienced a major node merger within a time \(\Delta t=\)[timeRecent] of the analysis time will be excluded from the analysis.
[particleCountMinimum] — The minimum particle count to assume when computing N-body errors on spins.
[massParticle] — The mass of the particle used in the N-body simulation from which spins were measured.
[energyEstimateParticleCountMaximum] — The maximum number of particles used in estimating halo energies when measuring spins from the N-body simulation.
[targetLabel] — Label for the target dataset.
[functionValueTarget] — The target function for likelihood calculations.
[functionCovarianceTarget] — The target function covariance for likelihood calculations.
[label] — A label for the mass function.
[comment] — A descriptive comment for the mass function.
[masses] — The masses corresponding to bin centers.
[covarianceBinomialBinsPerDecade] (default 10) — The number of bins per decade of halo mass to use when constructing HI mass function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMinimum] (default 1.0d8) — The minimum halo mass to consider when constructing HI mass function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMaximum] (default 1.0d16) — The maximum halo mass to consider when constructing HI mass function covariance matrices for main branch galaxies.
[targetLabel] — Label for the target dataset.
[functionValueTarget] — The target function for likelihood calculations.
[functionCovarianceTarget] — The target function covariance for likelihood calculations.
[label] — A label for the luminosity function.
[comment] — A descriptive comment for the luminosity function.
[magnitudesAbsolute] — The absolute magnitudes corresponding to bin centers.
[covarianceBinomialBinsPerDecade] (default 10) — The number of bins per decade of halo mass to use when constructing luminosity function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMinimum] (default 1.0d8) — The minimum halo mass to consider when constructing luminosity function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMaximum] (default 1.0d16) — The maximum halo mass to consider when constructing luminosity function covariance matrices for main branch galaxies.
[targetLabel] — Label for the target dataset.
[functionValueTarget] — The target function for likelihood calculations.
[functionCovarianceTarget] — The target function covariance for likelihood calculations.
[label] — A label for the luminosity function.
[comment] — A descriptive comment for the luminosity function.
[luminosities] — The luminosities corresponding to bin centers.
[covarianceBinomialBinsPerDecade] (default 10) — The number of bins per decade of halo mass to use when constructing luminosity function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMinimum] (default 1.0d8) — The minimum halo mass to consider when constructing luminosity function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMaximum] (default 1.0d16) — The maximum halo mass to consider when constructing luminosity function covariance matrices for main branch galaxies.
[includeNitrogenII] (default .false.) — If true, include contamination by the [NII] (6548AA \(+\) 6584AA) doublet.
[depthOpticalISMCoefficient] (default 1.0d0) — Multiplicative coefficient for optical depth in the ISM.
[targetLabel] — Label for the target dataset.
[functionValueTarget] — The target function for likelihood calculations.
[functionCovarianceTarget] — The target function covariance for likelihood calculations.
[label] — A label for the mass function.
[comment] — A descriptive comment for the mass function.
[masses] — The masses corresponding to bin centers.
[covarianceBinomialBinsPerDecade] (default 10) — The number of bins per decade of halo mass to use when constructing stellar mass function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMinimum] (default 1.0d8) — The minimum halo mass to consider when constructing stellar mass function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMaximum] (default 1.0d16) — The maximum halo mass to consider when constructing stellar mass function covariance matrices for main branch galaxies.
[targetLabel] — Label for the target dataset.
[functionValueTarget] — The target function for likelihood calculations.
[functionCovarianceTarget] — The target function covariance for likelihood calculations.
[rootVarianceFractionalMinimum] (default 0.0d0) — The minimum fractional root variance (relative to the target dataset).
[fileName] — The name of the file from which to read concentration distribution function parameters.
[comment] — A comment describing this analysis.
[label] — A label for this analysis.
[label] — A label for the concentration distribution function.
[comment] — A descriptive comment for the concentration distribution function.
[redshift] — The redshift at which to compute the concentration distribution function.
[massMinimum] — Minimum halo mass for the concentration distribution function.
[massMaximum] — Maximum halo mass for the concentration distribution function.
[concentrationMinimum] — Minimum concentration for the concentration distribution function.
[concentrationMaximum] — Maximum concentration for the concentration distribution function.
[countConcentrationsPerDecade] — Number of concentrations per decade at which to compute the concentration distribution function.
[timeRecent] — Halos which experienced a major node merger within a time \(\Delta t=\)[timeRecent] of the analysis time will be excluded from the analysis.
[massParticle] — The particle mass in the source N-body simulation.
[targetLabel] — Label for the target dataset.
[functionValueTarget] — The target function for likelihood calculations.
[functionCovarianceTarget] — The target function covariance for likelihood calculations.
[fileName] — The name of the file from which to read star forming main sequence function parameters.
[comment] — A comment describing this analysis.
[label] — A label for this analysis.
[label] — A label for the star forming main sequence function.
[comment] — A descriptive comment for the star forming main sequence function.
[massMinimum] — Minimum stellar mass for the star forming main sequence function.
[massMaximum] — Maximum stellar mass for the star forming main sequence function.
[countMassesPerDecade] — Number of masses per decade at which to compute the star forming main sequence function.
[targetLabel] — Label for the target dataset.
[meanValueTarget] — The target function for likelihood calculations.
[meanCovarianceTarget] — The target function covariance for likelihood calculations.
[label] — A label for the mass function.
[comment] — A descriptive comment for the mass function.
[separations] — The separations corresponding to bin centers.
[massMinima] — The minimum mass of each mass sample.
[massMaxima] — The maximum mass of each mass sample.
[massHaloBinsPerDecade] (default 10) — The number of bins per decade of halo mass to use when constructing the mass function covariance matrix for main branch galaxies.
[massHaloMinimum] (default 1.0d8) — The minimum halo mass to consider when constructing the mass function covariance matrix for main branch galaxies.
[massHaloMaximum] (default 1.0d16) — The maximum halo mass to consider when constructing the mass function covariance matrix for main branch galaxies.
[wavenumberCount] (default 60_c_size_t) — The number of bins in wavenumber to use in computing the correlation function.
[wavenumberMinimum] (default 1.0d-3) — The minimum wavenumber to use when computing the correlation function.
[wavenumberMaximum] (default 1.0d4) — The maximum wavenumber to use when computing the correlation function.
[integralConstraint] — The integral constraint for these correlation functions.
[depthLineOfSight] — The line-of-sight depth over which the correlation function was projected.
[halfIntegral] — Set to true if the projection integrand should be over line-of-sight depths greater than zero.
[binnedProjectedCorrelationTarget] — The target function for likelihood calculations.
[binnedProjectedCorrelationCovarianceTarget] — The target function covariance for likelihood calculations.
[targetLabel] (default var_str('')) — A label for the target dataset in a plot of this analysis.
[starFormationRateSpecificQuiescentLogarithmic] — The base-10 logarithm specific star formation rate (in units of Gyr\(^{-1}\)) separating quiescent and star-forming galaxies.
[starFormationRateSpecificLogarithmicError] — The observational fractional error in specific star formation rate (in units of dex) of galaxies.
[fileName] — The name of the file from which to read quiescent fraction function parameters.
[comment] — A comment describing this analysis.
[label] — A label for this analysis.
[label] — A label for the star forming main sequence function.
[comment] — A descriptive comment for the star forming main sequence function.
[massMinimum] — Minimum stellar mass for the star forming main sequence function.
[massMaximum] — Maximum stellar mass for the star forming main sequence function.
[countMassesPerDecade] — Number of masses per decade at which to compute the star forming main sequence function.
[targetLabel] — Label for the target dataset.
[meanValueTarget] — The target function for likelihood calculations.
[meanCovarianceTarget] — The target function covariance for likelihood calculations.
[nonAnalyticSolver] (default var_str('fallThrough')) — Selects how solutions are computed when no analytic solution is available. If set to “fallThrough” then the solution ignoring heating is used, while if set to “numerical” then numerical solvers are used to find solutions.
[radiusFractionalTruncateMinimum] (default 2.0d0) — The minimum radius (in units of the virial radius) to begin truncating the density profile.
[radiusFractionalTruncateMaximum] (default 4.0d0) — The maximum radius (in units of the virial radius) to finish truncating the density profile.
[nonAnalyticSolver] (default var_str('fallThrough')) — Selects how solutions are computed when no analytic solution is available. If set to “fallThrough” then the solution ignoring heating is used, while if set to “numerical” then numerical solvers are used to find solutions.
[velocityDispersionApproximate] (default .true.) — If true, radial velocity dispersion is computed using an approximate method in which we assume that \(\sigma_\mathrm{r}^2(r) \rightarrow \sigma_\mathrm{r}^2(r) - (2/3) \epsilon(r)\), where \(\epsilon(r)\) is the specific heating energy. If false then radial velocity dispersion is computed by numerically solving the Jeans equation.
[tolerateEnclosedMassIntegrationFailure] (default .false.) — If true, tolerate failures to find the mass enclosed as a function of radius.
[tolerateVelocityDispersionFailure] (default .false.) — If true, tolerate failures to compute the velocity dispersion.
[tolerateVelocityMaximumFailure] (default .false.) — If true, tolerate failures to find the radius of the maximum circular velocity.
[toleratePotentialIntegrationFailure] (default .false.) — If true, tolerate numerical failures when computing the gravitational potential of a heated dark matter profile, allowing the calculation to continue with a fallback result rather than aborting.
[toleranceRelativeVelocityDispersion] (default 1.0d-6) — The relative tolerance to use in numerical solutions for the velocity dispersion in dark-matter-only density profiles.
[toleranceRelativeVelocityDispersionMaximum] (default 1.0d-3) — The maximum relative tolerance to use in numerical solutions for the velocity dispersion in dark-matter-only density profiles.
[fractionRadiusFinalSmall] (default 1.0d-3) — The initial radius is limited to be no smaller than this fraction of the final radius. This can help avoid problems in profiles that are extremely close to being disrupted.
[toleranceRelativePotential] (default 1.0d-3) — The maximum allowed relative tolerance to use in numerical solutions for the gravitational potential in dark-matter-only density profiles before aborting.
[tolerateVelocityMaximumFailure] (default .true.) — If true, tolerate failures to find the radius of the peak in the rotation curve.
[lengthResolution] — The gravitational softening length \(\Delta x\) (in Mpc) of the N-body simulation, which sets the minimum spatial scale below which the dark matter profile is smoothed to avoid artificial two-body effects.
[massResolution] — The mass resolution \(\Delta M\) (in \(\mathrm{M}_\odot\)) of the N-body simulation, representing the minimum halo mass that can be resolved; profiles of halos near this limit are softened to account for particle discreteness effects.
[resolutionIsComoving] — If true, the resolution length is assumed to be fixed in comoving coordinates, otherwise in physical coordinates.
[nonAnalyticSolver] (default var_str('fallThrough')) — Selects how solutions are computed when no analytic solution is available. If set to “fallThrough” then the solution ignoring heating is used, while if set to “numerical” then numerical solvers are used to find solutions.
[C] (default 400.0d0) — The parameter \(C\) appearing in the halo concentration algorithm of Ludlow et al. (2016).
[f] (default 0.02d0) — The parameter \(f\) appearing in the halo concentration algorithm of Ludlow et al. (2016).
[timeFormationSeekDelta] (default 0.0d0) — The parameter \(\Delta \log t\) by which the logarithm of the trial formation time is incremented when stepping through the formation history of a node to find the formation time. If set to zero (or a negative value) the cumulative mass histories of nodes are assumed to be monotonic functions of time, and the formation time is instead found by a root finding algorithm,
[massBoundIsInactive] (default .false.) — Specifies whether or not the bound mass of the satellite component is inactive (i.e. does not appear in any ODE being solved).
[useLastIsolatedTime] (default .false.) — If true, evaluate the halo virial radius using a the virial density definition at the last isolated time of the halo.
[filterName] — The filter to select.
[filterType] — The filter type (rest or observed) to select.
[redshiftBand] — The redshift of the band (if not the output redshift).
[postprocessChain] — The postprocessing chain to use.
[cloudyTableFileName] (default var_str('%DATASTATICPATH%/hiiRegions/emissionLineLuminosities_BC2003_highResolution_imfChabrier.hdf5')) — The file of emission line luminosities to use.
[lineNames] — The emission lines to extract.
[component] — The component from which to extract star formation rate.
[toleranceRelative] (default 1.0d-3) — The relative tolerance used in integration over stellar population spectra.
[component] — The component from which to extract star formation rate.
[radiusCore] — The soliton core radius (in Mpc) characterizing the size of the quantum pressure-supported central core of the fuzzy dark matter halo; the density profile flattens inside this scale.
[densitySolitonCentral] — The central density (in \(\mathrm{M}_\odot\)/Mpc\(^3\)) of the solitonic core at \(r=0\), which sets the overall normalization of the density profile \(\rho(r) = \rho_\mathrm{c} [1+(r/r_c)^2]^{-8}\).
[toleranceRelativePotential] (default 1.0d-3) — The relative tolerance used in numerical ODE solutions for the gravitational potential of the solitonic core profile.
[dimensionless] (default .true.) — If true the soliton profile is treated as dimensionless (scale-free), allowing its radial and density quantities to be specified in arbitrary units.
[componentType] (default var_str('unknown')) — The galactic structure component type (e.g.dark matter halo, disk, spheroid) represented by this mass distribution, used for component-specific queries.
[massType] (default var_str('unknown')) — The mass type (e.g.dark matter, baryonic, total) represented by this mass distribution, used for mass-type-specific queries.
[radiusTransition] — The transition radius (in Mpc) at which the density profile smoothly switches from the halo profile to the accretion flow, controlled by the fourth-order transition function \(f_\mathrm{trans}(r)\).
[nonAnalyticSolver] (default var_str('fallThrough')) — Selects how solutions are computed when no analytic solution is available. If set to “fallThrough” then the solution ignoring heating is used, while if set to “numerical” then numerical solvers are used to find solutions.
[componentType] (default var_str('unknown')) — The component type that this mass distribution represents.
[massType] (default var_str('unknown')) — The mass type that this mass distribution represents.
[timeAge] — The age of the halo (in Gyr) since its formation, determining the total time available for SIDM self-interactions to thermalize the inner halo and produce an isothermal core.
[velocityRelativeMean] — Mean relative velocity to calculate self interaction cross section.
[nonAnalyticSolver] (default var_str('fallThrough')) — Selects how solutions are computed when no analytic solution is available. If set to “fallThrough” then the solution ignoring heating is used, while if set to “numerical” then numerical solvers are used to find solutions.
[componentType] (default var_str('unknown')) — The component type that this mass distribution represents.
[massType] (default var_str('unknown')) — The mass type that this mass distribution represents.
[nonAnalyticSolver] (default var_str('fallThrough')) — Selects how solutions are computed when no analytic solution is available. If set to “fallThrough” then the solution ignoring heating is used, while if set to “numerical” then numerical solvers are used to find solutions.
[componentType] (default var_str('unknown')) — The component type that this mass distribution represents.
[massType] (default var_str('unknown')) — The mass type that this mass distribution represents.
[tolerateVelocityMaximumFailure] (default .false.) — If true, tolerate failures to find the radius of the peak in the rotation curve.
[tolerateEnclosedMassIntegrationFailure] (default .false.) — If true, tolerate failures to find the mass enclosed as a function of radius.
[toleratePotentialIntegrationFailure] (default .false.) — If true, tolerate failures to compute the potential.
[fractionRadiusFinalSmall] (default 1.0d-3) — The initial radius is limited to be no smaller than this fraction of the final radius. This can help avoid problems in profiles that are extremely close to being disrupted.
[toleranceRelativePotential] (default 1.0d-3) — The maximum allowed relative tolerance to use in numerical solutions for the gravitational potential in dark-matter-only density profiles before aborting.
[lengthResolution] — The spatial resolution length scale (in Mpc) below which the underlying density profile is softened to a flat core, mimicking the finite force resolution of an N-body simulation.
[nonAnalyticSolver] (default var_str('fallThrough')) — Selects how solutions are computed when no analytic solution is available. If set to “fallThrough” then the solution ignoring heating is used, while if set to “numerical” then numerical solvers are used to find solutions.
[componentType] (default var_str('unknown')) — The component type that this mass distribution represents.
[massType] (default var_str('unknown')) — The mass type that this mass distribution represents.
[massMinimum] — The minimum halo mass (in \(\mathrm{M}_\odot\)) below which halos are excluded from the mass function histogram.
[massMaximum] — The maximum halo mass (in \(\mathrm{M}_\odot\)) above which halos are excluded from the mass function histogram.
[massCountPerDecade] — The number of logarithmic bins per decade of halo mass used when constructing the halo mass function.
[description] — A human-readable description of this mass function dataset, stored as metadata in the output file.
[simulationReference] — A bibliographic reference for the N-body simulation from which this mass function is derived, stored as metadata.
[simulationURL] — A URL pointing to the publicly accessible dataset or documentation for the N-body simulation, stored as metadata.
[bootstrapSampleCount] (default 30_c_size_t) — The number of bootstrap resamples of the particles that should be used.
[representativeMinimumCount] (default 10_c_size_t) — Minimum number of representative particles used to compute the center of a halo.
[tolerance] (default 1.0d-2) — The tolerance in the summed weight of bound particles which must be attained to declare convergence.
[bootstrapSampleRate] (default 1.0d0) — The sampling rate for particles.
[representativeFraction] (default 0.05d0) — Fraction of bound particles used to compute the center of a halo.
[analyzeAllParticles] (default .true.) — If true, all particles are assumed to be self-bound at the beginning of the analysis. Unbound particles at previous times are allowed to become bound in the current snapshot. If false and the self-bound information from the previous snapshot is available, only the particles that are self-bound at the previous snapshot are assumed to be bound at the beginning of the analysis.
[useVelocityMostBound] (default .false.) — If true, the velocity of the most bound particle in velocity space is used as the representative velocity of the satellite. If false, use the mass weighted mean velocity (center-of-mass velocity) of self-bound particles instead.
[orderRotation] (default var_str('none')) — The order in which evaluation of likelihoods should be rotated as a function of process number.
[logLikelihoodAccept] (default huge(0.0d0)) — The log-likelihood which should be “accepted”—once the log-likelihood reaches this value (or larger) no further updates to the chain will be made.
[report] (default .false.) — If true, report on the log-likelihood obtained.
[means] — The mean of the multivariate normal distribution.
[covariance] — The covariance matrix for the of the multivariate normal distribution.
[countForestsMaximum] (default -1_c_size_t) — If set to a positive number, this is the maximum number of forests that will be evolved.
[walltimeMaximum] (default -1_kind_int8) — If set to a positive number, this is the maximum wall time for which forest evolution is allowed to proceed before the task gives up.
[tolerateFailures] (default .false.) — If true then failures to evolve a forest are tolerated. The forest is evolved no further, but evolution of other forests continues.
[evolveForestsInParallel] (default .true.) — If true then each forest is evolved by a separate OpenMP thread. Otherwise, a single thread evolves all forests.
[suspendToRAM] (default .true.) — Specifies whether trees should be suspended to RAM (otherwise they are suspend to file).
[suspendPath] — The path to which tree suspension files will be stored.
[timeIntervalCheckpoint] (default -1_kind_int8) — If positive, gives the time in seconds between storing of checkpoint files. If zero or negative, no checkpointing is performed..
[fileNameCheckpoint] — The path to which checkpoint data will be stored.
[logM0] (default 10.0d0) — The parameter \(\log_{10} M_0\) (with \(M_0\) in units of \(\mathrm{M}_\odot\)) appearing in the star formation rate threshold expression for the star formation rate galactic filter class.
[logSFR0] (default 9.0d0) — The parameter \(\alpha_0\) appearing in the star formation rate threshold expression for the star formation rate galactic filter class.
[logSFR1] (default 0.0d0) — The parameter \(\alpha_1\) appearing in the star formation rate threshold expression for the star formation rate galactic filter class.
[cW] (default 3.78062835d0) — The parameter \(c_\mathrm{W}\) in the Bohr et al. (2021) power spectrum window function.
[beta] (default 3.4638743d0) — The parameter \(\beta\) in the Bohr et al. (2021) power spectrum window function.
[transferFunctionType] (default var_str('darkMatter')) — Specifies whether to use the darkMatter or total transfer function.
[fileName] — The name of the file from which to read a tabulated transfer function.
[redshift] (default 0.0d0) — The redshift of the transfer function to read.
[factorWavenumberSmoothExtrapolation] (default 0.0d0) — If positive, and extrapolation is used at high wavenumbers, the slope for extrapolation will be set by averaging over wavenumbers from \(k_\mathrm{max}/f\) to \(k_\mathrm{max}\), where \(f=\)[factorWavenumberSmoothExtrapolation] and \(k_\mathrm{max}\) is the highest wavenumber tabulated. This avoids spurious extrapolation for highly oscillatory transfer functions.
[acceptNegativeValues] (default .false.) — If true, negative values in the transfer function are allowed (and the absolute value is taken prior to interpolation). Otherwise, negative values result in an error.
[fractionalTimeStep] (default 0.01d0) — The fractional time step used when computing barrier crossing rates (i.e. the step used in finite difference calculations).
[fileName] (default var_str('none')) — The name of the file to/from which tabulations of barrier first crossing probabilities should be written/read. If set to “none” tables will not be stored.
[fractionalTimeStep] (default 0.01d0) — The fractional time step used when computing barrier crossing rates (i.e. the step used in finite difference calculations).
[varianceNumberPerUnitProbability] (default 1000) — The number of points to tabulate per unit variance for first crossing probabilities.
[varianceNumberPerUnit] (default 40) — The number of tabulation points per unit of \(\sigma^2\) used when building the rate look-up table for the Farahi excursion-set first-crossing distribution; higher values improve interpolation accuracy at the cost of memory and initialization time.
[varianceNumberPerDecade] (default 400) — The number of points to tabulate per decade of progenitor variance for first crossing rates.
[varianceNumberPerDecadeNonCrossing] (default 40) — The number of points to tabulate per decade of progenitor variance for non-crossing rates.
[timeNumberPerDecade] (default 10) — The number of tabulation points per decade of cosmic time used when building the first-crossing rate look-up table as a function of time; higher values improve temporal interpolation accuracy for rapidly evolving cosmologies.
[varianceIsUnlimited] (default .false.) — If true, the variance is assumed to have no upper limit (e.g. as in the case of CDM). This allows the tabulated solutions to be extended arbitrarily. Otherwise, tables are extended to encompass just the range of variance requested.
[linkingLength] (default 0.2d0) — The friends-of-friends linking length to use in computing virial density contrasts with the percolation analysis of More et al. (2011).

`nbodyOperatorMassFunctionCorrelation`¶

An N-body data operator which computes correlations between mass functions at different redshifts.

Methods

correlationLength — Return the current correlation length in the chains.
correlationLengthCompute — Compute correlation lengths in the chains.
postConvergenceCorrelationCount — Return the number of post-convergence correlation lengths that have accrued.

Parameters

[acceptedStateCount] (default 100) — The number of states to use in acceptance rate statistics.
[missingRootHalosAreFatal] (default .true.) — If true, if a hosted root halo is not found then a fatal error occurs. Otherwise, such cases are ignored and will not contribute to the halo mass function.
[massMinimum] — The minimum mass to consider counts.
[massMaximum] — The maximum mass to consider.
[massCountPerDecade] — The number of bins per decade of mass.
[countBootstraps] — The number of bootstrap samples to perform.
[description] — A description of this mass function.
[simulationReference] — A reference for the simulation.
[simulationURL] — A URL for the simulation.

`nbodyOperatorMassTotal`¶

An N-body data operator which computes the total mass of particles.

`nbodyOperatorMeanPosition`¶

An N-body data operator which computes the mass-weighted mean position and mean velocity of particles within an N-body halo. The selfBoundOnly parameter restricts the calculation to self-bound particles, and bootstrapSampleCount controls the number of bootstrap resamples for error estimation.

Parameters

[selfBoundParticlesOnly] — If true, the mean position and velocity are computed only for self-bound particles
[bootstrapSampleCount] (default 30_c_size_t) — The number of bootstrap resamples of the particles that should be used.

`nbodyOperatorMergerRates`¶

An N-body data operator which computes halo–halo merger rates from N-body simulation snapshots. Parameters control the snapshot index, primary and secondary halo mass ranges, whether to consider only always-isolated halos, and an optional suffix for labeling multiple merger rate outputs.

Parameters

[indexSnapshot] — The snapshot index of the descendant epoch at which halo mergers are identified; progenitors from the immediately preceding snapshot are matched to descendants at this epoch.
[massMinimum] — The minimum mass (of the secondary halo) for which to accumulate merging statistics.
[massMaximum] — The maximum mass (of the secondary halo) for which to accumulate merging statistics.
[massHostMinimum] — The minimum mass (of the primary halo) for which to accumulate merging statistics.
[massHostMaximum] — The maximum mass (of the primary halo) for which to accumulate merging statistics.
[missingHostIsFatal] (default .true.) — If true, missing host halos are cause for a fatal error. Otherwise they are ignored.
[alwaysIsolatedOnly] (default .true.) — If true, only mergers of halos which have been always isolated are considered. Otherwise, all halos are considered.
[suffix] (default var_str('')) — A suffix to append to the output merger rate attribute. Useful if you want to write output multiple merger rates.

`nbodyOperatorNull`¶

A null (no-op) N-body data operator that performs no transformation on the simulation data, used as a default or placeholder in operator pipelines.

(Default implementation)

Parameters

[dimensionless] (default .true.) — If true the null profile is considered to be dimensionless.

`nbodyOperatorPairCounts`¶

An N-body data operator which counts particle pairs in logarithmic bins of physical separation, enabling correlation function measurements. Parameters control the separation range and bin count, particle sampling rate, bootstrap resample count, and whether to compute cross-simulation or within-simulation pair counts.

Parameters

[crossCount] (default .false.) — If true, compute cross-simulation pair counts between the first and all simulations. Otherwise, compute pair counts within each simulation.
[bootstrapSampleCount] (default 30_c_size_t) — The number of bootstrap resamples of the particles that should be used.
[bootstrapSampleRate] (default 1.0d0) — The sampling rate for particles.
[separationMinimum] — The minimum pairwise separation (in the same units as particle positions) of pairs to include when counting and binning particle pairs.
[separationMaximum] — The maximum pairwise separation (in the same units as particle positions) of pairs to include when counting and binning particle pairs.
[separationCount] — The number of bins in separation for pair counts.
[includeUnbootstrapped] (default .true.) — If true, include results for the unbootstrapped (i.e. original) sample.

`nbodyOperatorPairwiseVelocityStatistics`¶

An N-body data operator which computes statistics of pairwise relative velocities between particles in bins of physical separation, useful for constraining galaxy peculiar velocity fields. Parameters control the separation range, redshift, Hubble flow correction, sampling rate, and bootstrap resample count.

Parameters

[crossCount] (default .false.) — If true, compute cross-simulation pairwise velocity statistics between the first and all simulations. Otherwise, compute pairwise velocity statistics within each simulation.
[addHubbleFlow] (default .false.) — If true, add the Hubble flow contribution \(H(z)\,\mathbf{r}\) to each particle velocity before computing pairwise statistics, converting peculiar velocities to total line-of-sight velocities.
[redshift] (default 0.0d0) — The redshift at which to evaluate the Hubble flow when adding the Hubble flow correction to particle velocities.
[bootstrapSampleCount] (default 30_c_size_t) — The number of bootstrap resamples of the particles that should be used.
[bootstrapSampleRate] (default 1.0d0) — The fraction of particles to sample randomly at each bootstrap iteration, between 0 and 1; values less than 1 reduce computational cost at the expense of statistical precision.
[separationMinimum] — The minimum separation to consider for pairwise velocity statistics.
[separationMaximum] — The maximum separation to consider for pairwise velocity statistics.
[separationCount] — The number of bins in separation for pairwise velocity statistics.
[includeUnbootstrapped] (default .true.) — If true, include results for the unbootstrapped (i.e. original) sample.

`nbodyOperatorPhysicalToComoving`¶

An N-body data operator which converts physical to comoving coordinates.

`nbodyOperatorPotentialEnergy`¶

An N-body data operator which computes the gravitational potential energy of each particle in an N-body halo using a tree-based algorithm. Parameters control whether only self-bound particles contribute, the tree opening-angle criterion, particle sampling rate, and bootstrap resample count.

Parameters

[selfBoundParticlesOnly] — If true, the gravitational potential is computed only from self-bound particles.
[bootstrapSampleCount] (default 30_c_size_t) — The number of bootstrap resamples of the particles that should be used.
[bootstrapSampleRate] (default 1.0d0) — The sampling rate for particles.
[thetaTolerance] (default 0.5d0) — The Barnes-Hut opening-angle tolerance \(\theta\); cells subtending an angle smaller than \(\theta\) at a particle are treated as point masses, trading accuracy for speed.

`nbodyOperatorProgenitorMassFunction`¶

An N-body data operator which computes the progenitor mass function by tabulating mass ratios of progenitor halos relative to their parent halos in logarithmic bins. Parameters specify the parent and progenitor mass ranges, bin counts per decade, the snapshot indices for parent and progenitor selection, and simulation metadata.

Parameters

[massParentMinimum] — The minimum halo mass (in \(\mathrm{M}_\odot\)) of parent halos to include in the progenitor mass function calculation.
[massParentMaximum] — The maximum halo mass (in \(\mathrm{M}_\odot\)) of parent halos to include in the progenitor mass function calculation.
[massParentCountPerDecade] — The number of logarithmic bins per decade of parent halo mass used when tabulating the progenitor mass function.
[massRatioProgenitorMinimum] — The minimum progenitor-to-parent mass ratio to include when binning the progenitor mass function.
[massRatioProgenitorMaximum] — The maximum progenitor-to-parent mass ratio to include when binning the progenitor mass function.
[massRatioProgenitorCountPerDecade] — The number of logarithmic bins per decade of progenitor-to-parent mass ratio used when tabulating the progenitor mass function.
[snapshotParents] — The snapshot index identifying the epoch at which parent halos are selected for the progenitor mass function.
[snapshotsProgenitors] — An array of snapshot indices identifying the epochs at which progenitor halos are selected and matched to parent halos in the progenitor mass function.
[description] — A description of this mass function.
[simulationReference] — A reference for the simulation.
[simulationURL] — A URL for the simulation.

`nbodyOperatorProgenitorOrderStatistics`¶

An N-body data operator which computes progenitor halo order statistics.

Parameters

[massParentMinimum] — The minimum parent mass to consider.
[massParentMaximum] — The maximum parent mass to consider.
[massParentCountPerDecade] — The number of bins per decade of parent mass.
[massRatioProgenitorMinimum] — The minimum mass ratio to consider.
[massRatioProgenitorMaximum] — The maximum mass ratio to consider.
[massRatioProgenitorCountPerDecade] — The number of bins per decade of mass ratio.
[snapshotParents] — The snapshot at which to select parent halos.
[orderMaximum] — The maximum order for which to compute the distribution.
[snapshotsProgenitors] — The snapshots at which to select progenitor halos.
[description] — A description of this mass function.
[simulationReference] — A reference for the simulation.
[simulationURL] — A URL for the simulation.

`nbodyOperatorRotationCurve`¶

An N-body data operator which computes the circular velocity rotation curve of an N-body halo by summing enclosed particle mass at a set of user-specified radii. Parameters specify the radii at which the rotation curve is evaluated, whether to restrict to self-bound particles, and the bootstrap resample count.

Parameters

[radiusSpecifiers] — A list of radius specifiers at which to output the rotation curve.
[includeRadii] (default .false.) — Specifies whether or not the radii at which rotation curve data are output should also be included in the output file.
[selfBoundParticlesOnly] — If true, the mean position and velocity are computed only for self-bound particles.
[radius] — An array of galactocentric radii (in the same units as particle positions) at which the circular velocity rotation curve will be evaluated.
[bootstrapSampleCount] (default 30_c_size_t) — The number of bootstrap resamples of the particles that should be used.

`nbodyOperatorSelectProperties`¶

An N-body data operator which filters N-body particles by retaining only those whose named integer property matches one of a user-supplied list of allowed values, enabling selection of subsets such as particles belonging to specific halos or simulation runs.

Parameters

[propertyName] — The name of the integer particle property whose values will be compared against the allowed list to determine which particles to retain.
[selectedValues] — A list of integer values of the named property that a particle must match to be retained; particles whose property value is not in this list are removed.

`nbodyOperatorSelfBound`¶

An N-body data operator which iteratively determines the gravitationally self-bound subset of particles within an N-body halo. Parameters control convergence tolerance, the minimum number of representative particles, the sampling rate, the bound fraction used to compute the halo center, and how prior-snapshot binding information is used.

Methods

tabulate — Tabulate the virial density contrast as a function of mass and time.
restoreTable — Restore a tabulated solution from file.
storeTable — Store a tabulated solution to file.

Parameters

[velocityCharacteristic] (default 250.0d0) — The velocity scale at which the SNe-driven outflow rate equals the star formation rate in disks.
[exponent] (default 3.5d0) — The velocity scaling of the SNe-driven outflow rate in disks.
[fraction] (default 0.01d0) — The normalization \(f\) of the outflow rate relative to the star formation rate at a reference halo velocity of 200 km/s and expansion factor of 1, setting the overall mass-loading amplitude of the halo-scaling feedback model.
[exponentVelocity] (default -2.0d0) — The exponent of virial velocity in the outflow rate in disks.
[exponentRedshift] (default 0.0d0) — The power-law exponent of the cosmological expansion factor \((1+z)\) in the halo-scaling outflow rate, allowing the mass-loading factor to evolve with redshift; a value of zero gives no redshift evolution.
[toleranceRelativeVelocityDispersion] (default 1.0d-6) — The relative tolerance to use in numerical solutions for the velocity dispersion in dark-matter-only density profiles.
[toleranceRelativeVelocityDispersionMaximum] (default 1.0d-3) — The maximum relative tolerance to use in numerical solutions for the velocity dispersion in dark-matter-only density profiles.
[radiusNormalization] (default 3.3d-6) — The initial value appearing in the radius-mass relation
[toleranceAbsoluteMass] (default 1.0d-6) — The mass tolerance used to judge whether the nuclear star cluster is physically plausible.
[toleranceRelativeMetallicity] (default 1.0d-4) — The metallicity tolerance for ODE solution.
[inactiveLuminositiesStellar] (default .false.) — Specifies whether or not nuclear star cluster stellar luminosities are inactive properties (i.e. do not appear in any ODE being solved).
[scaleRelativeMass] (default 1.0d-2) — The mass scale, relative to the total mass of the node, below which calculations in the delayed very simple hot halo component are allowed to become inaccurate.
[starveSatellites] (default .false.) — Specifies whether or not the hot halo should be removed (“starved”) when a node becomes a satellite.
[starveSatellitesOutflowed] (default .false.) — Specifies whether or not the outflowed hot halo should be removed (“starved”) when a node becomes a satellite.
[outflowReturnOnFormation] (default .false.) — Specifies whether or not outflowed gas should be returned to the hot reservoir on halo formation events.
[angularMomentumAlwaysGrows] (default .false.) — Specifies whether or not negative rates of accretion of angular momentum into the hot halo will be treated as positive for the purposes of computing the hot halo angular momentum.
[fractionBaryonLimitInNodeMerger] (default .false.) — Controls whether the hot gas content of nodes should be limited to not exceed the universal baryon fraction at node merger events. If set to true, hot gas (and angular momentum, abundances, and chemicals proportionally) will be removed from the merged halo to the unaccreted gas reservoir to limit the baryonic mass to the universal baryon fraction where possible.
[scaleAbsoluteMass] (default 100.0d0) — The absolute mass scale below which calculations in the very simple disk component are allowed to become inaccurate.
[toleranceAbsoluteMass] (default 1.0d-6) — The mass tolerance used to judge whether the disk is physically plausible.
[toleranceAbsoluteMass] (default 1.0d-6) — The mass tolerance used to judge whether the disk is physically plausible.
[toleranceRelativeMetallicity] (default 1.0d-4) — The metallicity tolerance for ODE solution.
[radiusStructureSolver] (default 1.0d0) — The radius (in units of the standard scale length) to use in solving for the size of the disk.
[structureSolverUseCole2000Method] (default .false.) — If true, use the method described in Cole et al. (2000) to correct for difference between thin disk and spherical mass distributions when solving for disk radii.
[diskNegativeAngularMomentumAllowed] (default .true.) — Specifies whether or not negative angular momentum is allowed for the disk.
[inactiveLuminositiesStellar] (default .false.) — Specifies whether or not disk stellar luminosities are inactive properties (i.e. do not appear in any ODE being solved).
[postStepZeroNegativeMasses] (default .true.) — If true, negative masses will be zeroed after each ODE step. Note that this can lead to non-conservation of mass.
[ratioAngularMomentumSolverRadius] (default ratioAngularMomentumSolverRadiusDefault) — The assumed ratio of the specific angular momentum at the structure solver radius to the mean specific angular momentum of the standard disk component.
[scaleAbsoluteMass] (default 100.0d0) — The absolute mass scale below which calculations in the very simple spheroid component are allowed to become inaccurate.
[toleranceAbsoluteMass] (default 1.0d-6) — The mass tolerance used to judge whether the spheroid is physically plausible.
[efficiencyEnergeticOutflow] (default 1.0d-2) — The proportionality factor relating mass outflow rate from the spheroid to the energy input rate divided by \(V_\mathrm{spheroid}^2\).
[toleranceRelativeMetallicity] (default 1.0d-4) — The metallicity tolerance for ODE solution.
[toleranceAbsoluteMass] (default 1.0d-6) — The mass tolerance used to judge whether the spheroid is physically plausible.
[inactiveLuminositiesStellar] (default .false.) — Specifies whether or not spheroid stellar luminosities are inactive properties (i.e. do not appear in any ODE being solved).
[postStepZeroNegativeMasses] (default .true.) — If true, negative masses will be zeroed after each ODE step. Note that this can lead to non-conservation of mass.
[ratioAngularMomentumScaleRadius] (default ratioAngularMomentumScaleRadiusDefault) — The assumed ratio of the specific angular momentum at the scale radius to the mean specific angular momentum of the standard spheroid component.
[outputMergers] (default .false.) — Determines whether or not properties of black hole mergers will be output.
[fileNames] — The name of the file(s) from which merger tree data should be read when using the [mergerTreeConstruct]\(=\)read tree construction method.
[forestSizeMaximum] (default 0_c_size_t) — The maximum number of nodes allowed in a forest before it will be broken up into trees and processed individually. A value of 0 implies that forests should never be split.
[presetMergerTimes] (default .true.) — Specifies whether merging times for subhalos should be preset when reading merger trees from a file.
[presetMergerNodes] (default .true.) — Specifies whether the target nodes for mergers should be preset (i.e. determined from descendant nodes). If they are not, merging will be with each satellite’s host node.
[presetSubhaloMasses] (default .true.) — Specifies whether subhalo mass should be preset when reading merger trees from a file.
[subhaloAngularMomentaMethod] (default var_str('summation')) — Specifies how to account for subhalo angular momentum when adding subhalo mass to host halo mass.
[presetSubhaloIndices] (default .true.) — Specifies whether subhalo indices should be preset when reading merger trees from a file.
[presetPositions] (default .true.) — Specifies whether node positions should be preset when reading merger trees from a file.
[presetScaleRadii] (default .true.) — Specifies whether node scale radii should be preset when reading merger trees from a file.
[scaleRadiiFailureIsFatal] (default .true.) — Specifies whether failure to set a node scale radii should be regarded as a fatal error. (If not, a fallback method to set scale radius is used in such cases.)
[presetScaleRadiiConcentrationMinimum] (default 3.0d0) — The lowest concentration (\(c=r_\mathrm{vir}/r_\mathrm{s}\)) allowed when setting scale radii, \(r_\mathrm{s}\).
[presetScaleRadiiConcentrationMaximum] (default 60.0d0) — The largest concentration (\(c=r_\mathrm{vir}/r_\mathrm{s}\)) allowed when setting scale radii, \(r_\mathrm{s}\).
[presetScaleRadiiMinimumMass] (default 0.0d0) — The minimum halo mass for which scale radii should be preset (if [presetScaleRadii]\(=\)true).
[presetUnphysicalAngularMomenta] (default .false.) — When reading merger trees from file and presetting halo angular momenta, detect unphysical (<=0) angular momenta and preset them using the selected halo spin method.
[presetAngularMomenta] (default .true.) — Specifies whether node angular momenta should be preset when reading merger trees from a file.
[presetAngularMomenta3D] (default .false.) — Specifies whether node 3-D angular momenta vectors should be preset when reading merger trees from a file.
[presetOrbits] (default .true.) — Specifies whether node orbits should be preset when reading merger trees from a file.
[presetOrbitsSetAll] (default .true.) — Forces all orbits to be set. If the computed orbit does not cross the virial radius, then select one at random instead.
[presetOrbitsAssertAllSet] (default .true.) — Asserts that all virial orbits must be preset. If any can not be set, Galacticus will stop.
[presetOrbitsBoundOnly] (default .true.) — Specifies whether only bound node orbits should be set.
[beginAt] (default -1_kind_int8) — Specifies the index of the tree to begin at. (Use -1 to always begin with the first tree.)
[outputTimeSnapTolerance] (default 0.0d0) — The relative tolerance required to “snap” a node time to the closest output time.
[missingHostsAreFatal] (default .true.) — Specifies whether nodes with missing host nodes should be considered to be fatal—see the discussion of missing host nodes in the class description above.
[treeIndexToRootNodeIndex] (default .false.) — Specifies whether tree indices should always be set to the index of their root node.
[allowBranchJumps] (default .true.) — Specifies whether nodes are allowed to jump between branches.
[allowSubhaloPromotions] (default .true.) — Specifies whether subhalos are permitted to be promoted to being isolated halos.
[alwaysPromoteMostMassive] (default .false.) — If true, the most massive progenitor is always promoted to be the primary progenitor even if it is a subhalo. Otherwise, isolated progenitors are given priority over subhalo progenitors, even if they are less massive.
[presetNamedReals] — Names of real datasets to be additionally read and stored in the nodes of the merger tree when using the [mergerTreeConstruct]\(=\)read tree construction method.
[presetNamedIntegers] — Names of integer datasets to be additionally read and stored in the nodes of the merger tree when using the [mergerTreeConstruct]\(=\)read tree construction method.
[fatalMismatches] (default .true.) — Specifies whether mismatches in cosmological parameter values between Galacticus and “Sussing Merger Trees” format (Srisawat et al., 2013) merger tree files should be considered fatal.
[fatalNonTreeNode] (default .true.) — Specifies whether nodes in snapshot files but not in the merger tree file should be considered fatal when importing from the “Sussing Merger Trees” format (Srisawat et al., 2013).
[subvolumeCount] (default 1) — Specifies the number of subvolumes along each axis into which a “Sussing Merger Trees” format (Srisawat et al., 2013) merger tree files should be split for processing through Galacticus.
[subvolumeBuffer] (default 0.0d0) — Specifies the buffer region (in units of Mpc\(/h\) to follow the format convention) around subvolumes of a “Sussing Merger Trees” format (Srisawat et al., 2013) merger tree file which should be read in to ensure that no halos are missed from trees.
[subvolumeIndex] (default [0,0,0]) — Specifies the index (in each dimension) of the subvolume of a “Sussing Merger Trees” format (Srisawat et al., 2013) merger tree file to process. Indices range from 0 to [subvolumeCount]\(-1\).
[badValue] (default -0.5d0) — Use for bad value detection in “Sussing” merger trees. Values for scale radius and halo spin which exceed this threshold are assumed to be bad.
[badValueTest] (default var_str('lessThan')) — Use for bad value detection in “Sussing” merger trees. Values which exceed the threshold in ths specified direction are assumed to be bad.
[treeSampleRate] (default 1.0d0) — Specify the probability that any given tree should processed (to permit subsampling).
[massOptions] (default var_str('default')) — Mass option for Sussing merger trees.
[mergeProbability] (default 0.1d0) — The largest probability of branching allowed in a timestep in merger trees built by the Cole et al. (2000) method.
[accretionLimit] (default 0.1d0) — The largest fractional mass change due to subresolution accretion allowed in a timestep in merger trees built by the Cole et al. (2000) method.
[redshiftMaximum] (default 1.0d5) — The highest redshift to which merger trees will be built in the Cole et al. (2000) method.
[toleranceTimeEarliest] (default 2.0d-6) — The fractional tolerance used to judge if a branch is at the earliest allowed time in the tree.
[branchIntervalStep] (default .true.) — If false use the original Cole et al. (2000) method to determine whether branching occurs in a timestep. If true draw branching intervals from a negative exponential distribution.
[toleranceResolutionSelf] (default 1.0d-6) — The fractional tolerance in node mass at the resolution limit below which branch mis-orderings will be ignored.
[toleranceResolutionParent] (default 1.0d-3) — The fractional tolerance in parent node mass at the resolution limit below which branch mis-orderings will be ignored.
[ignoreNoProgress] (default .false.) — If true, failure to make progress on a branch will be ignored (and the branch terminated).
[ignoreWellOrdering] (default .false.) — If true, non-well-ordered tree branches are pruned away instead of causing errors..
[redshiftBase] (default 0.0d0) — The redshift at which to plant the base node when building merger trees.
[timeSnapTolerance] (default 1.0d-6) — The fractional tolerance within which the tree base time will be snapped to a nearby output time.
[treeBeginAt] (default 0) — The index (in order of increasing base halo mass) of the tree at which to begin when building merger trees. A value of “0” means to begin with tree number 1 (if processing trees in ascending order), or equal to the number of trees (otherwise).
[processDescending] (default .true.) — If true, causes merger trees to be processed in order of decreasing mass.
[splitTrees] (default .false.) — If true, prune away any nodes of the tree that are not needed to determine evolution up to the latest time at which a node is present inside the lightcone. This typically leads to a tree splitting into a forest of trees.
[label] — A label for the mass function.
[comment] — A descriptive comment for the mass function.
[starFormationRates] — The star formation rates corresponding to bin centers.
[covarianceBinomialBinsPerDecade] (default 10) — The number of bins per decade of star formation rate to use when constructing star formation rate function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMinimum] (default 1.0d10) — The star formation rate to consider when constructing star formation rate function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMaximum] (default 1.0d12) — The maximum star formation rate to consider when constructing star formation rate function covariance matrices for main branch galaxies.
[targetLabel] — Label for the target dataset.
[functionValueTarget] — The target function for likelihood calculations.
[functionCovarianceTarget] — The target function covariance for likelihood calculations.
[likelihoodBins] — Controls which bins in the stellar mass–halo mass relation will be used in computing the likelihood: * not present: all bins are included in the likelihood calculation; * list of integers: use only the mass bin(s) given in this list in the likelihood calculation; * auto: use only bins which have a non-zero number of halos contributing to them in the likelihood calculation.
[fileNameTarget] — The name of the file containing the target data.
[redshiftInterval] (default 1) — The redshift interval to use.
[likelihoodNormalize] (default .false.) — If true, then normalize the likelihood to make it a probability density.
[computeScatter] (default .false.) — If true, the scatter in log10(stellar mass) is computed. Otherwise, the mean is computed.
[systematicErrorPolynomialCoefficient] (default [0.0d0]) — The coefficients of the systematic error polynomial for stellar mass in the stellar vs halo mass relation.
[systematicErrorMassHaloPolynomialCoefficient] (default [0.0d0]) — The coefficients of the systematic error polynomial for halo mass in the stellar vs halo mass relation.
[errorTolerant] (default .false.) — Error tolerance for the N-body spin distribution operator.
[logNormalRange] (default 100.0d0) — The multiplicative range of the log-normal distribution used to model the distribution of the mass and energy terms in the spin parameter. Specifically, the lognormal distribution is truncated outside the range \((\lambda_\mathrm{m}/R,\lambda_\mathrm{m} R\), where \(\lambda_\mathrm{m}\) is the measured spin, and \(R=\)[logNormalRange]
[fileName] — The name of the file from which to read spin distribution function parameters.
[comment] — A comment describing this analysis.
[label] — A label for this analysis.
[label] — A label for the spin distribution function.
[comment] — A descriptive comment for the spin distribution function.
[redshift] — The redshift at which to compute the spin distribution function.
[massMinimum] — Minimum halo mass for the spin distribution function.
[massMaximum] — Maximum halo mass for the spin distribution function.
[spinMinimum] — Minimum spin for the spin distribution function.
[spinMaximum] — Maximum spin for the spin distribution function.
[countSpinsPerDecade] — Number of spins per decade at which to compute the spin distribution function.
[timeRecent] — Halos which experienced a major node merger within a time \(\Delta t=\)[timeRecent] of the analysis time will be excluded from the analysis.
[particleCountMinimum] — The minimum particle count to assume when computing N-body errors on spins.
[massParticle] — The mass of the particle used in the N-body simulation from which spins were measured.
[energyEstimateParticleCountMaximum] — The maximum number of particles used in estimating halo energies when measuring spins from the N-body simulation.
[targetLabel] — Label for the target dataset.
[functionValueTarget] — The target function for likelihood calculations.
[functionCovarianceTarget] — The target function covariance for likelihood calculations.
[label] — A label for the mass function.
[comment] — A descriptive comment for the mass function.
[masses] — The masses corresponding to bin centers.
[covarianceBinomialBinsPerDecade] (default 10) — The number of bins per decade of halo mass to use when constructing HI mass function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMinimum] (default 1.0d8) — The minimum halo mass to consider when constructing HI mass function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMaximum] (default 1.0d16) — The maximum halo mass to consider when constructing HI mass function covariance matrices for main branch galaxies.
[targetLabel] — Label for the target dataset.
[functionValueTarget] — The target function for likelihood calculations.
[functionCovarianceTarget] — The target function covariance for likelihood calculations.
[label] — A label for the luminosity function.
[comment] — A descriptive comment for the luminosity function.
[magnitudesAbsolute] — The absolute magnitudes corresponding to bin centers.
[covarianceBinomialBinsPerDecade] (default 10) — The number of bins per decade of halo mass to use when constructing luminosity function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMinimum] (default 1.0d8) — The minimum halo mass to consider when constructing luminosity function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMaximum] (default 1.0d16) — The maximum halo mass to consider when constructing luminosity function covariance matrices for main branch galaxies.
[targetLabel] — Label for the target dataset.
[functionValueTarget] — The target function for likelihood calculations.
[functionCovarianceTarget] — The target function covariance for likelihood calculations.
[label] — A label for the luminosity function.
[comment] — A descriptive comment for the luminosity function.
[luminosities] — The luminosities corresponding to bin centers.
[covarianceBinomialBinsPerDecade] (default 10) — The number of bins per decade of halo mass to use when constructing luminosity function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMinimum] (default 1.0d8) — The minimum halo mass to consider when constructing luminosity function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMaximum] (default 1.0d16) — The maximum halo mass to consider when constructing luminosity function covariance matrices for main branch galaxies.
[includeNitrogenII] (default .false.) — If true, include contamination by the [NII] (6548AA \(+\) 6584AA) doublet.
[depthOpticalISMCoefficient] (default 1.0d0) — Multiplicative coefficient for optical depth in the ISM.
[targetLabel] — Label for the target dataset.
[functionValueTarget] — The target function for likelihood calculations.
[functionCovarianceTarget] — The target function covariance for likelihood calculations.
[label] — A label for the mass function.
[comment] — A descriptive comment for the mass function.
[masses] — The masses corresponding to bin centers.
[covarianceBinomialBinsPerDecade] (default 10) — The number of bins per decade of halo mass to use when constructing stellar mass function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMinimum] (default 1.0d8) — The minimum halo mass to consider when constructing stellar mass function covariance matrices for main branch galaxies.
[covarianceBinomialMassHaloMaximum] (default 1.0d16) — The maximum halo mass to consider when constructing stellar mass function covariance matrices for main branch galaxies.
[targetLabel] — Label for the target dataset.
[functionValueTarget] — The target function for likelihood calculations.
[functionCovarianceTarget] — The target function covariance for likelihood calculations.
[rootVarianceFractionalMinimum] (default 0.0d0) — The minimum fractional root variance (relative to the target dataset).
[fileName] — The name of the file from which to read concentration distribution function parameters.
[comment] — A comment describing this analysis.
[label] — A label for this analysis.
[label] — A label for the concentration distribution function.
[comment] — A descriptive comment for the concentration distribution function.
[redshift] — The redshift at which to compute the concentration distribution function.
[massMinimum] — Minimum halo mass for the concentration distribution function.
[massMaximum] — Maximum halo mass for the concentration distribution function.
[concentrationMinimum] — Minimum concentration for the concentration distribution function.
[concentrationMaximum] — Maximum concentration for the concentration distribution function.
[countConcentrationsPerDecade] — Number of concentrations per decade at which to compute the concentration distribution function.
[timeRecent] — Halos which experienced a major node merger within a time \(\Delta t=\)[timeRecent] of the analysis time will be excluded from the analysis.
[massParticle] — The particle mass in the source N-body simulation.
[targetLabel] — Label for the target dataset.
[functionValueTarget] — The target function for likelihood calculations.
[functionCovarianceTarget] — The target function covariance for likelihood calculations.
[fileName] — The name of the file from which to read star forming main sequence function parameters.
[comment] — A comment describing this analysis.
[label] — A label for this analysis.
[label] — A label for the star forming main sequence function.
[comment] — A descriptive comment for the star forming main sequence function.
[massMinimum] — Minimum stellar mass for the star forming main sequence function.
[massMaximum] — Maximum stellar mass for the star forming main sequence function.
[countMassesPerDecade] — Number of masses per decade at which to compute the star forming main sequence function.
[targetLabel] — Label for the target dataset.
[meanValueTarget] — The target function for likelihood calculations.
[meanCovarianceTarget] — The target function covariance for likelihood calculations.
[label] — A label for the mass function.
[comment] — A descriptive comment for the mass function.
[separations] — The separations corresponding to bin centers.
[massMinima] — The minimum mass of each mass sample.
[massMaxima] — The maximum mass of each mass sample.
[massHaloBinsPerDecade] (default 10) — The number of bins per decade of halo mass to use when constructing the mass function covariance matrix for main branch galaxies.
[massHaloMinimum] (default 1.0d8) — The minimum halo mass to consider when constructing the mass function covariance matrix for main branch galaxies.
[massHaloMaximum] (default 1.0d16) — The maximum halo mass to consider when constructing the mass function covariance matrix for main branch galaxies.
[wavenumberCount] (default 60_c_size_t) — The number of bins in wavenumber to use in computing the correlation function.
[wavenumberMinimum] (default 1.0d-3) — The minimum wavenumber to use when computing the correlation function.
[wavenumberMaximum] (default 1.0d4) — The maximum wavenumber to use when computing the correlation function.
[integralConstraint] — The integral constraint for these correlation functions.
[depthLineOfSight] — The line-of-sight depth over which the correlation function was projected.
[halfIntegral] — Set to true if the projection integrand should be over line-of-sight depths greater than zero.
[binnedProjectedCorrelationTarget] — The target function for likelihood calculations.
[binnedProjectedCorrelationCovarianceTarget] — The target function covariance for likelihood calculations.
[targetLabel] (default var_str('')) — A label for the target dataset in a plot of this analysis.
[starFormationRateSpecificQuiescentLogarithmic] — The base-10 logarithm specific star formation rate (in units of Gyr\(^{-1}\)) separating quiescent and star-forming galaxies.
[starFormationRateSpecificLogarithmicError] — The observational fractional error in specific star formation rate (in units of dex) of galaxies.
[fileName] — The name of the file from which to read quiescent fraction function parameters.
[comment] — A comment describing this analysis.
[label] — A label for this analysis.
[label] — A label for the star forming main sequence function.
[comment] — A descriptive comment for the star forming main sequence function.
[massMinimum] — Minimum stellar mass for the star forming main sequence function.
[massMaximum] — Maximum stellar mass for the star forming main sequence function.
[countMassesPerDecade] — Number of masses per decade at which to compute the star forming main sequence function.
[targetLabel] — Label for the target dataset.
[meanValueTarget] — The target function for likelihood calculations.
[meanCovarianceTarget] — The target function covariance for likelihood calculations.
[nonAnalyticSolver] (default var_str('fallThrough')) — Selects how solutions are computed when no analytic solution is available. If set to “fallThrough” then the solution ignoring heating is used, while if set to “numerical” then numerical solvers are used to find solutions.
[radiusFractionalTruncateMinimum] (default 2.0d0) — The minimum radius (in units of the virial radius) to begin truncating the density profile.
[radiusFractionalTruncateMaximum] (default 4.0d0) — The maximum radius (in units of the virial radius) to finish truncating the density profile.
[nonAnalyticSolver] (default var_str('fallThrough')) — Selects how solutions are computed when no analytic solution is available. If set to “fallThrough” then the solution ignoring heating is used, while if set to “numerical” then numerical solvers are used to find solutions.
[velocityDispersionApproximate] (default .true.) — If true, radial velocity dispersion is computed using an approximate method in which we assume that \(\sigma_\mathrm{r}^2(r) \rightarrow \sigma_\mathrm{r}^2(r) - (2/3) \epsilon(r)\), where \(\epsilon(r)\) is the specific heating energy. If false then radial velocity dispersion is computed by numerically solving the Jeans equation.
[tolerateEnclosedMassIntegrationFailure] (default .false.) — If true, tolerate failures to find the mass enclosed as a function of radius.
[tolerateVelocityDispersionFailure] (default .false.) — If true, tolerate failures to compute the velocity dispersion.
[tolerateVelocityMaximumFailure] (default .false.) — If true, tolerate failures to find the radius of the maximum circular velocity.
[toleratePotentialIntegrationFailure] (default .false.) — If true, tolerate numerical failures when computing the gravitational potential of a heated dark matter profile, allowing the calculation to continue with a fallback result rather than aborting.
[toleranceRelativeVelocityDispersion] (default 1.0d-6) — The relative tolerance to use in numerical solutions for the velocity dispersion in dark-matter-only density profiles.
[toleranceRelativeVelocityDispersionMaximum] (default 1.0d-3) — The maximum relative tolerance to use in numerical solutions for the velocity dispersion in dark-matter-only density profiles.
[fractionRadiusFinalSmall] (default 1.0d-3) — The initial radius is limited to be no smaller than this fraction of the final radius. This can help avoid problems in profiles that are extremely close to being disrupted.
[toleranceRelativePotential] (default 1.0d-3) — The maximum allowed relative tolerance to use in numerical solutions for the gravitational potential in dark-matter-only density profiles before aborting.
[tolerateVelocityMaximumFailure] (default .true.) — If true, tolerate failures to find the radius of the peak in the rotation curve.
[lengthResolution] — The gravitational softening length \(\Delta x\) (in Mpc) of the N-body simulation, which sets the minimum spatial scale below which the dark matter profile is smoothed to avoid artificial two-body effects.
[massResolution] — The mass resolution \(\Delta M\) (in \(\mathrm{M}_\odot\)) of the N-body simulation, representing the minimum halo mass that can be resolved; profiles of halos near this limit are softened to account for particle discreteness effects.
[resolutionIsComoving] — If true, the resolution length is assumed to be fixed in comoving coordinates, otherwise in physical coordinates.
[nonAnalyticSolver] (default var_str('fallThrough')) — Selects how solutions are computed when no analytic solution is available. If set to “fallThrough” then the solution ignoring heating is used, while if set to “numerical” then numerical solvers are used to find solutions.
[C] (default 400.0d0) — The parameter \(C\) appearing in the halo concentration algorithm of Ludlow et al. (2016).
[f] (default 0.02d0) — The parameter \(f\) appearing in the halo concentration algorithm of Ludlow et al. (2016).
[timeFormationSeekDelta] (default 0.0d0) — The parameter \(\Delta \log t\) by which the logarithm of the trial formation time is incremented when stepping through the formation history of a node to find the formation time. If set to zero (or a negative value) the cumulative mass histories of nodes are assumed to be monotonic functions of time, and the formation time is instead found by a root finding algorithm,
[massBoundIsInactive] (default .false.) — Specifies whether or not the bound mass of the satellite component is inactive (i.e. does not appear in any ODE being solved).
[useLastIsolatedTime] (default .false.) — If true, evaluate the halo virial radius using a the virial density definition at the last isolated time of the halo.
[filterName] — The filter to select.
[filterType] — The filter type (rest or observed) to select.
[redshiftBand] — The redshift of the band (if not the output redshift).
[postprocessChain] — The postprocessing chain to use.
[cloudyTableFileName] (default var_str('%DATASTATICPATH%/hiiRegions/emissionLineLuminosities_BC2003_highResolution_imfChabrier.hdf5')) — The file of emission line luminosities to use.
[lineNames] — The emission lines to extract.
[component] — The component from which to extract star formation rate.
[toleranceRelative] (default 1.0d-3) — The relative tolerance used in integration over stellar population spectra.
[component] — The component from which to extract star formation rate.
[radiusCore] — The soliton core radius (in Mpc) characterizing the size of the quantum pressure-supported central core of the fuzzy dark matter halo; the density profile flattens inside this scale.
[densitySolitonCentral] — The central density (in \(\mathrm{M}_\odot\)/Mpc\(^3\)) of the solitonic core at \(r=0\), which sets the overall normalization of the density profile \(\rho(r) = \rho_\mathrm{c} [1+(r/r_c)^2]^{-8}\).
[toleranceRelativePotential] (default 1.0d-3) — The relative tolerance used in numerical ODE solutions for the gravitational potential of the solitonic core profile.
[dimensionless] (default .true.) — If true the soliton profile is treated as dimensionless (scale-free), allowing its radial and density quantities to be specified in arbitrary units.
[componentType] (default var_str('unknown')) — The galactic structure component type (e.g.dark matter halo, disk, spheroid) represented by this mass distribution, used for component-specific queries.
[massType] (default var_str('unknown')) — The mass type (e.g.dark matter, baryonic, total) represented by this mass distribution, used for mass-type-specific queries.
[radiusTransition] — The transition radius (in Mpc) at which the density profile smoothly switches from the halo profile to the accretion flow, controlled by the fourth-order transition function \(f_\mathrm{trans}(r)\).
[nonAnalyticSolver] (default var_str('fallThrough')) — Selects how solutions are computed when no analytic solution is available. If set to “fallThrough” then the solution ignoring heating is used, while if set to “numerical” then numerical solvers are used to find solutions.
[componentType] (default var_str('unknown')) — The component type that this mass distribution represents.
[massType] (default var_str('unknown')) — The mass type that this mass distribution represents.
[timeAge] — The age of the halo (in Gyr) since its formation, determining the total time available for SIDM self-interactions to thermalize the inner halo and produce an isothermal core.
[velocityRelativeMean] — Mean relative velocity to calculate self interaction cross section.
[nonAnalyticSolver] (default var_str('fallThrough')) — Selects how solutions are computed when no analytic solution is available. If set to “fallThrough” then the solution ignoring heating is used, while if set to “numerical” then numerical solvers are used to find solutions.
[componentType] (default var_str('unknown')) — The component type that this mass distribution represents.
[massType] (default var_str('unknown')) — The mass type that this mass distribution represents.
[nonAnalyticSolver] (default var_str('fallThrough')) — Selects how solutions are computed when no analytic solution is available. If set to “fallThrough” then the solution ignoring heating is used, while if set to “numerical” then numerical solvers are used to find solutions.
[componentType] (default var_str('unknown')) — The component type that this mass distribution represents.
[massType] (default var_str('unknown')) — The mass type that this mass distribution represents.
[tolerateVelocityMaximumFailure] (default .false.) — If true, tolerate failures to find the radius of the peak in the rotation curve.
[tolerateEnclosedMassIntegrationFailure] (default .false.) — If true, tolerate failures to find the mass enclosed as a function of radius.
[toleratePotentialIntegrationFailure] (default .false.) — If true, tolerate failures to compute the potential.
[fractionRadiusFinalSmall] (default 1.0d-3) — The initial radius is limited to be no smaller than this fraction of the final radius. This can help avoid problems in profiles that are extremely close to being disrupted.
[toleranceRelativePotential] (default 1.0d-3) — The maximum allowed relative tolerance to use in numerical solutions for the gravitational potential in dark-matter-only density profiles before aborting.
[lengthResolution] — The spatial resolution length scale (in Mpc) below which the underlying density profile is softened to a flat core, mimicking the finite force resolution of an N-body simulation.
[nonAnalyticSolver] (default var_str('fallThrough')) — Selects how solutions are computed when no analytic solution is available. If set to “fallThrough” then the solution ignoring heating is used, while if set to “numerical” then numerical solvers are used to find solutions.
[componentType] (default var_str('unknown')) — The component type that this mass distribution represents.
[massType] (default var_str('unknown')) — The mass type that this mass distribution represents.
[massMinimum] — The minimum halo mass (in \(\mathrm{M}_\odot\)) below which halos are excluded from the mass function histogram.
[massMaximum] — The maximum halo mass (in \(\mathrm{M}_\odot\)) above which halos are excluded from the mass function histogram.
[massCountPerDecade] — The number of logarithmic bins per decade of halo mass used when constructing the halo mass function.
[description] — A human-readable description of this mass function dataset, stored as metadata in the output file.
[simulationReference] — A bibliographic reference for the N-body simulation from which this mass function is derived, stored as metadata.
[simulationURL] — A URL pointing to the publicly accessible dataset or documentation for the N-body simulation, stored as metadata.
[bootstrapSampleCount] (default 30_c_size_t) — The number of bootstrap resamples of the particles that should be used.
[representativeMinimumCount] (default 10_c_size_t) — Minimum number of representative particles used to compute the center of a halo.
[tolerance] (default 1.0d-2) — The tolerance in the summed weight of bound particles which must be attained to declare convergence.
[bootstrapSampleRate] (default 1.0d0) — The sampling rate for particles.
[representativeFraction] (default 0.05d0) — Fraction of bound particles used to compute the center of a halo.
[analyzeAllParticles] (default .true.) — If true, all particles are assumed to be self-bound at the beginning of the analysis. Unbound particles at previous times are allowed to become bound in the current snapshot. If false and the self-bound information from the previous snapshot is available, only the particles that are self-bound at the previous snapshot are assumed to be bound at the beginning of the analysis.
[useVelocityMostBound] (default .false.) — If true, the velocity of the most bound particle in velocity space is used as the representative velocity of the satellite. If false, use the mass weighted mean velocity (center-of-mass velocity) of self-bound particles instead.
[orderRotation] (default var_str('none')) — The order in which evaluation of likelihoods should be rotated as a function of process number.
[logLikelihoodAccept] (default huge(0.0d0)) — The log-likelihood which should be “accepted”—once the log-likelihood reaches this value (or larger) no further updates to the chain will be made.
[report] (default .false.) — If true, report on the log-likelihood obtained.
[means] — The mean of the multivariate normal distribution.
[covariance] — The covariance matrix for the of the multivariate normal distribution.
[countForestsMaximum] (default -1_c_size_t) — If set to a positive number, this is the maximum number of forests that will be evolved.
[walltimeMaximum] (default -1_kind_int8) — If set to a positive number, this is the maximum wall time for which forest evolution is allowed to proceed before the task gives up.
[tolerateFailures] (default .false.) — If true then failures to evolve a forest are tolerated. The forest is evolved no further, but evolution of other forests continues.
[evolveForestsInParallel] (default .true.) — If true then each forest is evolved by a separate OpenMP thread. Otherwise, a single thread evolves all forests.
[suspendToRAM] (default .true.) — Specifies whether trees should be suspended to RAM (otherwise they are suspend to file).
[suspendPath] — The path to which tree suspension files will be stored.
[timeIntervalCheckpoint] (default -1_kind_int8) — If positive, gives the time in seconds between storing of checkpoint files. If zero or negative, no checkpointing is performed..
[fileNameCheckpoint] — The path to which checkpoint data will be stored.
[logM0] (default 10.0d0) — The parameter \(\log_{10} M_0\) (with \(M_0\) in units of \(\mathrm{M}_\odot\)) appearing in the star formation rate threshold expression for the star formation rate galactic filter class.
[logSFR0] (default 9.0d0) — The parameter \(\alpha_0\) appearing in the star formation rate threshold expression for the star formation rate galactic filter class.
[logSFR1] (default 0.0d0) — The parameter \(\alpha_1\) appearing in the star formation rate threshold expression for the star formation rate galactic filter class.
[cW] (default 3.78062835d0) — The parameter \(c_\mathrm{W}\) in the Bohr et al. (2021) power spectrum window function.
[beta] (default 3.4638743d0) — The parameter \(\beta\) in the Bohr et al. (2021) power spectrum window function.
[transferFunctionType] (default var_str('darkMatter')) — Specifies whether to use the darkMatter or total transfer function.
[fileName] — The name of the file from which to read a tabulated transfer function.
[redshift] (default 0.0d0) — The redshift of the transfer function to read.
[factorWavenumberSmoothExtrapolation] (default 0.0d0) — If positive, and extrapolation is used at high wavenumbers, the slope for extrapolation will be set by averaging over wavenumbers from \(k_\mathrm{max}/f\) to \(k_\mathrm{max}\), where \(f=\)[factorWavenumberSmoothExtrapolation] and \(k_\mathrm{max}\) is the highest wavenumber tabulated. This avoids spurious extrapolation for highly oscillatory transfer functions.
[acceptNegativeValues] (default .false.) — If true, negative values in the transfer function are allowed (and the absolute value is taken prior to interpolation). Otherwise, negative values result in an error.
[fractionalTimeStep] (default 0.01d0) — The fractional time step used when computing barrier crossing rates (i.e. the step used in finite difference calculations).
[fileName] (default var_str('none')) — The name of the file to/from which tabulations of barrier first crossing probabilities should be written/read. If set to “none” tables will not be stored.
[fractionalTimeStep] (default 0.01d0) — The fractional time step used when computing barrier crossing rates (i.e. the step used in finite difference calculations).
[varianceNumberPerUnitProbability] (default 1000) — The number of points to tabulate per unit variance for first crossing probabilities.
[varianceNumberPerUnit] (default 40) — The number of tabulation points per unit of \(\sigma^2\) used when building the rate look-up table for the Farahi excursion-set first-crossing distribution; higher values improve interpolation accuracy at the cost of memory and initialization time.
[varianceNumberPerDecade] (default 400) — The number of points to tabulate per decade of progenitor variance for first crossing rates.
[varianceNumberPerDecadeNonCrossing] (default 40) — The number of points to tabulate per decade of progenitor variance for non-crossing rates.
[timeNumberPerDecade] (default 10) — The number of tabulation points per decade of cosmic time used when building the first-crossing rate look-up table as a function of time; higher values improve temporal interpolation accuracy for rapidly evolving cosmologies.
[varianceIsUnlimited] (default .false.) — If true, the variance is assumed to have no upper limit (e.g. as in the case of CDM). This allows the tabulated solutions to be extended arbitrarily. Otherwise, tables are extended to encompass just the range of variance requested.
[linkingLength] (default 0.2d0) — The friends-of-friends linking length to use in computing virial density contrasts with the percolation analysis of More et al. (2011).

`nbodyOperatorSelfBoundBarnesHut`¶

An N-body data operator which determines the subset of particles that are self-bound. The potential is computed using a tree method following Barnes and Hut (1986).

Parameters

[thetaTolerance] (default 0.5d0) — The Barnes-Hut opening-angle tolerance \(\theta\); cells subtending an angle smaller than \(\theta\) at a particle are treated as point masses, trading accuracy for speed.

`nbodyOperatorSelfFrictionAcceleration`¶

An N-body data operator which determines the acceleration of self-bound particles from unbound ones. The interaction between particles is computed using a tree method following Barnes and Hut (1986).

Parameters

[bootstrapSampleCount] (default 30_c_size_t) — The number of bootstrap resamples of the particles that should be used.
[thetaTolerance] (default 0.5d0) — The Barnes-Hut opening-angle tolerance \(\theta\); cells subtending an angle smaller than \(\theta\) at a particle are treated as point masses, trading accuracy for speed.

`nbodyOperatorSequence`¶

An N-body data operator which applies a sequence of other operators.

Methods

prepend — Prepend an operator to a sequence of weight operators.

`nbodyOperatorSetBoxSize`¶

An N-body data operator which assigns or overrides the periodic simulation box size for an N-body data set. This is useful when the box size metadata is absent or incorrect in the original data. The boxSize parameter specifies the side length to apply.

Parameters

[boxSize] — The side length of the periodic simulation box to assign to the dataset, in the same units as particle positions.

`nbodyOperatorShiftProperty`¶

An N-body data operator which shifts values of a property by an integer amount.

Parameters

[propertyName] — The name of the integer particle property whose values will be shifted by adding the constant shiftBy.
[shiftBy] — The integer constant added to every value of the named property, enabling renumbering or offset corrections for IDs and indices.

`nbodyOperatorSimulationSelector`¶

An N-body data operator which applies a delegate nbodyOperatorClass to a single element of the simulation array identified by [indexSimulation], leaving all other simulations unchanged.

Parameters

[indexSimulation] — The 1-based index into the simulation array selecting the single simulation dataset to which the delegate operator will be applied.

`nbodyOperatorSpinDistributionFunction`¶

An N-body data operator which computes the halo spin distribution function by binning halos as a function of dimensionless spin parameter within a specified mass and spin range. Mass limits and binning are set by [massMinimum], [massMaximum], and [massCountPerDecade], spin limits and binning by [spinMinimum], [spinMaximum], and [spinCountPerDecade].

Parameters

[massMinimum] — The minimum mass to consider.
[massMaximum] — The maximum mass to consider.
[massCountPerDecade] — The number of logarithmic bins per decade of mass used when constructing the spin distribution function.
[spinMinimum] — The minimum dimensionless spin parameter below which halos are excluded from the spin distribution function histogram.
[spinMaximum] — The maximum dimensionless spin parameter above which halos are excluded from the spin distribution function histogram.
[spinCountPerDecade] — The number of logarithmic bins per decade of spin parameter used when constructing the spin distribution function.
[description] — A human-readable description of this spin distribution function dataset, stored as metadata in the output file.
[simulationReference] — A bibliographic reference for the N-body simulation from which this spin distribution is derived, stored as output metadata.
[simulationURL] — A URL pointing to the publicly accessible dataset or documentation for the N-body simulation, stored as output metadata.

`nbodyOperatorSubhaloMassFunction`¶

An N-body data operator which computes subhalo mass functions.

Methods

finalizeAnalysis — Finalize analysis.

Parameters

[fileName] — The name of the file from which to read the target dataset.
[massRatioMinimum] (default 1.0d-4) — The minimum mass ratio to consider.
[massRatioMaximum] (default 1.0d0) — The maximum mass ratio to consider.
[countMassRatios] (default 10_c_size_t) — The number of bins in mass ratio to use.
[redshift] (default 0.0d0) — The redshift at which to compute the subhalo mass function.
[ignoreEmptyModelBins] (default .false.) — If true, model bins containing no subhalos are ignored. Otherwise, if the target data is non-zero, an impossible likelihood is returned.
[negativeBinomialScatterFractional] (default 0.18d0) — The fractional scatter (relative to the Poisson scatter) in the negative binomial distribution used in likelihood calculations.
[massHost] — The mass of the host halo.
[massRatioMinimum] — The minimum mass ratio to consider.
[massRatioMaximum] — The maximum mass ratio to consider.
[massCountPerDecade] — The number of bins per decade of mass ratio.
[description] — A description of this subhalo mass function.
[simulationReference] — A reference for the simulation.
[simulationURL] — A URL for the simulation.

`nbodyOperatorSubhaloRadiusFunction`¶

An N-body data operator which computes subhalo radial distribution functions.

Parameters

[massMinimum] — The minimum subhalo mass to include.
[radiusVirialHost] — The virial radius of the host halo.
[radiusRatioMinimum] — The minimum radius ratio to consider.
[radiusRatioMaximum] — The maximum radius ratio to consider.
[radiusCountPerDecade] — The number of bins per decade of radius ratio.
[description] — A description of this subhalo radial distribution function.
[simulationReference] — A reference for the simulation.
[simulationURL] — A URL for the simulation.

`nbodyOperatorSubhaloVelocityMaximumMeanFunction`¶

An N-body data operator which computes the mean subhalo \(V_\mathrm{max}\) as a function of subhalo mass.

Parameters

[massMinimum] — The minimum mass to consider.
[massMaximum] — The maximum mass to consider.
[massCountPerDecade] — The number of bins per decade of mass.
[description] — A description of this subhalo mass function.
[simulationReference] — A reference for the simulation.
[simulationURL] — A URL for the simulation.

`nbodyOperatorSubsample`¶

An N-body data operator which filters out particles based on a property range.

Parameters

[massThreshold] — The mass threshold, \(M_0\), below which subsampling is applied.
[fractionMassThreshold] — The fractional (relative to the tree mass) mass threshold, \(f_0\), below which subsampling is applied.
[subsamplingRateAtThreshold] — The subsampling rate at the mass threshold, \(P_0\).
[exponent] — The exponent, \(\alpha\), of mass in the subsampling probability, i.e. \(P(M) = P_0 (M/M_0)^\alpha\) for \(M < M_0\).
[factorMassGrowthConsolidate] (default 0.0d0) — The maximum factor by which the mass is allowed to grow between child and parent when consolidating nodes. A non-positive value prevents consolidation.
[destroyStubs] (default var_str('always')) — Parameter controlling when to destroy stub branches. Options are always, never, and sideBranchesOnly.
[rate] — The probability (between 0 and 1) that each particle is retained; particles are drawn independently from a uniform distribution to achieve the target subsampling fraction.

`nbodyOperatorTimeFormation`¶

An N-body data operator which computes progenitor halo order statistics.

Parameters

[massParentMinimum] — The minimum parent mass to consider.
[massParentMaximum] — The maximum parent mass to consider.
[massParentCountPerDecade] — The number of bins per decade of parent mass.
[redshiftFormationMinimum] — The minimum formation redshift to consider.
[redshiftFormationMaximum] — The maximum formation redshift to consider.
[redshiftFormationCountPerUnit] — The number of bins per unit of formation redshift.
[snapshotParents] — The snapshot at which to select parent halos.
[fractionFormation] — The mass fraction used to define the formation epoch.
[description] — A description of this mass function.
[simulationReference] — A reference for the simulation.
[simulationURL] — A URL for the simulation.
[timeRecent] — The parameter \(\Delta t\) (in units of Gyr) appearing in the formation time galactic filter class.

`nbodyOperatorTimeLastMajorMerger`¶

An N-body data operator which computes the time of the last major merger.

Parameters

[massParentMinimum] — The minimum parent mass to consider.
[massParentMaximum] — The maximum parent mass to consider.
[massParentCountPerDecade] — The number of bins per decade of parent mass.
[redshiftMergerMinimum] — The minimum merger redshift to consider.
[redshiftMergerMaximum] — The maximum merger redshift to consider.
[redshiftMergerCountPerUnit] — The number of bins per unit of merger redshift.
[snapshotParents] — The snapshot at which to select parent halos.
[ratioMajor] — The mass ratio used to define a major merger.
[description] — A description of this mass function.
[simulationReference] — A reference for the simulation.
[simulationURL] — A URL for the simulation.

`nbodyOperatorTimeOfFormation`¶

An N-body data operator which computes formation times of halos.

Parameters

[fractionFormation] — The mass fraction used to define the formation epoch.

`nbodyOperatorTimeSinceFormationFractional`¶

An N-body data operator which computes and stores the time since formation as a fraction of crossing time for halos.

`nbodyOperatorVelocityDispersion`¶

An N-body data operator which computes the rotation curve at a set of given radii.

Parameters

[radiusSpecifiers] — A list of radius specifiers at which to output the velocity dispersion.
[includeRadii] (default .false.) — Specifies whether or not the radii at which velocity dispersion data are output should also be included in the output file.
[integrationFailureIsFatal] (default .true.) — If true, failure of line-of-sight integrals is fatal. Otherwise, such errors are tolerated.
[toleranceRelative] (default 1.0d-3) — The relative tolerance to use in integrals.
[selfBoundParticlesOnly] — If true, the velocity dispersion is computed only for self-bound particles.
[radiusInner] — Inner radii of spherical shells within which the velocity dispersion should be computed.
[radiusOuter] — Outer radii of spherical shells within which the velocity dispersion should be computed.
[bootstrapSampleCount] (default 30_c_size_t) — The number of bootstrap resamples of the particles that should be used.

`nbodyOperatorVirialCrossingOrbitStatistics`¶

An N-body data operator which computes statistics of particle pairs on virial-crossing orbits, binning pairs as a function of radial velocity (in virial units) and pairwise separation (in virial radii). This is used to characterize the infall kinematics around dark matter halos. Parameters control the velocity and separation ranges, redshift, Hubble flow correction, velocity outlier cut, sampling rate, and bootstrap resample count.

Parameters

[crossCount] (default .false.) — If true, compute cross-simulation virial crossing orbit statistics between the first and all simulations. Otherwise, compute virial crossing orbit statistics within each simulation.
[addHubbleFlow] (default .false.) — If true, add the Hubble flow contribution \(H(z)\,\mathbf{r}\) to each particle velocity before computing virial crossing orbit statistics, converting peculiar velocities to total line-of-sight velocities.
[redshift] (default 0.0d0) — The redshift at which to evaluate the Hubble flow when adding the Hubble flow correction to particle velocities.
[bootstrapSampleCount] (default 30_c_size_t) — The number of bootstrap resamples of the particles that should be used.
[bootstrapSampleRate] (default 1.0d0) — The fraction of particles to sample randomly at each bootstrap iteration, between 0 and 1; values less than 1 reduce computational cost at the expense of statistical precision.
[velocityMinimum] — The minimum velocity to consider for virial crossing orbit statistics.
[velocityMaximum] — The maximum velocity to consider for virial crossing orbit statistics.
[velocityCount] — The number of bins in separation for virial crossing orbit statistics.
[separationMinimum] — The minimum separation (in units of virial radii) to consider for virial crossing orbit statistics.
[separationMaximum] — The maximum separation (in units of virial radii) to consider for virial crossing orbit statistics.
[includeUnbootstrapped] (default .true.) — If true, include results for the unbootstrapped (i.e. original) sample.
[velocityCut] (default huge(0.0d0)) — Pairs with (radial or tangential) velocities in excess of this value (in virial units) are cut from the sample.

N-Body Simulation Data Operators¶

Methods¶

nbodyOperatorAddAttributes¶

nbodyOperatorAngularMomentum¶

nbodyOperatorConcentrationDistributionFunction¶

nbodyOperatorConvexHull¶

nbodyOperatorConvexHullOverdensity¶

nbodyOperatorConvexHullVolume¶

nbodyOperatorDeleteProperties¶

nbodyOperatorDistanceFromPoint¶

nbodyOperatorEnergyTensors¶

nbodyOperatorEnvironmentalOverdensity¶

nbodyOperatorExportIRATE¶

nbodyOperatorFilterBox¶

nbodyOperatorFilterConvexHull¶

nbodyOperatorFilterID¶

nbodyOperatorFilterProperties¶

nbodyOperatorFilterUncontaminatedSphere¶

nbodyOperatorFlagAlwaysIsolated¶

nbodyOperatorHaloCrossingTime¶

nbodyOperatorHostedRootID¶

nbodyOperatorIdentifyFlybysMansfieldKravtsov2020¶

nbodyOperatorInertiaTensor¶

nbodyOperatorMassFunction¶

nbodyOperatorMassFunctionCorrelation¶

nbodyOperatorMassTotal¶

nbodyOperatorMeanPosition¶

nbodyOperatorMergerRates¶

nbodyOperatorNull¶

nbodyOperatorPairCounts¶

nbodyOperatorPairwiseVelocityStatistics¶

nbodyOperatorPhysicalToComoving¶

nbodyOperatorPotentialEnergy¶

nbodyOperatorProgenitorMassFunction¶

nbodyOperatorProgenitorOrderStatistics¶

nbodyOperatorRotationCurve¶

nbodyOperatorSelectProperties¶

nbodyOperatorSelfBound¶

nbodyOperatorSelfBoundBarnesHut¶

nbodyOperatorSelfFrictionAcceleration¶

nbodyOperatorSequence¶

nbodyOperatorSetBoxSize¶

nbodyOperatorShiftProperty¶

nbodyOperatorSimulationSelector¶

nbodyOperatorSpinDistributionFunction¶

nbodyOperatorSubhaloMassFunction¶

nbodyOperatorSubhaloRadiusFunction¶

nbodyOperatorSubhaloVelocityMaximumMeanFunction¶

nbodyOperatorSubsample¶

nbodyOperatorTimeFormation¶

nbodyOperatorTimeLastMajorMerger¶

nbodyOperatorTimeOfFormation¶

nbodyOperatorTimeSinceFormationFractional¶

nbodyOperatorVelocityDispersion¶

nbodyOperatorVirialCrossingOrbitStatistics¶

`nbodyOperatorAddAttributes`¶

`nbodyOperatorAngularMomentum`¶

`nbodyOperatorConcentrationDistributionFunction`¶

`nbodyOperatorConvexHull`¶

`nbodyOperatorConvexHullOverdensity`¶

`nbodyOperatorConvexHullVolume`¶

`nbodyOperatorDeleteProperties`¶

`nbodyOperatorDistanceFromPoint`¶

`nbodyOperatorEnergyTensors`¶

`nbodyOperatorEnvironmentalOverdensity`¶

`nbodyOperatorExportIRATE`¶

`nbodyOperatorFilterBox`¶

`nbodyOperatorFilterConvexHull`¶

`nbodyOperatorFilterID`¶

`nbodyOperatorFilterProperties`¶

`nbodyOperatorFilterUncontaminatedSphere`¶

`nbodyOperatorFlagAlwaysIsolated`¶

`nbodyOperatorHaloCrossingTime`¶

`nbodyOperatorHostedRootID`¶

`nbodyOperatorIdentifyFlybysMansfieldKravtsov2020`¶

`nbodyOperatorInertiaTensor`¶

`nbodyOperatorMassFunction`¶

`nbodyOperatorMassFunctionCorrelation`¶

`nbodyOperatorMassTotal`¶

`nbodyOperatorMeanPosition`¶

`nbodyOperatorMergerRates`¶

`nbodyOperatorNull`¶

`nbodyOperatorPairCounts`¶

`nbodyOperatorPairwiseVelocityStatistics`¶

`nbodyOperatorPhysicalToComoving`¶

`nbodyOperatorPotentialEnergy`¶

`nbodyOperatorProgenitorMassFunction`¶

`nbodyOperatorProgenitorOrderStatistics`¶

`nbodyOperatorRotationCurve`¶

`nbodyOperatorSelectProperties`¶

`nbodyOperatorSelfBound`¶

`nbodyOperatorSelfBoundBarnesHut`¶

`nbodyOperatorSelfFrictionAcceleration`¶

`nbodyOperatorSequence`¶

`nbodyOperatorSetBoxSize`¶

`nbodyOperatorShiftProperty`¶

`nbodyOperatorSimulationSelector`¶

`nbodyOperatorSpinDistributionFunction`¶

`nbodyOperatorSubhaloMassFunction`¶

`nbodyOperatorSubhaloRadiusFunction`¶

`nbodyOperatorSubhaloVelocityMaximumMeanFunction`¶

`nbodyOperatorSubsample`¶

`nbodyOperatorTimeFormation`¶

`nbodyOperatorTimeLastMajorMerger`¶

`nbodyOperatorTimeOfFormation`¶

`nbodyOperatorTimeSinceFormationFractional`¶

`nbodyOperatorVelocityDispersion`¶

`nbodyOperatorVirialCrossingOrbitStatistics`¶