About half of nearby galaxies have a central surface brightness >1 magnitude below that of the sky. The overall properties of these low-surface-brightness galaxies (LSBGs) remain understudied, and in particular we know very little about their massive black hole population. This gap must be closed to determine the frequency of massive black holes at z=0 as well as to understand their role in regulating galaxy evolution. Here we investigate the incidence and intensity of nuclear, accretion-powered X-ray emission in a sample of 32 nearby LSBGs with the Chandra X-ray Observatory. A nuclear X-ray source is detected in 4 galaxies (12.5%). Based on an X-ray binary contamination assessment technique developed for normal galaxies, we conclude that the detected X-ray nuclei indicate low-level accretion from massive black holes. The active fraction is consistent with that expected from the stellar mass distribution of the LSBGs, but not their total baryonic mass, when using a scaling relation from an unbiased X-ray survey of normal galaxies. This suggests that their black holes co-evolved with their stellar population. In addition, the apparent agreement nearly doubles the number of galaxies available within ~100 Mpc for which a measurement of nuclear activity can efficiently constrain the frequency of black holes as a function of stellar mass. We conclude by discussing the feasibility of measuring this occupation fraction to a few percent precision below 1e10 solar masses with high-resolution, wide-field X-ray missions currently under consideration.

We present high resolution spectroscopy taken with the Keck Echellete Spectrograph and Imager to measure stellar velocity dispersions for eight active dwarf galaxies ($M_{\ast}<3\times10^{9}~M_{\odot}$) with virial black hole masses. We double the number of systems in this stellar mass regime with measurements of both black hole mass ($M_{\rm BH}$) and stellar velocity dispersion ($\sigma_{\ast}$), and place them on the $M_{\rm BH}-\sigma_{\ast}$ relation. The tight relation between $M_{\rm BH}$ and $\sigma_{\ast}$ for higher mass galaxies is a strong piece of evidence for the co-evolution of BHs and their host galaxies, but it has been unclear whether this relation holds in the dwarf galaxy regime. Our sample is in good agreement with the extrapolation of the $M_{\rm BH}-\sigma_{\ast}$ relation to low BH/galaxy masses, suggesting that the processes which produce $M_{\rm BH}-\sigma_{\ast}$ can also operate in dwarf galaxies. These results provide important constraints for massive black hole seed formation models and models exploring the impact of AGN feedback in dwarf galaxies.

As part of an on-going study of radio transients in Epoch 1 (2017-2019) of the Very Large Array Sky Survey (VLASS), we have discovered a sample of 0.2<z<3.2 active galactic nuclei (AGN) selected in the optical/infrared that have recently brightened dramatically in the radio. These sources would have previously been classified as radio-quiet based on upper limits from the Faint Images of the Radio Sky at Twenty-centimeters (FIRST; 1993-2011) survey; however, they are now consistent with radio-loud quasars. We present a quasi-simultaneous, multi-band (1-18 GHz) VLA follow-up campaign of our sample of AGN with extreme radio variability. We conclude that the radio properties are most consistent with AGN that have recently launched jets within the past few decades, potentially making them among the youngest radio AGN known.

We report on small-amplitude optical variability and recent dissipation of the unusually persistent broad emission lines in the blue compact dwarf galaxy PHL 293B. The galaxy's unusual spectral features (P Cygni-like profiles with $\sim$800 km s$^{-1}$ blueshifted absorption lines) have resulted in conflicting interpretations of the nature of this source in the literature. However, analysis of new Gemini spectroscopy reveals the broad emission has begun to fade after being persistent for over a decade prior. Precise difference imaging light curves constructed with the Sloan Digital Sky Survey and the Dark Energy Survey reveal small-amplitude optical variability of $\sim$0.1 mag in the g band offset by $100\pm21$ pc from the brightest pixel of the host. The light curve is well-described by an active galactic nuclei (AGN)-like damped random walk process. However, we conclude that the origin of the optical variability and spectral features of PHL 293B is due to a long-lived stellar transient, likely a Type IIn supernova or non-terminal outburst, mimicking long-term AGN-like variability. This work highlights the challenges of discriminating between scenarios in such extreme environments, relevant to searches for AGNs in dwarf galaxies. This is the second long-lived transient discovered in a blue compact dwarf, after SDSS1133. Our result implies such long-lived stellar transients may be more common in metal-deficient galaxies. Systematic searches for low-level variability in dwarf galaxies will be possible with the upcoming Legacy Survey of Space and Time at Vera C. Rubin Observatory.

