This site uses cookies. By continuing to use this site you agree to our use of cookies. To find out more, see our Privacy and Cookies policy.

Vanadium Transitions in the Spectrum of Arcturus

, , , , , and

Published 2018 February 1 © 2018. The American Astronomical Society. All rights reserved.
, , Citation M. P. Wood et al 2018 ApJS 234 25 DOI 10.3847/1538-4365/aa9a41

Download Article PDF
DownloadArticle ePub

You need an eReader or compatible software to experience the benefits of the ePub3 file format.

0067-0049/234/2/25

Abstract

We derive a new abundance for vanadium in the bright, mildly metal-poor red giant Arcturus. This star has an excellent high-resolution spectral atlas and well-understood atmospheric parameters, and it displays a rich set of neutral vanadium lines that are available for abundance extraction. We employ a newly recorded set of laboratory FTS spectra to investigate any potential discrepancies in previously reported V i log(gf) values near 900 nm. These new spectra support our earlier laboratory transition data and the calibration method utilized in that study. We then perform a synthetic spectrum analysis of weak V i features in Arcturus, deriving log ε(V) = 3.54 ± 0.01 (σ = 0.04) from 55 lines. There are no significant abundance trends with wavelength, line strength, or lower excitation energy.

Export citation and abstract BibTeX RIS

1. Introduction

A fundamental astrophysical goal is to derive detailed chemical compositions of Milky Way populations. What are these compositions as functions of time and spatial distribution in the Galaxy? How many different variables are needed to define major populations, and how many heterogeneous sub-groups really exist in each of these populations? Answers to these questions depend heavily on determining accurate elemental abundances for as many stars as possible.

In the past few decades there have been substantial improvements in both the quantity and quality of stellar abundances, and these will only improve in the future through the accumulation of data from many large-sample surveys (e.g., APOGEE, Majewski et al. 2017; Gaia-ESO, Gilmore et al. 2012; GALAH, Freeman et al. 2012). Better spectrographs, larger telescopes, more detailed model stellar atmospheres and line formation physics, and more robust atomic and molecular data have all contributed to substantial improvements in abundance accuracy. More recently, the development of HgCdTe (HAWAII) detector arrays has opened new infrared (IR) wavelength regions for astrophysics (e.g., APOGEE, Holtzman et al. 2015; IGRINS, Yuk et al. 2010). These regions are of particular importance in stellar astrophysics, where IR lines are generally better isolated and the continuum level better understood and easier to determine. This blossoming of IR astronomy puts pressure on the laboratory astrophysics community to keep pace.

Our group has concentrated on producing laboratory atomic data for as many lines as possible of more than 25 species of neutron-capture elements (atomic number Z > 30) and iron-group elements (21 ≤ Z ≤ 30). In our papers we have applied these new lab data (transition probabilities, hyperfine and isotopic substructures) to derive improved abundances of the solar photosphere and selected very metal-poor stars, those with [Fe/H] < −2.6 The choices of these stellar targets are well motivated astrophysically but they do not exploit the new lab data fully. In stellar spectra the number and strength of atomic absorption lines increase toward shorter wavelengths, and for those targets most abundance information is obtained at wavelengths shortward of Hα (656.3 nm). For example, the Sun has mostly weak atomic lines in the red and IR spectral regions, and many very metal-poor stars have almost no spectroscopic information in these spectral regions.

However, our lab data include many lines in the red and IR regions. To fully utilize these data we have begun to explore their application to the spectrum of the cool, mildly metal-poor red giant Arcturus, and in this paper we concentrate on neutral vanadium. Lawler et al. (2014; hereafter La14) reported new transition probabilities for 836 V i lines. Recently, Holmes et al. (2016; hereafter Ho16) also determined transition probabilities for 208 lines of this species. For most lines in common between these two studies the results are in good agreement, but Ho16 suggest that the La14 log(gf) values are systematically too small for six lines in common near 900 nm. Lines in this region are often used to construct calibration bridges between optical and IR spectra, so any issue with these lines needs to be resolved before laboratory IR work on V i beyond 1 μm can proceed confidently. As a result, this disagreement warrants further investigation.

In Section 2 we investigate the discrepancy noted by Ho16 using a newly recorded set of laboratory spectra to independently check our earlier calibration method, finding no cause for revision of the La14 data. In Section 3 the transition probabilities from La14 are applied to the spectrum of Arcturus. Finally, in Sections 45 we discuss the new vanadium abundance in Arcturus in relation to previous work and suggest fruitful directions for future laboratory and stellar spectroscopy in this area.

2. Laboratory Branching Fraction (BF) Data Revisited

Ho16 recently reported transition probabilities for 208 lines of V i, of which 188 are common to and in generally good agreement with our earlier work (La14). However, Ho16 noted a discrepancy for six lines in common around 900 nm, indicating that the La14 results for these transitions are too small by 0.18 dex on average and should be corrected. The noted discrepancy is limited to the z4D°–a4P multiplet near 900 nm, and the relative strengths of the lines in the multiplet are not in dispute, since these lines are in a nearly pure LS or Russell–Saunders multiplet. The multiplet BFs are weak, BF < 0.04, compared to the dominant branches near 480 nm from the upper z4D° levels. The total BF for the multiplet is found to be approximately 0.032 in La14 and approximately 0.048 in Ho16. This overall 0.016 BF discordance has very little effect on the log(gf) values of other lines from the upper z4D° levels.

