Category Archives: Publications

Searching for magnetic fields in pulsating A-type stars

Numerous δ Sct and γ Dor pulsators are identified in the region of the Hertzsprung-Russell diagram that is occupied by chemically peculiar magnetic Ap stars. The connection between δ Sct and γ Dor pulsations and the magnetic field in Ap stars is however not clear: theory suggests for magnetic Ap stars some critical field strengths for pulsation mode suppression by computing the magnetic damping effect for selected p and g modes. To test these theoretical considerations, we obtained PEPSI spectropolarimetric snapshots of the typical Ap star HD 340577, for which δ Sct- like pulsations were recently detected in TESS data, and the γ Dor pulsator HR 8799, which is a remarkable system with multiple planets and two debris disks. Our measurements reveal the presence of a magnetic field with a strength of several hundred Gauss in HD 340577. The measured mean longitudinal field would be the strongest field measured so far in a δ Sct star if the pulsational character of HD 340577 is confirmed spectroscopically. No magnetic field is detected in HR 8799.

LSD Stokes I, V, and diagnostic null N spectra (from bottom to top) calculated for HD340577, for the mask with all lines and the masks for individual ions. The numbers in brackets relate to the number of lines used in each individual mask. The grey horizontal lines in the Stokes V and N spectra indicate the ±1σ error bars. Stokes V spectra are highlighted in red.

Read more: Hubrig et al. 2023, MNRAS 526, L83; arXiv 2308.09441

Zeeman Doppler Imaging of ksi Boo A and B

A magnetic-field surface map for both stellar components of the young visual binary ksi Boo AB (A: G8V, B: K5V) is presented. Employed are high resolution Stokes-V spectra obtained with the Potsdam Echelle Polarimetric and Spectroscopic Instrument (PEPSI) at the Large Binocular Telescope (LBT). Stokes V line profiles are inverted with our iMAP software and compared to previous inversions. We employed an iterative regularization scheme without the need of a penalty function and incorporated a three-component description of the surface magnetic-field vector. The spectral resolution of our data is 130,000 (0.040-0.055A) and have signal-to-noise ratios (S/N) of up to three thousand per pixel depending on wavelength. A singular-value decomposition (SVD) of a total of 1811 spectral lines is employed for averaging Stokes-V profiles. Our mapping is accompanied by a residual bootstrap error analysis. Magnetic flux densities of the radial field component of up to plus/minus 115 ± 5 G were reconstructed for ksi Boo A while up to plus/minus 55 ± 3G were reconstructed for ksi Boo B. ksi Boo A’s magnetic morphology is characterized by a very high latitude, nearly polar, spot of negative polarity and three low-to-mid latitude spots of positive polarity while ksi Boo B’s morphology is characterized by four low-to-mid latitude spots of mixed polarity. No polar magnetic field is reconstructed for the cooler ksi Boo B star. Both our maps are dominated by the radial field component, containing 86 and 89 percent of the magnetic energy of ksi Boo A and B, respectively. We found only weak azimuthal and meridional field densities on both stars (plus/minus 15-30 G), about a factor two weaker than what was seen previously for ksi Boo A. The phase averaged longitudinal field component and dispersion is +4.5 ± 1.5G for ksi Boo A and -5.0 ± 3.0 G for ksi Boo B.

Representative PEPSI spectra used in the study (panel a). Zeeman-Doppler images of ξ Boo A (panel b) and ξ Boo B (panel c) in orthographic projection. φ is the rotational phase. Magnetic field strength is color coded and identified in the right bars in units of Gauss. Positive polarity is depicted in red, negative polarity in blue. The length of the surface dashes is proportional to field strength.

Read more: K. G. Strassmeier, T. A. Carroll, & I. V. Ilyin 2023, A&A, 674, 118; (arXiv 2305.07470)

Constraints on Magnetic Braking from the G8 Dwarf Stars 61 UMa and Tau Cet

During the first half of their main-sequence lifetimes, stars rapidly lose angular momentum to their magnetized winds, a process known as magnetic braking. Recent observations suggest a substantial decrease in the magnetic braking efficiency when stars reach a critical value of the Rossby number, the stellar rotation period normalized by the convective overturn timescale. Cooler stars have deeper convection zones with longer overturn times, reaching this critical Rossby number at slower rotation rates. The nature and timing of the transition to weakened magnetic braking has previously been constrained by several solar analogs and two slightly hotter stars. In this Letter, we derive the first direct constraints from stars cooler than the Sun. We present new spectropolarimetry of the old G8 dwarf τ Cet from the Large Binocular Telescope, and we reanalyze a published Zeeman Doppler image of the younger G8 star 61 UMa, yielding the large-scale magnetic field strengths and morphologies. We estimate mass-loss rates using archival X-ray observations and inferences from Lyα measurements, and we adopt other stellar properties from asteroseismology and spectral energy distribution fitting. The resulting calculations of the wind braking torque demonstrate that the rate of angular momentum loss drops by a factor of 300 between the ages of these two stars (1.4–9 Gyr), well above theoretical expectations. We summarize the available data to help constrain the value of the critical Rossby number, and we identify a new signature of the long-period detection edge in recent measurements from the Kepler mission.