We present an analysis of the long-term optical variability for $\sim50,000$ nearby (z<0.055) galaxies from the NASA-Sloan Atlas, $35,000$ of which are low-mass ($M_{\ast}<10^{10}~M_{\odot}$). We use difference imaging of Palomar Transient Factory (PTF) R-band observations to construct light curves with typical baselines of several years. We then search for subtle variations in the nuclear light output. We determine whether detected variability is AGN-like by assessing the fit quality to a damped random walk model. We identify 424 variability-selected AGN, including 244 with stellar masses between $10^{7}$ and $10^{10}~M_{\odot}$. 75% of low-mass galaxies with AGN-like variability have narrow emission lines dominated by star formation. After controlling for nucleus magnitude, the fraction of variable AGN is constant down to $M_{\ast}=10^{9}~M_{\odot}$, suggesting no drastic decline in the BH occupation fraction down to this stellar mass regime. Combining our NASA-Sloan Atlas sample with samples of nearby galaxies with broad H$\alpha$ emission, we find no dependence of variability properties with black hole mass. However, we caution that the variable AGN fraction is strongly dependent on baseline. For baselines less than two years, the variable fraction for the full sample is 0.25%, compared to 1.0% for baselines longer than two years. Finally, comparing Stripe 82 light curves (Baldassare et al. 2018) to PTF light curves, we find populations of changing-look AGN: 8 galaxies that are variable in Stripe 82, but quiescent in PTF, and 15 galaxies where the reverse is true. Our PTF work demonstrates the promise of long-term optical variability searches in low-mass galaxies for finding AGNs missed by other selection techniques.

Active galactic nuclei (AGNs) show a correlation between the size of the broad line region (BLR) and the monochromatic continuum luminosity at 5100 \AA, allowing black hole mass estimation based on single-epoch spectra. However, the validity of the correlation is yet to be clearly tested for high-luminosity AGNs. We present the first reverberation-mapping results of the Seoul National University AGN monitoring program (SAMP), which is designed to focus on luminous AGNs for probing the high end of the size-luminosity relation. We report time lag measurements of two AGNs, namely, 2MASS J10261389+5237510 and SDSS J161911.24+501109.2, using the light curves obtained over a $\sim$1000 day period with an average cadence of $\sim$10 and $\sim$20 days, respectively for photometry and spectroscopy monitoring. Based on a cross-correlation analysis and H$\beta$ line width measurements, we determine the H$\beta$ lag as $41.8^{+4.9}_{-6.0}$ and $52.6^{+17.6}_{-14.7}$ days in the observed-frame, and black hole mass as $3.65^{+0.49}_{-0.57} \times 10^7 M_{\odot}$ and $23.02^{+7.81}_{-6.56} \times 10^7 M_{\odot}$, respectively for 2MASS J1026 and SDSS J1619.