To check the discordant multiplet, several new spectra (detailed in Table 1) of hollow cathode (HC) lamps were recorded using the National Institute of Standards and Technology (NIST) 2 m Fourier transform spectrometer (FTS). Two of these spectra were recorded with sealed HC lamps from Heraeus7 run at a current of 20 mA, while the rest result from a custom water-cooled high-current HC lamp. A tungsten standard lamp provided a relative radiometric calibration of these spectra, similar to the method employed by Ho16. One spectrum of a sealed V/Ar lamp (V032417.001) was recorded that provided independent confirmation of the calibration using Ar branching ratios (BRs) as described below. However, the vanadium lines of interest in this spectrum were relatively weak and had uncertainties of 15%–20%.

Table 1.  Newly Recorded Fourier Transform Spectra of Hollow Cathode (HC) Lamps Recorded Using the NIST 2 m FTS

Index Spectruma Serial Number HC Lamp Typeb Buffer Gas Lamp Current (mA) Limit of Resolution (cm−1) Coadds
A Ti041715b 001 Water-cooled Ar 500 0.015 25
B V032417 001 Sealed Ar 20 0.08 64
C V062117a 001 Sealed Ne 20 0.025 32
D V062317a 001 Water-cooled Ne 200 0.03 100
E V062317a 002 Water-cooled Ne 700 0.025 115
F V062617b 001 Water-cooled Ar 250 0.02 95
G V062617b 002 Water-cooled Ar 500 0.02 60

Notes. All spectra were recorded using Si detectors, silver optics, and an aluminum beamsplitter, and have a usable wavenumber range of approximately 8500–24,000 cm−1.

aThe date can be determined from the numbers in the spectrum name (MM/DD/YY). bLamp types include commercially available sealed HC lamps from Heraeus and a high-current water-cooled HC lamp with a demountable window.

Download table as:  ASCIITypeset image

The log(gf) agreement between the two earlier studies for shorter wavelength lines from the four z4D° levels is very good, 〈log(gf)Ho16–log(gf)La14〉 = 0.00 ± 0.03,8 and it is therefore reasonable to use BRs to make a relative comparison between the multiplet lines in question and the other lines from these levels. We choose specifically to calculate the BR of lines from the z4D°–a4P multiplet around 900 nm compared to reference lines from the z4D°–a4D multiplet around 820 nm. BRs determined from spectra C through G in Table 1 are shown in Figure 1, where each panel represents a separate upper level and the different symbols represent BRs from Ho16, La14, and the present study. The error bars on the Ho16 results are derived using their quoted BF uncertainties. Error bars on the La14 results are derived in a similar way, with the 0.001%/cm−1 calibration uncertainty with respect to the dominant line (see La14 for details) removed from the BF uncertainties. The error bars on the BRs from this study include uncertainties in the integrated line intensities and in the relative radiometric calibration of the spectra. Each panel includes a single line for which all three studies are in perfect agreement at 1.0; this represents the reference line against which BRs are calculated for each upper level, therefore this line does not include an error bar. These new spectra were analyzed using different software from that in La14, which, when combined with the different instrumentation and calibration method, results in an independent test of the La14 results. As expected, the BRs from La14, Ho16, and the present study are in good agreement for the lines near 820 nm. However, the Ho16 BRs for the 900 nm multiplet are systematically above the BRs from both La14 and this study, which are still in good agreement. This independent check of the La14 results suggests that the discordance for this multiplet is confined to the Ho16 results.

Figure 1.

Figure 1. Branching ratios for lines from each of the upper z4D° levels to the ${a}^{4}{D}$ and ${a}^{4}{P}$ lower levels. The orange squares, black crosses, and green circles are results from Ho16, La14, and the present study, respectively. The associated error bars are included for all three studies. The single line for which all three symbols perfectly overlap at 1.0 is the reference line against which all BRs for a given upper level (specified in the upper right corner) are calculated.

Standard image High-resolution image

Uncertainties in BF studies covering a large wavelength range are often dominated by systematic uncertainties in the relative radiometric calibration. Ho16 used a tungsten standard lamp to determine their calibration, while the relative calibration of the La14 log(gf)s is based on Ar i and Ar ii BRs from Whaling et al. (1993) and earlier work by Hashiguchi & Hasikuni (1985), Danzmann & Kock (1982), and Adams & Whaling (1981). Ho16 claim the 0.18 dex discrepancy is the result of an error in this Ar BR calibration method, so it is worthwhile to compare the advantages and disadvantages of each method. The tungsten lamp is advantageous since it provides a continuum calibration. The response of the NIST 2 m FTS changes rapidly between 800 and 900 nm due to a dip in the reflectivity of aluminum in this region and Ar BRs may not have sufficient density to keep track of the changing response (see Figure 2). However, the tungsten lamp does not take into account the transmission of the window of the HC lamp. A separate measurement of the window transmission can be made if the lamp can be dismantled, but this is not possible with the sealed lamps used in the present study. The Ar BR method is advantageous in this case since it is recorded as part of the data spectrum and therefore the light is guaranteed to follow the same optical path, including through the HC lamp window. Additionally, any possible wavelength-dependent reflection/scattering of light in the lamp is automatically captured in the Ar BR calibration (see Lawler & Den Hartog 2018). Hence, the most reliable calibration can be obtained by combining the two methods, provided that the lines of interest can be excited in a lamp containing Ar.

