II Zw 096: Luminous Infrared Galaxy Merger

• II Zw 096 is a luminous infrared galaxy merger characterized by a buried, off-nuclear source (D1) that dominates the infrared emission.
• Multiwavelength analyses with Spitzer, JWST, and MUSE reveal complex merger dynamics, intricate gas kinematics, and distinct tidal features.
• Observations indicate the coexistence of a compact, obscured AGN with intense star formation, offering insights into merger-driven compaction and nuclear growth.

Searching arXiv for recent and foundational papers on II Zw 096 to ground the article. {"query":"II Zw 096 arXiv", "max_results": 10, "sort_by": "relevance"} Searching for the JWST and MUSE studies of II Zw 096, plus earlier Spitzer work. {"query":"\"II Zw 096\" JWST MUSE Spitzer arXiv", "max_results": 10, "sort_by": "relevance"} II Zw 096 is a luminous infrared galaxy merger, also known as CGCG 448-020 and IRAS 20550+1655, whose defining characteristic is that the dominant power source is an exceptionally obscured off-nuclear region rather than either optical nucleus. Across Spitzer, HST, Chandra, AKARI, JWST, ALMA, and VLT/MUSE studies, the system is described through the main galaxies A and B, the eastern C+D complex, and the faint E region; within D, the compact source D1 is the most enshrouded component and dominates a large fraction of the infrared luminosity (Inami et al., 2010, García-Bernete et al., 2024, Riesco et al., 8 Jul 2025).

1. System identification and large-scale structure

The system has been classified as a merging LIRG with $\log(L_{\rm IR}/L_\odot)=11.94$ and was described in earlier GOALS work as a Class III merger with two identifiable nuclei and well-developed tidal structure (Inami et al., 2010). A later optical IFU study gives $L_{\rm IR}=8.7\times10^{11}\,L_\odot$ and classifies the object as a stage c merger, emphasizing strong tidal tails, amorphous or disturbed disks, and clear pre-coalescence interaction signatures (Riesco et al., 8 Jul 2025). Recent literature reports $z=0.036$, $z=0.0362$, and $z=0.0365$, with luminosity distances of 161 Mpc and 159 Mpc in the Spitzer-era and JWST studies, respectively (Inami et al., 2010, García-Bernete et al., 2024).

Morphologically, the traditional nomenclature identifies source A as the southern nucleus and source B as the northern nucleus, with the red off-nuclear complex split into sources C and D (Inami et al., 2010). Higher-resolution HST/NICMOS imaging resolved D into D0 and D1, and JWST work treats D1 as the most enshrouded source in the interacting system (García-Bernete et al., 2024). The 2025 MUSE analysis argues that II Zw 096 is more complex than a simple binary merger and contains three or more distinct galaxies, with the western side dominated by the still recognizable spiral galaxies A and B, and the eastern side containing the highly disrupted C+D complex and the faint E region (Riesco et al., 8 Jul 2025).

This shift in description is central to the modern interpretation. Earlier work foregrounded the contrast between two optical nuclei and one buried extranuclear infrared source (Inami et al., 2010). Later optical IFU work instead interprets the full system as a compact group-like multiple merger or collapsing small galaxy group, with different subcomponents at different dynamical stages (Riesco et al., 8 Jul 2025). This suggests that II Zw 096 is unusual not only because of its buried power source, but also because of the multiplicity and asymmetry of the interaction.

2. The buried infrared powerhouse: source D and the D1 nucleus

The defining empirical result of the 2010 multiwavelength analysis is that roughly 80% of the total infrared output is generated by a compact source identified as source D, part of the red C+D complex east of the main merging disks (Inami et al., 2010). Using the Spitzer 24 and 70 $\mu$m data, that work estimated $$L_{\rm IR}(8-1000\,\mu{\rm m}) \approx 6.87\times10^{11}\,L_\odot$$ for source D and found that source D contributes $67.2\pm10.6\%$ of the total galaxy flux at 24 $\mu$m and $87.9\pm13.8\%$ at 70 $L_{\rm IR}=8.7\times10^{11}\,L_\odot$0m, with measured fluxes $L_{\rm IR}=8.7\times10^{11}\,L_\odot$1 Jy and $L_{\rm IR}=8.7\times10^{11}\,L_\odot$2 Jy, respectively (Inami et al., 2010). At the spatial resolution of those data, source D was already recognized as extraordinary because it is off-nuclear yet system-dominating.