The origin and properties of black hole seeds that grow to produce the detected population of supermassive black holes are unconstrained at present. Despite the existence of several potentially feasible channels for the production of initial seeds in the high redshift universe, since even actively growing seeds are not directly observable at these epochs, discriminating between models remains challenging. Several new observables that encapsulate information about seeding have been proposed in recent years, and these offer exciting prospects for truly unraveling the nature of black hole seeds in the coming years. One of the key challenges for this task lies in the complexity of the problem, the required disentangling of the confounding effects of accretion physics and mergers, as mergers and accretion events over cosmic time stand to erase these initial conditions. Nevertheless, some unique signatures of seeding do survive and still exist in: local scaling relations between black holes and their galaxy hosts at low-masses; in high-redshift luminosity functions of accreting black holes; and in the total number and mass functions of gravitational wave coalescence events from merging binary black holes. One of the clearest discriminants for seed models are these high redshift gravitational wave detections of mergers from space detectable in the milliHertz range. These predicted event rates offer the most direct constraints on the properties of initial black hole seeds. Improving our theoretical understanding of black hole dynamics and accretion will also be pivotal in constraining seeding models in combination with the wide range of multi-messenger data.

In this white paper we explore the capabilities required to identify and study supermassive black holes formed from heavy seeds ($\mathrm{M_{\bullet}} \sim 10^4 - 10^6 \, \mathrm{M_{\odot}}$) in the early Universe. To obtain an unequivocal detection of heavy seeds we need to probe mass scales of $\sim 10^{5-6} \, \mathrm{M_{\odot}}$ at redshift $z \gtrsim 10$. From this theoretical perspective, we review the observational requirements and how they compare with planned/proposed instruments, in the infrared, X-ray and gravitational waves realms. In conclusion, detecting heavy black hole seeds at $z \gtrsim 10$ in the next decade will be challenging but, according to current theoretical models, feasible with upcoming/proposed facilities. Their detection will be fundamental to understand the early history of the Universe, as well as its evolution until now. Shedding light on the dawn of black holes will certainly be one of the key tasks that the astronomical community will focus on in the next decade.

This document illustrates the feasibility of a few per cent level measurement of the local black hole occupation fraction in low mass galaxies through wide-field, high angular resolution X-ray imaging observations of local volume galaxies. The occupation fraction, particularly at the low end of the galaxy luminosity function, is a key benchmark for any model which aims to reproduce the formation and growth of super-massive black holes and their host galaxies. Our proposed measurement will complement orthogonal efforts that are planned in X-rays at high red-shifts, as well as in the local Universe with ground-based facilities.

We characterize the incidence and intensity of low-level super-massive black hole activity within the Fornax cluster, through X-ray observations of the nuclei of 29 quiescent early-type galaxies. Using the \textit{Chandra X-ray Telescope}, we target 17 galaxies from the HST Fornax Cluster Survey, down to a uniform (3$\sigma$) limiting X-ray luminosity threshold of $5\cdot10^{38}$ergs$^{-1}$, which we combine with deeper, archival observations for an additional 12 galaxies. A nuclear X-ray point-source is detected in 11 out of 29 targets. After accounting for the low mass X-ray binary contamination to the nuclear X-ray signal, the X-ray active fraction is measured at $26.6\% \pm 9.6\%$. The results from this analysis are compared to similar investigations targeting quiescent early types in the Virgo cluster, as well as the field. After correcting for the different mass distributions, the measured Fornax active fraction is less than the field fraction, at more than 3$\sigma$, confirming that the funneling of gas to the nuclear regions of cluster members is inhibited compared to those galaxies in the field. At the same time, we find no statistically significant difference between Fornax and Virgo galaxies, with only marginal evidence for a lower active fraction in Fornax (1 $\sigma$); if real, owing to Fornax's higher galaxy number density, this could indicate that galaxy-galaxy interactions are more effective at gas removal than galaxy-gas effects.

