Category Archives: Publications

Does the solar oxygen abundance change over the solar cycle? An investigation into activity-induced variations of the O I infrared triplet

The determination of the solar oxygen abundance remains a central problem in astrophysics, as its accuracy is limited not only by models but also by systematics. While many of these factors have been thoroughly characterized, the effect of the solar activity cycle has so far remained unexplored. Due to its relative strength and accessibility, the O I infrared triplet is typically the primary choice for abundance studies. However, previous investigations have shown that abundances inferred from this triplet tend to be higher than expected on active stars, whereas such an overabundance effect is not observed for the much weaker forbidden O I 6300 Å line. This raises the question of whether a similar trend can be found for the Sun. To address this question, we analyze two decades’ worth of synoptic disk-integrated Sun-as-a-star datasets from the FEROS, HARPS-N, PEPSI, and NEID spectrographs, focusing on the infrared triplet (7772, 7774, 7775 Å) and the forbidden O I 6300 Å line. The excellent signal-to-noise ratio of the PEPSI observations allows us to detect a weak but significant variation in the equivalent widths of the infrared triplet, corresponding to about 0.01 dex difference in abundance between activity minimum and maximum. This value is significantly smaller than the typical uncertainties on the solar oxygen abundance. Due to higher scatter, no comparable trend is found in the other data sets. Based on these results, we conclude that within the typical uncertainties presented in other works, we can assume the inferred solar oxygen abundance to be stable across the solar cycle, but that this effect may be significant for other, more active stars.

Measured EW of the O I 7772 Å line with FEROS (blue), PEPSI (orange), and NEID (green), including corresponding uncertainties. The gray curve shows the sunspot number as an activity reference, while the black horizontal line at 50 distinguishes active from quiet times.

Read more: Pietrow et al. 2026, A&A, in press. arXiv:2601.15919

Surface image and activity-corrected orbit of the RS CVn binary HR7275. Disentangling activity tracers.

Quantifying stellar parameters and magnetic activity for cool stars in double-lined spectroscopic binaries (SB2) is not straightforward, as both stars contribute to the observed composite spectra and are likely variable. Disentangled component spectra allow a detailed analysis of a component’s magnetic activity. We aim at separating the spectra of the two stellar components of the HR7275 SB2 system. Our further aim is a more accurate orbital solution by cleaning the observed radial velocities (RV) from activity perturbations of the spotted primary (“RV jitter”) and obtain a surface image of this component. We provide time-series high- and ultra-high resolution optical spectra and apply two different disentangling methods. RV residuals are modeled with three-sine function fits. The primary’s spectral-line profiles are modeled with the Doppler imaging code iMAP. Magnetic fields are measured for the primary based on least-square deconvolved Stokes-V line profiles. Chromospheric emission is determined from the line-cores of Ca II H&K, Ca II IRT 8542 Å, and Balmer Hα. Before applying those analyses, we provide a disentangling technique to determine the system properties more accurately. The Doppler image of the primary shows two large cool spots of size ≈20% of the visible hemisphere plus three smaller spots, each still ≈13% in size. In total, HR7275a exhibited an impressive spottedness of ≈40% of its entire surface in May-June 2022. The RV is modulated by the rotation of the primary with maximum amplitudes of 320 ms−1 and 650 ms−1 for two different modulation behaviors during the 250 d of our observations. This jitter is primarily caused by the varying asymmetries of the apparent disk brightness due to the cool spots. Its removal resulted in roughly ten times higher precision of the orbital elements. Our snapshot magnetic-field measurements reveal phase dependent (large-scale) surface fields between +0.6±2.0G at phase 0.1 and −15.2±2.7G at phase 0.6, indicating a complex magnetic morphology related to the location of the photospheric spots. We also obtain a logarithmic lithium abundance of 0.58±0.1 for HR7275a, indicating considerable mixing, and 0.16 for HR7275b, which is an extremely low value.

 

Spectral disentangling by median subtraction. Panel a: The mean spectrum of the secondary star is plotted as a black line (top). The composite spectrum (blue) and a spectrum of the primary star (red) are plotted with a vertical shift of 0.15 at orbital phase 0.82. Panel b: Time series composite spectra phase folded with the orbital period.

Read more: Adebali, Weber, Strassmeier, et al. 2026, A&A, in press. arXiv:2512.09521

Untangling the Sources of Abundance Dispersion in Low-metallicity Stars II:Neutron Capture Elements