Figure 2.

Figure 2. A relative radiometric calibration of the 2 m FTS at NIST based on 107 Ar i and Ar ii lines (blue open circles) to a W strip standard lamp calibration (black line) in the near-IR and optical. The gray shaded region represents the uncertainty in the standard lamp calibration combined in quadrature with the HC window transmission uncertainty. The error bars represent individual Ar line intensity uncertainties combined in quadrature with the gain correction uncertainty.

Standard image High-resolution image

Unfortunately, no standard lamp spectra were recorded in conjunction with the V spectra from the National Solar Observatory digital archive analyzed by La14, but an independent check of the Ar BRs has recently been made. Figure 2 presents a comparison between a relative radiometric calibration based on 107 Ar i and Ar ii lines (recorded simultaneously with a Ti HC lamp spectrum) and that of a tungsten standard lamp for spectrum A in Table 1. These earlier data were recorded using the NIST 2 m FTS as part of a separate study to investigate the primary tungsten standard for the FTS using the Ar BRs as a secondary calibration. Care was taken to ensure the FTS setup was similar between the standard lamp spectrum and Ar spectrum, with the only changes being the orientation of the entrance mirror to select each lamp, the spectral resolution, and the gain of the FTS detectors. A change in the spectral resolution has no effect on the response of the spectrometer, since it varies slowly with wavenumber. Changing the gain of the detectors can affect the response of the spectrometer if the amplifier gain settings have different bandwidths, but this is necessary due to the large difference in intensity between our standard lamp and the HC lamps. We have therefore measured the bandwidth of our detectors at different gain settings and corrected our response curve accordingly. In addition, two other corrections are incorporated in Figure 2 to ensure a fair comparison between the two calibrations: (1) the standard lamp includes a correction to account for the transmission of the HC window, and (2) the standard lamp includes an aging correction determined by comparison to a similar calibrated lamp with fewer hours of usage. There is generally good agreement between the Ar BR and tungsten standard lamp calibrations over the full calibrated range of the spectrum from 8500 to 24,000 cm−1. While there exist a few outliers, the large number of Ar calibration lines makes a comparison relatively straightforward. In particular, Figure 2 does not show a large systematic offset of the Ar BRs near 900 nm (≈11,000 cm−1) and does not indicate the need for a correction to the Ar BR calibration as suggested by Ho16. A similar comparison of Ar BRs from spectrum B in Table 1 and shown in Figure 3 indicates that the calibration with the tungsten standard lamp agrees with that derived from the Ar BRs to within ±10% in the region from 750 to 1100 nm (9000–13,000 cm−1). Ar BRs derived from this spectrum are compared to the literature values detailed above, and a value of 1.0 represents perfect agreement between derived and literature BRs. Due to the use of a sealed lamp with this spectrum, a window transmission correction could not be determined and applied to Figure 3, but the good agreement indicates that the change in transmission of the window across this region is small. The BRs of the vanadium lines from 800 to 900 nm in this spectrum also agree with La14 within their uncertainty of 15%–20%.

Figure 3.

Figure 3. Ratio of Ar BRs derived from NIST 2 m FTS spectra calibrated using a tungsten standard lamp to Ar BRs in the literature (see the text for details). The different colors represent BRs from different upper levels in Ar i and Ar ii. A ratio of 1.0 indicates that the BRs derived using the standard lamp calibrated spectra are in perfect agreement with the literature. The symbols and colors represent different line groups, with numbers taken from Whaling et al. (1993).

Standard image High-resolution image

Based on the independent verification of the La14 results and multiple checks of the Ar BR calibration method employed in that study, the oscillator strengths from La14 do not require any adjustment as suggested by Ho16.

3. The Vanadium Abundance of Arcturus

Arcturus is a very bright red giant star (V  = −0.05, V–K = 2.96) that has become an abundance standard for cool star spectroscopy. For example, Arcturus and the Sun are used as reference stars by the APOGEE team to test and improve infrared (H-Band) spectroscopic data (Shetrone et al. 2015). The optical spectrum of this star was studied extensively with the early photographic atlas of Griffin (1968). Arcturus now has an excellent public high-resolution, high signal-to-noise electronic spectrum9 covering the ultraviolet (115–380 nm; Hinkle et al. 2005), optical (373–930 nm; Hinkle et al. 2000), and IR (0.9–5 μm; Hinkle & Wallace 2005) spectral regions. Hereafter, we refer collectively to these publications as the Arcturus Atlas. This spectrum is so line-rich that it is difficult to do effective abundance work at wavelengths below 500 nm. The line density is not as high in the yellow-red spectral region (λ > 500 nm), and inspection of the Arcturus Atlas reveals many strong and relatively unblended V i absorption features in this region.

