MIST

MISTChemistry()

Returns a singleton struct representing the MIST chemical mixture model. MIST assumes the protostellar Asplund et al. [5] solar abundances. Sum of protostellar hydrogen, helium, metal mass fractions from last row of Table 4 sums to 0.9999, not 1 as it should. To keep calculations consistent, the protostellar values are normalized to sum to 1 here.

julia> using BolometricCorrections.MIST: MISTChemistry, X, Y, Z, X_phot, Y_phot, Z_phot,
                                         MH;

julia> chem = MISTChemistry();

julia> X(chem) + Y(chem) + Z(chem) ≈ 1 # solar protostellar values
true

julia> X_phot(chem) + Y_phot(chem) + Z_phot(chem) ≈ 1 # solar photospheric values
true

julia> MH(chem, Z(chem) * 0.1) ≈ -1.0189881255814277
true

julia> Z(chem, -1.0189881255814277) ≈ Z(chem) * 0.1
true

MISTBCGrid(grid::AbstractString)

Load and return the MIST bolometric corrections for the given photometric system grid. This type is used to create instances of MISTBCTable that have fixed dependent grid variables ([Fe/H], Av, Rv). This can be done either by calling an instance of MISTBCGrid with (feh, Av) arguments or by using the appropriate constructor for MISTBCTable.

julia> grid = MISTBCGrid("JWST")
MIST bolometric correction grid for photometric system MIST_JWST

julia> grid(-1.01, 0.11) # Can be called to construct table with interpolated [Fe/H], Av
MIST bolometric correction table with [Fe/H] -1.01 and V-band extinction 0.11

MISTBCTable(grid::MISTBCGrid, feh::Real, Av::Real)

Interpolates the MIST bolometric corrections in grid to a fixed value of [Fe/H] (feh) and V-band extinction (Av), leaving only Teff and logg as dependent variables (the MIST BCs have only one Rv value). Returns an instance that is callable with arguments (Teff, logg) to interpolate the bolometric corrections as a function of temperature and surface gravity.

julia> grid = MISTBCGrid("JWST")
MIST bolometric correction grid for photometric system MIST_JWST

julia> table = MISTBCTable(grid, -1.01, 0.011)
MIST bolometric correction table with [Fe/H] -1.01 and V-band extinction 0.011

julia> length(table(2755, 0.01)) == 29 # Returns BC in each filter
true

julia> size(table([2755, 2756], [0.01, 0.02])) # `table(array, array)` is also supported
(29, 2)

julia> using TypedTables: Table # `table(Table, array, array)` will return result as a Table

julia> table(Table, [2755, 2756], [0.01, 0.02]) isa Table
true

This constant is an instance of MISTZeropoints. See the docs for more informations on supported operations. This constant is returned when calling zeropoints on instances of MIST.MISTBCGrid and MIST.MISTBCTable.

Each set of bolometric corrections is specified on either the Vega or AB magnitude system. For ease of conversion amongst the AB, ST, and Vega systems this table is provided containing one line for each filter in our collection. The contents of the file can be used to convert amongst AB, ST, and Vega magnitude systems as follows.

Above is the information for WISE W1, which is tabulated by default in Vega mags (as noted in column "system"). To convert WISE mags from Vega to AB is a simple operation: mag(AB) = mag(Vega) + "mag(Vega/AB)" [column 4 in the file]. For example, a star with WISE_W1(Vega) = 0.0 would have WISE_W1(AB) = 2.66.

To ease compatibility between Symbol and String representations, the mag(Vega/ST) and mag(Vega/AB) columns have been simplified to VegaST and VegaAB.

Struct to contain the MIST zeropoint information. A constant instance is available as BolometricCorrections.MIST.zpt. Instances are callable with a filtername::AbstractString which will return the zeropoint entry for the relevant filter. The full table can be retrieved with Table(zpt) and the list of available filters can be retrieved with filternames(zpt).