Far-infrared (FIR) images and photometry are presented for 201 Luminous and Ultraluminous Infrared Galaxies [LIRGs: log$(L_{\rm IR}/L_\odot) = 11.00 - 11.99$, ULIRGs: log$(L_{\rm IR}/L_\odot) = 12.00 - 12.99$], in the Great Observatories All-Sky LIRG Survey (GOALS) based on observations with the $Herschel$ $Space$ $Observatory$ Photodetector Array Camera and Spectrometer (PACS) and the Spectral and Photometric Imaging Receiver (SPIRE) instruments. The image atlas displays each GOALS target in the three PACS bands (70, 100, and 160 $\mu$m) and the three SPIRE bands (250, 350, and 500 $\mu$m), optimized to reveal structures at both high and low surface brightness levels, with images scaled to simplify comparison of structures in the same physical areas of $\sim$$100\times100$ kpc$^2$. Flux densities of companion galaxies in merging systems are provided where possible, depending on their angular separation and the spatial resolution in each passband, along with integrated system fluxes (sum of components). This dataset constitutes the imaging and photometric component of the GOALS Herschel OT1 observing program, and is complementary to atlases presented for the Hubble Space Telescope (Evans et al. 2017, in prep.), Spitzer Space Telescope (Mazzarella et al. 2017, in prep.), and Chandra X-ray Observatory (Iwasawa et al. 2011, 2017, in prep.). Collectively these data will enable a wide range of detailed studies of AGN and starburst activity within the most luminous infrared galaxies in the local Universe.

We present an analysis of the morphological classification of 89 luminous infrared galaxies (LIRGs) from the Great Observatories All-sky LIRG Survey (GOALS) sample using non-parametric coefficients and compare their morphology as a function of wavelength. We rely on images obtained in the optical (B- and I-band) as well as in the infrared (H-band and 5.8$\mu$m). Our classification is based on the calculation of $Gini$ and the second order of light ($M_{20}$) non-parametric coefficients which we explore as a function of stellar mass ($M_\star$), infrared luminosity ($L_{IR}$) and star formation rate (SFR). We investigate the relation between $M_{20}$, the specific SFR (sSFR) and the dust temperature ($T_{dust}$) in our galaxy sample. We find that $M_{20}$ is a better morphological tracer than $Gini$, as it allows to distinguish systems formed by double systems from isolated and post-merger LIRGs. The multi-wavelength analysis allows us to identify a region in the $Gini$-$M_{20}$ parameter space where ongoing mergers reside, regardless of the band used to calculate the coefficients. In particular when measured in the H-band, this region can be used to identify ongoing mergers, with a minimal contamination from LIRGs in other stages. We also find that while the sSFR is positively correlated with $M_{20}$ when measured in the mid-infrared, i.e. star-bursting galaxies show more compact emission, it is anti-correlated with the B-band based $M_{20}$. We interpret this as the spatial decoupling between obscured and un-obscured star formation, whereby the ultraviolet/optical size of a LIRGs experience an intense dust enshrouded central starburst is larger than in the one in the mid-infrared since the contrast between the nuclear to the extended disk emission is smaller in the mid-infrared. This has important implications for high redshift surveys of dusty sources. [abridged]

We present Australia Telescope Compact Array radio data of the dwarf irregular galaxy ESO 324-G024 which is seen in projection against the giant, northern lobe of the radio galaxy Centaurus A (Cen A, NGC 5128). The distorted morphology and kinematics of ESO 324-G024, as observed in the 21 cm spectral line emission of neutral hydrogen, indicate disruptions by external forces. We investigate whether tidal interactions and/or ram pressure stripping are responsible for the formation of the HI tail stretching to the northeast of ESO 324-G024 with the latter being most probable. Furthermore, we closely analyze the sub-structure of Cen A's polarized radio lobes to ascertain whether ESO 324-G024 is located in front, within or behind the northern lobe. Our multi-wavelength, multi-component approach allows us to determine that ESO 324-G024 is most likely behind the northern radio lobe of Cen A. This result helps to constrain the orientation of the lobe, which is likely inclined to our line of sight by approximately 60 degrees if NGC 5128 and ESO 324-G024 are at the same distance.

We present ALMA Cycle-0 observations of the CO (6-5) line emission and of the 435um dust continuum emission in the central kpc of NGC 1614, a local luminous infrared galaxy (LIRG) at a distance of 67.8 Mpc (1 arcsec = 329 pc). The CO emission is well resolved by the ALMA beam (0".26 x 0".20) into a circum-nuclear ring, with an integrated flux of f_{CO(6-5)} = 898 (+-153) Jy km/s, which is 63(+-12)% of the total CO(6-5) flux measured by Herschel. The molecular ring, located between 100pc < r < 350pc from the nucleus, looks clumpy and includes seven unresolved (or marginally resolved) knots with median velocity dispersion of 40 km/s. These knots are associated with strong star formation regions with \Sigma_{SFR} 100 M_\sun/yr/kpc^{2} and \Sigma_{Gas} 1.0E4 M_\sun/pc^{2}. The non-detections of the nucleus in both the CO (6-5) line emission and the 435um continuum rule out, with relatively high confidence, a Compton-thick AGN in NGC 1614. Comparisons with radio continuum emission show a strong deviation from an expected local correlation between \Sigma_{Gas} and \Sigma_{SFR}, indicating a breakdown of the Kennicutt-Schmidt law on the linear scale of 100 pc.

