A cooler look at the environment of Cygnus X-1: Searching for dynamical interactions within cold molecular gas

Technical Summary: A Cooler Look at the Environment of Cygnus X-1

Problem and Motivation
Black hole X-ray binaries (BHXBs) significantly influence their surrounding interstellar medium (ISM) through relativistic jets and, in the case of high-mass companions, powerful stellar winds. While the feedback mechanisms of jets in active galactic nuclei and low-mass X-ray binaries are relatively well-studied, the interaction of multiple simultaneous outflows—specifically the coupling of a relativistic jet and a massive stellar wind—remains complex and difficult to disentangle using standard continuum imaging. Cygnus X-1 (Cyg X-1), a high-mass X-ray binary (HMXB) consisting of a black hole and an O-type supergiant, presents a unique laboratory for this study. It possesses both persistent relativistic jets and a strong stellar wind ($\dot{M}_w \sim 10^{-6} M_\odot \text{yr}^{-1}$, $v_w \sim 1800\text{--}2600 \text{ km s}^{-1}$). Previous observations have identified a radio/optical bow shock to the north, attributed to the approaching jet or wind, but the molecular gas environment, particularly the potential signatures of dynamical interactions in cold gas, had not been fully characterized.

Methodology
The authors utilized observations from the IRAM–30m telescope conducted on December 12, 2018, to map the molecular gas environment of Cyg X-1. The study targeted the $^{12}\text{CO}(J=2-1)$ and $^{13}\text{CO}(J=2-1)$ transitions at central frequencies of 230.538 GHz and 220.399 GHz, respectively, achieving an angular resolution of approximately $11''$. Data reduction was performed using the GILDAS software package, applying standard calibration, baseline removal, and masking techniques to generate spectral cubes covering the local standard of rest (LSR) velocity range of $[-100, +100] \text{ km s}^{-1}$.

The analysis focused on constructing integrated intensity maps and performing spectral and kinematic analyses on five specific regions of interest surrounding the source:

1. The Spider: An irregular structure immediately surrounding the BHXB.
2. The Southeastern Cloud: A compact region off the receding jet axis.
3. The Approaching Jet Direction: Along the known jet axis toward the radio bow shock.
4. The Handlebar: An elongated feature northeast of the source.
5. The Southwestern Cloud: An extended feature to the southwest.

Position-Velocity (P-V) diagrams were generated to analyze gas kinematics, and spectral profiles were examined for line widths, asymmetries, and multi-peaked structures indicative of shocks or outflows.

Key Results

• Absence of Molecular Gas in the Bow Shock: Contrary to expectations, the region corresponding to the known radio/optical bright bow shock to the north is largely devoid of molecular gas within the observed velocity range.
• Discovery of the "Spider": The authors identified a previously unseen, irregular molecular structure surrounding the BHXB, termed the "spider." Spectral analysis of this region reveals complex, multi-peaked line profiles with significant velocity gradients (up to $\sim 3 \text{ km s}^{-1} \text{ arcmin}^{-1}$). The P-V diagrams show an undulating velocity pattern that is consistently blueshifted relative to the systemic velocity but does not align with a simple jet-axis expansion or the system's proper motion.
• The Southeastern Cloud: A second distinct structure was found in the receding jet direction. It exhibits a clear redshifted velocity gradient increasing with distance from the source, consistent with gas being pushed away, though it is spatially offset by $\sim 25^\circ$ from the known jet axis.
• Other Features: The "handlebar" and "southwestern cloud" structures showed kinematics inconsistent with direct jet interaction (e.g., redshifted motion in the approaching jet direction) or appeared linked to external structures like the Tulip Nebula, suggesting they are not primary products of Cyg X-1 feedback.

Physical Interpretation and Models
The paper proposes two primary scenarios for the origin of the "spider" structure:

1. High Relative Motion (Favored): The authors argue that Cyg X-1 moves through the ambient ISM with a peculiar velocity of $\sim 20 \text{ km s}^{-1}$. This motion, combined with the massive stellar wind, creates a bow shock. As the system moves, the wind sweeps up ISM material which cools radiatively ($t_{\text{cool}} \sim 10^1\text{--}10^3$ years) to form dense molecular gas. The relativistic jets, inflating cocoons, confine this wind-driven shell to the orbital plane, creating the anisotropic, belt-like "spider" morphology. Energy budget calculations suggest the molecular gas luminosity ($L_{\text{CO}} \sim 3 \times 10^{31} \text{ erg s}^{-1}$) is consistent with a momentum-driven interaction where the wind supplies the energy.
2. Low Relative Motion: A scenario where the ambient ISM is co-moving with the source, leading to an energy-driven bubble. However, the authors find this requires fine-tuned conditions and fails to explain the observed low velocities without invoking a foreground cloud interaction, making it less likely.

Regarding the southeastern cloud, while its redshifted kinematics suggest interaction with a receding jet, the significant spatial offset from the jet axis makes a direct causal link uncertain. The authors suggest it may lie within a receding jet cocoon but requires higher resolution to confirm.

Significance and Claims
The paper claims that molecular line imaging is a critical tool for disentangling complex feedback mechanisms in environments with multiple outflows. Specifically, the study demonstrates that:

• Molecular gas tracers can reveal interaction sites (like the spider) that are invisible or distinct in continuum radio or optical observations.
• The "spider" structure provides evidence that stellar winds and relativistic jets can simultaneously shape the ISM, with the wind driving the accumulation of cold gas and the jets sculpting its geometry.
• The lack of molecular gas in the immediate vicinity of the radio bow shock suggests that the bow shock is likely a hot, ionized structure driven by the jet, while the cold molecular component is driven by the stellar wind and system motion.

The authors conclude that while the association between Cyg X-1 and the spider is statistically probable given the unique velocity structure, future high-resolution observations (e.g., with NOEMA) are necessary to fully resolve the 3D geometry and confirm the hydrodynamical interaction models. The work validates the utility of "astrochemical" approaches in constraining jet and wind energetics in HMXBs.