We present the abundances of 23 elements, including 11 heavy elements (Cu, Zn, Sr, Y, Zr, Ba, La, Ce, Nd, Sm, Eu) for up to 86 metal-poor (−2 ≲ [Fe/H] ≲ −1) subgiants. We use KORG, a state of the art spectral synthesis package, to derive 1D-LTE abundances from high-SNR and high-resolution spectra taken by the Large Binocular Telescope with the Potsdam Echelle Polarimetric and Spectroscopic Instrument. These precise spectra and abundance measurements minimize the impact of photon-noise (≲ 0.06 dex), allowing us to robustly measure the intrinsic abundance scatter in [X/Fe]. After removing two stars with exceptional s-process enhancement, we find that the intrinsic scatter among the s- and r-process elements tends to be larger than for the lighter elements, with heavy element scatter ranging from 0.11 (Zn) to 0.27 (Eu) dex. Intrinsic abundance scatter could have multiple origins, including starto-star variations in the ratios of nucleosynthetic sources as well as stochastic sampling of the progenitor supernovae properties, such as mass, rotation, and magnetic field strength. We explore the expected abundance scatter signature caused by stochastic sampling, finding that a fraction of both rapidly rotating CCSN and magnetorotationally driven SN are needed to reach the observed abundances and intrinsic scatter. This analysis is limited by the restrictive parameter spaces spanned by existing yield sets. A diverse, finely sampled grid of supernovae yields is needed to robustly model stochastic abundance scatter.

Heavy element line windows for three stars with [Fe/H] of −2.12 (green, 2MASS
J04315411-0632100), −1.58 (blue, 2MASS J15581861+0203059), and −1.04 (purple, 2MASS J17140534+1407170). Spectra have been offset by 0.1 for clarity. We show one window for each heavy element, with line centers labeled, including blended lines for Zr and Eu.

Read more: Griffith, E. J. et al. 2026, ApJ, in press. arXiv:2512.02122

Elemental abundance pattern and temperature inversion on thedayside of HAT-P-70b observed with CARMENES and PEPSI

We observed the dayside thermal emission spectrum of UHJ HAT-P-70b using the high-resolution spectrographs CARMENES and PEPSI. Through our cross-correlation analysis, we detected emission signals for Al I, AlH, Ca II, Cr I, Fe I, Fe II, Mg I, Mn I, and Ti I, marking the first detection of Al I and AlH in an exoplanetary atmosphere. Tentative signals of C I, Ca I, Na I, NaH, and Ni I were also identified. Based on those detections, we were able to perform atmospheric retrievals to constrain the thermal profile and elemental abundances of the planet’s dayside hemisphere. The retrieved temperature-pressure profile reveals a strong temperature inversion layer. The chemical free retrieval yielded a metallicity of [Fe/H] = 0.38, while the chemical equilibrium retrieval resulted in [Fe/H] = 0.23, with both values consistent with the solar metallicity. We also tentatively found an enriched abundance of Ni, which could result from the accretion of Ni-rich planetesimals during the planet’s formation.

Ca II triplet emission lines observed with PEPSI. They were combined over two nights and shifted to the planetary rest frame using the best-fit velocity from the Ca II CCF signal. The dashed blue lines indicate the expected positions of the Ca II triplet lines. Among them, only the lines near 8500.36 Å
and 8664.52 Å (vacuum wavelengths) were detected, with their
Gaussian fits in red.

Read more: B. Guo, F. Yan, Th. Henning, et al. 2026, A&A, in press. arXiv:2512.21470

Surface activity of a Rossby sequence of cool Hyades stars

We present quantitative surface-activity information for a sequence of 21 Hyades dwarf stars with effective temperatures all cooler than the red edge of the lithium dip and Rossby (Ro) numbers between 0.14 to 0.54 with respect to the Sun (Ro(Sun)=1). High-resolution high-S/N PEPSI Stokes-IV spectra and least-squares deconvolution of thousands of spectral lines per spectrum are employed for measuring surface magnetic fields, rotational velocities, lithium abundances, and chromospheric CaII IRT fluxes. Lithium abundances A(Li) range from 95 times solar on the warm end of the sample to 1/25 solar on the cool end. We confirm the tight relation with T(eff)  and extent it to K-M stars. A formal relation with rotational period and velocity in the sense higher A(Li) for faster rotators is present. Targets rotating faster than vsini of 6 km/s appear Li saturated. CaII IRT fluxes also show a relation with T(eff), P(rot) and vsini, but opposite to A(Li), in an inverse sense with higher radiative losses for the slower=cooler rotators. Disk-integrated, unsigned, magnetic-field strengths of 15.4±3.6(rms)G are measured for targets warmer than 5000K and 91±61(rms)G for targets cooler than this. These field strengths relate to P(rot), vsini, and Ro, but in a bi-modal fashion. We conclude that the Rossby-number dependency of the surface activity tracers on our Hyades dwarf sequence primarily originates from convective motions, expressed by its turnover time, and only to a smaller and sometimes inverse extent from surface rotation and its related extra mixing.