We present an analysis of the nuclear variability of $\sim28,000$ nearby ($z<0.15$) galaxies with Sloan Digital Sky Survey (SDSS) spectroscopy in Stripe 82. We construct light curves using difference imaging of SDSS g-band images, which allows us to detect subtle variations in the central light output. We select variable AGN by assessing whether detected variability is well-described by a damped random walk model. We find 135 galaxies with AGN-like nuclear variability. While most of the variability-selected AGNs have narrow emission lines consistent with the presence of an AGN, a small fraction have narrow emission lines dominated by star formation. The star-forming systems with nuclear AGN-like variability tend to be low-mass ($M_{\ast}<10^{10}~M_{\odot}$), and may be AGNs missed by other selection techniques due to star formation dilution or low-metallicities. We explore the AGN fraction as a function of stellar mass, and find that the fraction of variable AGN increases with stellar mass, even after taking into account the fact that lower mass systems are fainter. There are several possible explanations for an observed decline in the fraction of variable AGN with decreasing stellar mass, including a drop in the supermassive black hole occupation fraction, a decrease in the ratio of black hole mass to galaxy stellar mass, or a change in the variability properties of lower-mass AGNs. We demonstrate that optical photometric variability is a promising avenue for detecting AGNs in low-mass, star formation-dominated galaxies, which has implications for the upcoming Large Synoptic Survey Telescope.

RGG 118 (SDSS 1523+1145) is a nearby ($z=0.0243$), dwarf disk galaxy ($M_{\ast}\approx2\times10^{9} M_{\odot}$) found to host an active $\sim50,000$ solar mass black hole at its core (Baldassare et al. 2015). RGG 118 is one of a growing collective sample of dwarf galaxies known to contain active galactic nuclei -- a group which, until recently, contained only a handful of objects. Here, we report on new \textit{Hubble Space Telescope} Wide Field Camera 3 UVIS and IR imaging of RGG 118, with the main goal of analyzing its structure. Using 2-D parametric modeling, we find that the morphology of RGG 118 is best described by an outer spiral disk, inner component consistent with a pseudobulge, and central PSF. The luminosity of the PSF is consistent with the central point source being dominated by the AGN. We measure the luminosity and mass of the "pseudobulge" and confirm that the central black hole in RGG 118 is under-massive with respect to the $M_{BH}-M_{\rm bulge}$ and $M_{BH}-L_{\rm bulge}$ relations. This result is consistent with a picture in which black holes in disk-dominated galaxies grow primarily through secular processes.

Motivated by theoretical expectations that Nuclear Star Clusters (NSCs) in galactic centers may provide a favorable environment for super-massive black holes to form and/or efficiently grow, we set out to measure the fraction of nearby nucleated galaxies that also host an Active Galactic Nucleus (AGN). We targeted a distance-limited sample of 98 objects with the Chandra X-ray Telescope, down to a uniform X-ray luminosity threshold of $\sim$10$^{38}$ erg s$^{-1}$. The sample is composed of 47 late-types and 51 early-types, enabling us to further investigate the active fraction as a function of galactic morphology. After correcting for contamination to the nuclear X-ray signal from bright X-ray binaries, we measure an active fraction $f$=11.2$\%^{+7.4}_{-4.9}$ (1$\sigma$ C.L.) across the whole sample, in agreement with previous estimates based on an heterogeneous combination of optical, X-ray and radio diagnostics, by Seth et al. (2008). After accounting for the different stellar mass distributions in our samples, we find no statistically significant difference in the active fraction of early- vs. late-type nucleated galaxies, with $f$=10.6$\%^{+11.9}_{-4.9}$ and 10.8$\%^{+11.3}_{-6.3}$, respectively. For the early-type nucleated galaxies, we are able to carry out a controlled comparison with a parent sample of non-nucleated galaxies covering the same stellar mass range, finding again no statistically significant difference in the active fraction. Taken at face value, our findings suggest that the presence of a NSC does not facilitate nor enhance accretion-powered emission from a nuclear super-massive black hole. This is true even for late-type nucleated galaxies, home to bluer NSCs and arguably larger gas reservoirs.

