Benchmark Brown Dwarf Systems I: Chemical Abundance Analysis of FGK Stars with Wide-Separation Brown Dwarf Companions Using PEPS

Technical Summary: Benchmark Brown Dwarf Systems I: Chemical Abundance Analysis of FGK Stars with Wide-Separation Brown Dwarf Companions Using PEPSI

Problem Statement
Brown dwarfs present a significant challenge in astrophysics due to a mass-temperature-age degeneracy, where their cooling histories make it difficult to independently determine their formation mechanisms and evolutionary states. While "benchmark" systems—where age, mass, or chemistry can be constrained independently—offer a path forward, the chemical composition of brown dwarf companions is often assumed to be solar or identical to their host stars without rigorous verification. Furthermore, most brown dwarfs are found in isolation, but wide-separation brown dwarf companions (>100 AU) to FGK stars offer a unique opportunity to link substellar chemistry directly to the host star's natal environment. However, previous studies have often relied on heterogeneous datasets or assumed solar abundances for host stars, potentially introducing biases in atmospheric retrievals and formation interpretations. There is a critical need for a uniform, high-precision chemical abundance analysis of host stars to serve as a compositional anchor for interpreting the atmospheres of their brown dwarf companions, particularly in the era of JWST.

Methodology
The authors present a spectroscopic survey of 32 FGK stars hosting wide-separation brown dwarfs (with the exception of two closer systems: HD 89744 at 63 AU and BD+01 299 at 35–45 AU).

• Data Acquisition: High-resolution optical spectra were obtained using the PEPSI spectrograph on the Large Binocular Telescope (LBT) between February 2023 and January 2024. The observations utilized two resolution modes: $R = 130,000$ (200 $\mu$m fiber) for the majority of the sample and $R = 50,000$ (300 $\mu$m fiber) for five targets. The spectral coverage focused on cross-dispersers CD III (4800–5441 Å) and CD VI (7419–9140 Å) to access key atomic lines, including the O I triplet.
• Reduction and Calibration: Spectra were reduced using the SDS4PEPSI pipeline, followed by a secondary continuum normalization using low-order polynomial fits to correct for residual continuum variations.
• Stellar Parameter Determination: Initial photometric parameters were derived using EXOFASTv2 with MIST isochrones and multi-band photometry (Gaia, 2MASS, WISE). However, the final spectroscopic parameters ($T_{\text{eff}}$, $\log g$, $[M/H]$, $v_{\text{mic}}$) were determined using the BACCHUS (Brussels Automatic Code for Characterizing High accUracy Spectra) framework. This code employs spectral synthesis within the 1D LTE Turbospectrum radiative transfer code and MARCS model atmospheres, utilizing an excitation-ionization balance method based on Fe I and Fe II lines. For cooler K-dwarfs where Fe II lines were weak, $\log g$ was fixed to photometric values to ensure stability.
• Abundance Analysis: Elemental abundances for 11 elements (C, O, Mg, Si, Ca, Al, Ti, Fe, Y, S, N) were derived using four independent methods within BACCHUS ($\chi^2$, synthesis, equivalent width, and line core matching). Carbon abundances for cool stars were determined using C$_2$ Swan bands, while Nitrogen was derived from CN molecular features. Non-LTE (NLTE) corrections were applied to Oxygen abundances using lookup tables.
• Validation: The methodology was validated against a solar spectrum observed by PEPSI and compared against the Brewer Catalog (J. M. Brewer & D. A. Fischer 2016) and other literature values for overlapping targets.

Key Contributions and Results

• Uniform Abundance Catalog: The paper provides a uniform abundance analysis for 32 host stars, achieving typical errors of $\sim 42$ K in $T_{\text{eff}}$ and $\sim 0.03$ dex in $[Fe/H]$ for high signal-to-noise data ($S/N > 200$).
• Chemical Diversity: The sample exhibits a significant dispersion in chemical ratios compared to the solar neighborhood. While many stars are near-solar, several metal-poor stars were identified (e.g., HIP 63506, NLTT 1011, BD+06 2986), complicating abundance measurements due to weak absorption features.
• Cloud Formation Predictions: Using the host star's Mg/Si ratios, the authors predict the dominant silicate cloud species in the brown dwarf companions. Most systems have Mg/Si ratios $\gtrsim 0.9$, predicting Enstatite + Forsterite clouds. However, a subset (e.g., BD+60 1417, GJ 417) has Mg/Si $\lesssim 0.9$, suggesting the potential formation of Quartz (SiO$_2$) clouds.
• C/O Ratios and Formation Pathways: The study finds a wide range of C/O ratios in the host stars, challenging the assumption of a uniform solar C/O ratio. The authors note that while some brown dwarfs previously reported with super-solar C/O ratios have been re-evaluated (e.g., Gliese 229B), the diversity in host star C/O ratios implies that brown dwarf formation environments are chemically heterogeneous.
• Oxygen Sinks: The authors calculate the fraction of oxygen sequestered in silicate clouds (the "oxygen sink"), finding an average of $21.40 \pm 5.39\%$. This suggests that atmospheric retrievals must account for oxygen depletion to accurately infer bulk C/O ratios.
• Age Constraints via [Y/Mg]: The study evaluates the [Y/Mg] stellar clock for the sample. While the method provides tighter constraints for some systems, it yields implausible ages (exceeding the age of the universe) for four targets, highlighting limitations in applying this clock to all FGK stars without independent age priors.
• Host-Companion Connection: For the system HD 125141, the host star's Mg/Si ratio ($0.95 \pm 0.07$) aligns closely with recent atmospheric retrievals of the brown dwarf companion ($1.00$), supporting the hypothesis that wide-separation brown dwarfs retain the bulk elemental abundance patterns of their host stars.

Significance and Claims
The authors assert that this work provides a timely and uniform abundance analysis essential for the "benchmark brown dwarf" framework. By establishing precise host star abundances, the study:

1. Reduces Degeneracies: Offers a physically meaningful reference point for atmospheric retrievals of brown dwarfs, reducing biases caused by assuming solar abundances.
2. Informs Cloud Physics: Provides the necessary chemical inputs (Mg/Si, Ca/Al) to predict cloud properties and spectral diversity in brown dwarf atmospheres, directly supporting upcoming JWST observations.
3. Tests Formation Models: Enables rigorous tests of whether wide-separation brown dwarfs form via gravitational collapse (sharing host chemistry) or disk processes, by comparing host and companion compositions.
4. Supports JWST Science: The uniform dataset is explicitly designed to support extended wavelength observations of brown dwarf companions by JWST, allowing for more accurate interpretation of their formation histories and atmospheric compositions.

The paper concludes that while the host star abundances generally serve as a reliable proxy for the companion's bulk composition, the observed chemical diversity necessitates case-by-case analysis rather than broad assumptions of solar composition.