JWST spectroscopy refined this picture by isolating D1 inside the D region and showing that D1 alone contributes 40–70% of the total $L_{\rm IR}=8.7\times10^{11}\,L_\odot$3 emission (Riesco et al., 8 Jul 2025). The JWST NIRSpec+MIRI/MRS study detected gas-phase H$L_{\rm IR}=8.7\times10^{11}\,L_\odot$4O $L_{\rm IR}=8.7\times10^{11}\,L_\odot$5 at $L_{\rm IR}=8.7\times10^{11}\,L_\odot$6–$L_{\rm IR}=8.7\times10^{11}\,L_\odot$7m and $L_{\rm IR}=8.7\times10^{11}\,L_\odot$8CO $L_{\rm IR}=8.7\times10^{11}\,L_\odot$9 at $z=0.036$0–$z=0.036$1m in D1 (García-Bernete et al., 2024). The CO band shows the R- and P-branches in emission, while the H$z=0.036$2O band shows the P-branch in emission and the R-branch in absorption (García-Bernete et al., 2024). The H$z=0.036$3O R-branch absorption is interpreted as the signature of a compact IR-bright embedded source, whereas the broad CO component is interpreted as tracing a highly turbulent environment (García-Bernete et al., 2024).

The same JWST analysis modeled D1 with four components: a torus, an outflow or turbulent gas component, an SF-extended warm component, and a cold-extended component (García-Bernete et al., 2024). The preferred torus solution gives a diameter of $z=0.036$4 pc at $z=0.036$5m, a diameter of $z=0.036$6 pc at 1 mm, $z=0.036$7 K, $z=0.036$8, inclination $z=0.036$9, half-opening angle $z=0.0362$0, and bolometric luminosity $z=0.0362$1 (García-Bernete et al., 2024). Using

$z=0.0362$2

the study derives $z=0.0362$3 for the compact dusty structure (García-Bernete et al., 2024).

These results transformed D1 from a compact red knot into a resolved physical scenario: a very small, warm, optically thick dusty nucleus embedded within a broader star-forming and turbulent molecular environment (García-Bernete et al., 2024). The earlier statement that source D dominates the system and the later statement that D1 alone contributes 40–70% of the total infrared emission are therefore not competing descriptions; they correspond to different spatial decompositions of the same buried eastern power source (Inami et al., 2010, Riesco et al., 8 Jul 2025).

3. Starburst versus obscured AGN

The earliest detailed interpretation of the buried source favored a starburst. Spitzer IRS spectroscopy showed strong PAH emission, deep silicate absorption, and strong low- to intermediate-ionization fine-structure lines, but no high-ionization AGN lines (Inami et al., 2010). The paper quotes a $z=0.0362$4 upper limit $z=0.0362$5, $z=0.0362$6, and therefore $z=0.0362$7; it also gives $z=0.0362$8 and $z=0.0362$9 (Inami et al., 2010). Strong 3.3 $z=0.0365$0m and 6.2 $z=0.0365$1m PAH emission, with 3.3 $z=0.0365$2m EQW $z=0.0365$3m and 6.2 $z=0.0365$4m EQW $z=0.0365$5m, were taken to imply a starburst-dominated energy source (Inami et al., 2010). The same work also found $z=0.0365$6, consistent with star formation rather than AGN radio emission (Inami et al., 2010).

Later studies do not remove that starburst component, but they substantially strengthen the obscured-AGN interpretation. The JWST molecular-band analysis argues that the compact warm dusty core, the very high $z=0.0365$7, the H$z=0.0365$8O R-branch absorption, HCN and C$z=0.0365$9H$\mu$0 absorption, HCN-vib emission, and OH megamasers are difficult to reconcile with a pure starburst, and therefore favor a dust-obscured AGN with an exceptionally high covering factor (García-Bernete et al., 2024). The paper notes that a pure starburst would require $\mu$1 O stars within $\mu$2 pc and that the derived surface brightness is $\mu$3 times above the theoretical maximum warm starburst scale of $\mu$4 (García-Bernete et al., 2024).