Lithium-activity-rotation relation versus normalized Rossby number (Ro_n). From left to right: A(Li) logarithmic lithium abundance, R‘(IRT) logarithmic radiative loss in the three CaII infrared-triplet lines, <|B|> disk-integrated, unsigned magnetic field strength in Gauss.

 

Read more: Strassmeier et al. 2025, A&A, in press; e-print arXiv:2510.21456

PEPSI Investigation, Retrieval, and Atlas of Numerous Giant Atmospheres (PIRANGA). IV.High-Resolution Phased-Resolved Spectroscopy of The Ultra Hot Jupiter KELT-20 b

We present five datasets of high-resolution optical emission spectra of the ultra-hot Jupiter KELT-20 b with the PEPSI spectrograph. Using a Bayesian retrieval framework, we constrain its dayside pressure-temperature profile and abundances of Fe, Ni, and Ca, providing the first measurements for Ni and Ca for KELT-20 b in emission. We retrieve the pre- and post-eclipse datasets separately (corresponding to the evening and morning sides, respectively), and compare the constraints on their thermal structures and chemical abundances. We constrain lower abundances in the pre-eclipse datasets compared to the post-eclipse datasets. We interpret these results with an equilibrium chemistry model which suggests ∼ 10 − 30× supersolar refractory abundances. Due to the well-known degeneracy between absolute abundances and continuum opacities, the abundance ratios are more precise probes of the planetary abundances. Therefore we measure the abundance ratios [Ni/Fe] and [Ca/Fe] across these datasets and find they agree within 1σ. We constrain [Ni/Fe] to be consistent with solar within 2σ, and [Ca/Fe] to be 0.001-0.01× solar, not accounting for ionization.

Top-down view of KELT-20 b’s orbit as it passes behind its host star, moving counterclockwise from pre-eclipse to post-eclipse. The orange side depicts the hot day side, while the blue side corresponds to the cooler night side of the tidally-locked planet. The translucent gray side depicts the side of the planet that is outside of the view of the observer, who is positioned off the bottom of the page.

Read more: Bonidie et al. 2025, AJ, in press.  eprint arXiv:2511.09720

PEPSI Investigation, Retrieval, and Atlas of Numerous Giant Atmospheres (PIRANGA). II.Phase-Resolved Cross-Correlation Transmission Spectroscopy of KELT-20b

ELT-20b is a well-studied exoplanet within the highly observable ultra hot Jupiter (UHJ) regime, yet its multidimensional atmospheric structure remains largely unconstrained. Recent advances in instrumentation and increased general circulation model (GCM) complexity have enabled observers to resolve the imprints of more intricate physical mechanisms in time-resolved data. We performed high-resolution cross-correlation transmission spectroscopy (HRCCTS) on a single transit time series of KELT-20b, observed with PEPSI on the LBT. We detect Fe I (11.9σ) and Fe II (23.7σ) and tentatively detect Na I (3.4σ) and Cr I (3.3σ) upon combining nineteen in-transit exposures.

Observability score of each neutral species considered in this study. Colored elements are more likely to be detected since their absorption lines are covered by the PEPSI bandpass ranges we use. Dark cells have fewer and/or shallower lines. Light gray cells were skipped because they don’t have significant lines in the PEPSI bandpass or did not have available opacities. Species with observability score of ∼ 0.2 or greater have been detected in at least one study, except for Ni I.

Read more: Lenhart et al. 2025, AJ, in press. eprint arXiv:2503.07719

PEPSI investigation, retrieval, and atlas of numerous giant atmospheres (PIRANGA) –III. Composition and winds in the atmosphere of TOI-1518b

Ultra-hot Jupiters (UHJs) orbit close to their host stars and experience extreme conditions, making them important laboratories to explore atmospheric composition and dynamics. Transmission spectroscopy is a useful tool to reveal chemical species and their vertical and longitudinal distribution in the atmosphere. We use transmission spectra from the PEPSI (Potsdam Echelle Polarimetric and Spectroscopic Instrument) spectrograph on the Large Binocular Telescope to search for species and measure their time-resolved wind velocities in the atmosphere of TOI-1518 b. We detect Fe I at 7.8 σ and Fe II at 8.9 σ, and tentatively detect Cr I at 4.4 σ and Ni I at 4.0 σ. The time-resolved wind velocities of Fe I show a velocity pattern that is consistent with the velocity pattern of Fe II . TOI-1518 b joins a small sample of UHJs for which time-resolved wind velocities have been measured.