La14 determined transition probabilities for 219 V i lines in the 500–920 nm wavelength domain. We searched the Arcturus Atlas for all of these lines, ignoring the line selection of La14 for their study of the Sun and the metal-poor star HD 84937. We first eliminated transitions that are severely blended and/or vanishingly weak. As in the past papers of this series, to aid this effort we used the solar line compendium of Moore et al. (1966) and the comprehensive Kurucz (2011) atomic and molecular line database.10 This process yielded about 100 V i lines that were judged worthy of further inspection. All of these remaining lines have λ > 550 nm, since line blending by atomic and molecular (C2, MgH) features block potentially useful lines blueward of this wavelength.

We generated synthetic spectra of small spectral regions around each of the remaining V i lines with the current version of the local thermodynamic equilibrium (LTE) plane-parallel line analysis code MOOG (Sneden 1973).11 Arcturus spectral analyses generally yield similar model parameters. In the PASTEL stellar parameter catalog (Soubiran et al. 2016)12 the simple means of the 42 entries for Arcturus yield effective temperature $\langle {T}_{\mathrm{eff}}\rangle $ = 4290 ± 100 K, surface gravity $\langle \mathrm{log}\,g\rangle $ = 1.6 ± 0.3, and metallicity $\langle [\mathrm{Fe}/{\rm{H}}]\rangle $ = −0.5 ± 0.1 (microturbulent velocities vt are not tabulated in this database). Recent detailed analyses include those of Ramírez & Allende Prieto (2011), who derived Teff = 4286 K, log g = 1.66, vt = 1.74 km ${{\rm{s}}}^{-1}$, and [Fe/H] = −0.52, and Peterson et al. (2017), who recommended Teff = 4275 K, log g = 1.3, vt = 1.6 km ${{\rm{s}}}^{-1}$, and [Fe/H] = −0.55. We adopted the latter model for our abundance computations, and comment on this choice below. We interpolated a model atmosphere with the Peterson et al. (2017) parameters from the Kurucz (2011) model atmosphere grid13 using software kindly provided by Andy McWilliam and Inese Ivans.

Line lists for the synthetic spectra were generated as described in our previous papers (Sneden et al. 2016, and references therein). To summarize, we began with the atomic and molecular lines of the Kurucz (2011) database in 0.4 nm spectral intervals centered on the V i transition of interest. We then substituted or added to these lists the atomic transitions from our group's previous publications on iron-group and neutron-capture element species, CN from Sneden et al. (2014), C2 from Ram et al. (2014), and MgH from Hinkle et al. (2013). As in La14, we expanded the V i lines to include hyperfine structure (HFS) components when these data were available. For this we searched the atomic physics literature anew, and in the Appendix we discuss the HFS data and provide a table that can be useful for synthetic spectrum computations.

Comparisons between these syntheses and the Arcturus Atlas resulted in the elimination of many candidate V i lines. Rejection of a line would result from one or more issues that made it unsuitable for abundance determinations: (1) the observed line was too weak (reduced equivalent width log(RW) = log(EW/λ) < −6.2); (2) the observed line was too blended with other atomic and/or molecular species (contaminant feature >20% of the total feature); (3) the synthesized line had a predicted strong blend in the Kurucz (2011) database that could not justifiably be neglected, even if it was not apparent in the observed line; (4) the Arcturus Atlas feature was too contaminated with very strong telluric O2 or H2O features; or (5) the observed line was far too strong. The first four reasons for line rejection are standard for this kind of abundance work. Our justification for neglecting some of the strongest V i lines was to ensure that uncertainties in the assumed microturbulent velocity would not play a significant role in derived abundances, as discussed below.

After eliminating these problematic lines, we subjected the remaining V i features to more detailed comparisons between synthetic and observed spectra. In Figure 4 we show examples of these comparisons. Each panel shows points from the Arcturus Atlas as well as syntheses assuming various V i abundances (the derived log ε abundance for the feature is noted in the panel). A red line shows the synthesis without any V i contribution, and a black line (often masked by the points) shows the best-fitting synthesis/observation match. The blue, green, and orange lines are for syntheses that assume V abundances that are factors of 1/4, 1/2, and 2 times the best-fit abundance, respectively (e.g., Δlog ε values that are −0.6, −0.3, and +0.3 shifted from the derived best-fit abundance for the feature). Panels (a)–(c) show representative lines in the yellow-red spectral region, and panels (d)–(f) show lines near 900 nm for which Ho16 derived substantially different transition probabilities than La14. The gaps in the Arcturus Atlas points in panels (c) and (f) are small wavelength intervals where telluric absorptions are so deep that the underlying stellar spectrum cannot be reliably detected.

Figure 4.

Figure 4. Representative Arcturus Atlas spectra (open circles) of V i transitions and synthetic spectra (lines). In each panel the black line represents the best-fit synthesis to the observation, and its assumed V abundance is written in the lower right corner of the panel. The orange line is the synthesis with the V abundance increased by 0.3 dex, the green line is for a decrease of 0.3 dex, the blue lines are for a decrease of 0.6 dex, and the red line is for no V contribution to the feature.

Standard image High-resolution image