The Great Observatories All-Sky LIRG Survey (GOALS) is a comprehensive, multiwavelength study of luminous infrared galaxies (LIRGs) in the local universe. Here we present the results of a multi-component, spectral decomposition analysis of the low resolution mid-IR Spitzer IRS spectra from 5-38um of 244 LIRG nuclei. The detailed fits and high quality spectra allow for characterization of the individual PAH features, warm molecular hydrogen emission, and optical depths for silicate dust grains and water ices. We find that starbursting LIRGs, which make up the majority of GOALS, are very consistent in their MIR properties (i.e. tau_9.7um, tau_ice, neon line and PAH feature ratios). However, as their PAH EQW decreases, usually an indicator of an increasingly dominant AGN, LIRGs cover a larger spread in these MIR parameters. The contribution from PAHs to the total L(IR) in LIRGs varies from 2-29% and LIRGs prior to their first encounter show higher L(PAH)/L(IR) ratios on average. We observe a correlation between the strength of the starburst (IR8) and the PAH fraction at 8um but not with the 7.7 to 11.3 PAH ratio, suggesting the fractional PDR emission, and not the overall grain properties, is associated with the rise in IR8 for galaxies off the starburst main sequence. We detect crystalline silicate features in 6% of the sample but only in the most obscured sources (s_9.7um < -1.24). Ice absorption features are observed in 11% (56%) of GOALS LIRGs (ULIRGs). Most GOALS LIRGs have L(H2)/L(PAH) ratios elevated above those observed for normal star-forming galaxies and exhibit a trend for increasing L(H2)/L(PAH) ratio with increasing L(H2). While star formation appears to be the dominant process responsible for exciting the H2 in most of the GOALS galaxies, a subset of LIRGs (10%) show excess H2 emission that is inconsistent with PDR models and may be excited by shocks or AGN-induced outflows.

While most ram pressure studies have focused on ram pressure stripping in galaxy clusters, we devise a novel approach based on a kinematic measurement of ram pressure perturbations in HI velocity fields for intergalactic material (IGM) densities and relative velocities that are one to two orders of magnitude lower than in galaxies showing ram pressure stripping. Our model evaluates ram pressure induced kinematic terms in gas disks with constant inclination as well as those with a warped geometry. Ram pressure perturbations are characterized by kinematic modes of even order, m=0 and m=2, corresponding to a ram wind perpendicular and parallel to the gas disk, respectively. Long-term consequences of ram pressure, such as warped disks as well as uncertainties in the disk geometry typically generate uneven modes (m=1 and m=3), that are clearly distinguishable from the kinematic ram pressure terms. We have applied our models to three nearby isolated galaxies, utilizing Markov Chain Monte Carlo fitting routines to determine ram pressure perturbations in the velocity fields of NGC 6946 and NGC 3621 of ~30km s$^{-1}$ (effective line-of-sight velocity change) at HI column densities below (4-10)$\times$10$^{20}$cm$^{-2}$ (at radial scales greater than ~15kpc). In contrast, NGC 628 is dominated by a strongly warped disk. Our model fits reveal the three-dimensional vector of the galaxies' movement with respect to the IGM rest-frame and provide constraints on the product of speed with IGM density, opening a new window for extragalactic velocity measurements and studies of the intergalactic medium.

We present Herschel/PACS observations of extended [CII]157.7{\mu}m line emission detected on ~ 1 - 10 kpc scales in 60 local luminous infrared galaxies (LIRGs) from the Great Observatories All-sky LIRG Survey (GOALS). We find that most of the extra-nuclear emission show [CII]/FIR ratios >~ 4 x 10^-3, larger than the mean ratio seen in the nuclei, and similar to those found in the extended disks of normal star-forming galaxies and the diffuse inter-stellar medium (ISM) of our Galaxy. The [CII] "deficits" found in the most luminous local LIRGs are therefore restricted to their nuclei. There is a trend for LIRGs with warmer nuclei to show larger differences between their nuclear and extra-nuclear [CII]/FIR ratios. We find an anti-correlation between [CII]/FIR and the luminosity surface density, {\Sigma}_IR, for the extended emission in the spatially-resolved galaxies. However, there is an offset between this trend and that found for the LIRG nuclei. We use this offset to derive a beam filling-factor for the star-forming regions within the LIRG disks of ~ 6 % relative to their nuclei. We confront the observed trend to photo-dissociation region (PDR) models and find that the slope of the correlation is much shallower than the model predictions. Finally, we compare the correlation found between [CII]/FIR and {\Sigma}_IR with measurements of high-redshift starbursting IR-luminous galaxies.