Stokes V polarization profile for τ Cet from LBT observations on 2022 September 18. The mean profile is shown as a black line with uncertainties indicated by the gray shaded area. The dashed blue line is an axisymmetric model profile assuming dipole morphology with the inclination fixed at i = 20 deg.

Read more: Metcalfe, T., Strassmeier, K., Ilyin, I., et al. 2023, ApJ Letters, in press (arXiv:2304.09896)

The PEPSI Exoplanet Transit Survey. III: The detection of Fe I, Cr I and Ti I in the atmosphere of MASCARA-1 b through high-resolution emission spectroscopy

Hot giant planets like MASCARA-1 b are expected to have thermally inverted atmospheres, that makes them perfect laboratory for the atmospheric characterization through high-resolution spectroscopy. Nonetheless, previous attempts of detecting the atmosphere of MASCARA-1 b in transmission have led to negative results. In this paper we aim at the detection of the optical emission spectrum of MASCARA-1 b. We used the high-resolution spectrograph PEPSI to observe MASCARA-1 (spectral type A8) near the secondary eclipse of the planet. We cross-correlated the spectra with synthetic templates computed for several atomic and molecular species. We obtained the detection of FeI, Cr I and Ti I in the atmosphere of MASCARA-1 b with a S/N ~7, 4 and 5 respectively, and confirmed the expected systemic velocity of ~13 km/s and the radial velocity semi-amplitude of MASCARA-1 b of ~200 km/s. The detection of Ti is of particular importance in the context of the recently proposed Ti cold-trapping below a certain planetary equilibrium temperature. We confirm the presence of an the atmosphere around MASCARA-1 b through emission spectroscopy. We conclude that the atmospheric non detection in transmission spectroscopy is due to the high gravity of the planet and/or to the overlap between the planetary track and its Doppler shadow.

Atmospheric planetary signal obtained with the Ti I template. The left column shows the Kp–vsys bi-dimensional map, with the color bar encoding the significance of the detection. The cross in white dashes marks the position where the planetary signal is expected. The central (right) column shows the cut along the vertical (horizontal) direction in the map and corresponding to the peak detection, together with its gaussian best fit (blue dashed line).

Read more: G. Scandariato, F. Borsa, A.S. Bonomo, et al. 2023, A&A, in press (eprint arXiv:2304.03328)

VPNEP: Detailed characterization of TESS targets around the Northern Ecliptic Pole

We embarked on a high-resolution optical spectroscopic survey of bright Transiting Exoplanet Survey Satellite (TESS) stars around the Northern Ecliptic Pole (NEP), dubbed the Vatican-Potsdam-NEP (VPNEP) survey. Our NEP coverage comprises 1067 stars, of which 352 are bona fide dwarf stars and 715 are giant stars, all cooler than spectral type F0 and brighter than V=8.5. Our aim is to characterize these stars for the benefit of future studies in the community. We analyzed the spectra via comparisons with synthetic spectra. Particular line profiles were analyzed by means of eigen-profiles, equivalent widths, and relative emission-line fluxes (when applicable). Two R=200 000 spectra were obtained for each of the dwarf stars with the Vatican Advanced Technology Telescope (VATT) and the Potsdam Echelle Polarimetric and Spectroscopic Instrument (PEPSI), with typically three R=55 000 spectra obtained for the giant stars with STELLA and the STELLA Echelle Spectrograph (SES). Combined with V-band magnitudes, Gaia eDR3 parallaxes, and isochrones from the Padova and Trieste Stellar Evolutionary Code, the spectra can be used to obtain radial velocities, effective temperatures, gravities, rotational and turbulence broadenings, stellar masses and ages, and abundances for 27 chemical elements, as well as isotope ratios for lithium and carbon, line bisector spans, convective blue-shifts (when feasible), and levels of magnetic activity from Halpha, Hbeta, and the Ca II infrared triplet. In this initial paper, we discuss our analysis tools and biases, presenting our first results from a pilot sub-sample of 54 stars (27 bona-fide dwarf stars observed with VATT+PEPSI and 27 bona-fide giant stars observed with STELLA+SES) and making all reduced spectra available to the community.