We present new Chandra X-ray Observatory and Hubble Space Telescope observations of eight optically selected broad-line AGN candidates in nearby dwarf galaxies ($z<0.055$). Including archival Chandra observations of three additional sources, our sample contains all ten galaxies from Reines et al. (2013) with both broad H$\alpha$ emission and narrow-line AGN ratios (6 AGNs, 4 Composites), as well as one low-metallicity dwarf galaxy with broad H$\alpha$ and narrow-line ratios characteristic of star formation. All eleven galaxies are detected in X-rays. Nuclear X-ray luminosities range from $L_{0.5-7 \rm{keV}}\approx5\times10^{39}$ to $1\times10^{42}$ $\rm{erg}\rm{s^{-1}}$. In all cases except for the star forming galaxy, the nuclear X-ray luminosities are significantly higher than would be expected from X-ray binaries, providing strong confirmation that AGN and composite dwarf galaxies do indeed host actively accreting BHs. Using our estimated BH masses (which range from $\sim7\times10^{4}-1\times10^{6}~M_{\odot}$), we find inferred Eddington fractions ranging from $\sim0.1-50\%$, i.e. comparable to massive broad-line quasars at higher redshift. We use the HST imaging to determine the ratio of ultraviolet to X-ray emission for these AGN, finding that they appear to be less X-ray luminous with respect to their UV emission than more massive quasars (i.e. $\alpha_{\rm OX}$ values an average of 0.36 lower than expected based on the relation between $\alpha_{\rm OX}$ and $2500{\rm \AA}$ luminosity). Finally, we discuss our results in the context of different accretion models onto nuclear BHs.

We use time-domain optical spectroscopy to distinguish between broad emission lines powered by accreting black holes (BHs) or stellar processes (i.e., supernovae) for 16 galaxies identified as AGN candidates by Reines \etal (2013). Our study is primarily focused on those objects with narrow emission-line ratios dominated by star formation. Based on follow-up spectra taken with the Magellan Echellette Spectrograph (MagE), the Dual Imaging Spectrograph, and the Ohio State Multi-Object Spectrograph, we find that the broad H$\alpha$ emission has faded or was ambiguous for all of the star-forming objects (14/16) over baselines ranging from 5 to 14 years. For the two objects in our follow-up sample with narrow-line AGN signatures (RGG 9 and RGG 119), we find persistent broad H$\alpha$ emission consistent with an AGN origin. Additionally, we use our MagE observations to measure stellar velocity dispersions for 15 objects in the Reines et al. (2013) sample, all with narrow-line ratios indicating the presence of an AGN. Stellar masses range from $\sim5\times10^{8}$ to $3\times10^{9}$~\msun, and we measure $\sigma_{\ast}$ ranging from $28-71~{\rm km~s^{-1}}$. These $\sigma_{\ast}$ correspond to some of the lowest-mass galaxies with optical signatures of AGN activity. We show that RGG 119, the one object which has both a measured $\sigma_{\ast}$ and persistent broad H$\alpha$ emission, falls near the extrapolation of the $\rm M_{BH}-\sigma_{\star}$ relation to the low-mass end.

We examine the hypothesis that the red near-infrared colors of some L dwarfs could be explained by a "dust haze" of small particles in their upper atmospheres. This dust haze would exist in conjunction with the clouds found in dwarfs with more typical colors. We developed a model which uses Mie theory and the Hansen particle size distributions to reproduce the extinction due to the proposed dust haze. We apply our method to 23 young L dwarfs and 23 red field L dwarfs. We constrain the properties of the dust haze including particle size distribution and column density using Markov-Chain Monte Carlo methods. We find that sub-micron range silicate grains reproduce the observed reddening. Current brown dwarf atmosphere models include large grain (1--100~$\mu m$) dust clouds but not submicron dust grains. Our results provide a strong proof of concept and motivate a combination of large and small dust grains in brown dwarf atmosphere models.

We present a kinematic analysis of 152 low surface gravity M7-L8 dwarfs by adding 8 parallaxes, 38 radial velocities, and 19 proper motions. We find 39 objects to be high-likelihood or bona fide members of nearby moving groups, 92 objects to be ambiguous members and 21 objects that are non-members. We find that gravity classification and photometric color separate 5-150 Myr sources from > 3 Gyr field objects, but they do not correlate one-to-one with the narrower 5 -150 Myr age range. The absolute magnitudes of low-gravity sources from J band through W3 show a flux redistribution when compared to equivalent field sources that is correlated with spectral subtype. Clouds, which are a far more dominant opacity source for L dwarfs, are the likely cause. On CMDs, the latest-type low-gravity L dwarfs drive the elbow of the L/T transition up to 1 mag redder and 1 mag fainter than field dwarfs at M_J but are consistent with or brighter than the elbow at M_W1 and M_W2. Furthermore, there is an indication on CMD's (such as M_J versus (J-W2) of increasingly redder sequences separated by gravity classification. Examining bolometric luminosities for planets and low-gravity objects, we confirm that young M dwarfs are overluminous while young L dwarfs are normal compared to the field. This translates into warmer M dwarf temperatures compared to the field sequence while lower temperatures for L dwarfs.

Scaling relations between black hole (BH) masses and their host galaxy properties have been studied extensively over the last two decades, and point towards co-evolution of central massive BHs and their hosts. However, these relations remain poorly constrained for BH masses below $\sim10^{6}$ M_sun. Here we present optical and X-ray observations of the dwarf galaxy RGG 118 taken with the Magellan Echellette Spectrograph on the 6.5m Clay Telescope and Chandra X-ray Observatory. Based on Sloan Digital Sky Survey spectroscopy, RGG 118 was identified as possessing narrow emission line ratios indicative of photoionization partly due to an active galactic nucleus. Our higher resolution spectroscopy clearly reveals broad H$\alpha$ emission in the spectrum of RGG 118. Using virial BH mass estimate techniques, we calculate a BH mass of $\sim50,000$ \msun. We detect a nuclear X-ray point source in RGG 118, suggesting a total accretion powered luminosity of $L=4\times10^{40}~{\rm erg~s^{-1}}$, and an Eddington fraction of $\sim1$ per cent. The BH in RGG 118 is the smallest ever reported in a galaxy nucleus and we find that it lies on the extrapolation of the $M_{\rm BH}-\sigma_{\ast}$ relation to the lowest masses yet.

Distinct seed formation mechanisms are imprinted upon the fraction of dwarf galaxies currently containing a central supermassive black hole. Seeding by Pop III remnants is expected to produce a higher occupation fraction than is generated with direct gas collapse precursors. Chandra observations of nearby early-type galaxies can directly detect even low-level supermassive black hole activity, and the active fraction immediately provides a firm lower limit to the occupation fraction. Here, we use the volume-limited AMUSE surveys of ~200 optically-selected early-type galaxies to characterize simultaneously, for the first time, the occupation fraction and the scaling of nuclear X-ray luminosity with stellar mass, accounting for intrinsic scatter, measurement uncertainties, and X-ray limits. For early-type galaxies with log(M_star/M_sun)<10, we obtain a lower limit to the occupation fraction of >20% (at 95% confidence), but full occupation cannot be excluded. The preferred dependence of log(L_X) upon log(M_star) has a slope of about 0.7-0.8, consistent with the "downsizing" trend previously identified from the AMUSE dataset, and a uniform Eddington efficiency is disfavored at ~2 sigma. We provide guidelines for the future precision with which these parameters may be refined with larger or more sensitive samples.