We present ALMA Cycle-0 observations of the CO (6-5) line emission (rest-frame frequency = 691.473 GHz) and of the 435$\mu m$ dust continuum emission in the nuclear region of NGC 34, a local luminous infrared galaxy (LIRG) at a distance of 84 Mpc (1" = 407 pc) which contains a Seyfert 2 active galactic nucleus (AGN) and a nuclear starburst. The CO emission is well resolved by the ALMA beam ($\rm 0."26\times 0."23$), with an integrated flux of $\rm f_{CO~(6-5)} = 1004\; (\pm 151) \; Jy\; km\; s^{-1}$. Both the morphology and kinematics of the CO (6-5) emission are rather regular, consistent with a compact rotating disk with a size of 200 pc. A significant emission feature is detected on the red-shifted wing of the line profile at the frequency of the $\rm H^{13}CN\; (8-7)$ line, with an integrated flux of $\rm 17.7 \pm 2.1 (random) \pm 2.7 (sysmatic)\; Jy\;km\; s^{-1}$. However, it cannot be ruled out that the feature is due to an outflow of warm dense gas with a mean velocity of $\rm 400\; km\; s^{-1}$. The continuum is resolved into an elongated configuration, and the observed flux corresponds to a dust mass of $\rm M_{dust} = 10^{6.97\pm 0.13}\; M_{sun}$. An unresolved central core ($\rm radius \simeq 50\; pc$) contributes $28\%$ of the continuum flux and $19\%$ of the CO (6-5) flux, consistent with insignificant contributions of the AGN to both emissions. Both the CO (6-5) and continuum spatial distributions suggest a very high gas column density ($\rm >= 10^4\; M_{sun}\; pc^{-2}$) in the nuclear region at $\rm radius <= 100\; pc$.

We consider the most commonly occurring circumstances which apply to galaxies, namely membership in galaxy groups of about $10^{13}h^{-1} M_\odot$ total mass, and estimate the accompanying physical conditions of intergalactic medium (IGM) density and the relative galaxy-IGM space velocity. We then investigate the dynamical consequences of such a typical galaxy-IGM interaction on a rotating gaseous disk within the galaxy potential. We find that the rotating outer disk is systematically distorted into a characteristic "warp" morphology, of the type that has been well-documented in the majority of well-studied nearby systems. The distortion is established rapidly, within two rotation periods, and is long-lived, surviving for at least ten. A second consequence of the interaction is the formation of a one arm retrograde spiral wave pattern that propagates in the disk. We suggest that the ubiquity of the warp phenomenon might be used to reconstruct both the IGM density profile and individual member orbits within galaxy groups.

We present the data and our analysis of MIR fine-structure emission lines detected in Spitzer/IRS high-res spectra of 202 local LIRGs observed as part of the GOALS project. We detect emission lines of [SIV], [NeII], [NeV], [NeIII], [SIII]18.7, [OIV], [FeII], [SIII]33.5, and [SiII]. Over 75% of our galaxies are classified as starburst (SB) sources in the MIR. We compare ratios of the emission line fluxes to stellar photo- and shock-ionization models to constrain the gas properties in the SB nuclei. Comparing the [SIV]/[NeII] and [NeIII]/[NeII] ratios to the Starburst99-Mappings III models with an instantaneous burst history, the line ratios suggest that the SB in our LIRGs have ages of 1-4.5Myr, metallicities of 1-2Z_sun, and ionization parameters of 2-8e7cm/s. Based on the [SIII]/[SIII] ratios, the electron density in LIRG nuclei has a median electron density of ~300cm-3 for sources above the low density limit. We also find that strong shocks are likely present in 10 SB sources. A significant fraction of the GOALS sources have resolved neon lines and 5 show velocity differences of >200km/s in [NeIII] or [NeV] relative to [NeII]. Furthermore, 6 SB and 5 AGN LIRGs show a trend of increasing line width with ionization potential, suggesting the possibility of a compact energy source and stratified ISM in their nuclei. We confirm a strong correlation between the [NeII]+[NeIII] emission, as well as [SIII]33.5, with both the IR luminosity and the 24um warm dust emission measured from the spectra. Finally, we find no correlation between the hardness of the radiation field or the line width and the ratio of the total IR to 8um emission (IR8). This may be because the IR luminosity and the MIR fine-structure lines are sensitive to different timescales over the SB, or that IR8 is more sensitive to the geometry of the warm dust region than the radiation field producing the HII region emission.

We present the first results of a survey of the [CII]158um emission line in 241 luminous infrared galaxies (LIRGs) comprising the Great Observatories All-sky Survey (GOALS) sample, obtained with the PACS instrument on board Herschel. The [CII] luminosities of the LIRGs in GOALS range from ~10^7 to 2x10^9 Lsun. We find that LIRGs show a tight correlation of [CII]/FIR with far-IR flux density ratios, with a strong negative trend spanning from ~10^-2 to 10^-4, as the average temperature of dust increases. We find correlations between the [CII]/FIR ratio and the strength of the 9.7um silicate absorption feature as well as with the luminosity surface density of the mid-IR emitting region (Sigma_MIR), suggesting that warmer, more compact starbursts have substantially smaller [CII]/FIR ratios. Pure star-forming (SF) LIRGs have a mean [CII]/FIR ~ 4x10^-3, while galaxies with low 6.2um PAH equivalent widths (EWs), indicative of the presence of active galactic nuclei (AGN), span the full range in [CII]/FIR. However, we show that even when only pure SF galaxies are considered, the [CII]/FIR ratio drops by an order of magnitude, from 10^-2 to 10^-3, with Sigma_MIR and Sigma_IR, implying that the [CII] luminosity is not a good indicator of the star formation rate (SFR) for most LIRGs, for it does not scale linearly with the warm dust emission. Moreover, even in LIRGs in which we detect an AGN in the mid-IR, the majority (2/3) of galaxies show [CII]/FIR >= 10^-3 typical of high 6.2um PAH EW sources, suggesting that most AGNs do not contribute significantly to the far-IR emission. We provide an empirical relation between the [CII]/FIR and the specific SFR (SSFR) for SF LIRGs. Finally, we present predictions for the starburst size based on the observed [CII] and far-IR luminosities which should be useful for comparing with results from future surveys of high-redshift galaxies with ALMA and CCAT.

The distribution of gas in the central kiloparsec of a galaxy has a dynamically rapid evolution. Nonaxisymmetries in the gravitational potential of the galactic disk, such as a large scale stellar bar or spiral, can lead to significant radial motion of gaseous material from larger radii to the central region. The large influx of gas and the subsequent star formation keep the central region constantly changing. However, the ability of gas to reach the nucleus proper to fuel an AGN phase is not guaranteed. Gas inflow can be halted at a circumnuclear star forming ring several hundred parsec away. The nearby galaxy NGC5248 is especially interesting in this sense since it is said to host 2 circumnuclear star forming rings at 100pc and 370pc from its quiescent nucleus. Here we present new subarcsecond PdBI+30m CO(2-1) emission line observations of the central region. For the first time the molecular gas distribution at the smallest stellar ring is resolved into a gas ring, consistent with the presence of a quiescent nucleus. However, the molecular gas shows no ring structure at the larger ring. We combine analyses of the gaseous and stellar content in the central kiloparsec of this galaxy to understand the gas distribution and dynamics of this star forming central region. We discuss the probability of two scenarios leading to the current observations, given our full understanding of this system, and discuss whether there are really two circumnuclear star forming rings in this galaxy.

Nuclear stellar cusps are defined as central excess light component in the stellar light profiles of galaxies and are suggested to be stellar relics of intense compact starbursts in the central ~100-500 pc region of gas-rich major mergers. Here we probe the build-up of nuclear cusps during the actual starburst phase for a complete sample of Luminous Infrared Galaxy systems (85 LIRGs, with 11.4<log[LIR/L_sun]<12.5) in the GOALS sample. Cusp properties are derived via 2-dimensional fitting of the nuclear stellar light imaged in the near-infrared by the Hubble Space Telescope and have been combined with mid-IR diagnostics for AGN/starburst characterization. We find that nuclear stellar cusps are resolved in 76% of LIRGs (merger and non-interacting galaxies). The cusp strength and luminosity increases with far-IR luminosity (excluding AGN) and merger stage, confirming theoretical models that starburst activity is associated with the build-up of nuclear stellar cusps. Evidence for ultra compact nuclear starbursts is found in ~13% of LIRGs, which have a strong unresolved central NIR light component but no significant contribution of an AGN. The nuclear near-IR surface density (measured within 1 kpc radius) increases by a factor of ~5 towards late merger stages. A careful comparison to local early-type galaxies with comparable masses reveals (a) that local (U)LIRGs have a significantly larger cusp fraction and (b) that the majority of the cusp LIRGs have host galaxy luminosities (H-band) similar to core ellipticals which is roughly one order in magnitude larger than for cusp ellipticals.

A {\it Hubble Space Telescope} ({\it HST}) / Advanced Camera for Surveys (ACS) study of the structural properties of 85 luminous and ultraluminous ($L_{\rm IR} > 10^{11.4}$ L$_\odot$) infrared galaxies (LIRGs and ULIRGs) in the Great Observatories All-sky LIRG Survey (GOALS) sample is presented. Two-dimensional GALFIT analysis has been performed on F814W "{\it I}-band" images to decompose each galaxy, as appropriate, into bulge, disk, central PSF and stellar bar components. The fraction of bulge-less disk systems is observed to be higher in LIRGs (35%) than in ULIRGs (20%), with the disk+bulge systems making up the dominant fraction of both LIRGs (55%) and ULIRGs (45%). Further, bulge+disk systems are the dominant late-stage merger galaxy type and are the dominant type for LIRGs and ULIRGs at almost every stage of galaxy-galaxy nuclear separation. The mean {\it I}-band host absolute magnitude of the GOALS galaxies is $-22.64\pm$0.62 mag (1.8$^{+1.4}_{-0.4}$ L$^*_I$), and the mean bulge absolute magnitude in GOALS galaxies is about 1.1 magnitude fainter than the mean host magnitude. Almost all ULIRGs have bulge magnitudes at the high end (-20.6 to - 23.5 mag) of the GOALS bulge magnitude range. Mass ratios in the GOALS binary systems are consistent with most of the galaxies being the result of major mergers, and an examination of the residual-to-host intensity ratios in GOALS binary systems suggests that smaller companions suffer more tidal distortion than the larger companions. We find approximately twice as many bars in GOALS disk+bulge systems (32.8%) than in pure-disk mergers (15.9%) but most of the disk+bulge systems that contain bars are disk-dominated with small bulges. The bar-to-host intensity ratio, bar half-light radius, and bar ellipticity in GOALS galaxies are similar to those found in nearby spiral galaxies.(abridged)

The Great Observatories All-Sky LIRG Survey (GOALS) is a multiwavelength study of luminous infrared galaxies (LIRGs) in the local universe. Here we present low resolution Spitzer spectra covering 5-38um and provide a basic analysis of the mid-IR spectral properties for nearby LIRGs. In a companion paper, we discuss detailed fits to the spectra. The GOALS sample of 244 nuclei in 180 luminous and 22 ultraluminous IR galaxies represents a complete subset of the IRAS RBGS and covers a range of merger stages, morphologies and spectral types. The majority (>60%) of GOALS LIRGs have high 6.2um PAH equivalent widths (EQW > 0.4um) and low levels of silicate absorption (s_9.7um >-1.0). There is a general trend among the U/LIRGs for silicate depth and MIR slope to increase with LIR. U/LIRGs in the late stages of a merger also have on average steeper MIR slopes and higher levels of dust obscuration. Together these trends suggest that as gas & dust is funneled towards the center of a coalescing merger, the nuclei become more compact and obscured. The sources that depart from these correlations have very low PAH EQW (EQW < 0.1um) consistent with their MIR emission being dominated by an AGN. The most heavily dust obscured sources are the most compact in their MIR emission, suggesting that the obscuring (cool) dust is associated with the outer regions of the starburst. As the merger progresses a marked decline is seen for the fraction of high EQW (star formation dominated) sources while the fraction of composite sources increases but the fraction of AGN-dominated sources remains low. When compared to the MIR spectra of submillimeter galaxies (SMGs) at z~2, the average GOALS LIRG is more absorbed at 9.7um and has more PAH emission. However, when the AGN contributions to both the local LIRGs and the high-z SMGs are removed, the average local starbursting LIRG closely resembles the starbursting SMGs.

Luminous and ultraluminous infrared galaxies ((U)LIRGs) are the most extreme star forming galaxies in the universe. The local (U)LIRGs provide a unique opportunity to study their multi-wavelength properties in detail for comparison to their more numerous counterparts at high redshifts. We present common large aperture photometry at radio through X-ray wavelengths, and spectral energy distributions (SEDs) for a sample of 53 nearby LIRGs and 11 ULIRGs spanning log (LIR/Lsun) = 11.14-12.57 from the flux-limited Great Observatories All-sky LIRG Survey (GOALS). The SEDs for all objects are similar in that they show a broad, thermal stellar peak and a dominant FIR thermal dust peak, where nuLnu(60um) / nuLnu(V) increases from ~2-30 with increasing LIR. When normalized at IRAS-60um, the largest range in the luminosity ratio, R(lambda)=log[nuLnu(lambda)/nuLnu(60um)] observed over the full sample is seen in the Hard X-rays (HX=2-10 keV). A small range is found in the Radio (1.4GHz), where the mean ratio is largest. Total infrared luminosities, LIR(8-1000um), dust temperatures, and dust masses were computed from fitting thermal dust emission modified blackbodies to the mid-infrared (MIR) through submillimeter SEDs. The new results reflect an overall ~0.02 dex lower luminosity than the original IRAS values. Total stellar masses were computed by fitting stellar population synthesis models to the observed near-infrared (NIR) through ultraviolet (UV) SEDs. Mean stellar masses are found to be log(M/Msun) = 10.79+/-0.40. Star formation rates have been determined from the infrared (SFR_IR~45Msun/yr) and from the monochromatic UV luminosities (SFR_UV~1.3Msun/yr), respectively. Multiwavelength AGN indicators have be used to select putative AGN: about 60% of the ULIRGs would have been classified as an AGN by at least one of the selection criteria.

Results of observations with the Spitzer, Hubble, GALEX, Chandra, and XMM-Newton space telescopes are presented for the Luminous Infrared Galaxy (LIRG) merger Mrk 266. The SW (Seyfert 2) and NE (LINER) nuclei reside in galaxies with Hubble types SBb (pec) and S0/a (pec), respectively. Both galaxies have L > L*, and they are inferred to each contain a ~2.5x10^8 M_sun black hole. Mrk 266 SW is likely the primary source of a bright Fe K-alpha line detected from the system, consistent with the reflection-dominated X-ray spectrum of a heavily obscured AGN. Optical knots embedded in an arc with aligned radio continuum radiation, combined with luminous H_2 line emission, provide evidence for a radiative bow shock in an AGN-driven outflow surrounding the NE nucleus. Soft X-ray emission modeled as shock-heated plasma is co-spatial with radio continuum emission between the galaxies. Mid-IR diagnostics indicate roughly equal contributions of AGN and starburst radiation powering the bolometric luminosity. Approximately 120 star clusters have been detected, with most having estimated ages < 50 Myr. Detection of 24 micron emission extending ~34 arcsec (20 kpc) north of the galaxies is interpreted as ~2x10^7 M_sun of dust entrained in an outflowing superwind; at optical wavelengths this region is resolved into a fragmented morphology indicative of Rayleigh-Taylor instabilities in an expanding shell of ionized gas. Mrk 266 demonstrates that the dust "blow-out" phase can begin in a LIRG well before the galaxies fully coalesce during a subsequent ULIRG phase, and rapid gas consumption in luminous dual AGNs (kpc scale separations) early in the merger process may explain the paucity of detected binary QSOs (sub-pc scale orbits) in large surveys. An evolutionary sequence is proposed representing a progression from dual to binary AGNs, accompanied by an increase in observed L_x/L_ir ratios by 10^4 or more.

We have mapped the key mid-IR diagnostics in eight major merger systems of the Toomre Sequence (NGC4676, NGC7592, NGC6621, NGC2623, NGC6240, NGC520, NGC3921, and NGC7252) using the Spitzer Infrared Spectrograph (IRS). With these maps, we explore the variation of the ionized-gas, PAH, and warm-gas (H_2) properties across the sequence and within the galaxies. While the global PAH interband strength and ionized gas flux ratios ([Ne III]/[Ne II]) are similar to those of normal star forming galaxies, the distribution of the spatially resolved PAH and fine structure line flux ratios is significant different from one system to the other. Rather than a constant H_2/PAH flux ratio, we find that the relation between the H_2 and PAH fluxes is characterized by a power law with a roughly constant exponent (0.61+/-0.05) over all merger components and spatial scales. While following the same power law on local scales, three galaxies have a factor of ten larger integrated (i.e. global) H_2/PAH flux ratio than the rest of the sample, even larger than what it is in most nearby AGNs. These findings suggest a common dominant excitation mechanism for H_2 emission over a large range of global H_2/PAH flux ratios in major mergers. Early merger systems show a different distribution between the cold (CO J=1-0) and warm (H_2) molecular gas component, which is likely due to the merger interaction. Strong evidence for buried star formation in the overlap region of the merging galaxies is found in two merger systems (NGC6621 and NGC7592) as seen in the PAH, [Ne II], [Ne III], and warm gas line emission, but with no apparent corresponding CO (J=1-0) emission. Our findings also demonstrate that the variations of the physical conditions within a merger are much larger than any systematic trends along the Toomre Sequence.

We present results from the second part of our analysis of the extended mid-infrared (MIR) emission of the Great Observatories All-Sky LIRG Survey (GOALS) sample based on 5-14 micron low-resolution spectra obtained with the IRS on Spitzer. We calculate the fraction of extended emission as a function of wavelength for all galaxies in the sample, FEE_lambda, and spatially separate the MIR spectrum of galaxies into their nuclear and extended components. We find that the [NeII] emission line is as compact as the hot dust MIR continuum, while the polycyclic aromatic hydrocarbon (PAH) emission is more extended. The 6.2 and 7.7 micron PAH emission is more compact than that of the 11.3 micron PAH, which is consistent with the formers being enhanced in a more ionized medium. The presence of an AGN or a powerful nuclear starburst increases the compactness of the hot dust MIR continuum, but has a negligible effect on the spatial extent of the PAH emission on kpc-scales. Globally, the spectra of the extended emission component are homogeneous for all galaxies in GOALS. This suggests that the physical properties of star formation taking place at distances farther than 1.5 kpc from the nuclei of (U)LIRGs are very similar, resembling local star-forming galaxies with L_IR < 10^11 Lsun, as well as star formation-dominated ULIRGs at z~2. In contrast, the MIR spectra of the nuclear component of local (U)LIRGs are very diverse. This implies that the observed variety of their integrated MIR properties arise, on average, only from the processes that are taking place in their cores.

We use the Expanded Very Large Array to image radio continuum emission from local luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) in 1 GHz windows centered at 4.7, 7.2, 29, and 36 GHz. This allows us to probe the integrated radio spectral energy distribution (SED) of the most energetic galaxies in the local universe. The 4-8 GHz flux densities agree well with previous measurements. They yield spectral indices \alpha \approx -0.67 (where F_\nu \propto \nu^\alpha) with \pm 0.15 (1\sigma) scatter, typical of nonthermal (synchrotron) emission from star-forming galaxies. The contrast of our 4-8 GHz data with literature 1.5 and 8.4 GHz flux densities gives further evidence for curvature of the radio SED of U/LIRGs. The SED appears flatter near \sim 1 GHz than near \sim 6 GHz, suggesting significant optical depth effects at the lower frequencies. The high frequency (28-37 GHz) flux densities are low compared to extrapolations from the 4-8 GHz data. We confirm and extend to higher frequency a previously observed deficit of high frequency radio emission for luminous starburst galaxies.