In our analysis we dealt only with relatively weak V i lines, ones with sensitivity to abundance but not microturbulent velocity vt. We computed Arcturus V i curves of growth with our model atmosphere; these predict equivalent width (EW) as a function of assumed abundance for lines with no HFS substructure. The tests suggest that the very weak-line portion of the curve of growth extends up to residual width log(RW) = log(EW/λ) ≈ −5.5. Residual widths show gradually decreasing sensitivity to abundance and increased dependence on assumed vt as line strength increases, and abundances for lines with log(RW) > −4.75 are greatly dependent on the assumed microturbulent velocity. We measured approximate EWs for all potential transitions, and then conservatively chose to neglect V i lines stronger than log(RW) > −4.8.

The wide HFS of many V i lines complicates statements about line strengths. HFS varies substantially from transition to transition, as illustrated in Figure 4. It is nearly negligible for the 554.5 nm line (panel a), but is about 0.05 nm for the 675.3 and 893.3 nm lines (panels b and d). For lines on the linear portion of the curve of growth, HFS broadens their profiles but does not change their total absorptions (and thus their EWs). However, HFS desaturates strong lines, increasing their total absorption while moving them toward the linear part of the curve of growth. Many of the stronger lines included in the figure have so much HFS broadening that they effectively act as unsaturated transitions. However, lines with log(RW) > −4.8 systematically yield abundances about 0.05–0.15 dex lower than the weaker transitions. This trend can be canceled by decreasing the microturbulent velocity to vt ≈ 1.5 km ${{\rm{s}}}^{-1}$, but this would be a purely empirical correction. Exploration of this issue would involve a more thorough investigation of line formation in the Arcturus atmosphere and is beyond the scope of this paper.

In the end we employed 55 V i lines in the 550–920 nm spectral region to determine a new V abundance in Arcturus. Individual line abundances are listed in Table 2, derived from both the log(gf) values in La14 (columns 3 and 4) and in Ho16 (columns 5 and 6). From these lines we compute a mean value of $\langle \mathrm{log}\,\varepsilon \rangle $ = 3.54 ± 0.01 (sample standard deviation σ = 0.04). La14 derived a solar photospheric abundance of $\langle \mathrm{log}\,{\varepsilon }_{\odot }\rangle \,=3.96\pm 0.01$ (σ = 0.04), which yields [V/H] = −0.42, or [V/Fe] = +0.13. In Figure 5 we plot the individual line abundances as a function of wavelength. The red symbols are used to indicate the lines near 900 nm for which Ho16 derived log(gf) values that are on average 0.18 dex larger than those of La14. In this figure we write the abundance statistics separately for transitions below and above 890 nm. There is excellent agreement between the 49 lines at shorter wavelengths and the 6 lines near 900 nm. In short, there is no abundance trend with wavelength in Arcturus when the La14 transition probabilities are employed, as expected from the small σ value of the mean abundance statistics. This agreement of derived V abundances near 900 nm and below this wavelength cannot be duplicated if the Ho16 log(gf) values are employed.

Figure 5.

Figure 5. Abundances from V i lines derived using La14 transition data, plotted as a function of wavelength. Abundances below and above 890 nm are color-coded, and their means and standard deviations are written in the figure. The solid horizontal line is drawn at the mean value for the lines with λ < 826 nm, and the dotted lines are drawn at ±1σ.

Standard image High-resolution image

Table 2.  Vanadium Abundances from V i lines in Arcturus

λ Elower log(gf) log ε log(gf) log ε
(nm) (eV) La14 La14 Ho16 Ho16
554.4858 1.051 −2.56 3.60 −2.57 3.61
555.7450 0.017 −3.45 3.52 −3.47 3.54
556.0547 0.040 −3.63 3.53 ... ...
557.6502 1.064 −2.45 3.57 −2.49 3.61
557.8373 1.051 −2.58 3.55 ... ...
558.4497 1.064 −1.18 3.48 ... ...
558.8457 1.081 −2.25 3.55 −2.30 3.60
559.2960 0.040 −3.21 3.46 −3.24 3.49
559.7797 1.064 −2.42 3.61 −2.57 3.76
560.4935 1.043 −1.26 3.48 −1.18 3.40
563.2454 0.069 −3.23 3.51 −3.29 3.57
568.7753 2.379 −1.14 3.58 ... ...
572.5639 2.365 −0.17 3.55 ... ...
575.0642 1.955 −1.29 3.56 ... ...
577.2412 1.931 −0.55 3.53 ... ...
577.6687 1.081 −1.54 3.53 ... ...
600.2303 1.218 −1.79 3.51 ... ...
601.7920 1.195 −2.39 3.59 −2.34 3.54
605.8142 1.043 −1.40 3.48 −1.37 3.45
608.7496 1.051 −2.51 3.53 ... ...
618.1862 1.051 −2.45 3.54 ... ...
619.0506 1.064 −2.53 3.53 ... ...
624.3107 0.301 −0.94 3.48 −0.95 3.49
628.5160 0.275 −1.54 3.47 −1.55 3.48
632.6836 1.868 −0.76 3.53 ... ...
634.9473 1.853 −0.95 3.57 ... ...
635.5572 2.122 −1.42 3.58 ... ...
635.7291 1.849 −1.06 3.58 ... ...
645.2344 1.195 −1.22 3.53 −1.18 3.49
653.1421 1.218 −0.85 3.48 −0.82 3.45
654.3504 1.195 −1.71 3.51 −1.66 3.46
656.5883 1.183 −2.05 3.51 −2.05 3.51
660.5974 1.195 −1.34 3.51 −1.31 3.48
660.7828 1.350 −1.98 3.63 ... ...
662.4845 1.218 −1.30 3.48 −1.26 3.44
675.3010 1.081 −1.50 3.53 ... ...
678.5000 1.051 −1.88 3.55 ... ...
681.2400 1.043 −2.14 3.55 ... ...
683.2442 1.064 −2.43 3.53 ... ...
684.1877 1.051 −2.53 3.56 ... ...
733.8931 2.138 −0.74 3.56 ... ...
736.3154 2.122 −1.01 3.53 ... ...
802.7366 1.064 −1.85 3.55 −1.86 3.56
809.3468 1.051 −1.76 3.60 −1.76 3.60
811.6789 1.081 −1.07 3.57 −1.03 3.53
814.4560 1.043 −1.90 3.56 −1.87 3.53
824.1599 1.051 −1.90 3.58 −1.87 3.55
825.3506 1.081 −1.85 3.58 −1.84 3.57
825.5896 1.064 −1.75 3.58 −1.70 3.53
891.9847 1.218 −1.49 3.58 −1.33 3.42
893.2947 1.195 −1.74 3.51 −1.58 3.35
897.1673 1.183 −2.13 3.57 −1.98 3.42
903.7613 1.183 −2.12 3.55 −1.94 3.37
904.6693 1.195 −2.02 3.53 −1.79 3.30
908.5231 1.218 −2.08 3.55 −1.91 3.38