Data delivery example: 8000-Å CN region of HR 5844. Diamond symbols are the observed spectrum and the lines are fits with three different 12C/13C isotope ratios, using the CN line list of Plez. The best average from eight fits is 12C/13C=17.8 using a relative nitrogen abundance of 0.50±0.05 and a relative iron abundance of –0.25±0.03. Average S/N per pixel of this PEPSI spectrum in this wavelength range is 400.

Read more: Strassmeier et al. 2023, A&A, 671, A7 (arXiv:2302.01794)

Press release: To new worlds with quantitative spectroscopy

The PEPSI Exoplanet Transit Survey (PETS). II. A Deep Search for Thermal Inversion Agents in KELT-20 b/MASCARA-2 b with Emission and Transmission Spectroscopy

Recent observations have shown that the atmospheres of ultra hot Jupiters (UHJs) commonly possess temperature inversions, where the temperature increases with increasing altitude. Nonetheless, which opacity sources are responsible for the presence of these inversions remains largely observationally unconstrained. We used LBT/PEPSI to observe the atmosphere of the UHJ KELT-20 b in both transmission and emission in order to search for molecular agents which could be responsible for the temperature inversion. We validate our methodology by confirming previous detections of Fe I in emission at 16.9 σ. Our search for the inversion agents TiO, VO, FeH, and CaH results in nondetections. Using injection-recovery testing we set 4σ upper limits upon the volume mixing ratios for these constituents as low as ∼1×10−9 for TiO. For TiO, VO, and CaH, our limits are much lower than expectations from an equilibrium chemical model, while we cannot set constraining limits on FeH with our data. We thus rule out TiO and CaH as the source of the temperature inversion in KELT-20 b, and VO only if the line lists are sufficiently accurate.

Shifted and combined CCFs for emission from Fe I. We detect the Fe I emission at a significance of 16.9σ, reproducing previous results. We were able to recover the signal near the expected planetary radial velocity. Yan et al. (2022b) obtained a 7.7σ detection, and Borsa et al. (2022) at 7.1σ with the Calar Alto 3.5m and the TNG 3.6m, respectively.

Read more: Marshall et al. 2023, AJ, in press (arXiv:2205.12162)

Untangling the Sources of Abundance Dispersion in Low-Metallicity Stars

We measured abundances of 12 elements (Na, Mg, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni) in a sample of 86 metal-poor subgiant stars in the solar neighborhood. Abundances are derived from high-resolution spectra taken with the Potsdam Echelle Polarimetric and Spectroscopic Instrument on the Large Binocular Telescope, modeled using iSpec and MOOG. By carefully quantifying the impact of photon-noise (< 0.05 dex for all elements) we robustly measure the intrinsic scatter of abundance ratios. At fixed [Fe/H] the RMS intrinsic scatter in [X/Fe] ranges from 0.04 dex (Cr) to 0.16 dex (Na), with a median of 0.08 dex. Scatter in [X/Mg] is similar, and accounting for [alpha/Fe] only reduces the overall scatter moderately. We consider several possible origins of the intrinsic scatter with particular attention to fluctuations in the relative enrichment by core-collapse supernovae (CCSN) and Type Ia supernovae (SNIa) and stochastic sampling of the CCSN progenitor mass distribution. The stochastic sampling scenario provides a good quantitative explanation of our data if the effective number of CCSN contributing to the enrichment of a typical sample star is N approx 50. At the median metallicity of our sample, this interpretation implies that the CCSN ejecta are mixed over a gas mass 10^5 MSun before forming stars. The scatter of elemental abundance ratios is a powerful diagnostic test for simulations of star formation, feedback, and gas mixing in the early phases of the Galaxy.

Spectral section from 5495−5535 A of three stars with [Fe/H] near −2 (green), −1.5 (blue), and −1.0 (purple). Line features for Fe I, Ti I, Ca II, Sc II, and Mg I are labeled. Fluxes are normalized, with the [Fe/H] ≈ −1.5 spectrum offset by 0.5 and the [Fe/H] ≈ −2 spectrum offset by 1.

Read more:  Griffith et al. 2022, ApJ, in press  (arXiv:2210.01821)

Io’s Optical Aurorae in Jupiter’s Shadow

Recent studies have established that the majority of Io’s molecular atmosphere, SO2 and SO, condenses during its passage through Jupiter’s shadow. The eclipse response of Io’s atomic atmosphere is less certain, having been characterized solely by ultraviolet aurorae. Here we explore the response of optical aurorae for the first time. We find oxygen to be indifferent to the changing illumination with [O I] brightness merely tracking the plasma density at Io’s position in the torus. In shadow, line ratios confirm sparse SO2 coverage relative to O since their collisions would otherwise quench the emission. Io’s sodium aurora mostly disappears in eclipse and e-folding timescales for decline and recovery differ sharply: ~10 minutes at ingress and nearly 2 hours at egress. Only ion chemistry can produce such a disparity; Io’s molecular ionosphere is much weaker at egress due to rapid recombination. Auroral emission is also evident from potassium, confirming K as the major source of far red emissions seen in situ. In all cases, direct electron impact on atomic gas is sufficient to explain the brightness without invoking significant dissociative excitation of molecules. The non-response of O and rapid depletion of Na during Io’s eclipse phase is surprisingly inverted from the eclipse phase behavior of the SO2 and NaCl parent molecules.