The optical IFU evidence points in the same direction. In the MUSE NFM data, integrated D1 is classified as Seyfert in the [N II] BPT plane, while in the [S II] and [O I] planes the nominal point falls in SF but the error bars allow Seyfert; some pixels within D1 itself are directly classified as Seyfert (Riesco et al., 8 Jul 2025). The study further states that starburst plus shock models cannot reproduce the observed D1 line ratios and that an additional ionization source is required (Riesco et al., 8 Jul 2025). It quotes an extinction-corrected $\mu$5, $\mu$6 for D1, and an upper-limit-style inference $\mu$7 if the 100 GHz emission is AGN dominated (Riesco et al., 8 Jul 2025).

The resulting picture is not a replacement of “starburst” by “AGN,” but a change from “buried compact starburst” to “heavily obscured AGN embedded in a dense compact starburst” (Inami et al., 2010, García-Bernete et al., 2024, Riesco et al., 8 Jul 2025). This suggests that the main disagreement among studies is driven less by contradictory data than by improved spatial resolution and access to diagnostics—especially mid-IR molecular bands and AO-assisted optical IFU spectroscopy—that were unavailable in the earlier work.

4. Gas kinematics, shocks, and merger dynamics

The 2025 MUSE study provides the most detailed dynamical description of II Zw 096 (Riesco et al., 8 Jul 2025). It combines Wide Field Mode observations for system-wide context with Narrow Field Mode adaptive-optics observations for the compact inner regions. The WFM seeing FWHM is 0.76–0.94 arcsec, corresponding to 0.55–0.74 kpc, whereas the NFM analysis assumes $\mu$8, corresponding to $\mu$9 pc (Riesco et al., 8 Jul 2025). Emission-line fitting used H$$L_{\rm IR}(8-1000\,\mu{\rm m}) \approx 6.87\times10^{11}\,L_\odot$$0, [O III], [O I], H$$L_{\rm IR}(8-1000\,\mu{\rm m}) \approx 6.87\times10^{11}\,L_\odot$$1, [N II], and [S II], with one- and two-component Gaussian models and fixed [O III] and [N II] doublet ratios of 1:3 (Riesco et al., 8 Jul 2025).

At least two rotating components are identified. For galaxy A, the fitted rotating disk has $$L_{\rm IR}(8-1000\,\mu{\rm m}) \approx 6.87\times10^{11}\,L_\odot$$2 km s$$L_{\rm IR}(8-1000\,\mu{\rm m}) \approx 6.87\times10^{11}\,L_\odot$$3, $$L_{\rm IR}(8-1000\,\mu{\rm m}) \approx 6.87\times10^{11}\,L_\odot$$4, $$L_{\rm IR}(8-1000\,\mu{\rm m}) \approx 6.87\times10^{11}\,L_\odot$$5, radial extent 2 kpc, and maximum velocity 57 km s$$L_{\rm IR}(8-1000\,\mu{\rm m}) \approx 6.87\times10^{11}\,L_\odot$$6 (Riesco et al., 8 Jul 2025). For galaxy B, the fitted parameters are $$L_{\rm IR}(8-1000\,\mu{\rm m}) \approx 6.87\times10^{11}\,L_\odot$$7 km s$$L_{\rm IR}(8-1000\,\mu{\rm m}) \approx 6.87\times10^{11}\,L_\odot$$8, $$L_{\rm IR}(8-1000\,\mu{\rm m}) \approx 6.87\times10^{11}\,L_\odot$$9, $67.2\pm10.6\%$0, 9 kpc, and 165 km s$67.2\pm10.6\%$1, respectively (Riesco et al., 8 Jul 2025). Galaxy B is therefore the clearest large rotating disk, while galaxy A shows stronger tidal distortion and non-circular motion (Riesco et al., 8 Jul 2025).