Download table as:  ASCIITypeset image

In Figure 6 we plot the line abundances as a function of transition strength. As in previous papers of this (e.g., La14)series we define relative strengths of lines of a single species as products of their transition probabilities and Boltzmann factors: log(gf)–θχ, where θ = 5040/T. In our previous papers we have assumed θ = 1.0 as a compromise between the solar value θ ≈ 0.9 and those of typical very metal-poor giants θ ≈ 1.1. Arcturus is much cooler, θ = 1.18, and we have used this value for Figure 6. Inspection of Figure 6 indicates that there are no secular abundance trends with line strength, and the 900 nm lines span nearly 1 dex in transition strength. In Figure 7 we plot the V abundances as a function of lower level energy (excitation potential), and as with the previous figures there are no significant abundance trends. We note that all of the 900 nm lines arise from closely spaced lower energy levels near 1.2 eV, and these levels give rise to lines at shorter wavelengths that we have included in our abundance computations. Taken together, these figures show good agreement in derived V abundance for lines below λ < 826 nm and the lines near 900 nm that were questioned by Ho16. That there are no obvious trends with wavelength, transition strength, or excitation potential further supports our assertion that the data from La14 are consistent in both wavelength regions.

Figure 6.

Figure 6. Abundances from V i lines derived using La14 transition data, plotted as a function of transition strength. The symbols and lines are the same as in Figure 5. See the text for further explanation of the transition strength definition adopted here.

Standard image High-resolution image
Figure 7.

Figure 7. Abundances from V i lines derived using La14 transition data, plotted as a function of lower level energy (excitation potential). The symbols and lines are the same as in Figure 5.

Standard image High-resolution image

Ramírez & Allende Prieto (2011) reported V abundances in Arcturus through an EW analysis of eight weak V i lines, log(RW) = −6.0 to −5.3, treating them as single absorbers. They reported [V/Fe] = +0.20 ± 0.05, which is consistent with our value given that they used different V i line parameters and a different model atmosphere than we employed. However, their preferred vt = 1.74 km ${{\rm{s}}}^{-1}$ leads to more abundance inconsistencies between weak and strong V i lines than we derived with the Peterson et al. value of vt = 1.6 km ${{\rm{s}}}^{-1}$. Further investigation of microturbulence in the Arcturus atmosphere is beyond the scope of our work.

4. Neutral Vanadium Transitions in the Infrared

The longest-wavelength V i line reported by La14 is at 915.66 nm, but Hinkle et al. (1995) listed 18 identifications of this species in their first Arcturus infrared atlas. Ho16 reported new log(gf) values for nine IR lines, but none of these appeared in the Hinkle et al. list. We searched without success for these lines in the Arcturus Atlas. Nearly all of them lie in the thicket of telluric absorption that defines the division between the photometric H and K bands. Only one of the Ho16 lines, 1798.75 nm, is not overly contaminated by neighboring telluric or stellar features. However, the V i absorption is less than 1% deep at most in the observed Arcturus Atlas spectrum and thus can yield no reliable abundance value.

We explored this wavelength region with V i transitions listed in the Kurucz (2011) line compendium. We were able to detect many of the lines suggested by Hinkle et al. (1995), but the observed lines are weak. Having no lab data to guide this work, we adopted the Kurucz log(gf) values and made synthetic/observed spectrum matches of all clean V i lines that we identified. This exercise yielded V abundances in qualitative agreement with those determined from optical-region lines, but with large (±0.15) line-to-line scatter. This spectral region will be studied in the future with support from the new V i lab spectra acquired as part of this study.

5. Vanadium Trends with Metallicity