Io’s [O I] 6300 Å line emission amidst Jupiter scattered light on 24 April 2019 with LBT/PEPSI. The red-to-blue color scheme shows elapsed time in shadow. Scattered light levels increase as Io approaches the jovian limb and dotted lines show the fitting of a background spectrum to this component. Lower panel shows the residual after scattered light subtraction. Gaussians are fit to Io’s emission line, which Doppler shifts slightly in time. The observing geometry for the LBT/PEPSI observations is shown on the left panel.

Read more: Schmidt et al. 2022, PSJ, in press (AAS Planetary Science Journal)

On the lithium abundance of the visual binary components ksi Boo A (G8V) and ksi Boo B (K5V)

A spectroscopic investigation of the lithium resonance doublet in ξ Boo A and ξ Boo B in terms of both abundance and isotopic ratio is presented. We obtained new R=130 000 spectra with a signal-to-noise ratio (S/N) per pixel of up to 3200 using the 11.8m LBT and PEPSI. From fits with synthetic line profiles based on 1D-LTE MARCS model atmospheres and 3D-NLTE corrections, we determine the abundances of both isotopes. For ξ Boo A, we find A(Li) = 2.40±0.03 dex and 6Li/7Li < 1.5±1.0 % in 1D LTE, which increases to ≈2.45 for the 3D-NLTE case. For ξ Boo B we obtain A(Li) = 0.37±0.09 dex in 1D-LTE with an unspecified 6Li/7Li level. Therefore, no 6Li is seen on any of the two stars. We consider a spot model for the Li fit for ξ Boo B and find A(Li) = 0.45±0.09 dex. The 7Li abundance is 23 times higher for ξ Boo A than the Sun’s, but three times lower than the Sun’s for ξ Boo B while both fit the trend of single stars in the similar-aged M35 open cluster. Effective temperatures are redetermined from the TiO band head strength. We note that the best-fit global metallicities are –0.13±0.01 dex for ξ Boo A but +0.13±0.02 dex for ξ Boo B. Lithium abundance for the K5V benchmark star 61 Cyg A was obtained to A(Li)≈0.53 dex when including a spot model but to ≈0.15 dex without a spot model.

Lithium 6708 Å of ξ Boo A (black dots) and ξ Boo B (red dots). The wavelengths of the two doublets from the two lithium isotopes are marked as vertical lines. Also indicated are the blending features from the line list of Meléndez et al. (2012). Additional molecular features from the sunspot umbral spectrum of Wallace et al. (1999) are indicated as well.

Read more: Strassmeier & Steffen 2022, AN, 343, article id. e20220036

The recurrent nova V3890 Sgr: a near-infrared and optical study of the red giant component and its environment

We present an analysis of the red giant component of the recurrent nova V3890 Sgr, using data obtained before and after its 2019 eruption. Its effective temperature is Teff = 3050 ± 200 K for log g = 0.7, although there are modest changes in Teff. There is an overabundance of both carbon (0.20 ± 0.05 dex) and sodium (1.0 ± 0.3 dex) relative to their solar values, possibly the result of ejecta from the 1990 nova eruption being entrained into the red giant photosphere. We find 12C/13C =25 ± 2, a value similar to that found in red giants in other recurrent novae. The spectrum in the region of the Na I D lines is complex, and includes at least six interstellar components, together with likely evidence for interaction between ejecta from the 2019 eruption and material accumulated in the plane of the binary.

PEPSI (R = 130, 000) observed velocity profiles of various strong emission lines in V3890 Sgr post the 2019 RN eruption at a phase of 0.709. Top: Hα and Hβ. Middle: Permitted lines of He I 5875Å, 6678Å, and the forbidden O III 5006Å line. Bottom: The region of the Na I D lines, for which the velocities are calculated with respect to the D2 line (5889.995Å air). The two narrow features marked by the blue circles are night sky sodium lines. Velocity components of various absorption features are indicated. Note that there are only three distinct Na I D absorption systems: the features at −110.21 km/s and −190.47 km/s are from the same systems as those at +112 km/s and +193 km/s, but for D1.

Read more: Kaminsky et al. 2022, MNRAS, 517, 6064