Within the NUclei of GAlaxies project we have obtained IRAM PdBI and 30m 12CO(1-0) and 12CO(2-1) observations of the spiral galaxy NGC 6951. Previous work shows that there is indirect evidence of gas inflow from 3 kpc down to small radii: a large-scale stellar bar, a prominent starburst ring (r~580 pc) and a LINER/Seyfert 2 nucleus. In this paper we study the gas kinematics as traced by the CO line emission in detail. We quantify the influence of the large-scale stellar bar by constructing an analytical model of the evolution of gas particles in a barred potential. From this model gravitational torques and mass accumulation rates are computed. We compare our model-based gravitational torque results with previous observationally-based ones. The model also shows that the large-scale stellar bar is indeed the dominant force for driving the gas inward, to the starburst ring. Inside the ring itself a nuclear stellar oval might play an important role. Detailed analysis of the CO gas kinematics there shows that emission arises from two co-spatial, but kinematically distinct components at several locations. The main emission component can always be related to the overall bar-driven gas kinematics. The second component exhibits velocities that are larger than expected for gas on stable orbits, has a molecular gas mass of 1.8x10^6Msun, is very likely connected to the nuclear stellar oval, and is consistent with inflowing motion towards the very center. This may form the last link in the chain of gas inflow towards the active galactic nucleus in NGC 6951.

