The onset of warps in Spitzer observations of edge-on spiral galaxies

We analyze warps in the nearby edge-on spiral galaxies observed in the {\em Spitzer/IRAC} 4.5 micron band. In our sample of 24 galaxies we find evidence of warp in 14 galaxies. We estimate the observed onset radii for the warps in a subsample of 10 galaxies. The dark matter distribution in each of these galaxies are calculated using the mass distribution derived from the observed light distribution and the observed rotation curves. The theoretical predictions of the onset radii for the warps are then derived by applying a self-consistent linear response theory to the obtained mass models for 6 galaxies with rotation curves in the literature. By comparing the observed onset radii to the theoretical ones, we find that discs with constant thickness can not explain the observations; moderately flaring discs are needed. The required flaring is consistent with the observations. Our analysis shows that the onset of warp is not symmetric in our sample of galaxies. We define a new quantity called the onset-asymmetry index and study its dependence on galaxy properties. The onset asymmetries in warps tend to be larger in galaxies with smaller disc scale lengths. We also define and quantify the global asymmetry in the stellar light distribution, that we call the edge-on asymmetry in edge-on galaxies. It is shown that in most cases the onset asymmetry in warp is actually anti-correlated with the measured edge-on asymmetry in our sample of edge-on galaxies and this could plausibly indicate that the surrounding dark matter distribution is asymmetric.

💡 Research Summary

This paper presents a systematic investigation of warp phenomena in a sample of nearby edge‑on spiral galaxies using Spitzer/IRAC 4.5 µm imaging. The authors selected 24 galaxies with clear edge‑on orientations and identified warp signatures in 14 of them through visual inspection of the vertical light profiles. For a subsample of 10 galaxies, the radial location where the warp first becomes apparent—the “onset radius” (R_onset)—was measured by locating the radius at which the vertical displacement of the mid‑plane deviates significantly from a straight line on either side of the galaxy.

To interpret these measurements, the authors constructed mass models for each galaxy. The 4.5 µm surface‑brightness profiles were converted into stellar mass distributions using appropriate mass‑to‑light ratios (M/L), while published rotation curves supplied the total dynamical mass, allowing the authors to infer the contribution of dark matter. The resulting axisymmetric mass models provide the surface‑density Σ(R) and the vertical scale‑height h(R) of the stellar disc, as well as a parametrized dark‑matter halo (typically a pseudo‑isothermal sphere).

With these mass models in hand, the authors applied a self‑consistent linear response theory. In this framework, an external torque—originating from a misaligned dark‑matter halo, satellite galaxies, or large‑scale tidal fields—is treated as a small perturbation. The linear response of a thin, rotating disc to such a torque can be expressed analytically in terms of Σ(R), h(R), and the epicyclic frequency κ(R). The theory predicts a characteristic radius at which the disc’s vertical displacement begins to grow, i.e., the theoretical warp onset radius (R_theory).

Two families of disc structure were examined: (i) a constant‑thickness disc (h = const.) and (ii) a modestly flaring disc in which the scale‑height increases with radius by roughly 10–30 % over the optical extent. The constant‑thickness models systematically under‑predicted the observed R_onset, often by a factor of two, indicating that a simple flat disc cannot reproduce the data. By contrast, the flaring models yielded R_theory values that matched the observed onset radii within the uncertainties for the six galaxies for which reliable rotation curves were available. This agreement is reinforced by independent measurements of disc flaring from the 4.5 µm light profiles, suggesting that the vertical structure of real stellar discs is indeed mildly flared.

Beyond the radial location, the authors introduced a quantitative “onset‑asymmetry index” (OAI) to capture the difference between the warp onset radii on the two sides of a galaxy. OAI = (R_onset,East – R_onset,West) / (R_onset,East + R_onset,West). An OAI close to zero denotes a symmetric warp, while larger absolute values indicate stronger asymmetry. They found a clear trend: galaxies with smaller exponential disc scale lengths (h_R) tend to have larger |OAI|, implying that compact discs are more susceptible to asymmetric torques.

In parallel, a new metric called “edge‑on asymmetry” (EOA) was defined to quantify the overall left‑right asymmetry of the stellar light distribution in an edge‑on view. EOA is computed as the normalized difference between the total flux on the two sides of the galaxy’s minor axis. Strikingly, for the majority of the sample, OAI and EOA are anti‑correlated: galaxies whose stellar light is brighter on one side often exhibit the warp onset first on the opposite side. This anti‑correlation suggests that the source of the torque is not directly linked to the luminous disc but rather to an asymmetric distribution of dark matter surrounding the galaxy.

The paper’s conclusions can be summarized as follows: (1) Spitzer 4.5 µm imaging provides a robust tracer of the stellar mass distribution and enables precise measurement of warp onset radii; (2) a modestly flaring stellar disc is required to reconcile observed warp onsets with predictions from linear response theory, confirming that real discs are not perfectly flat; (3) warp onset asymmetries are more pronounced in galaxies with short disc scale lengths, indicating a heightened sensitivity to external perturbations; and (4) the observed anti‑correlation between warp onset asymmetry and stellar light asymmetry points to an underlying asymmetry in the dark‑matter halo as a plausible driver of warps.

The study opens several avenues for future work. High‑resolution HI and CO kinematic maps could refine the external torque estimates and test whether gas warps follow the same onset behavior as stellar warps. Numerical simulations that incorporate realistic halo triaxiality, satellite accretion, and disc flaring would help to quantify the relative importance of each mechanism. Finally, extending the analysis to a larger, statistically complete sample—including galaxies with a range of inclinations—would allow the community to assess how universal the flaring requirement and the OAI–EOA anti‑correlation truly are.