Raw CCF for Fe I. Dark (black) regions are anticorrelated and light (white) regions are correlated. The horizontal blue lines show the ingress and egress of the transit, calculated using the transit duration determined by Duck et al.. The red dotted line shows the expected velocity of the atmospheric signal. , φ is the orbital phase, and Vsys = −11.17 ± 0.035 km s−1 is the systemic velocity offset in the LBT PEPSI frame (S. Petz et al. 2025 ). The dark shaded region is the Doppler shadow.

Read more: Basinger et al. (2025), MNRAS, 543, 4136

A search for Maunder-minimum candidate stars

Stars with very low levels of magnetic activity provide an opportunity for a more quantitative comparison with the Sun during its Maunder minimum. We employ spectra from the RAVE survey in a search for particularly low-activity stars with the goal of identifying candidates for so-called Maunder-minimum stars. Spectra were used to measure the relative flux in the cores of the Ca II infrared-triplet (IRT) lines. Those were converted to absolute emission-line fluxes and were corrected with target fluxes from high-resolution STELLA and ultra-high-resolution PEPSI spectra. Absolute Ca II IRT fluxes for a total of 78 111 RAVE dwarf stars are presented and compared with fluxes of the 123 stars from our high-resolution STELLA+PEPSI sample. RAVE fluxes appear higher than the STELLA and PEPSI fluxes by on average 19% for IRT-1, 21% for IRT-2, and 25% for IRT-3 due to their lower spectral resolution. Our sample also spans a metallicity [Fe/H] range relative to the Sun of −1.5 to +0.5 dex. We confirm the strong dependency of IRT fluxes on metallicity and quantify it to be at most ±14% in the B–V range 0.53–0.73. Without a metallicity correction, practically all very-low-activity RAVE dwarfs show a super-solar metallicity. After correcting for spectral resolution and for metallicity, we find 13 RAVE stars out of 13 326 (0.1%) that fall well below our empirical lower flux bound from high-resolution versus B–V. For solar B–V, this relates to a photospheric uncorrected radiative loss in the IRT lines of log RIRT = −4.13 (≈20% below the solar-minimum value in late 2016). However, 11 targets turned out to be evolved stars based on their Gaia DR3 parallaxes. Two stars, TIC 352227373 (G2V) and TYC 7560-477-1 (G7V), are our only Maunder-minimum candidates from the present search. Contrary to the initial suggestion from the Mount-Wilson H&K Survey, we conclude that such stars are very rare.

Metallicity-corrected RAVE Ca II IRT-1 fluxes (13 326 targets). Crosses are the possible candidate targets. Boldfaced crosses are the two MM candidates, TIC 352227373 (RAVE J191213.1- 760732) and TYC 7560-477-1 (RAVE J025410.6-383603).

Read more: Järvinen & Strassmeier 2025, A&A, 698, A93;  arXiv

PEPSI Investigation, Retrieval, and Atlas of Numerous Giant Atmospheres (PIRANGA). I. The Ubiquity of Fe I Emission and Inversions in Ultra Hot Jupiter Atmospheres

We present high-resolution optical emission spectroscopy observations of the ultra hot Jupiters (UHJs) TOI-1431b and TOI-1518b using the PEPSI spectrograph on the LBT. We detect emission lines from Fe I with a significance of 5.68σ and 7.68σ for TOI 1431b and TOI-1518b, respectively. We also tentatively detect Cr I emission from TOI-1431b at 4.32σ. For TOI-1518 b, we tentatively detect Ni I, Fe II, and Mg I at significance levels ranging from 3−4σ. Detection of emission lines indicates that both planets possess temperature inversions in their atmospheres, providing further evidence of the ubiquity of stratospheres among UHJs. By analyzing the population of hot Jupiters, we compare models that predict the distribution of planets in the temperature-gravity space, and find a recent global circulation model suite from Roth et al. (2024) provides a reasonable match to the observed onset of inversions at Teq∼2000 K. The ubiquity of strong Fe I emission lines among UHJs, together with the paucity of detections of TiO, suggest that atomic iron is the dominant optical opacity source in their atmospheres and can be responsible for the inversions.

Left: Phase coverage of the observations used in this work. The solid inner circle shows the stellar surface, while the middle and outer circles show the orbits of TOI-1518 b and TOI-1431 b, respectively, to scale. The observer is off the bottom of the page, and the vertical dashed lines show the line of sight, such that transit occurs in the lower intersection of the dashed lines and planetary orbits, and the secondary eclipse in the upper intersection. The planets orbit counter clockwise. The colored points show the portions of the orbit where we obtained data. Right: signal-to-noise ratio of our observations as a function of orbital phase. The PEPSI red and blue arms are shown as the solid and dotted lines, respectively. The SNR values shown are, for each spectrum, the 95th quantile per-pixel signal to-noise ratios.

Read more: Petz et al. 2025, AJ, in press; arXiv