We present results of Hubble Space Telescope NICMOS H-band imaging of 73 of most luminous (i.e., log[L_IR/L_0]>11.4) Infrared Galaxies (LIRGs) in the Great Observatories All-sky LIRG Survey (GOALS). This dataset combines multi-wavelength imaging and spectroscopic data from space (Spitzer, HST, GALEX, and Chandra) and ground-based telescopes. In this paper we use the high-resolution near-infrared data to recover nuclear structure that is obscured by dust at optical wavelengths and measure the evolution in this structure along the merger sequence. A large fraction of all galaxies in our sample possess double nuclei (~63%) or show evidence for triple nuclei (~6%). Half of these double nuclei are not visible in the HST B-band images due to dust obscuration. The majority of interacting LIRGs have remaining merger timescales of 0.3 to 1.3 Gyrs, based on the projected nuclear separations and the mass ratio of nuclei. We find that the bulge luminosity surface density increases significantly along the merger sequence (primarily due to a decrease of the bulge radius), while the bulge luminosity shows a small increase towards late merger stages. No significant increase of the bulge Sersic index is found. LIRGs that show no interaction features have on average a significantly larger bulge luminosity, suggesting that non merging LIRGs have larger bulge masses than merging LIRGs. This may be related to the flux limited nature of the sample and the fact that mergers can significantly boost the IR luminosity of otherwise low luminosity galaxies. We find that the projected nuclear separation is significantly smaller for ULIRGs (median value of 1.2 kpc) than for LIRGs (mean value of 6.7 kpc), suggesting that the LIRG phase appears earlier in mergers than the ULIRG phase.

