ATLAS9

This submodule (BolometricCorrections.YBC.ATLAS9) enables interaction with the "YBC" bolometric correction (BC) grid calculated from the ATLAS9 model atmospheres of Castelli and Kurucz (Jan 2003), specifically the "ODFNEW" atmospheres. Details on the BC calculations are given in section 3.1 of Chen et al. (2019). A webpage hosting and describing the atmospheres is maintained by StSCI here and SVO here. Robert Kurucz's webpage hosting the grid is here. These atmosphere models have been widely used since their release as they cover a wide range of stellar parameters and provide good fits to observed spectra of hotter stars (see, e.g., Bertone et al. (2004)). Other atmospheres like PHOENIX and MARCS may be preferred for lower temperature stars (for example, Chen et al. (2019) use PHOENIX for $T_e < 5500$ K and ATLAS9 for $T_e > 6500$ K, interpolating between them in the overlapping region). The full range of dependent variables covered by this grid is given in the table below.

minmax
Teff3,467 K50,118 K
logg0.05.0
[Fe/H]-2.5 dex0.5 dex
Av0.0 mag20.0 mag
Rv3.13.1

The full grid of unique values for the dependent variables is available in BolometricCorrections.YBC.ATLAS9.gridinfo.

using BolometricCorrections
keys(BolometricCorrections.YBC.ATLAS9.gridinfo)
(:logTeff, :logg, :MH, :Av, :Rv)
Example block output

Types

BolometricCorrections.YBC.ATLAS9.ATLAS9YBCGridType
ATLAS9YBCGrid(grid::AbstractString)

Load and return the YBC ATLAS9 bolometric corrections for the given photometric system grid, which must be a valid entry in BolometricCorrections.YBC.systems. This type is used to create instances of ATLAS9YBCTable that have fixed dependent grid variables ([M/H], Av). This can be done either by calling an instance of ATLAS9YBCGrid with (mh, Av) arguments or by using the appropriate constructor for ATLAS9YBCTable.

julia> grid = ATLAS9YBCGrid("acs_wfc")
YBC ATLAS9 bolometric correction grid for photometric system YBC/acs_wfc.

julia> grid(-1.01, 0.11) # Can be called to construct table with interpolated [M/H], Av
YBC ATLAS9 bolometric correction table for system YBC/acs_wfc with [M/H] -1.01 and V-band extinction 0.11
source
BolometricCorrections.YBC.ATLAS9.ATLAS9YBCTableType
ATLAS9YBCTable(grid::ATLAS9YBCGrid, mh::Real, Av::Real)

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

ATLAS9YBCTable(grid::AbstractString, mh::Real, Av::Real)

Loads the data necessary to construct the BC table for the provided grid (e.g., "acs_wfc") at [M/H] = mh and V-band extinction Av. This method does not support interpolation in metallicity or extinction, so the arguments mh and Av must be among the values provided by ATLAS9 (see BolometricCorrections.YBC.ATLAS9.gridinfo.MH).

julia> grid = ATLAS9YBCGrid("acs_wfc")
YBC ATLAS9 bolometric correction grid for photometric system YBC/acs_wfc.

julia> table = ATLAS9YBCTable(grid, -1.01, 0.011) # Interpolate table from full grid
YBC ATLAS9 bolometric correction table for system YBC/acs_wfc with [M/H] -1.01 and V-band extinction 0.011

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

julia> size(table([4025, 4225], [0.01, 0.02])) # `table(array, array)` is also supported
(12, 2)

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

julia> table(Table, [4025, 4225], [0.01, 0.02]) isa Table
true

julia> ATLAS9YBCTable("acs_wfc", -2.0, 0.5) isa ATLAS9YBCTable # Can construct without a ATLAS9YBCGrid
true
source

See here for a list of the photometric systems available that can be used as input arguments to ATLAS9YBCGrid and ATLAS9YBCTable.

Chemistry

These models use solar chemical abundances from Grevesse and Sauval (1998). These results do not consider metal diffusion, so the photospheric solar abundances are the same as the primordial solar abundances (i.e., X(ATLAS9Chemistry) == X_phot(ATLAS9Chemistry)). All metallicities have scaled-solar abundance patterns.

BolometricCorrections.YBC.ATLAS9.ATLAS9ChemistryType
ATLAS9Chemistry()

Returns a singleton struct representing the ATLAS9 chemical mixture model. ATLAS9 assumes the Grevesse and Sauval (1998) solar abundances for which photospheric abundances are equal to protostellar abundances – this is what they mean when they write "the effects of element migration at the bottom of the convective zone ... are not observed." In ATLAS9, the helium abundance is not scaled with the metallicity, so only the hydrogen mass fraction and metal mass fraction change as a function of metal abundance (i.e., dY/dZ = 0). See also Girardi et al. (2007).

julia> using BolometricCorrections.YBC.ATLAS9: ATLAS9Chemistry, X, Y, Z, X_phot, Y_phot, Z_phot, MH;

julia> chem = ATLAS9Chemistry();

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> X(chem) == X_phot(chem) # photospheric and protostellar abundances equal (assume no diffusion)
true

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

julia> Z(chem, -1.00894760736597) ≈ Z(chem) * 0.1
true
source
chem = chemistry(ATLAS9YBCTable)
(X = X(chem), Y = Y(chem), Z = Z(chem))
(X = 0.735, Y = 0.248, Z = 0.017)

ATLAS9 References

This page cites the following references: