We present the results of a {\it Hubble Space Telescope} ACS/HRC FUV, ACS/WFC optical study into the cluster populations of a sample of 22 Luminous Infrared Galaxies in the Great Observatories All-Sky LIRG Survey. Through integrated broadband photometry we have derived ages and masses for a total of 484 star clusters contained within these systems. This allows us to examine the properties of star clusters found in the extreme environments of LIRGs relative to lower luminosity star-forming galaxies in the local Universe. We find that by adopting a Bruzual \& Charlot simple stellar population (SSP) model and Salpeter initial mass function, the age distribution of clusters declines as $dN/d\tau = \tau^{-0.9 +/- 0.3}$, consistent with the age distribution derived for the Antennae Galaxies, and interpreted as evidence for rapid cluster disruption occuring in the strong tidal fields of merging galaxies. The large number of $10^{6} M_{\odot}$ young clusters identified in the sample also suggests that LIRGs are capable of producing more high-mass clusters than what is observed to date in any lower luminosity star-forming galaxy in the local Universe. The observed cluster mass distribution of $dN/dM = M^{-1.95 +/- 0.11}$ is consistent with the canonical -2 power law used to describe the underlying initial cluster mass function (ICMF) for a wide range of galactic environments. We interpret this as evidence against mass-dependent cluster disruption, which would flatten the observed CMF relative to the underlying ICMF distribution.

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.

New optical HST, Spitzer, GALEX, and Chandra observations of the single-nucleus, luminous infrared galaxy (LIRG) merger IC 883 are presented. The galaxy is a member of the Great Observatories All-sky LIRG Survey (GOALS), and is of particular interest for a detailed examination of a luminous late-stage merger due to the richness of the optically-visible star clusters and the extended nature of the nuclear X-ray, mid-IR, CO and radio emission. In the HST ACS images, the galaxy is shown to contain 156 optically visible star clusters distributed throughout the nuclear regions and tidal tails of the merger, with a majority of visible clusters residing in an arc ~ 3-7 kpc from the position of the mid-infrared core of the galaxy. The luminosity functions of the clusters have an alpha_F435W ~ -2.17+/-0.22 and alpha_F814W ~ -2.01+/-0.21. Further, the colors and absolute magnitudes of the majority of the clusters are consistent with instantaneous burst population synthesis model ages in the range of a few x10^7 - 10^8 yrs (for 10^5 M_sun clusters), but may be as low as few x10^6 yrs with extinction factored in. The X-ray and mid-IR spectroscopy are indicative of predominantly starburst-produced nuclear emission, and the star formation rate is ~ 80 M_sun / yr. The kinematics of the CO emission and the morphology of both the CO and radio emission are consistent with the nuclear starburst being situated in a highly inclined disk 2 kpc in diameter with an infrared surface brightness mu_IR ~ 2x10^11 L_sun kpc^-2, a factor of 10 less than that of the Orion star-forming region. Finally, the detection of the [Ne V] 14.32 um emission line is evidence that an AGN is present. The faintness of the line (i.e., [Ne V] / [Ne II] um ~ 0.01) and the small equivalent width of the 6.2 um PAH feature ($= 0.39\mu$m) are both indicative of a relatively weak AGN. (abridged)

We present X-ray data for a complete sample of 44 luminous infrared galaxies (LIRGs), obtained with the Chandra X-ray Observatory. These are the X-ray observations of the high luminosity portion of the Great Observatory All-sky LIRG Survey (GOALS), which includes the most luminous infrared selected galaxies, log (Lir/Lsun) > 11.73, in the local universe, z < 0.088. X-rays were detected from 43 out of 44 objects, and their arcsec-resolution images, spectra, and radial brightness distributions are presented. With a selection by hard X-ray colour and the 6.4 keV iron line, AGN are found in 37% of the objects, with higher luminosity sources more likely to contain an AGN. These AGN also tend to be found in late-stage mergers. The AGN fraction would increase to 48% if objects with mid-IR [Ne V] detection are included. Double AGN are clearly detected only in NGC 6240 among 24 double/triple systems. Other AGN are found either in single nucleus objects or in one of the double nuclei at similar rates. Objects without conventional X-ray signatures of AGN appear to be hard X-ray quiet, relative to the X-ray to far-IR correlation for starburst galaxies, as discussed elsewhere. Most objects also show extended soft X-ray emission, which is likely related to an outflow from the nuclear region, with a metal abundance pattern suggesting enrichment by core collapse supernovae, as expected for a starburst.

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

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.

We present a study of the nearby Seyfert galaxy NGC 1068 using mid- and far- infrared data acquired with the IRAC, IRS, and MIPS instruments aboard the Spitzer Space Telescope. The images show extensive 8 um and 24 um emission coinciding with star formation in the inner spiral approximately 15" (1 kpc) from the nucleus, and a bright complex of star formation 47" (3 kpc) SW of the nucleus. The brightest 8 um PAH emission regions coincide remarkably well with knots observed in an Halpha image. Strong PAH features at 6.2, 7.7, 8.6, and 11.3 um are detected in IRS spectra measured at numerous locations inside, within, and outside the inner spiral. The IRAC colors and IRS spectra of these regions rule out dust heated by the AGN as the primary emission source; the SEDs are dominated by starlight and PAH emission. The equivalent widths and flux ratios of the PAH features in the inner spiral are generally consistent with conditions in a typical spiral galaxy ISM. Interior to the inner spiral, the influence of the AGN on the ISM is evident via PAH flux ratios indicative of a higher ionization parameter and a significantly smaller mean equivalent width than observed in the inner spiral. The brightest 8 and 24 um emission peaks in the disk of the galaxy, even at distances beyond the inner spiral, are located within the ionization cones traced by [O III]/Hbeta, and they are also remarkably well aligned with the axis of the radio jets. Although it is possible that radiation from the AGN may directly enhance PAH excitation or trigger the formation of OB stars that subsequently excite PAH emission at these locations in the inner spiral, the orientation of collimated radiation from the AGN and star formation knots in the inner spiral could be coincidental. (abridged)

New millimeter-wave CO and HCN observations of the host galaxies of infrared-excess Palomar Green quasi-stellar objects (PG QSOs) previously detected in CO are presented. These observations are designed to assess the validity of using the infrared luminosity to estimate star formation rates of luminous AGN by determining the relative significance of dust-heating by young, massive stars and active galactic nuclei (AGN) in QSO hosts and IRAS galaxies with warm, AGN-like infrared colors. The HCN data show the PG QSO host IZw1 and most of the warm IRAS galaxies to have high L_IR / L'_HCN (>1600) relative to the cool IRAS galaxy population for which the median L_IR / L'_HCN ~ 890(+440,-470). If the assumption is made that the infrared emission from cool IRAS galaxies is reprocessed light from embedded star-forming regions, then high values of L_IR / L'_HCN are likely the result of dust heating by the AGN. Further, if the median ratio of L'_HCN / L'_CO ~ 0.06 observed for Seyfert galaxies and IZw1 is applied to the PG QSOs not detected in HCN, then the derived L_IR / L'_HCN correspond to a stellar contribution to the production of L_IR of ~ 7-39%, and star formation rates ~ 2-37 M_sun/yr are derived for the QSO hosts. Alternatively, if the far-infrared is adopted as the star formation component of the total infrared in cool galaxies, the stellar contributions in QSO hosts to their L_FIR are up to 35% higher than the percentages derived for L_IR. This raises the possibility that the L_FIR in several of the PG QSO hosts, including IZw1, could be due entirely to dust heated by young, massive stars. Finally, there is no evidence that the global HCN emission is enhanced relative to CO in galaxies hosting luminous AGN.