We present a statistical analysis of the mid-infrared (MIR) spectra of 248 luminous infrared (IR) galaxies (LIRGs) which comprise the Great Observatories All-sky LIRG Survey (GOALS) observed with the Infrared Spectrograph (IRS) on-board the Spitzer Space Telescope. The GOALS sample enables a direct measurement of the relative contributions of star-formation and active galactic nuclei (AGN) to the total IR emission from a large sample of local LIRGs. The AGN contribution to the MIR emission (f-AGN) is estimated by employing several diagnostics based on the properties of the [NeV], [OIV] and [NeII] fine structure gas emission lines, the 6.2 microns PAH and the shape of the MIR continuum. We find that 18% of all LIRGs contain an AGN and that in 10% of all sources the AGN contributes more than 50% of the total IR luminosity. Summing up the total IR luminosity contributed by AGN in all our sources suggests that AGN supply ~12% of the total energy emitted by LIRGs. The average spectrum of sources with an AGN looks similar to the average spectrum of sources without an AGN, but it has lower PAH emission and a flatter MIR continuum. AGN dominated LIRGs have higher IR luminosities, warmer MIR colors and are found in interacting systems more often than pure starbursts LIRGs. However we find no linear correlations between these properties and f-AGN. We used the IRAC colors of LIRGs to confirm that finding AGN on the basis of their MIR colors may miss ~40% of AGN dominated (U)LIRGs

