2026ApJ...997..118Y

We present atmospheric retrievals of the benchmark brown dwarf binary Luhman 16AB using high-resolution VLT/CRIRES spectra and the differentiable framework ExoJAX. We derive elemental abundances and temperature-pressure (T-P) profiles while explicitly testing the robustness of the results against major sources of systematic uncertainty. We first perform retrievals with a power-law T-P profile and assess the sensitivity of inferred molecular abundances and C/O ratios to different CO line lists (ExoMol, HITEMP with air- and H2-broadening). We then introduce a flexible Gaussian process-based T-P profile, allowing a nonparametric characterization of the thermal structure and a more conservative treatment of uncertainties. For both components, we infer C/O ratios of about 0.67, slightly above solar, with line list systematics at the 7% level emerging as the dominant source of uncertainty, whereas assumptions about T-P parameterization or photometric variability play a lesser role. The retrieved T-P profiles and molecular abundances are broadly consistent with atmospheric models and equilibrium chemistry. Our results establish Luhman 16AB as a key anchor for substellar C/O measurements, demonstrate the utility of flexible T-P modeling in high-resolution retrievals, and highlight the importance of systematic tests-particularly line list uncertainties-for robust comparisons between brown dwarfs and giant exoplanets.