The optical light profiles of nearby early type galaxies are known to exhibit a smooth transition from nuclear light deficits to nuclear light excesses with decreasing galaxy mass, with as much as 80 per cent of the galaxies with stellar masses below 10^10 Msun hosting a massive nuclear star cluster. At the same time, while all massive galaxies are thought to harbor nuclear super-massive black holes (SMBHs), observational evidence for SMBHs is slim at the low end of the mass function. Here, we explore the environmental dependence of the nucleation fraction by comparing two homogeneous samples of nearby field vs. cluster early type galaxies with uniform Hubble Space Telescope (HST) coverage. Existing Chandra X-ray Telescope data for both samples yield complementary information on low-level accretion onto nuclear SMBHs. Specifically, we report on dual-band (F475W & F850LP) Advanced Camera for Surveys (ACS) imaging data for 28 out of the 103 field early type galaxies that compose the AMUSE-Field Chandra survey, and compare our results against the companion HST and Chandra surveys for a sample of 100 Virgo cluster early types (ACS Virgo Cluster and AMUSE-Virgo surveys, respectively). We model the two-dimensional light profiles of the field targets to identify and characterize NSCs, and find a field nucleation fraction of 26% (+17%, -11%; at the 1-sigma level), consistent with the measured Virgo nucleation fraction of 30% (+17%, -12%), across a comparable mass distribution. Coupled with the Chandra result that SMBH activity is higher for the field, our findings indicate that, since the last epoch of star formation, the funneling of gas to the nuclear regions has been inhibited more effectively for Virgo galaxies, arguably via ram pressure stripping.

Ultraluminous X-ray sources (ULXs) are some of the brightest phenomena found outside of a galaxy's nucleus, and their explanation typically invokes accretion of material onto a black hole. Here, we perform the largest population study to date of ULXs in early-type galaxies, focusing on whether a galaxy's large scale environment can affect its ULX content. Using the AMUSE survey, which includes homogeneous X-ray coverage of 100 elliptical galaxies in the Virgo cluster and a similar number of elliptical galaxies in the field (spanning stellar masses of 1e8-1e12 M_sun), we identify 37.9 +/- 10.1 ULXs in Virgo and 28.1 +/- 8.7 ULXs in the field. Across both samples, we constrain the number of ULXs per unit stellar mass, i.e., the ULX specific frequency, to be 0.062 +/- 0.013 ULXs per 10^10 M_sun (or about 1 ULX per 1.6e11 M_sun of galaxy stellar mass). We find that the number of ULXs, the specific frequency of ULXs, and the average ULX spectral properties are all similar in both cluster and field environments. Contrary to late-type galaxies, we do not see any trend between specific ULX frequency and host galaxy stellar mass, and we show that dwarf ellipticals host fewer ULXs than later-type dwarf galaxies at a statistically meaningful level. Our results are consistent with ULXs in early-type galaxies probing the luminous tail of the low-mass X-ray binary population, and are briefly discussed in context of the influence of gravitational interactions on the long-term evolution of a galaxy's (older) stellar population.

We analyze the morphological properties of a large sample of 1503 70 micron selected galaxies in the COSMOS field spanning the redshift range 0.01<z< 3.5 with a median redshift of 0.5 and an infrared luminosity range of 10^8<L_IR<10^14L_sun with a median luminosity of 10^11.4 L_sun. In general these galaxies are massive, with a stellar mass range of 10^10-10^12 M_sun, and luminous, with -25<M_K<-20. We find a strong correlation between the fraction of major mergers and L_IR, with the fraction at the highest luminosity being up to 50%. We also find that the fraction of spirals drops dramatically with L_IR. Minor mergers likely play a role in boosting the infrared luminosity for sources with low luminosities. The precise fraction of mergers in any given L_IR bin varies by redshift due to sources at z>1 being difficult to classify and subject to the effects of band pass shifting, therefore, these numbers can only be considered lower limits. At z<1, where the morphological classifications are most robust, major mergers clearly dominate the ULIRG population (50-80%) and are important for the LIRG population (25-40%). At z>1 the fraction of major mergers is at least 30-40% for ULIRGs. Although the general morphological trends agree with what has been observed for local (U)LIRGs, the fraction of major mergers is slightly lower than seen locally. This is in part due to the difficulty of identifying merger signatures at high redshift. We argue that given the number of major gas-rich mergers observed and the relatively short timescale that they would be observable in the (U)LIRG phase that it is plausible for the observed red sequence of massive ellipticals (<10^12 M_sun) to have been formed entirely by gas-rich major mergers.