We present an analysis of the extended mid-infrared (MIR) emission of the Great Observatories All-Sky LIRG Survey (GOALS) sample based on 5-15um low resolution spectra obtained with the IRS on Spitzer. We calculate the fraction of extended emission as a function of wavelength for the galaxies in the sample, FEE_lambda. We can identify 3 general types of FEE_lambda: one where it is constant, one where features due to emission lines and PAHs appear more extended than the continuum, and a third which is characteristic of sources with deep silicate absorption at 9.7um. More than 30% of the galaxies have a median FEE_lambda larger than 0.5 implying that at least half of their MIR emission is extended. Luminous Infrared Galaxies (LIRGs) display a wide range of FEE in their warm dust continuum (0<=FEE_13.2um<=0.85). The large values of FEE_13.2um that we find in many LIRGs suggest that their extended MIR continuum emission originates in scales up to 10kpc. The mean size of the LIRG cores at 13.2um is 2.6kpc. However, once the LIR of the systems reaches the threshold of ~10^11.8Lsun, all sources become clearly more compact, with FEE_13.2um<=0.2, and their cores are unresolved. Our estimated upper limit for the core size of ULIRGs is less than 1.5kpc. The analysis indicates that the compactness of systems with LIR>~10^11.25Lsun strongly increases in those classified as mergers in their final stage of interaction. The FEE_13.2um is also related to the contribution of an active galactic nucleus (AGN) to the MIR. Galaxies which are more AGN-dominated are less extended, independently of their LIR. We finally find that the extent of the MIR continuum emission is correlated with the far-IR IRAS log(f_60um/f_100um) color. This enables us to place a lower limit to the area in a galaxy from where the cold dust emission may originate, a prediction which can be tested soon with the Herschel Space Telescope.

An analysis of data from the Spitzer Space Telescope, Hubble Space Telescope, Chandra X-ray Observatory, and AKARI Infrared Astronomy Satellite is presented for the z=0.036 merging galaxy system II Zw 096 (CGCG 448-020). Because II Zw 096 has an infrared luminosity of log(L_IR/L_sun) = 11.94, it is classified as a Luminous Infrared Galaxy (LIRG), and was observed as part of the Great Observatories All-sky LIRG Survey (GOALS). The Spitzer data suggest that 80% of the total infrared luminosity comes from an extremely compact, red source not associated with the nuclei of the merging galaxies. The Spitzer mid-infrared spectra indicate no high-ionization lines from a buried active galactic nucleus in this source. The strong detection of the 3.3 micron and 6.2 micron PAH emission features in the AKARI and Spitzer spectra also implies that the energy source of II Zw 096 is a starburst. Based on Spitzer infrared imaging and AKARI near-infrared spectroscopy, the star formation rate is estimated to be 120 M_sun/yr and > 45 M_sun/yr, respectively. Finally, the high-resolution B, I, and H-band images show many star clusters in the interacting system. The colors of these clusters suggest at least two populations - one with an age of 1-5 Myr and one with an age of 20-500 Myr, reddened by 0-2 magnitudes of visual extinction. The masses of these clusters span a range between 10^6-10^8 M_sun. This starburst source is reminiscent of the extra-nuclear starburst seen in NGC 4038/9 (the Antennae Galaxies) and Arp 299 but approximately an order of magnitude more luminous than the Antennae. The source is remarkable in that the off-nuclear infrared luminosity dominates the enitre system.

