2.1. Lines BRASS Data Interface
An important source of uncertainty in stellar spectrum synthesis calculations is the accuracy of atomic data of permitted transitions. It is crucial to constrain atomic data uncertainties for reliable measurements of the thermal conditions and chemical composition of stellar atmospheres. For BRASS we retrieve ∼400,000 transition entries from various online atomic databases: VALD-3, NIST, Spectr-W, TIPbase, TOPbase, CHIANTI, and SpectroWeb. We collect the atomic transition data of neutral species and ions up to 5 for wavelengths between 420 and 680 nm. The datasets are homogenized and cross-matched against the BRASS atomic line list compilation. The BRASS list is composed of Kurucz and NIST V4.0 lines containing for each transition the species (element and ionization stage), line rest-wavelength, log(), upper and lower electronic configurations and energy levels, J-values, and the corresponding literature references. Our cross-matching is performed in two different ways: the parametric cross-match method is based on wavelength- and level energy-values for finding the same transition of a given species. On the other hand, the non-parametric cross-match method is based on detailed electronic configuration information for finding transitions that are physically identical between the datasets. The cross-matching accounts for atomic fine structure, the provided isotopic information, and the type of transition. It however does not account for currently missing hyperfine structure information.
The BRASS compilation was initially tested with theoretical spectrum calculations of the solar flux spectrum [
1] and using Mercator-HERMES [
2] spectra of selected B-, A-, F-, G-, and K-type stars (see [
3]). The BRASS list has been also cleaned from numerous un-observed lines, spurious atomic & molecular background features, and duplicated lines have been excluded. Note that the SpectroWeb atomic lines list was previously compiled from VALD-2 and NIST data (V2.0 through V4.0), and was also extensively tested similar to the BRASS list with theoretical spectrum calculations of high-quality hot and cool star spectra [
4].
Table 1 lists the number of retrieved lines, source databases, dates of retrieval, and various atomic data values collected from each database. We have made extensive use of the online VAMDC portal offering homogenized datasets which has expedited the comparison and cross-matching of the datasets we have retrieved for BRASS. We partly incorporate data from TIPbase and TOPbase and include some of our expansions into fine-structure transitions [
5]. We also calculate line log(
)-values for Spectr-W
using the
-values they offer.
For the BRASS project [
6] have used the non-parametric cross-match method to explore differences between multiple occurrences of identical transitions in the retrieved datasets. Detailed comparisons of
vs.
,
E vs.
E,
vs.
E, and
vs.
log(
)-values mainly reveal the presence of small-scale conversion precision differences. Large-scale systematic correlations are detected for a few cases only. However, the comparison of the line log(
)-values reveals differences in excess of 2 dex, which has important implications for quantitative stellar spectroscopy.
An investigation of duplicated transitions (also accounting for hyperfine-, isotopic-, and E2-M1 forbidden-transitions) in the retrieved datasets show a significant number of almost 2% in VALD-3 lists. These duplicates could be sourced back to the original work in 99% of cases, hence they were not produced by the databases from which the BRASS datasets are retrieved. The duplicated transitions for example have not been detected in the line datasets retrieved from NIST.
The cross-matched atomic datasets, including the BRASS atomic lines compilation, have been incorporated in the online Lines BRASS Data Interface (LBDI) at
brass.sdf.org. Lists of duplicated lines are also offered there for a variety of data formats (HTML, ASCII, PDF). The left-hand panel of
Figure 1 shows the LBDI that can be queried for a given element in a user-defined wavelength interval. In case a cross-match listing for every element is requested the users can set the Element input field to
all. The query results can be sorted by increasing rest-wavelengths or
log(
) (standard deviation)-values marked in blue in
Figure 2. The query results can be exported and saved in extensive line lists or per user-selected line to machine-readable (tab-separated ASCII) tables.
Figure 2 shows for example the cross-matched atomic data of a S
ii line retrieved from seven atomic data sources providing five different log(
)-values ranging from −0.341 dex to −0.059 dex. The literature references of the log(
)-values are offered together with the upper and lower electronic configurations and level energies. The right-hand panel of
Figure 1 also shows a subset of LBDI dynamic plots of the BRASS compilation log(
)-values vs. log(
)-difference values for VALD3-BRASS and NIST-BRASS of cross-matched Fe
i lines and of Fe
ii lines. These dynamic plots can be interactively zoomed and the data of individual lines marked and displayed by mouse interaction. The log(
)-difference plots are provided per query for all atomic data sources and are ordered by neutral, singly, and multiply ionized species (from left to right). This provides users with an interactive and comprehensive overview of all cross-matched log(
) datasets offered in BRASS. Note that the BRASS Data Interface also offers comprehensive Help pages (under the main green tab) for a number of BRASS usecases and corresponding tutorial videos.
2.2. Spectra BRASS Data Interface
For the BRASS project we observe benchmark spectra of a variety of bright stars (
V < 7
) with the HERMES and ESO-VLT-UVES high-resolution spectrographs. We investigate HERMES benchmark spectra of 6 dwarf stars of F-, G-, and K spectral types observed with very high signal-to-noise ratios (SNR) of ∼800–1000: 51 Peg, 70 Oph, 70 Vir, 10 Tau,
Eri, and
Com (
between 5000 K and 6000 K). The spectra are modelled in detail with advanced LTE synthesis calculations using 1-D hydrostatic atmosphere models (see [
3]). The detailed spectrum modeling determines
-, log(
g)-, [M/H]-,
-,
vsin
i-, and [
/Fe]-values we also compare to published stellar parameters measured with high-resolution spectra. The BRASS benchmark stars exclude binaries and are selected for non-variability and non-peculiarity. They are normal dwarf stars with narrow absorption lines having small rotational velocities below 6
and metallicities very close to solar values. Metal-poor stars are excluded to avoid non-LTE effects in our theoretical spectrum calculations.
The spectra of 11 BRASS benchmark stars (including the solar FTS spectrum) and the theoretical spectra are incorporated in the Spectra BRASS Data Interface (SBDI) shown in
Figure 3. Users can interactively display up to four spectral Regions in two benchmark stars selected from the left-hand menus. The wavelengths of identified atomic lines are marked with red and blue labels. The red label numbers mark investigated lines. The red (and blue) labels can be clicked for displaying BRASS atomic data (log(
)- and
-values in a red ’graded list’ in the right-hand panels), together with measured line properties such as the observed line equivalent widths (
and
-error), and the type of quality assessment we performed for the accuracy of the line log(
)- and rest-wavelength -values (see
Section 2.3). By double-clicking the red line labels users can build up lists of BRASS line data values (marked in green) for saving to their local computer disc. Clicking the ‘View data quality’ link in the red (or green) tables populates the ‘Atomic Data Quality’ tab in the central SBDI panel for a complete overview of the atomic line data quality assessment results BRASS offers for the investigated line.
2.3. Atomic Line Data Quality Assessments
A large-scale homogeneous selection of atomic lines was performed by [
7] for BRASS by calculating the theoretical spectra of the 6 FGK benchmark HERMES spectra and the solar FTS spectrum. A selection of 1091 theoretically deep and sufficiently un-blended lines in the wavelength range 420 nm to 680 nm proved to be suitable for advanced quality assessments of the accuracy of the atomic datasets collected in BRASS. We determine astrophysical (semi-empiric) log(
)-values for these 1091 transitions using two commonly employed analysis methods. The agreement of the measured log(
)-values is used for selecting well-behaving lines for the quality assessment work. A total of 845 atomic lines are found to be suitable for quality assessment, of which 408 are robust against any systematic differences between both analysis methods. Around 53% of the quality-assessed lines are found to have at least one literature log(
)-value in agreement with the calculated values, although the remaining values can disagree by as much as 0.5 dex (see
Section 2.4).
For selecting atomic lines we calculate the amounts of blending in the 82337 BRASS lines of the solar and 51 Peg benchmark spectra. To reduce the impact of the line blending amounts on the atomic data quality assessment work a cut-off for blending of ≤10% is used, selected as a good balance between blending of the line core and the number of investigated lines. An additional cut-off on the central line core depth ≥0.02 is also used to ensure the observed line profiles can be measured with sufficient accuracy. A total of 1515 atomic lines is initially selected as ‘un-blended’ lines in both stars. The lines selection procedure does not place limits on the atomic species. The equivalent line widths of the 1515 un-blended lines are automatically measured in all seven benchmark spectra using a single Gaussian fit profile. The line fit procedure is optimized using Gauss-Newton non-linear regression, or Nelder-Mead minimization in the case of slow convergence [
8]. The best fit to the observed line fluxes is limited to the wavelength interval between two local flux maxima in both line wings exceeding 2% of the normalized continuum flux level. Beyond the local flux maxima the
integration is extended for Gaussian line wings. A goodness-of-fit value of
≤ 0.95 is used to remove poorly fitted (in addition to visual inspection), non-existent, too blended, or Earth line contaminated absorption features. The SBDI also offers an interactive
measurement tool under the Gauss line fit tab of the central panel in
Figure 3. Users can select lines in observed BRASS spectra and display the single Gaussian best fit result together with the list of measured line properties for saving to their computer’s disc.
The astrophysical log()-values are determined with two commonly employed methods. The measured line equivalent widths are converted into log()-values using the theoretical curve-of-growth calculated for the line in each benchmark star (called COG method). The other method varies the log()-values in detailed radiative transfer calculations for determining the best-fit value to the observed line profile. The latter method is called GRID because it involves an iterative line modeling procedure for which a grid of spectra is calculated and the best fitting spectrum for a range of log()- and -values is obtained with -minimization by interpolating in steps of (0.005 Å) and log() (0.01 dex) using a bivariate cubic spline fit. Both methods introduce assessable uncertainties resulting from the accuracy of the best fit procedures to the observed -value and the continuum normalized line flux distribution. The uncertainties can be attributed to the spectral SNR, specific atmosphere modeling assumptions, the continuum flux level normalization procedure, and blending with the observed line unaccounted for in our theoretical spectrum calculations. For the seven BRASS benchmark spectra we measure an intrinsic scatter between both methods of ±0.04 dex (1 standard deviation) for line blending levels below 3–4%. The value of ±0.04 dex is therefore used as a constraint on the lines selection for limiting the impact of systematic differences between both methods on the atomic data quality assessment results. In our analysis method the close agreement between the COG and GRID log()-values is required for quality assessing the literature log()-values retrieved for BRASS. The COG log()-value is calculated for a given transition with the observed -value of an absorption feature we can attribute to the line, while the GRID log()-value results from complete theoretical spectrum calculations that fit the observed spectrum incorporating the (sufficiently un-blended) line profile.
The SBDI offers atomic data quality assessment pages showing plots and data values for the 1091 investigated lines.
Figure 4 shows a screencopy of the SBDI Atomic Data Quality tab for the Ni
i 6598 line observed in the BRASS benchmark spectra (solid black line with dots) over-plotted with the theoretical profiles we calculate for the atomic data values retrieved from four atomic databases. The line profiles calculated for log(
)-values we determine from the GRID and COG analysis methods are over-plotted in blue and green colors, respectively. Users can interactively zoom-in, pan, and reset these line profile plots for each benchmark star. By clicking the check boxes above the plots the theoretical line profiles calculated with the log(
)-values in the atomic databases are also over-plotted for user inspection. The Quality assessment table shown below the line profile plots lists the GRID and COG line log(
)- and rest-wavelength (
)-values, together with the differences (
log(
) and
) with respect to the GRID values. The last column of this table offers a Yes/No flag indicating if the
log(
)-value is within the errors of the GRID log(
)-value. The flags and
-values are useful for determining if the log(
)-values retrieved from the databases for BRASS are sufficiently accurate for detailed spectrum synthesis calculations. For example, for Ni
i 6598 we determine GRID log(
)- and COG log(
)-values within errors of each other (hence having quality-assessable atomic data), but not within ±0.04 dex of each other, signaling the line is not robust against the analysis method. The bottom table with Equivalent widths offers the observed (Measured) and theoretical
-values (in mÅ) we calculate for the investigated line per database in all the benchmark stars. Note that we also add small corrections listed for
to the observed line equivalent width values in case the line saturates on the curve-of-growth and Voigt profile corrections are introduced in our best Gaussian line fit procedure.
2.4. Comparison of Atomic Data Quality Assessment Results
We find 845 of the 1091 investigated lines to be quality-assessable, and 408 are also analysis-independent lines. Nearly half of the investigated and quality-assessable lines are of Fe
i, while another ∼10% belong to singly ionized species. The retrieved literature log(
)-values of a quality-assessable line are considered in agreement with our results and can be recommended in theoretical spectrum calculations only in case they agree within the errors of the mean (averaged over all benchmarks) GRID log(
)-value and its standard deviation. We do not consider any literature errorbars because they are not available for the vast majority of investigated lines. In most cases we adopt the mean GRID log(
)-value as the BRASS reference value because the GRID method yields smaller
-values than the COG method. About 53% of the quality-assessable lines have literature log(
)-values in agreement with the mean GRID log(
)-values. A similar percentage of the 408 analysis-independent lines have sufficiently accurate atomic data. The majority of Fe-group species (V
i, Cr
i, Mn
i, Co
i, Ni
i, Ti
i, and Sc
ii, Ti
ii, Fe
ii) have a good number of lines with accurate atomic data for 70–75% of the lines. The Fe
i lines, however, have only ∼38% with sufficiently accurate atomic data (see
Section 2.5).
The right-hand panel of
Figure 5 shows mean GRID log(
)- (blue dots) and mean COG log(
)-values (black dots) compared to the log(
)-values in the BRASS (input) dataset for the 408 analysis-independent lines (where both COG and GRID astrophysical values agree within ±0.04 dex). We find sizable differences with the BRASS log(
)-values for a considerable number of lines. Difference log(
)-values in excess of ±0.5 dex are observed. The inset panel shows lines with smaller log(
) differences (≤0.2 dex), although many are not in agreement within the derived errorbars.
It is important to point out that the large differences between the literature log()-values we retrieve for BRASS and the mean astrophysical log()-values calculated with FGK BRASS benchmark spectra are also detected using a linear approximation method. Absorption lines on the linear part of the curve-of-growth follow a linear relationship between - and log()-values. For these lines the difference between observed and theoretical log()-values equals log(/), where is the line equivalent width we calculate with the theoretical log()-value.
The left-hand panel of
Figure 5 shows the mean of the
log(
)-values we calculate for the seven BRASS benchmark stars against the retrieved BRASS log(
)-values. The largest mean
log(
)-values can also exceed 0.5 dex, although the standard deviations are ≤0.02 dex for the majority of investigated lines (blue errorbars) (see [
8]). Note however that the mean COG and GRID errors are 0.065 dex and 0.05 dex, respectively, or about 3 times larger. Similar to the GRID vs. COG quality assessment method the mean log(
)-differences we calculate with this linear approximation method remain typically below ±1 dex and are chiefly observed for the medium-strong lines having −3 ≤ log(
) ≤ −0.5. The lines with negative log(
)-differences were also found in a separate analysis of Fe-group element lines in the solar FTS spectrum and in HERMES and UVES spectra of Procyon and
Eri [
4]. For these lines the literature log(
)-values are overestimated yielding theoretical
-values that exceed observed values. Similar to the full-fledged GRID vs. COG analysis method smaller
log(
)-values are also found towards the weakest (log(
) < −3.5) and strongest (log(
) > 0) investigated lines.
2.5. Multiplet Analysis of Fe i Transitions in BRASS
The rather small percentage of only ∼38% of sufficiently accurate atomic data for the Fe
i lines in BRASS calls for an investigation of its origin we briefly discuss. The left-hand panel of
Figure 6 shows the curve-of-growth for Fe
i lines we observe in the solar benchmark spectrum. The black dots show observed (reduced)
-values against log(
) (co-added with other terms), for log(
)-values in the (input) BRASS compilation of
Table 1. We find considerable scatter for the transitions on the linear part of the curve or mainly for the weak and medium-strong Fe
i lines. The large scatter is due to the limited accuracy of the literature log(
)-values for these lines. We find that this scatter across the curve substantially reduces after replacing the literature log(
)-values with the ones we calculate from the linear approximation method in
Section 2.4 shown with red dots. By replacing the log(
)-values with the ones we calculate from the COG vs. GRID method the scatter nearly vanishes and the curve assumes the smooth (and narrow) shape required for atomic lines belonging to the same species in stellar spectra. The large percentage we find of over 60% of literature Fe
i atomic data with limited quality mainly results from medium-strong (and weak) lines having −3 ≤ log(
) ≤ −0.5 in
Figure 5.
A more extensive analysis of the Fe
i fine structure data we retrieve for BRASS reveals that the lines with limited/poor log(
) quality mostly have
> 4 eV. We compile 25 electric dipole multiplets of 69 Fe
i lines in BRASS, also shown in the right-hand panel of
Figure 6. For each line of these multiplets we calculate the relative line strength ratios assuming single-configuration Russell–Saunders (LS) coupling and obeying the selection rules for these permitted transitions. The calculated multiplet line strengths (using Wigner 6j-symbol calculations) are normalized by scaling the LS-coupling log(
)-values to the strongest available principal line (marked with
or
in
Figure 6), or the largest log(
)-value we calculate with the GRID method using the benchmark spectra. For a number of multiplets we find reasonable to good agreement between the LS and GRID relative log(
)-values. For example multiplet
(
= 2.17–2.23 eV) shows very similar distributions across its principal (
x) and satellite (
y and
z) transitions. For multiplet
(
= 4.73–4.84 eV) the relative log(
) distributions agree less, but also show differences between the literature and GRID log(
)-values of ∼1.0 dex. This is also the case for
(
= 4.9–5.1 eV) and
(
= 4.1–4.23 eV) multiplets for which the relative LS and GRID log(
)-distributions across the
x,
y, and
z line series are dissimilar and a re-normalization cannot remove the large differences above 0.5–1.0 dex for individual lines. The right-hand panel of
Figure 6 shows that the differences between the normalized LS and the (best) BRASS GRID log(
)-values increase towards larger
-values. For lines having
> 4 eV in the 25 studied multiplets the differences can increase above 0.5–1.0 dex mainly for the satellite transitions (marked
y).
Using the interactive NIST Grotrian diagrams we find that both the and multiplet energy levels have level components that lie very close in energy to neighboring high energy levels of other atomic terms. For example the (J = 1) level at 4.217 eV falls next to another energy level at 4.220 eV corresponding to the (J = 5) level. In addition, the (J = 2) level at 5.085 eV lies near a (J = 2) level at 5.070 eV. The proximity of other nearby energy levels for these multiplet lower levels yields significant configuration interaction between the levels. The LS-coupling calculations cannot accurately predict the relative line strengths assuming single electronic configuration interaction. The increasing differences between the literature, LS and GRID log()-values towards larger results from inaccurate theoretical Fe i log()-values due to poorly constrained configuration mixing coefficients and inaccurate or incomplete theoretical energy levels for the large number of close energy levels above ∼4 eV in the neutral Fe atom.