One of the important goals of our studies is to understand the changing nature of detailed stellar abundance ratios throughout the history of the Galaxy, i.e., Galactic chemical evolution (GCE). As detailed in Section 3, the present study indicates that the vanadium abundance in Arcturus is [V/H] = −0.42, or [V/Fe] = +0.13.

In Figure 8 we compare this newly determined abundance with other abundance surveys of Galactic stars including Feltzing & Gustafsson (1998), Barklem et al. (2005), Roederer et al. (2014), and Battistini & Bensby (2015); see also Reddy et al. (2006) for additional observations. These compilations include both halo and disk stars and extend over a more metal-rich region including super-solar metallicity, i.e., [Fe/H] > 0. At metallicities near [Fe/H] = −0.5 there is a small overabundance seen in [V/Fe], and our new value for Arcturus is consistent with that trend. Previous studies have also indicated modest [V/Fe] overabundances at lower metallicities typical of halo stars, as illustrated in Figure 8 (see also Sneden et al. 2016).

Figure 8.

Figure 8. Newly determined V abundance in Arcturus compared to abundance values of other Galactic disk and halo stars. (See the text for details.)

Standard image High-resolution image

GCE models have attempted to match these abundance observations as a function of metallicity. In Figure 8 we illustrate two such models, one by Kobayashi et al. (2011) and a Kobayashi model with jet effects (Sneden et al. 2016). It is clear that these models mimic the general trends but are not yet a good fit to the observed V abundances. In particular, it is seen that the deviation between GCE models and the abundance data grows larger at higher metallicities. We note, however, that (arbitrarily) adding 0.3 [V/Fe] to the predictions of the GCE model with jet effects would result in a good fit to almost all of the data, including the new abundance ratio for Arcturus.

Some other iron-group elements show a similar lack of fit with current GCE models (see Sneden et al. 2016). At metallicities typical of Arcturus the iron-group elements are produced by both core-collapse supernovae (SNe) and Type Ia SNe. Therefore, new measured values for iron-group abundances like that of V at these metallicities help to constrain the frequency and models of both SNe. Furthermore, additional high-resolution stellar observations spanning a wide metallicity range combined with abundances derived using precise laboratory transition data will provide meaningful inputs for improving GCE models and help our understanding of the time evolution of the elemental abundances in the Galaxy.

6. Summary

The increasing availability of IR astrophysical spectra generates a significant need for accurate and reliable IR laboratory atomic transition data. A recent study of V i oscillator strengths by Ho16 found generally good agreement with an earlier study (La14). However, Ho16 noted disagreement with a small subset of lines in the near-IR around 900 nm, suggesting that the discordance indicated a calibration error by La14. In the present study, several new laboratory FTS spectra of V HC lamps were recorded at NIST to assess this discordance and compare the two calibration methods employed in the earlier studies. After analyzing these new data, we find no evidence for a problem with the La14 results or the Ar BR calibration technique near 900 nm. Furthermore, we derived a new V abundance in the very bright red giant star Arcturus as a secondary test of the reliability of the La14 log(gf) values. That there is good abundance agreement among all V i lines utilized and that no obvious trends are seen in the abundances further bolsters our confidence in the earlier La14 results. Our new Arcturus [V/Fe] ratio is in good agreement with those of other mildly metal-poor stars. Finally, a thorough survey of the literature has resulted in a greatly expanded set of hyperfine structure component patterns for V i lines that may be of interest to future stellar abundance studies.

The authors thank the referee for the many helpful and detailed suggestions that improved the quality of the paper. Our lab/stellar research has been supported by grants from NASA and NSF, most recently by NASA grant NNX16AE96G (J.E.L.), by NSF grant AST-1516182 (J.E.L., E.A.D.H.) and NSF grant AST-1211585 and AST-161040 (C.S.). M.P.W. and J.E.L. served as guest researchers at NIST while recording the vanadium spectra analyzed in the present study. The earlier titanium spectrum was recorded while M.P.W. served as an NIST NRC Postdoctoral Fellow. J.J.C. acknowledges support by the National Science Foundation under grant No. PHY-1430152 (JINA Center for the Evolution of the Elements).

Appendix:

There has been substantial new work on the HFS of lines of V i since La14. Complete HFS line component patterns for a total of 94 lines of V i are provided in La14 based on both FTS measurements of lines connecting to 25 upper levels and published HFS data from Childs & Goodman (1967), Gough et al. (1985), Cochrane et al. (1998), Palmeri et al. (1995, 1997), Lefèbvre et al. (2002), and Güzelçimen et al. (2011). The 94 lines selected in La14 were lines found to be useful in abundance determinations for the Sun and HD 84937. An article by Güzelçimen et al. (2014) appeared while La14 was in the galley proof stage that confirmed all but one of the FTS measurements in La14 to a fraction of an error bar. There is substantial variation in the accuracy and precision of HFS measurements depending on the technique used. The radio frequency technique as used by Childs & Goodman (1967), applicable to the ground and low metastable levels, yields results to better than 1 kHz. The Doppler-free laser technique as used by Gough et al. (1985) yields HFS constants to better than 1 MHz (30 MHz = 0.001 cm−1). Nonlinear least-square fitting to FTS line profiles as used by La14 yields results accurate to around 0.001-0.002 cm−1. For elemental abundance studies on the Sun and other stars, HFS component patterns to FTS accuracy are satisfactory. The additional publications employed herein to generate complete HFS line component patterns are Güzelçimen et al. (2014, 2015), Unkel et al. (1989), and Ho16 for the ${x}^{6}{{\rm{D}}}_{3/2}$ level at 28368.753 cm−1. The Casimir formula is from the text by Woodgate (1980):