The Great Observatories All-sky LIRG Survey (GOALS) consists of a complete sample of 202 Luminous Infrared Galaxies (LIRGs) selected from the IRAS Revised Bright Galaxy Sample (RBGS). The galaxies span the full range of interaction stages, from isolated galaxies to interacting pairs to late stage mergers. We present a comparison of the UV and infrared properties of 135 galaxies in GOALS observed by GALEX and Spitzer. For interacting galaxies with separations greater than the resolution of GALEX and Spitzer (2-6"), we assess the UV and IR properties of each galaxy individually. The contribution of the FUV to the measured SFR ranges from 0.2% to 17.9%, with a median of 2.8% and a mean of 4.0 +/- 0.4%. The specific star formation rate of the GOALS sample is extremely high, with a median value (3.9*10^{-10} yr^{-1}) that is comparable to the highest specific star formation rates seen in the Spitzer Infrared Nearby Galaxies Survey sample. We examine the position of each galaxy on the IR excess-UV slope (IRX-beta) diagram as a function of galaxy properties, including IR luminosity and interaction stage. The LIRGs on average have greater IR excesses than would be expected based on their UV colors if they obeyed the same relations as starbursts with L_IR < 10^{11}L_0 or normal late-type galaxies. The ratio of L_IR to the value one would estimate from the IRXg-beta relation published for lower luminosity starburst galaxies ranges from 0.2 to 68, with a median value of 2.7. A minimum of 19% of the total IR luminosity in the RBGS is produced in LIRGs and ULIRGs with red UV colors (beta > 0). Among resolved interacting systems, 32% contain one galaxy which dominates the IR emission while the companion dominates the UV emission. Only 21% of the resolved systems contain a single galaxy which dominates both wavelengths.

The Great Observatories All-sky LIRG Survey (GOALS) combines data from NASA's Spitzer, Chandra, Hubble and GALEX observatories, together with ground-based data into a comprehensive imaging and spectroscopic survey of over 200 low redshift Luminous Infrared Galaxies (LIRGs). The LIRGs are a complete subset of the IRAS Revised Bright Galaxy Sample (RBGS). The LIRGs targeted in GOALS span the full range of nuclear spectral types defined via traditional optical line-ratio diagrams as well as interaction stages. They provide an unbiased picture of the processes responsible for enhanced infrared emission in galaxies in the local Universe. As an example of the analytic power of the multi-wavelength GOALS dataset, we present data for the interacting system VV 340 (IRAS F14547+2449). Between 80-95% of the total far-infrared emission (or about 5E11 solar luminosities) originates in VV 340 North. While the IRAC colors of VV 340 North and South are consistent with star-forming galaxies, both the Spitzer IRS and Chandra ACIS data indicate the presence of a buried AGN in VV 340 North. The GALEX far and near-UV fluxes imply a extremely large infrared "excess" (IRX) for the system (IR/FUV = 81) which is well above the correlation seen in starburst galaxies. Most of this excess is driven by VV 340 N, which alone has an IR excess of nearly 400. The VV 340 system seems to be comprised of two very different galaxies - an infrared luminous edge-on galaxy (VV 340 North) that dominates the long-wavelength emission from the system and which hosts a buried AGN, and a face-on starburst (VV 340 South) that dominates the short-wavelength emission.

To examine the role of the host galaxy structure in fueling nuclear activity, we estimated gas flow rates from several kpc down to the inner few 10 pc for seven nearby spiral galaxies, selected from the NUGA sample (NUclei of GAlaxies). We calculated gravitational torques from near-IR images and determined gas in/out-flow rates as a function of radius and location within the galactic disks, based on high angular resolution interferometric observations of molecular (CO using PdBI) and atomic (HI using the VLA) gas. The results are compared with kinematic evidence for radial gas flows and the dynamical state of the galaxies (via resonances) derived from several different methods. We show that gravitational torques are very efficient at transporting gas from the outer disk all the way into the galaxies centers at ~100 pc; previously assumed dynamical barriers to gas transport, such as the Corotation Resonance of stellar bars, seem to be overcome by gravitational torque induced gas flows from other non-axisymmmetric structures. The resulting rates of gas mass inflow range from 0.01 to 50 solar masses per year and are larger for the galaxy center than for the outer disk. Our gas flow maps show the action of nested bars within larger bars for 3 galaxies. Non-circular streaming motions found in the kinematic maps are larger in the center than in the outer disk and appear to correlate only loosely with the in/out-flow rates as a function of radius. We demonstrate that spiral gas disks are very dynamic systems that undergo strong radial evolution on timescales of a few rotation periods (e.g. 5 times 10^8 yrs at a radius of 5 kpc), due to the effectiveness of gravitational torques in redistributing the cold galactic gas.

We create a complete gravity torque map of the disk of the LINER/Seyfert 1.9 galaxy NGC4579. We quantify the efficiency of angular momentum transport and search for signatures of secular evolution in the fueling process from r~15kpc down to the inner r~50pc around the Active Galactic Nucleus (AGN). We use both the 1-0 and 2-1 line maps of CO obtained with the Plateau de Bure Interferometer (PdBI) as part of the NUclei of Galaxies-(NUGA)-project. We derive the stellar potential from a NIR (K band) wide field image of the galaxy. The K-band image, which reveals a stellar bar, together with a high resolution HI map of NGC4579 obtained with the Very Large Array (VLA), allow us to extend the gravity torque analysis to the outer disk. The bulk of the gas response traced by the CO PdBI maps follows the expected gas flow pattern induced by the bar potential in the presence of two Inner Lindblad Resonances (ILR). We also detect an oval distortion in the inner r~200pc of the K-band image. The oval is not aligned with the large-scale bar, a signature of dynamical decoupling. The morphology of the outer disk suggests that the neutral gas is currently piling up in a pseudo-ring formed by two winding spiral arms that are morphologically decoupled from the bar structure. In the outer disk, the decoupling of the spiral allows the gas to efficiently produce net gas inflow on intermediate scales. The corotation barrier seems to be overcome due to secular evolution processes. The gas in the inner disk is efficiently funneled by gravity torques down to r~300pc. Closer to the AGN, the two m=2 modes (bar and oval) act in concert to produce net gas inflow down to r~50pc, providing a clear smoking gun evidence of fueling with associated short dynamical time-scales.

We present a comprehensive spectroscopic imaging survey of the distribution and kinematics of atomic hydrogen (HI) in 16 nearby spiral galaxies hosting low luminosity AGN, observed with high spectral and spatial resolution (resolution: ~20 arcsec, 5 km/s) using the NRAO Very Large Array (VLA). The sample contains a range of nuclear types, ranging from Seyfert to star-forming nuclei and was originally selected for the NUclei of GAlaxies project (NUGA) - a spectrally and spatially resolved interferometric survey of gas dynamics in nearby galaxies designed to identify the fueling mechanisms of AGN and the relation to host galaxy evolution. Here we investigate the relationship between the HI properties of these galaxies, their environment, their stellar distribution and their AGN type. The large-scale HI morphology of each galaxy is classified as ringed, spiral, or centrally concentrated; comparison of the resulting morphological classification with AGN type reveals that ring structures are significantly more common in LINER than in Seyfert host galaxies, suggesting a time evolution of the AGN activity together with the redistribution of the neutral gas. Dynamically disturbed HI disks are also more prevalent in LINER host galaxies than in Seyfert host galaxies. While several galaxies are surrounded by companions (some with associated HI emission), there is no correlation between the presence of companions and the AGN type (Seyfert/LINER).

We are conducting the first systematic 3D spectroscopic imaging survey to quantify the properties of the atomic gas (HI) in a distance-limited sample of 28 Seyfert galaxies and a sample of 28 inactive control galaxies with well-matched optical properties (the VHIKINGS survey). This study aims to address the role of the host galaxy in nuclear activity and confront outstanding controversies in optical/IR imaging surveys. Early results show possible relationships between Seyfert activity and HI extent, content and the prevalence of small, nearby gas-rich dwarf galaxies (M(HI)~10^7 Msun); results will be tested via rigorous comparison with control galaxies. Initial results from our optical followup study of 15 of our galaxies using the SAURON integral field unit on the WHT suggest a possible difference between Seyfert and inactive stellar and gaseous kinematics that support the conclusion that internal kinematics of galaxies are the key to nuclear activity.