vice.milkyway.default_mass_loading

The default mass loading factor as a function of galactocentric radius in kpc.

Signature: vice.milkyway.default_mass_loading(rgal)

Parameters

rgalreal number

Galactocentric radius in kpc.

Returns

etareal number

The mass loading factor at that radius, defined by:

\[\eta(r) = (y_\text{O}^\text{CC}) / Z_\text{O}^\odot 10^{0.08(r - 4\text{ kpc}) - 0.3} - 0.6\]

where \(Z_\text{O}^\odot\) is the solar abundance by mass of oxygen and \(y_\text{O}^\text{CC}\) is the IMF-averaged CCSN yield of oxygen. While these values are customizable through vice.solar_z and vice.yields.ccsne.settings, this function assumes a value of \(Z_\text{O}^\odot\) = 0.00572 (Asplund et al. 2009 [1]) and \(y_\text{O}^\text{CC}\) = 0.015 (Johnson & Weinberg 2020 [2], Johnson et al. 2021 [3]).

See also

vice.milkyway.mass_loading

Example Code

>>> import vice
>>> vice.milkyway.default_mass_loading(0)
0.029064576665950193
>>> vice.milkyway.default_mass_loading(8)
2.1459664721614495