A major kinematic result is the prevalence of double-peaked emission lines in overlap regions and around C+D (Riesco et al., 8 Jul 2025). The authors interpret these profiles mainly as superposed merger structures at different velocities rather than classical bipolar outflows (Riesco et al., 8 Jul 2025). Regions with total fitted $67.2\pm10.6\%$2 km s$67.2\pm10.6\%$3 are frequently attributed to unresolved blending of multiple narrow components (Riesco et al., 8 Jul 2025). In the NFM-East data, D1 and D0 have $67.2\pm10.6\%$4 km s$67.2\pm10.6\%$5, while a region between C1, C2, and C3 reaches $67.2\pm10.6\%$6 km s$67.2\pm10.6\%$7 (Riesco et al., 8 Jul 2025).

One special case is the broad component near the center of galaxy A, labeled PO, with velocity offsets of about $67.2\pm10.6\%$8 to $67.2\pm10.6\%$9 km s$\mu$0, $\mu$1 km s$\mu$2 and locally $\mu$3 km s$\mu$4, and $\mu$5 (Riesco et al., 8 Jul 2025). The paper interprets this as a potential outflow but does not treat that interpretation as secure (Riesco et al., 8 Jul 2025). Another important structure is the redshifted stripe north of D that appears to converge toward D1 and D0; this is interpreted as possible inflowing material feeding the buried nucleus (Riesco et al., 8 Jul 2025).

Ionization diagnostics reinforce the dynamical picture. The MUSE analysis employs the standard BPT-style planes

$\mu$6

with corresponding [S II] and [O I] variants (Riesco et al., 8 Jul 2025). Galaxy A and most of galaxy B are SF-dominated, but the tidal tail from A toward C+D, the outskirts of C and D, the overlap regions, and the PO component show mixed excitation or strong shock signatures (Riesco et al., 8 Jul 2025). The tail from A toward C+D is especially notable: it is composite in [N II], SF to Seyfert in [S II], and SF to Seyfert/LINER-like in [O I], with inferred shock fractions of 40–100% (Riesco et al., 8 Jul 2025).

5. Star formation, stellar populations, and the interstellar medium

Star formation in II Zw 096 is both intense and temporally structured. The 2010 HST analysis identified 128 clusters, with 97 detected in FUV and 88 in $\mu$7, and found that the color distributions are consistent with at least two populations: one with ages of roughly 1–5 Myr and one with ages of roughly 20–500 Myr (Inami et al., 2010). The inferred cluster masses span $\mu$8–$\mu$9 (Inami et al., 2010). Source D itself was assigned a stellar mass of $87.9\pm13.8\%$0–$87.9\pm13.8\%$1, assuming an age of roughly 5–7 Myr and $87.9\pm13.8\%$2 mag (Inami et al., 2010).

Integrated star-formation indicators are correspondingly large. From the infrared luminosity of source D, the Spitzer-era work estimated $87.9\pm13.8\%$3, while AKARI Br$87.9\pm13.8\%$4 gave an uncorrected $87.9\pm13.8\%$5 and about $87.9\pm13.8\%$6 after adopting $87.9\pm13.8\%$7 mag (Inami et al., 2010). The same study derived $87.9\pm13.8\%$8 and therefore $87.9\pm13.8\%$9 mag toward source C+D (Inami et al., 2010). In the optical IFU work, extinction-corrected H$L_{\rm IR}=8.7\times10^{11}\,L_\odot$00 gives $L_{\rm IR}=8.7\times10^{11}\,L_\odot$01 for galaxy A, $L_{\rm IR}=8.7\times10^{11}\,L_\odot$02 for galaxy B, and $L_{\rm IR}=8.7\times10^{11}\,L_\odot$03 for region E, while the appendix uses $L_{\rm IR}=8.7\times10^{11}\,L_\odot$04 for D1 when constructing starburst models (Riesco et al., 8 Jul 2025).