where ΔE is the shift in energy (or wavenumber) of a HFS sub-level (F, J) from the center of gravity of the fine structure level (J),

F is the total atomic angular momentum, J is the total electronic angular momentum, and I = 3/2 is the nuclear spin for 51V. Although Childs & Goodman report HFS C (magnetic octupole) coefficients, their error bars overlap zero and they are not used here. Non-zero B (electric quadrupole) coefficients are included for some levels in this work, but only reliable HFS A (magnetic dipole) coefficients are necessary for abundance studies on stellar spectra.

Rather than add HFS component patterns for a small number of lines to the 94 published in La14, we decided to survey all 836 lines of V i with log(gf)s in La14 and generate complete HFS component patterns for every line with a known upper and lower level HFS A coefficient. In order to provide some added value in this exercise we compared the generated component patterns to our best FTS line profiles. Although nearly all of the hundreds of published HFS A and B coefficients were found to be reliable to FTS standards, several problems were identified and corrected in the published results. The HFS A = 73 MHz coefficient of the 2G9/2 level at 37361.95 cm−1 reported by Güzelçimen et al. (2014) is not reliable and should be 0.0100 ± 0.001 cm−1 (≈300 MHz) from our nonlinear least-square fits to FTS line profiles. Similarly, their HFS A = 118 MHz coefficient of the 4I11/2 level at 37315.93 cm−1 is not reliable and should be 0.0103 (≈309 MHz) for astrophysical studies. This upper 4I11/2 level lacks lines in our FTS data with resolved structure and we are not able to set a secure uncertainty on the new HFS A coefficient. Finally, a line from the ${x}^{2}{{\rm{F}}}_{3/2}$ level at 36766.041 cm−1 is useful in Arcturus and has some partially resolved HFS in our FTS data. Nonlinear least-square fits to FTS line profiles yield 0.0100 ± 0.001 cm−1 for the HFS A of the ${x}^{2}{{\rm{F}}}_{3/2}$ level, and this value is used with published HFS coefficients in the generation of line component patterns.

The machine-readable version of Table 3 has complete HFS line component patterns for more than 650 of the 839 lines with log(gf)s in La14. Although only a few new HFS A coefficients were used to generate the table, the complete patterns are convenient in astrophysical studies. This table represents a sevenfold expansion of the La14 table of HFS line component patterns and includes the majority of the lines in La14 with reliable log(gf) values.

Table 3.  HFS Line Component Patterns for 658 Transitions of 51V i

Transition Wavenumbera Wavelength in Aira Fupp Flow Component Position Component Position Normalized Strengthb
(cm−1) (nm)     (cm−1) (pm)  
41620.972 240.1904 7 6 0.01867 −0.1078 0.23437
41620.972 240.1904 6 6 −0.04559 0.2631 0.03385
41620.972 240.1904 6 5 0.01876 −0.1083 0.16927
41620.972 240.1904 5 6 −0.10067 0.5810 0.00260
41620.972 240.1904 5 5 −0.03632 0.2096 0.05320

Notes. The first five entries are shown to guide the reader as to the form and content of the table. As a result, the normalized strengths of the components shown in this table stub do not sum to unity. These component patterns are computed from the best published HFS coefficients (see the text) and energy levels from Thorne et al. (2011).

aCenter-of-gravity wavenumbers and wavelengths are given, with component positions relative to those values. Air wavelengths are computed using the standard index of air (Peck & Reeder 1972). bComponent strengths are normalized to 1.0 for each transition.

Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.

Download table as:  DataTypeset image

Footnotes

  • We use standard abundance notations. For elements X and Y, the relative abundances are written as $[{\rm{X}}/{\rm{Y}}]={\mathrm{log}}_{10}{({{\rm{N}}}_{{\rm{X}}}/{{\rm{N}}}_{{\rm{Y}}})}_{\mathrm{star}}-{\mathrm{log}}_{10}{({{\rm{N}}}_{{\rm{X}}}/{{\rm{N}}}_{{\rm{Y}}})}_{\odot }$. For element X, the "absolute" abundance is written as $\mathrm{log}\,\varepsilon ({\rm{X}})\,={\mathrm{log}}_{10}({{\rm{N}}}_{{\rm{X}}}/{{\rm{N}}}_{{\rm{H}}})+12$. Metallicity is defined as [Fe/H]. The unit dex stands for decimal exponent, such that x dex = 10x.

  • Any commercial equipment, instruments, or materials are identified in this paper in order to specify the experimental procedure adequately. Such identification is not intended to imply recommendation or endorsement by NIST, nor is it intended to imply that the materials or equipment identified are necessarily the best available for the purpose.

  • All uncertainties in this paper are one standard deviation.

  • 10 
  • 11 
  • 12 
  • 13 
Please wait… references are loading.
10.3847/1538-4365/aa9a41