The molecular and nebular ISM also show distinct physical components. Warm molecular hydrogen traced by H$L_{\rm IR}=8.7\times10^{11}\,L_\odot$05 S(3), S(2), and S(1) yields a temperature of 329 K and $L_{\rm IR}=8.7\times10^{11}\,L_\odot$06 (Inami et al., 2010). The [Ne III]/[Ne II] ratio is $L_{\rm IR}=8.7\times10^{11}\,L_\odot$07, which was interpreted as evidence for a hard radiation field associated with many massive stars (Inami et al., 2010). The [S III] $L_{\rm IR}=8.7\times10^{11}\,L_\odot$08m ratio of $L_{\rm IR}=8.7\times10^{11}\,L_\odot$09 implies $L_{\rm IR}=8.7\times10^{11}\,L_\odot$10 for $L_{\rm IR}=8.7\times10^{11}\,L_\odot$11 K over the large LH aperture (Inami et al., 2010).

JWST resolves the obscured eastern source into a compact buried nucleus plus extended PDR-like star formation (García-Bernete et al., 2024). The narrow CO emission is attributed to an SF-extended warm component, with $L_{\rm IR}=8.7\times10^{11}\,L_\odot$12, similar to the Orion Bar (García-Bernete et al., 2024). The 7.7 $L_{\rm IR}=8.7\times10^{11}\,L_\odot$13m to 11.3 $L_{\rm IR}=8.7\times10^{11}\,L_\odot$14m PAH ratio is $L_{\rm IR}=8.7\times10^{11}\,L_\odot$15, interpreted as evidence for a strong ionized-PAH component and intense circumnuclear PDR activity (García-Bernete et al., 2024). The appendix further derives $L_{\rm IR}=8.7\times10^{11}\,L_\odot$16 K from CO(3–2), consistent with widespread PDRs inside the JWST beam (García-Bernete et al., 2024).

6. Physical interpretation and broader significance

II Zw 096 occupies a distinctive position among nearby interacting infrared galaxies because the energetically dominant source is off-nuclear and deeply obscured. The 2010 study compared the buried source to the extranuclear starbursts in the Antennae and Arp 299, but emphasized that source D is more than an order of magnitude more luminous than the Antennae overlap-region source and contributes $L_{\rm IR}=8.7\times10^{11}\,L_\odot$17 of the system luminosity rather than a minor fraction (Inami et al., 2010). The later JWST and MUSE results retain the off-nuclear emphasis but reinterpret the buried source as a compact obscured nucleus hosting a likely AGN embedded in vigorous star formation (García-Bernete et al., 2024, Riesco et al., 8 Jul 2025).

Two common simplifications are therefore misleading. The first is that II Zw 096 is a standard two-galaxy merger. The optical IFU data instead support three or more distinct galaxies, a triple overlap region, and an eastern side that may be a disrupted companion or remnant undergoing a different merger stage from the western pair (Riesco et al., 8 Jul 2025). The second is that the D region is either a pure starburst or an unambiguous AGN. The observational record is more specific: early mid-IR diagnostics showed a starburst-dominated spectrum with no secure direct AGN lines, whereas later JWST molecular-band modeling and high-resolution optical IFU data strongly favor a heavily obscured AGN coexisting with a dense compact starburst (Inami et al., 2010, García-Bernete et al., 2024, Riesco et al., 8 Jul 2025).

This makes II Zw 096 a particularly valuable local laboratory for merger-driven compaction, obscured black-hole growth, and the coexistence of shocks, compact starbursts, and buried nuclear activity. The D1 nucleus links several normally separate observational regimes: pc-scale warm dust, rovibrational H$L_{\rm IR}=8.7\times10^{11}\,L_\odot$18O and CO bands, radio/mm compactness, extreme optical extinction, Seyfert-like line ratios in high-resolution optical data, and a merger geometry capable of driving both star formation and inflow (García-Bernete et al., 2024, Riesco et al., 8 Jul 2025). A plausible implication is that the system is an unusually clear case in which improved angular resolution progressively revealed that the same buried source can appear starburst-dominated in low-resolution global diagnostics while harboring a compact obscured AGN at its core.