Over the past several decades, advances in telescope/detector technologies and deep imaging techniques have pushed surface brightness limits to ever fainter levels. We can now both detect and measure the diffuse, extended star light that surrounds galaxies and permeates galaxy clusters, enabling the study of galaxy halos, tidal streams, diffuse galaxy populations, and the assembly history of galaxies and galaxy clusters. With successes come new challenges, however, and pushing even deeper will require careful attention to systematic sources of error. In this review I highlight recent advances in the study of diffuse starlight in galaxies, and discuss challenges faced by the next generation of deep imaging campaigns.

Deep observations of faint surface brightness stellar tidal streams in external galaxies with LSST are addressed in this White Paper contribution. We propose using the Wide--Fast--Deep survey that contains several nearby galaxies (at distances where the stars themselves are not resolved, i.e., beyond 20 Mpc). In the context of hierarchical galaxy formation, it is necessary to understand the prevalence and properties of tidal substructure around external galaxies based on integrated (i.e., unresolved) diffuse light. This requires collecting observations on much larger samples of galaxies than the Milky Way and M31. We will compare the observed structures to the predictions of cosmological models of galactic halo formation that inform us about the number and properties of streams around Milky Way-like galaxies. The insight gained from these comparisons will allow us to infer the properties of stream progenitors (masses, dynamics, metallicities, stellar populations). The changes in the host galaxies caused by the interactions with the dissolving companion galaxies will be another focus of our studies. We conclude by discussing synergies with WFIRST and Euclid, and also provide concrete suggestions for how the effects of scattered light could be minimized in LSST images to optimize the search for low surface brightness features, such as faint unresolved stellar tidal streams.

Using deep, high resolution optical imaging from the Next Generation Virgo Cluster Survey we study the properties of nuclear star clusters (NSCs) in a sample of nearly 400 quiescent galaxies in the core of Virgo with stellar masses $10^{5}\lesssim M_{*}/M_{\odot} \lesssim10^{12}$. The nucleation fraction reaches a peak value $f_{n}\approx90\%$ for $M_{*} \approx 10^{9} M_{\odot}$ galaxies and declines for both higher and lower masses, but nuclei populate galaxies as small as $M_{*} \approx5\times10^{5} M_{\odot}$. Comparison with literature data for nearby groups and clusters shows that at the low-mass end nucleation is more frequent in denser environments. The NSC mass function peaks at $M_{NSC}\approx7\times10^{5} M_{\odot}$, a factor 3-4 times larger than the turnover mass for globular clusters (GCs). We find a nonlinear relation between the stellar masses of NSCs and of their host galaxies, with a mean nucleus-to-galaxy mass ratio that drops to $M_{NSC}/M_{*}\approx3.6\times10^{-3}$ for $M_{*} \approx 5\times10^{9} M_{\odot}$ galaxies. Nuclei in both more and less massive galaxies are much more prominent: $M_{NSC}\propto M_{*}^{0.46}$ at the low-mass end, where nuclei are nearly 50% as massive as their hosts. We measure an intrinsic scatter in NSC masses at fixed galaxy stellar mass of 0.4 dex, which we interpret as evidence that the process of NSC growth is significantly stochastic. At low galaxy masses we find a close connection between NSCs and GC systems, including a very similar occupation distribution and comparable total masses. We discuss these results in the context of current dissipative and dissipationless models of NSC formation.

We report the detection of diffuse starlight in an extragalactic HI cloud in the nearby Leo~I galaxy group. We detect the source, dubbed BST1047+1156, in both broadband optical and GALEX ultraviolet light. Spanning ~ 2 kpc in radius, it has a peak surface brightness of $\mu_B$=28.8 mag/arcsec$^2$, making it the lowest surface brightness object ever detected via integrated light. Although the object is extremely gas-rich, with a gas fraction of $f_g$=0.99, its peak HI column density is well below levels where star formation is typically observed in galaxies. Nonetheless, BST1047+1156 shows evidence for young stellar populations: along with the detected UV emission, the object is extremely blue, with B-V=0.14$\pm$0.09. The object sports two tidal tails and is found embedded within diffuse gas connecting the spiral galaxy M96 to the group's extended HI Leo Ring. The nature of BST1047+1156 is unclear. It could be a disrupting tidal dwarf, recently spawned from star formation triggered in the Leo I Group's tidal debris. Alternatively, the object may have been a pre-existing galaxy --- the most extreme example of a gas-rich field LSB known to date --- which had a recent burst of star formation triggered by encounters in the group environment.

Intra-cluster (IC) populations are expected to be a natural result of the hierarchical assembly of clusters, yet their low space densities make them difficult to detect and study. We present the first definitive kinematic detection of an IC population of globular clusters (GCs) in the Virgo cluster, around the central galaxy, M87. This study focuses on the Virgo core for which the combination of NGVS photometry and follow-up spectroscopy allows us to reject foreground star contamination and explore GC kinematics over the full Virgo dynamical range. The GC kinematics changes gradually with galactocentric distance, decreasing in mean velocity and increasing in velocity dispersion, eventually becoming indistinguishable from the kinematics of Virgo dwarf galaxies at $\mathrm{R>320\, kpc}$. By kinematically tagging M87 halo and intra-cluster GCs we find that 1) the M87 halo has a smaller fraction ($52\pm3\%$) of blue clusters with respect to the IC counterpart ($77\pm10\%$), 2) the $(g'-r')_{0}$ vs $(i'-z')_{0}$ color-color diagrams reveal a galaxy population that is redder than the IC population that may be due to a different composition in chemical abundance and progenitor mass, and 3) the ICGC distribution is shallower and more extended than the M87 GCs, yet still centrally concentrated. The ICGC specific frequency, $S_{N,\mathrm{ICL}}=10.2\pm4.8$, is consistent with what is observed for the population of quenched, low-mass galaxies within 1~Mpc from the cluster's center. The IC population at Virgo's center is thus consistent with being an accreted component from low-mass galaxies tidally stripped or disrupted through interactions, with a total mass of $\mathrm{M_{ICL,tot}=10.8\pm0.1\times10^{11}M_{\odot}}$.

We use deep Hubble Space Telescope imaging in the outskirts of the nearby spiral M101 to study stellar populations in the galaxy's outer disk and halo. Our ACS field lies 17.6 arcmin (36 kpc) from the center of M101 and targets the blue "NE Plume" of M101's outer disk, while the parallel WFC3 field lies at a distance of 23.3 arcmin (47 kpc) to sample the galaxy's stellar halo. The WFC3 halo field shows a well-defined red giant branch characterized by low metallicity ([M/H]=-1.7 $\pm$ 0.2), with no evidence of young stellar populations. In contrast, the ACS disk field shows multiple stellar populations, including a young main sequence, blue and red helium burning stars, and old RGB and AGB populations. The mean metallicity of these disk stars is quite low: [M/H]=-1.3 $\pm$ 0.2 for the RGB population, and -1.15 $\pm$ 0.2 for the younger helium burning sequences. Of particular interest is a bunching of stars along the BHeB sequence, indicative of an evolving cohort of massive young stars. We show that the young stellar populations in this field are well-described by a decaying burst of star formation that peaked ~ 300-400 Myr ago, along with a more extended star formation history to produce the older RGB and AGB populations. These results confirm and extend the results from our previous deep surface photometry of M101's outer disk, providing an important cross-check on stellar population studies using resolved stellar populations versus integrated light photometry. We discuss our results in the context of halo formation models and the interaction history of M101 and its companions.

We present the discovery of a vast cloud of ionized gas 13$^{\prime}$ (32 kpc) north of the interacting system M51. We detected this cloud via deep narrow-band imaging with the Burrell Schmidt Telescope, where it appears as an extended, diffuse H$\alpha$-emitting feature with no embedded compact regions. The Cloud spans $\sim$10$^{\prime}\times$3$^{\prime}$ (25$\times$7.5 kpc) in size and has no stellar counterpart; comparisons with our previous deep broadband imaging show no detected continuum light to a limit of $\mu_{\rm lim, B} \sim$30 mag arcsec$^{-2}$. WIYN SparsePak observations confirm the cloud's kinematic association with M51, and the high NII/H$\alpha$, SII/H$\alpha$, and OI/H$\alpha$ line ratios we measure imply a hard ionization source such as AGN photoionization or shock heating rather than photoionization due to young stars. Given the strong NII emission, we infer roughly solar metallicity for the cloud, ruling out an origin due to infall of primordial gas. Instead we favor models where the gas has been expelled from the inner regions of the M51 system due to tidal stripping or starburst/AGN winds and has been subsequently ionized either by shocks or a fading AGN. This latter scenario raises the intriguing possibility that M51 may be the nearest example of an AGN fossil nebula or light echo, akin to the famous "Hanny's Voorwerp" in the IC 2497 system.

We present Keck/DEIMOS spectroscopy of globular clusters (GCs) around the ultra-diffuse galaxies (UDGs) VLSB-B, VLSB-D, and VCC615 located in the central regions of the Virgo cluster. We spectroscopically identify 4, 12, and 7 GC satellites of these UDGs, respectively. We find that the three UDGs have systemic velocities ($V_{sys}$) consistent with being in the Virgo cluster, and that they span a wide range of velocity dispersions, from $\sim 16$ to $\sim 47$ km/s, and high dynamical mass-to-light ratios within the radius that contains half the number of GCs ($ 407^{+916}_{-407}$, $21^{+15}_{-11}$, $60^{+65}_{-38}$, respectively). VLSB-D shows possible evidence for rotation along the stellar major axis and its $V_{sys}$ is consistent with that of the massive galaxy M84 and the center of the Virgo cluster itself. These findings, in addition to having a dynamically and spatially ($\sim 1$ kpc) off-centered nucleus and being extremely elongated, suggest that VLSB-D could be tidally perturbed. On the contrary, VLSB-B and VCC615 show no signals of tidal deformation. Whereas the dynamics of VLSB-D suggest that it has a less massive dark matter halo than expected for its stellar mass, VLSB-B and VCC615 are consistent with a $\sim 10^{12}$ M$_{\odot}$ dark matter halo. Although our samples of galaxies and GCs are small, these results suggest that UDGs may be a diverse population, with their low surface brightnesses being the result of very early formation, tidal disruption, or a combination of the two.

We present new Chandra observations of the outer halo of the giant elliptical galaxy NGC 4472 (M49) in the Virgo Cluster. The data extend to 130 kpc (28'), and have a combined exposure time of 150 ks. After identifying optical counterparts using the Next Generation Virgo Cluster Survey to remove background active galactic nuclei and globular cluster (GC) sources, and correcting for completeness, we find that the number of field low-mass X-ray binaries (LMXBs) per unit stellar V-band light increases significantly with galactocentric radius. Because the flux limit of the complete sample corresponds to the Eddington limit for neutron stars in NGC 4472, many of the ~90 field LMXBs in this sample could host black holes. The excess of field LMXBs at large galactocentric radii may be partially caused by natal kicks on black holes and neutron stars in binary systems in the inner part of the galaxy. Furthermore, since the metallicity in the halo of NGC 4472 strongly decreases towards larger galactocentric radii, the number of field LMXBs per unit stellar mass is anti-correlated with metallicity, opposite to what is observed in GCs. Another way to explain the spatial distribution of field LMXBs is therefore a reversed metallicity effect, although we have not identified a mechanism to explain this in terms of stellar and binary evolution.

We present a multiwavelength study of star formation within the nearby M101 Group, including new deep H$\alpha$ imaging of M101 and its two companions. We perform a statistical analysis of the H$\alpha$ to FUV flux ratios in HII regions located in three different environments: M101's inner disk, M101's outer disk, and M101's lower mass companion galaxy NGC~5474. We find that, once bulk radial trends in extinction are taken into account, both the median and scatter in F$_{\rm H\alpha}/$F$_{\rm FUV}$ in HII regions are invariant across all of these environments. Also, using Starburst99 models, we are able to qualitatively reproduce the distributions of F$_{\rm H\alpha}/$F$_{\rm FUV}$ throughout these different environments using a standard Kroupa IMF, hence we find no need to invoke truncations in the upper mass end of the IMF to explain the young star-forming regions in the M101 Group even at extremely low surface density. This implies that star formation in low density environments differs from star formation in high density environments only by intensity and not by cloud-to-cloud physics.

We present a deep imaging survey of the Virgo Cluster, designed to study the connection between cluster galaxies and Virgo's diffuse intracluster light (ICL). Our observations span roughly 16 square degrees and reach a 3-sigma depth of mu(B)=29.5 and mu(V)=28.5 mag/arcsec^2. At these depths, the limiting systematic uncertainties are astrophysical: scattered starlight from foreground Galactic dust, and variations in faint background sources. The dust-scattered starlight is well-traced by deep far-infrared imaging, making it distinguishable from true Virgo diffuse light. Our imaging maps the Virgo core around M87 and the adjacent M86/M84 region, in subcluster B around M49, and in the more distant W' cloud around NGC 4365. Most of the detected ICL is found in the Virgo core and within the W' cloud, with little evidence for extensive ICL in subcluster B. The large amount of diffuse light seen in the infalling W' cloud likely illustrates the importance of the group environment for generating ICL. The bulk of the detected ICL is fairly red (B-V=0.7-0.9), indicative of old stellar populations. We estimate a total Virgo ICL fraction of 7-15%, somewhat smaller than expected for massive, evolved clusters, suggesting that Virgo is still growing its ICL component. We trace M87's extremely boxy halo out to ~ 150 kpc, and show that the current stripping rate of low luminosity galaxies is insufficient to have built M87's outer halo over a Hubble time. Finally, we identify another large ultra-diffuse galaxy in Virgo, likely in the process of being shredded by the cluster tidal field.

The Next Generation Virgo Cluster Survey is a deep (with a $2\sigma$ detection limit $\mu_g$ = 29~mag~arcsec$^{-2}$ in the $g-$band) optical panchromatic survey targeting the Virgo cluster from its core to virial radius, for a total areal coverage of 104 square degrees. As such, the survey is well suited for the study of galaxies' outskirts, haloes and low surface brightness features that arise from dynamical interactions within the cluster environment. We report the discovery of extremely faint ($\mu_g$ $>$ 25 mag arcsec$^{-2}$) shells in three Virgo cluster early-type dwarf galaxies, VCC~1361, VCC~1447 and VCC~1668. Among them, VCC~1447 has an absolute magnitude M$_{g}$ = -11.71 mag and is {\it the least massive galaxy with a shell system discovered to date}. We present a detailed study of these low surface brightness features. We detect between three and four shells in each of our galaxies. Within the uncertainties, we find no evidence of a color difference between the galaxy main body and shell features. The observed arcs of the shells are located upto several effective radii of the galaxies. We further explore the origin of these low surface brightness features with the help of idealized numerical simulations. We find that a near equal mass merger is best able to reproduce the main properties of the shells, including their quite symmetric appearance and their alignment along the major axis of the galaxy. The simulations provide support for a formation scenario in which a recent merger, between two near-equal mass, gas-free dwarf galaxies forms the observed shell systems.

We use deep optical photometry from the Next Generation Virgo Cluster Survey [NGVS] to investigate the color-magnitude diagram for the galaxies inhabiting the core of this cluster. The sensitivity of the NGVS imaging allows us to continuously probe galaxy colors over a factor of $\sim 2 \times 10^5$ in luminosity, from brightest cluster galaxies to scales overlapping classical satellites of the Milky Way [$M_{g^{\prime}}$ $\sim$ $-$9; $M_{*}$ $\sim 10^6$ M$_{\odot}$], within a single environment. Remarkably, we find the first evidence that the RS flattens in all colors at the faint-magnitude end [starting between $-$14 $\le$ $M_{g^{\prime}}$ $\le$ $-$13, around $M_{*}$ $\sim 4 \times 10^7$ M$_{\odot}$], with the slope decreasing to $\sim$60% or less of its value at brighter magnitudes. This could indicate that the stellar populations of faint dwarfs in Virgo's core share similar characteristics [e.g. constant mean age] over $\sim$3 mags in luminosity, suggesting that these galaxies were quenched coevally, likely via pre-processing in smaller hosts. We also compare our results to galaxy formation models, finding that the RS in model clusters have slopes at intermediate magnitudes that are too shallow, and in the case of semi-analytic models, do not reproduce the flattening seen at both extremes [bright/faint] of the Virgo RS. Deficiencies in the chemical evolution of model galaxies likely contribute to the model-data discrepancies at all masses, while overly efficient quenching may also be a factor at dwarf scales. Deep UV and near-IR photometry are required to unambiguously diagnose the cause of the faint-end flattening.

We present an analysis of high-quality photometry for globular clusters (GCs) in the Virgo cluster core region, based on data from the Next Generation Virgo Cluster Survey (NGVS) pilot field, and in the Milky Way (MW) based on VLT/X-Shooter spectrophotometry. We find significant discrepancies in color-color diagrams between sub-samples from different environments, confirming that the environment has a strong influence on the integrated colors of GCs. GC color distributions along a single color are not sufficient to capture the differences we observe in color-color space. While the average photometric colors become bluer with increasing radial distance to the cD galaxy M87, we also find a relation between the environment and the slope and intercept of the color-color relations. A denser environment seems to produce a larger dynamic range in certain color indices. We argue that these results are not due solely to differential extinction, IMF variations, calibration uncertainties, or overall age/metallicity variations. We therefore suggest that the relation between the environment and GC colors is, at least in part, due to chemical abundance variations, which affect stellar spectra and stellar evolution tracks. Our results demonstrate that stellar population diagnostics derived from model predictions which are calibrated on one particular sample of GCs may not be appropriate for all extragalactic GCs. These results advocate a more complex model of the assembly history of GC systems in massive galaxies that goes beyond the simple bimodality found in previous decades.

We present deep, wide-field imaging of the M96 galaxy group (also known as the Leo I Group). Down to surface brightness limits of $\mu_{B}=30.1$ and $\mu_{V}=29.5$, we find no diffuse, large-scale optical counterpart to the "Leo Ring", an extended HI ring surrounding the central elliptical M105 (NGC 3379). However, we do find a number of extremely low surface-brightness ($\mu_{B} \gtrsim 29$) small-scale streamlike features, possibly tidal in origin, two of which may be associated with the Ring. In addition we present detailed surface photometry of each of the group's most massive members -- M105, NGC 3384, M96 (NGC 3368), and M95 (NGC 3351) -- out to large radius and low surface brightness, where we search for signatures of interaction and accretion events. We find that the outer isophotes of both M105 and M95 appear almost completely undisturbed, in contrast to NGC 3384 which shows a system of diffuse shells indicative of a recent minor merger. We also find photometric evidence that M96 is accreting gas from the HI ring, in agreement with HI data. In general, however, interaction signatures in the M96 Group are extremely subtle for a group environment, and provide some tension with interaction scenarios for the formation of the Leo HI Ring. The lack of a significant component of diffuse intragroup starlight in the M96 Group is consistent with its status as a loose galaxy group in which encounters are relatively mild and infrequent.

We present results from deep, wide-field surface photometry of three nearby (D=4--7 Mpc) spiral galaxies: M94 (NGC 4736), M64 (NGC 4826), and M106 (NGC 4258). Our imaging reaches limiting surface brightnesses of $\mu_{B} \sim$ 28 -- 30 mag arcsec$^{-2}$ and probes colors down to $\mu_{B} \sim$ 27.5 mag arcsec$^{-2}$. We compare our broadband optical data to available ultraviolet and high column-density HI data to better constrain the star forming history and stellar populations of the outermost parts of each galaxy's disk. Each galaxy has a well-defined radius beyond which little star formation occurs and the disk light appears both azimuthally smooth and red in color, suggestive of old, well-mixed stellar populations. Given the lack of ongoing star formation or blue stellar populations in these galaxies' outer disks, the most likely mechanisms for their formation are dynamical processes such as disk heating or radial migration, rather than inside-out growth of the disks. This is also implied by the similarity in outer disk properties despite each galaxy showing distinct levels of environmental influence, from a purely isolated galaxy (M94) to one experiencing weak tidal perturbations from its satellite galaxies (M106) to a galaxy recovering from a recent merger (M64), suggesting that a variety of evolutionary histories can yield similar outer disk structure. While this suggests a common secular mechanism for outer disk formation, the large extent of these smooth, red stellar populations---which reach several disk scalelengths beyond the galaxies' spiral structure---may challenge models of radial migration given the lack of any non-axisymmetric forcing at such large radii.

We present measurements of the galaxy luminosity and stellar mass function in a 3.71 deg$^2$ (0.3 Mpc$^2$) area in the core of the Virgo cluster, based on $ugriz$ data from the Next Generation Virgo Cluster Survey (NGVS). The galaxy sample consists of 352 objects brighter than $M_g=-9.13$ mag, the 50% completeness limit of the survey. Using a Bayesian analysis, we find a best-fit faint end slope of $\alpha=-1.33 \pm 0.02$ for the g-band luminosity function; consistent results are found for the stellar mass function as well as the luminosity function in the other four NGVS bandpasses. We discuss the implications for the faint-end slope of adding 92 ultra compact dwarfs galaxies (UCDs) -- previously compiled by the NGVS in this region -- to the galaxy sample, assuming that UCDs are the stripped remnants of nucleated dwarf galaxies. Under this assumption, the slope of the luminosity function (down to the UCD faint magnitude limit, $M_g = -9.6$ mag) increases dramatically, up to $\alpha = -1.60 \pm 0.06$ when correcting for the expected number of disrupted non-nucleated galaxies. We also calculate the total number of UCDs and globular clusters that may have been deposited in the core of Virgo due to the disruption of satellites, both nucleated and non-nucleated. We estimate that ~150 objects with $M_g\lesssim-9.6$ mag and that are currently classified as globular clusters, might, in fact, be the nuclei of disrupted galaxies. We further estimate that as many as 40% of the (mostly blue) globular clusters in the core of Virgo might once have belonged to such satellites; these same disrupted satellites might have contributed ~40% of the total luminosity in galaxies observed in the core region today. Finally, we use an updated Local Group galaxy catalog to provide a new measurement of the luminosity function of Local Group satellites, $\alpha=-1.21\pm0.05$.

The largest stellar halos in the universe are found in massive galaxy clusters, where interactions and mergers of galaxies, along with the cluster tidal field, all act to strip stars from their host galaxies and feed the diffuse intracluster light (ICL) and extended halos of brightest cluster galaxies (BCGs). Studies of the nearby Virgo Cluster reveal a variety of accretion signatures imprinted in the morphology and stellar populations of its ICL. While simulations suggest the ICL should grow with time, attempts to track this evolution across clusters spanning a range of mass and redshift have proved difficult due to a variety of observational and definitional issues. Meanwhile, studies of nearby galaxy groups reveal the earliest stages of ICL formation: the extremely diffuse tidal streams formed during interactions in the group environment.

We report on the properties of the most massive ultra-compact dwarf galaxy (UCD) in the nearby Virgo Cluster of galaxies using imaging from the Next Generation Virgo Cluster Survey (NGVS) and spectroscopy from Keck/DEIMOS. This object (M59-UCD3) appears to be associated with the massive Virgo galaxy M59 (NGC 4621), has an integrated velocity dispersion of 78 km/s, a dynamical mass of $3.7\times10^8 M_\odot$, and an effective radius ($R_e$) of 25 pc. With an effective surface mass density of $9.4\times10^{10} M_\odot/kpc^2$, it is the densest galaxy in the local Universe discovered to date, surpassing the density of the luminous Virgo UCD, M60-UCD1. M59-UCD3 has a total luminosity of $M_{g'}=-14.2$ mag, and a spectral energy distribution consistent with an old (14 Gyr) stellar population with [Fe/H]=0.0 and [$\alpha$/Fe]=+0.2. We also examine deep imaging around M59 and find a broad low surface brightness stream pointing towards M59-UCD3, which may represent a tidal remnant of the UCD progenitor. This UCD, along with similar objects like M60-UCD1 and M59cO, likely represents an extreme population of tidally stripped galaxies more akin to larger and more massive compact early-type galaxies than to nuclear star clusters in present-day dwarf galaxies.

We use imaging from the Next Generation Virgo cluster Survey (NGVS) to present a comparative study of ultra-compact dwarf (UCD) galaxies associated with three prominent Virgo sub-clusters: those centered on the massive, red-sequence galaxies M87, M49 and M60. We show how UCDs can be selected with high completeness using a combination of half-light radius and location in color-color diagrams ($u^*iK_s$ or $u^*gz$). Although the central galaxies in each of these sub-clusters have nearly identical luminosities and stellar masses, we find large differences in the sizes of their UCD populations, with M87 containing ~3.5 and 7.8 times more UCDs than M49 and M60, respectively. The relative abundance of UCDs in the three regions scales in proportion to sub-cluster mass, as traced by X-ray gas mass, total gravitating mass, number of globular clusters, and number of nearby galaxies. We find that the UCDs are predominantly blue in color, with ~85% of the UCDs having colors similar to blue GCs and stellar nuclei of dwarf galaxies. We present evidence that UCDs surrounding M87 and M49 may follow a morphological sequence ordered by the prominence of their outer, low surface brightness envelope, ultimately merging with the sequence of nucleated low-mass galaxies, and that envelope prominence correlates with distance from either galaxy. Our analysis provides evidence that tidal stripping of nucleated galaxies is an important process in the formation of UCDs.

We present kinematic and photometric evidence for an accretion event in the halo of the cD galaxy M87 in the last Gyr. Using velocities for ~300 planetary nebulas (PNs) in the M87 halo, we identify a chevron-like substructure in the PN phase-space. We implement a probabilistic Gaussian mixture model to identify the PNs that belong to the chevron. From analysis of deep V-band images of M87, we find that the region with the highest density of PNs associated to the chevron, is a crown-shaped substructure in the optical light. We assign a total of N_(PN,sub)=54 to the substructure, which extends over ~50 kpc along the major axis where we also observe radial variations of the ellipticity profile and a colour gradient. The substructure has highest surface brightness in a 20kpc x 60kpc region around 70 kpc in radius. In this region, it causes an increase in surface brightness by >60%. The accretion event is consistent with a progenitor galaxy with a V-band luminosity of L=2.8\pm1.0 x 10^9 L_(sun,V), a colour of (B-V)=0.76\pm0.05, and a stellar mass of M=6.4\pm2.3 x 10^9 M_sun. The accretion of this progenitor galaxy has caused an important modification of the outer halo of M87 in the last Gyr. By itself it is strong evidence that the galaxy's cD halo is growing through the accretion of smaller galaxies as predicted by hierarchical galaxy evolution models.

We present deep, wide-field imaging of the M51 system using CWRU's Burrell Schmidt telescope at KPNO to study the faint tidal features that constrain its interaction history. Our images trace M51's tidal morphology down to a limiting surface brightness of $\mu_{B,lim}\sim $30 mag arcsec$^{-2}$, and provide accurate colors ($\sigma_{B-V} < 0.1$) down to $\mu_B\sim 28$. We identify two new tidal streams in the system (the South and Northeast Plumes) with surface brightnesses of $\mu_B =29$ and luminosities of $\sim 10^6 L_{\odot,B}$. While the Northeast Plume may be a faint outer extension of the tidal "crown" north of NGC 5195 (M51b), the South Plume has no analogue in any existing M51 simulation and may represent a distinct tidal stream or disrupted dwarf galaxy. We also trace the extremely diffuse Northwest Plume out to a total extent of 20' (43 kpc) from NGC 5194(M51a), and show it to be physically distinct from the overlapping bright tidal streams from M51b. The Northwest Plume's morphology and red color ($B-V=0.8$) instead argue that it originated from tidal stripping of M51a's extreme outer disk. Finally, we confirm the strong segregation of gas and stars in the Southeast Tail, and do not detect any diffuse stellar component in the HI portion of the tail. Extant simulations of M51 have difficulty matching both the wealth of tidal structure in the system and the lack of stars in the HI tail, motivating new modeling campaigns to study the dynamical evolution of this classic interacting system.

We report on a large-scale study of the distribution of globular clusters (GCs) throughout the Virgo cluster, based on photometry from the Next Generation Virgo Cluster Survey, a large imaging survey covering Virgo's primary subclusters to their virial radii. Using the g', (g'-i') color-magnitude diagram of unresolved and marginally-resolved sources, we constructed 2-D maps of the GC distribution. We present the clearest evidence to date showing the difference in concentration between red and blue GCs over the extent of the cluster, where the red (metal-rich) GCs are largely located around the massive early-type galaxies, whilst the blue (metal-poor) GCs have a more extended spatial distribution, with significant populations present beyond 83' (215 kpc) along the major axes of M49 and M87. The GC distribution around M87 and M49 shows remarkable agreement with the shape, ellipticity and boxiness of the diffuse light surrounding both galaxies. We find evidence for spatial enhancements of GCs surrounding M87 that may be indicative of recent interactions or an ongoing merger history. We compare the GC map to the locations of Virgo galaxies and the intracluster X-ray gas, and find good agreement between these baryonic structures. The Virgo cluster contains a total population of 67300$\pm$14400 GCs, of which 35% are located in M87 and M49 alone. We compute a cluster-wide specific frequency S_N,CL=$2.8\pm0.7$, including Virgo's diffuse light. The GC-to-baryonic mass fraction is e_b=$5.7\pm1.1\times10^{-4} $and the GC-to-total cluster mass formation efficiency is e_t=$2.9\pm0.5\times10^{-5}$, values slightly lower than, but consistent with, those derived for individual galactic halos. Our results show that the production of the complex structures in the unrelaxed Virgo cluster core (including the diffuse intracluster light) is an ongoing process.(abridged)

We present an investigation into the origins of a series of interlaced narrow filamentary stellar structures, loops and plumes in the vicinity of the Virgo Cluster, edge-on spiral galaxy, NGC 4216 that were previously identified by the Blackbird Telescope. Using the deeper, higher-resolution and precisely calibrated optical CFHT/MegaCam images obtained as part of the Next Generation Virgo Cluster Survey (NGVS), we confirm the previously identified features and identify a few additional structures. The NGVS data allowed us to make a physical study of these low-surface brightness features and investigate their origin. The likely progenitors of the structures were identified as either already catalogued VCC dwarfs or newly discovered satellites caught in the act of being destroyed. They have the same g-i color index and likely contain similar stellar populations. The alignment of three dwarfs along an apparently single stream is intriguing, and we cannot totally exclude that these are second-generation dwarf galaxies being born inside the filament from the debris of an original dwarf. The observed complex structures, including in particular a stream apparently emanating from a satellite of a satellite, point to a high rate of ongoing dwarf destruction/accretion in the region of the Virgo Cluster where NGC 4216 is located. We discuss the age of the interactions and whether they occurred in a group that is just falling into the cluster and shows signs of so-called "pre-processing" before it gets affected by the cluster environment, or in a group which already ventured towards the central regions of Virgo Cluster.

We study two nearby, early-type galaxies, NGC4342 and NGC4291, that host unusually massive black holes relative to their low stellar mass. The observed black hole-to-bulge mass ratios of NGC4342 and NGC4291 are ~6.9% and ~1.9%, respectively, which significantly exceed the typical observed ratio of ~0.2%. As a consequence of the exceedingly large black hole-to-bulge mass ratios, NGC4342 and NGC4291 are ~5.1 sigma and ~3.4 sigma outliers from the M_BH - M_bulge scaling relation, respectively. In this paper, we explore the origin of the unusually high black hole-to-bulge mass ratio. Based on Chandra X-ray observations of the hot gas content of NGC4342 and NGC4291, we compute gravitating mass profiles, and conclude that both galaxies reside in massive dark matter halos, which extend well beyond the stellar light. The presence of dark matter halos around NGC4342 and NGC4291 and a deep optical image of the environment of NGC4342 indicate that tidal stripping, in which >90% of the stellar mass was lost, cannot explain the observed high black hole-to-bulge mass ratios. Therefore, we conclude that these galaxies formed with low stellar masses, implying that the bulge and black hole did not grow in tandem. We also find that the black hole mass correlates well with the properties of the dark matter halo, suggesting that dark matter halos may play a major role in regulating the growth of the supermassive black holes.

Context: We study the peculiar interacting galaxy system of VCC1249/M49 located in the core of the Virgo B subcluster. Owing to a recent interaction between the dwarf galaxy VCC1249 and the halo gas of the gE M49, neutral hydrogen has been displaced from the interstellar medium of this dwarf into the Virgo ICM. Observations also reveal multiple compact star-forming regions that are embedded in this HI cloud, with a projected separation up to 13 kpc from VCC1249 in the northwest direction. Aims: Motivated by recent NUV imaging from GUViCS of the VCC1249/M49 system that shows significant ongoing/recent star formation in the compact regions, we aim to constrain the origin of these outlying HII regions with a multi-wavelength approach. Methods: Using deep optical (u, g, i, z) imaging from NGVS and new Halpha imaging obtained at the San Pedro Martir observatory together with Keck long-slit spectroscopy, we characterize the SFR, ages, and metallicity of VCC1249 and its outlying compact regions. Moreover, we analyze the color and luminosity profile of the galaxy to investigate its recent interaction with M49. Results: Our new observations indicate that VCC1249 underwent a recent interaction with M49 in which both ram-pressure stripping and tidal interaction occured. The joint action of the two mechanisms led to the removal of the HI gas from the ISM of VCC1249, while the gravitational tides triggered the stellar tail and counter-tail of VCC1249. Our SED analysis reveals that the star formation in this galaxy was truncated around 200 Myr ago and that the outlying HII regions were born in situ about 10 Myr ago out of pre-enriched gas removed from the dwarf galaxy. These observations also reveal that interactions between central and satellite galaxies similar to the one between VCC1249/M49 may be an effective way of dispersing metals into the halos of massive galaxies.

The halos of galaxies preserve unique records of their formation histories. We carry out the first combined observational and theoretical study of phase-space halo substructure in an early-type galaxy: M87, the central galaxy in the Virgo cluster. We analyze an unprecedented wide-field, high-precision photometric and spectroscopic data set for 488 globular clusters (GCs), which includes new, large-radius Subaru/Suprime-Cam and Keck/DEIMOS observations. We find signatures of two substructures in position-velocity phase-space. One is a small, cold stream associated with a known stellar filament in the outer halo; the other is a large shell-like pattern in the inner halo that implies a massive, hitherto unrecognized accretion event. We perform extensive statistical tests and independent metallicity analyses to verify the presence and characterize the properties of these features, and to provide more general methodologies for future extragalactic studies of phase-space substructure. The cold outer stream is consistent with a dwarf galaxy accretion event, while for the inner shell there is tension between a low progenitor mass implied by the cold velocity dispersion, and a high mass from the large number of GCs, which might be resolved by a ~0.5 L* E/S0 progenitor. We also carry out proof-of-principle numerical simulations of the accretion of smaller galaxies in an M87-like gravitational potential. These produce analogous features to the observed substructures, which should have observable lifetimes of ~1 Gyr. The shell and stream GCs together support a scenario where the extended stellar envelope of M87 has been built up by a steady rain of material that continues until the present day. This phase-space method demonstrates unique potential for detailed tests of galaxy formation beyond the Local Group.

We present 3.6 and 4.5 micron Spitzer IRAC imaging over 0.77 square degrees at the Virgo cluster core for the purpose of understanding the formation mechanisms of the low surface brightness intracluster light features. Instrumental and astrophysical backgrounds that are hundreds of times higher than the signal were carefully characterized and removed. We examine both intracluster light plumes as well as the outer halo of the giant elliptical M87. For two intracluster light plumes, we use optical colors to constrain their ages to be greater than 3 & 5 Gyr, respectively. Upper limits on the IRAC fluxes constrain the upper limits to the masses, and optical detections constrain the lower limits to the masses. In this first measurement of mass of intracluster light plumes we find masses in the range of 5.5 x 10^8 - 4.5 x 10^9 and 2.1 x 10^8 - 1.5 x 10^9 solar masses for the two plumes for which we have coverage. Given their expected short lifetimes, and a constant production rate for these types of streams, integrated over Virgo's lifetime, they can account for the total ICL content of the cluster implying that we do not need to invoke ICL formation mechanisms other than gravitational mechanisms leading to bright plumes. We also examined the outer halo of the giant elliptical M87. The color profile from the inner to outer halo of M87 (160 Kpc) is consistent with either a flat or optically blue gradient, where a blue gradient could be due to younger or lower metallicity stars at larger radii. The similarity of the age predicted by both the infrared and optical colors (> few Gyr) indicates that the optical measurements are not strongly affected by dust extinction.

We present the initial results of a 40 night contiguous ground-based campaign of time series photometric observations of a 1.39 sq. deg field located within the NASA Kepler mission field of view. The goal of this pre-launch survey was to search for transiting extrasolar planets and to provide independent variability information of stellar sources. We have gathered a data set containing light curves of 54,687 stars from which we have created a statistical sub-sample of 13,786 stars between 14< r <18.5 and have statistically examined each light curve to test for variability. We present a summary of our preliminary photometric findings including the overall level and content of stellar variability in this portion of the Kepler field and give some examples of unusual variable stars found within. We present a preliminary catalog of 2,457 candidate variable stars, of which 776 show signs of periodicity. We also present three potential exoplanet candidates, all of which should be observable in detail by the Kepler mission.

We continue our deep optical imaging survey of the Virgo cluster using the CWRU Burrell Schmidt telescope by presenting B-band surface photometry of the core of the Virgo cluster in order to study the cluster's intracluster light (ICL). We find ICL features down to mu_b ~ 29 mag sq. arcsec, confirming the results of Mihos et al. (2005), who saw a vast web of low-surface brightness streams, arcs, plumes, and diffuse light in the Virgo cluster core using V-band imaging. By combining these two data sets, we are able to measure the optical colors of many of the cluster's low-surface brightness features. While much of our imaging area is contaminated by galactic cirrus, the cluster core near the cD galaxy, M87, is unobscured. We trace the color profile of M87 out to over 2000 arcsec, and find a blueing trend with radius, continuing out to the largest radii. Moreover, we have measured the colors of several ICL features which extend beyond M87's outermost reaches and find that they have similar colors to the M87's halo itself, B-V ~ 0.8. The common colors of these features suggests that the extended outer envelopes of cD galaxies, such as M87, may be formed from similar streams, created by tidal interactions within the cluster, that have since dissolved into a smooth background in the cluster potential.

Using a suite of N-body simulations of galaxy clusters specifically tailored to study the intracluster light (ICL) component, we measure the quantity of ICL using a number of different methods previously employed in the literature for both observational and simulation data sets. By measuring the ICL of the clusters using multiple techniques, we identify systematic differences in how each detection method identifies the ICL. We find that techniques which define the ICL solely based on the current position of the cluster luminosity, such as a surface brightness or local density threshold, tend to find less ICL than methods utilizing time or velocity information, including stellar particles' density history or binding energy. The range of ICL fractions (the fraction of the clusters' total luminosity found in the ICL component) we measure at z=0 across all our clusters using any definition span the range from 9-36%, and even within a single cluster different methods can change the measured ICL fraction by up to a factor of two. Separating the cluster's central galaxy from the surrounding ICL component is a challenge for all ICL techniques, and because the ICL is centrally concentrated within the cluster, the differences in the measured ICL quantity between techniques are largely a consequence of this central galaxy/ICL separation. We thoroughly explore the free parameters involved with each measurement method, and find that adjusting these parameters can change the measured ICL fraction by up to a factor of two. While for all definitions the quantity of ICL tends to increase with time, the ICL fraction does not grow at a uniform rate, nor even monotonically under some definitions. Thus, the ICL can be used as a rough indicator of dynamical age, where more dynamically advanced clusters will on average have higher ICL fractions.

We use deep V-band surface photometry of five of the brightest elliptical galaxies in the Virgo cluster to search for diffuse tidal streams, shells, and plumes in their outer halos (r > 50 kpc). We fit and subtract elliptical isophotal models from the galaxy images to reveal a variety of substructure, with surface brightnesses in the range mu_V= 26-29 mag/arcsec^2. M49 possesses an extended, interleaved shell system reminiscent of the radial accretion of a satellite companion, while M89's complex system of shells and plumes suggests a more complicated accretion history involving either multiple events or a major merger. M87 has a set of long streamers as might be expected from stripping of low luminosity dwarfs on radial orbits in Virgo. M86 also displays a number of small streams indicative of stripping of dwarf companions, but these comprise much less luminosity than those of M87. Only M84 lacks significant tidal features. We quantify the photometric properties of these structures, and discuss their origins in the context of each galaxy's environment and kinematics within the Virgo cluster.

Paper has been withdrawn due to a critical error in the numerical evaluation of the capture fraction.

We compare the distribution of diffuse intracluster light detected in the Virgo Cluster via broadband imaging with that inferred from searches for intracluster planetary nebulae (IPNe). We find a rough correspondence on large scales (~ 100 kpc) between the two, but with very large scatter (~ 1.3 mag/arcsec^2). On smaller scales (1 -- 10 kpc), the presence or absence of correlation is clearly dependent on the underlying surface brightness. On these scales, we find a correlation in regions of higher surface brightness (mu_V < ~27) which are dominated by the halos of large galaxies such as M87, M86, and M84. In those cases, we are likely tracing PNe associated with galaxies rather than true IPNe. In true intracluster fields, at lower surface brightness, the correlation between luminosity and IPN candidates is much weaker. While a correlation between broadband light and IPNe is expected based on stellar populations, a variety of statistical, physical, and methodological effects can act to wash out this correlation and explain the lack of a strong correlation at lower surface brightness found here. [abridged]

Using N-body simulations, we have modeled the production and evolution of substructures in the intracluster light (ICL) of a simulated galaxy cluster. We use a density-based definition of ICL, where ICL consists of luminous particles which are at low densities, to identify ICL particles and track their evolution. We have implemented a friends-of-friends-type clustering algorithm which finds groups of particles correlated in both position and velocity space to identify substructures in the ICL, hereafter referred to as ``streams''. We find that ~40% of the cluster's ICL is generated in the form of these massive (M > 7.0x10^8 Msun), dynamically cold streams. The fraction of the ICL generated in streams is greater early in the cluster's evolution, when galaxies are interacting in the group environment, than later in its evolution when the massive cluster potential has been assembled. The production of streams requires the strong tidal fields associated with close interactions between pairs of galaxies, and is usually associated with merging pairs of galaxies, or fast, close encounters with the cluster's central galaxy. Once streams are formed, they begin to decay as they are disrupted by the tidal field of the cluster. We find that streams have decay timescales which are ~1.5 times their dynamical time in the cluster.

We report the results from a comprehensive study of 74 ultraluminous infrared galaxies (ULIRGs) and 34 Palomar-Green (PG) quasars within z ~ 0.3$ observed with the Spitzer Infrared Spectrograph (IRS). The contribution of nuclear activity to the bolometric luminosity in these systems is quantified using six independent methods that span a range in wavelength and give consistent results within ~ +/-10-15% on average. The average derived AGN contribution in ULIRGs is ~35-40%, ranging from ~15-35% among "cool" (f_25/f_60 =< 0.2) optically classified HII-like and LINER ULIRGs to ~50 and ~75% among warm Seyfert 2 and Seyfert 1 ULIRGs, respectively. This number exceeds ~80% in PG QSOs. ULIRGs fall in one of three distinct AGN classes: (1) objects with small extinctions and large PAH equivalent widths are highly starburst-dominated; (2) systems with large extinctions and modest PAH equivalent widths have larger AGN contributions, but still tend to be starburst-dominated; and (3) ULIRGs with both small extinctions and small PAH equivalent widths host AGN that are at least as powerful as the starbursts. The AGN contributions in class 2 ULIRGs are more uncertain than in the other objects, and we cannot formally rule out the possibility that these objects represent a physically distinct type of ULIRGs. A morphological trend is seen along the sequence (1)-(2)-(3), in general agreement with the standard ULIRG - QSO evolution scenario and suggestive of a broad peak in extinction during the intermediate stages of merger evolution. However, the scatter in this sequence, implies that black hole accretion, in addition to depending on the merger phase, also has a strong chaotic/random component, as in local AGN. (abridged)

We present a detailed comparison of the CO(3-2) emitting molecular gas between a local sample of luminous infrared galaxies (U/LIRGs) and a high redshift sample that comprises submm selected galaxies (SMGs), quasars, and Lyman Break Galaxies (LBGs). The U/LIRG sample consists of our recent CO(3-2) survey using the Submillimeter Array while the CO(3-2) data for the high redshift population are obtained from the literature. We find that the L(CO(3-2)) and L(FIR) relation is correlated over five orders of magnitude, which suggests that the molecular gas traced in CO(3-2) emission is a robust tracer of dusty star formation activity. The near unity slope of 0.93 +/- 0.03 obtained from a fit to this relation suggests that the star formation efficiency is constant to within a factor of two across different types of galaxies residing in vastly different epochs. The CO(3-2) size measurements suggest that the molecular gas disks in local U/LIRGs (0.3 - 3.1 kpc) are much more compact than the SMGs (3 - 16 kpc), and that the size scales of SMGs are comparable to the nuclear separation (5 - 40 kpc) of the widely separated nuclei of U/LIRGs in our sample. We argue from these results that the SMGs studied here are predominantly intermediate stage mergers, and that the wider line-widths arise from the violent merger of two massive gas-rich galaxies taking place deep in a massive halo potential.

We model the NFW potential to determine if, and under what conditions, the NFW halo appears consistent with the observed velocity fields of low surface brightness (LSB) galaxies. We present mock DensePak IFU velocity fields and rotation curves of axisymmetric and non-axisymmetric potentials that are well-matched to the spatial resolution and velocity range of our sample galaxies. We find that the DensePak IFU can accurately reconstruct the velocity field produced by an axisymmetric NFW potential and that a tilted-ring fitting program can successfully recover the corresponding NFW rotation curve. We also find that non-axisymmetric potentials with fixed axis ratios change only the normalization of the mock velocity fields and rotation curves and not their shape. The shape of the modeled NFW rotation curves does not reproduce the data: these potentials are unable to simultaneously bring the mock data at both small and large radii into agreement with observations. Indeed, to match the slow rise of LSB galaxy rotation curves, a specific viewing angle of the non-axisymmetric potential is required. For each of the simulated LSB galaxies, the observer's line-of-sight must be along the minor axis of the potential, an arrangement which is inconsistent with a random distribution of halo orientations on the sky.

We present new data obtained with the Submillimeter Array for a sample of fourteen nearby luminous and ultraluminous infrared galaxies. The galaxies were selected to have luminosity distances D < 200 Mpc and far-infrared luminosities log(L_FIR) > 11.4. The galaxies were observed with spatial resolutions of order 1 kpc in the CO J=3-2, CO J=2-1, 13CO J=2-1, and HCO+ J=4-3 lines as well as the continuum at 880 microns and 1.3 mm. We have combined our CO and continuum data to measure an average gas-to-dust mass ratio of 120 +/- 28 (rms deviation 109) in the central regions of these galaxies, very similar to the value of 150 determined for the Milky Way. This similarity is interesting given the more intense heating from the starburst and possibly accretion activity in the luminous infrared galaxies compared to the Milky Way. We find that the peak H_2 surface density correlates with the far-infrared luminosity, which suggests that galaxies with higher gas surface densities inside the central kiloparsec have a higher star formation rate. The lack of a significant correlation between total H_2 mass and far-infrared luminosity in our sample suggests that the increased star formation rate is due to the increased availability of molecular gas as fuel for star formation in the central regions. In contrast to previous analyses by other authors, we do not find a significant correlation between central gas surface density and the star formation efficiency, as trace by the ratio of far-infrared luminosity to nuclear gas mass. Our data show that it is the star formation rate, not the star formation efficiency, that increases with increasing central gas surface density in these galaxies.

We discuss the near-infrared properties of the nuclei in the 11 merging galaxies of the Toomre sequence, based on high spatial resolution J, H, and K imaging data using NICMOS onboard the Hubble Space Telescope (HST). The observations are less affected by dust extinction than our previous HST/WFPC2 observations and offer higher spatial resolution than existing ground-based near-IR data. We see a marginal trend for the nuclei to become bluer with advancing merger stage, which we attribute to a dispersal of dust at late times in the merging process. Our data also indicate a statistically significant trend for the nuclei in the sequence to become more luminous, within an aperture of fixed physical size and after correcting for dust extinction, with advancing merger stage. We derive K-band surface brightness profiles for those nuclei for which the morphology allows a meaningful isophotal analysis, and fit the profiles with a ``Nuker law'' for comparison with other samples of galaxies observed with HST. The majority of the nuclei have steep profiles that can be characterized as power-law type. In general, the Toomre sequence galaxies tend to have steeper profiles and higher central luminosity surface densities than E/S0's. We derive V-K color profiles for the nuclei to further address this possibility, but find that the large amounts of dust extinction complicate their interpretation. Overall, our results are consistent with the generic predictions of N-body simulations of spiral galaxy mergers. If left to evolve and fade for several Gyrs, it is possible that the properties of the Toomre sequence nuclei would resemble theproperties of the nuclei of normal E/S0 galaxies. Our results therefore support the view that mergers of spiral galaxies can lead to the formation of early-type galaxies (Abridged).

Using N-body simulations, we have modeled the production and evolution of diffuse, low surface brightness intracluster light (ICL) in three simulated galaxy clusters. Using an observational definition of ICL to be luminosity at a surface brightness mu_V>26.5 mag/sq.arcsec, we have found that the fraction of cluster luminosity contained in ICL generally increases as clusters evolve, although there are large deviations from this trend over short timescales, including sustained periods of decreasing ICL luminosity. Most ICL luminosity increases come in short, discrete events which are highly correlated with group accretion events within the cluster. In evolved clusters we find that ~10-15% of the clusters' luminosity is at ICL surface brightness. The morphological structure of the ICL changes with time, evolving from a complex of filaments and small-scale, relatively high surface brightness features early in a cluster's history, to a more diffuse and amorphous cluster-scale ICL envelope at later times. Finally, we also see a correlation between the evolution of ICL at different surface brightnesses, including a time delay between the evolution of faint and extremely faint surface brightness features which is traced to the differing dynamical timescales in the group and cluster environment.

Using N-body simulations, we have modeled the production and evolution of low surface brightness, diffuse intra-cluster light (ICL) in galaxy clusters. By creating simulated observations of the clusters we have measured the evolution of the ICL luminosity throughout the dynamical history of the clusters. We find that ICL production tends to occur in short, discrete events, which correlate very strongly with strong, small-scale interactions and accretions between groups within the clusters.

Hypervelocity stars have been recently discovered in the outskirts of galaxies, such as the unbound star in the Milky Way halo, or the three anomalously fast intracluster planetary nebulae (ICPNe) in the Virgo Cluster. These may have been ejected by close 3-body interactions with a binary supermassive black hole (SMBBH), where a star which passes within the semimajor axis of the SMBBH can receive enough energy to eject it from the system. Stars ejected by SMBBHs may form a significant sub-population with very different kinematics and mean metallicity than the bulk of the intracluster stars. The number, kinematics, and orientation of the ejected stars may constrain the mass ratio, semimajor axis, and even the orbital plane of the SMBBH. We investigate the evolution of the ejected debris from a SMBBH within a clumpy and time-dependent cluster potential using a high resolution, self-consistent cosmological N-body simulation of a galaxy cluster. We show that the predicted number and kinematic signature of the fast Virgo ICPNe is consistent with 3-body scattering by a SMBBH with a mass ratio $10:1$ at the center of M87.

(abridged) We investigate the detailed response of gas to the formation of transient and long-lived dynamical structures induced in the early stages of a disk-disk collision, and identify observational signatures of radial gas inflow through a detailed examination of the collision simulation of an equal mass bulge dominated galaxy. Stars respond to the tidal interaction by forming both transient arms and long lived m=2 bars, but the gas response is more transient, flowing directly toward the central regions within about 10^8 years after the initial collision. The rate of inflow declines when more than half of the total gas supply reaches the inner few kpc, where the gas forms a dense nuclear ring inside the stellar bar. The average gas inflow rate to the central 1.8 kpc is \~7 Msun/yr with a peak rate of 17 Msun/yr. The evolution of gas in a bulgeless progenitor galaxy is also discussed, and a possible link to the ``chain galaxy'' population observed at high redshifts is inferred. The evolution of the structural parameters (the asymmetry and concentration) of both stars and gas are studied in detail. Further, a new structural parameter (the compactness parameter K) that traces the evolution of the size scale of the gas relative to the stellar disk is introduced. Non-circular gas kinematics driven by the perturbation of the non-axisymmetric structure can produce distinct emission features in the "forbidden velocity quadrants'' of the position-velocity diagram (PVD). The dynamical mass calculated using the rotation curve derived from fitting the emission envelope of the PVD can determine the true mass to within 20% to 40%. The evolution of the molecular fraction $M_H2/M_(H2 + HI) and the compactness (K) are potential tracers to quantitatively assign the age of the interaction.

We report the search for intracluster light in four Abell Type II/III (non-cD) galaxy clusters: Abell 801, 1234, 1553, & 1914. We find on average that these clusters contain $\sim$ 10% of their detected stellar luminosity in a diffuse component. We show that for two of the clusters the intracluster light closely follows the galaxy distribution, but in the other two cases, there are noticeable differences between the spatial distribution of the galaxies and the intracluster light. We report the results of a search for intracluster tidal debris in each cluster, and note that Abell 1914 in particular has a number of strong tidal features likely due to its status as a recent cluster merger. One of the Abell 1914 features appears to be spatially coincident with an extension seen in weak lensing maps, implying the feature traces a large amount of mass. We compare these results to numerical simulations of hierarchically-formed galaxy clusters, and find good general agreement between the observed and simulated images, although we also find that our observations sample only the brightest features of the intracluster light. Together, these results suggest that intracluster light can be a valuable tool in determining the evolutionary state of galaxy clusters.

We present HST/WFPC2 observations of the five bluest E+A galaxies (z~0.1) in the Zabludoff et al. sample to study whether their detailed morphologies are consistent with late-to-early type evolution and to determine what drives that evolution. The morphologies of four galaxies are disturbed, indicating that a galaxy-galaxy merger is at least one mechanism that leads to the E+A phase. Two-dimensional image fitting shows that the E+As are generally bulge-dominated systems, even though at least two E+As may have underlying disks. In the Fundamental Plane, E+As stand apart from the E/S0s mainly due to their high effective surface brightness. Fading of the young stellar population and the corresponding increase in their effective radii will cause these galaxies to migrate toward the locus of E/S0s. E+As have profiles qualitatively like those of normal power-law early-type galaxies, but have higher surface brightnesses. This result provides the first direct evidence supporting the hypothesis that power-law ellipticals form via gas-rich mergers. In total, at least four E+As are morphologically consistent with early-type galaxies. We detect compact sources, possibly young star clusters, associated with the galaxies. These sources are much brighter (M_R ~ -13) than Galactic globular clusters, have luminosities consistent with the brightest clusters in nearby starburst galaxies, and have blue colors consistent with the ages estimated from the E+A galaxy spectra (several 10^8 yr). Further study of such young star cluster candidates might provide the elusive chronometer needed to break the age/burst-strength degeneracy for these post-merger galaxies.

We present the results from multiple surveys for intracluster planetary nebulae (IPNe) in nearby galaxy clusters and groups. We find that in the case of clusters, our observations imply: 1) the amount of intracluster starlight is significant, up to 20% of the total starlight, 2) the Virgo Cluster is elongated along our line of sight, and 3) the intracluster light is clustered on the sky, implying ongoing tidal stripping. In contrast, searches for IPNe in groups have found little or no intra-group population, implying there may be something in the cluster environment that significantly enhances intracluster star production. From high-resolution N-body simulations, we find that the IPNe should create observable features in position-velocity space, and that these features may eventually allow us to place limits on the dynamics of galaxy clusters.

We give an update of our ongoing survey for intracluster light (ICL), in a sample of distant Abell clusters. We find that the amount of intracluster starlight is comparable to that seen in nearby clusters, and that tidal debris appears to be common.

We present results from an HST study of the morphology and kinematics of NGC 6240. This merging galaxy with a double nucleus is one of the nearest and best-studied ultraluminous infrared galaxies. HST resolves both nuclei into seperate components. The distance between the northern and southern optical/near-infrared components is greater than that observed in radio and X-ray studies, arguing that even in K-band we may not be seeing all the way through the dust to the true nuclei. The ionized gas does not display rotation around either of the nuclei, or equilibrium motion in general. There is a strong velocity gradient between the nuclei, similar to what is seen in CO data. There is no such gradient in our stellar kinematics. The velocity dispersion of the gas is larger than expected for a cold disk. We also map and model the emission-line velocity field at an off-nuclear position where a steep velocity gradient was previously detected in ground-based data. Overall, the data indicate that line-of-sight projection effects, dust absorption, non-equilibrium merger dynamics, and the possible influence of a wind may be playing an important role in the observed kinematics. Chandra observations of hard X-rays have shown that both of the nuclei contain an Active Galactic Nucleus (AGN). The HST data show no clear sign of the two AGNs: neither continuum nor narrow-band imaging shows evidence for unresolved components in the nuclei, and there are no increased emission line widths or rapid rotation near the nuclei. This underscores the importance of X-ray data for identifying AGNs in highly dust-enshrouded environments.

We report on the properties of nuclear regions in the Toomre Sequence of merging galaxies, based on imaging data gathered with the Hubble Space Telescope WFPC2 camera. We have imaged the 11 systems in the proposed evolutionary merger sequence in the F555W and F814W broad-band filters, and in H-alpha+[NII] narrow-band filters. The broad-band morphology of the nuclear regions varies from non-nucleated starburst clumps through dust-covered nuclei to a nucleated morphology. There is no unambiguous trend in the morphology with merger stage. The emission-line morphology is extended beyond the nucleus in most cases, but centrally concentrated (within 1 kpc) emission-line gas can be seen in the four latest-stage merger systems. We have quantified the intrinsic luminosity densities and colors within the inner 100 pc and 1 kpc of each identified nucleus. We find little evidence for a clear trend in nuclear properties along the merger sequence, other than a suggestive rise in the nuclear luminosity density in the most evolved members of the sequence. The lack of clear trends in nuclear properties is likely due both to the effects of obscuration and geometry, as well as the physical variety of galaxies included in the Toomre Sequence.

We present some early results from our deep imaging survey of galaxy clusters intended to detect and study intracluster light (ICL). From our observations to date, we find that ICL is common in galaxy clusters, and that substructure in the ICL also appears to be common as well. We also discuss some initial comparisons of our imaging results to high-resolution numerical simulations of galaxy clusters, and give avenues for future research.

We report the results of a radial velocity survey of planetary nebulae (PNe) located in the tidal features of the well-known interacting system NGC 5194/95 (M51). We find clear kinematic evidence that M51's northwestern tidal debris consists of two discrete structures which overlap in projection -- NGC 5195's own tidal tail, and diffuse material stripped from NGC 5194. We compare these kinematic data to a new numerical simulation of the M51 system, and show that the data are consistent with the classic ``single passage'' model for the encounter, with a parabolic satellite trajectory and a 2:1 mass ratio. We also comment on the spectra of two unusual objects: a high-velocity PN which may be associated with NGC 5194's halo, and a possible interloping high-redshift galaxy.

We report the initial results of a deep imaging survey of galaxy clusters. The primary goals of this survey are to quantify the amount of intracluster light as a function of cluster properties, and to quantify the frequency of tidal debris. We outline the techniques needed to perform such a survey, and we report findings for the first two galaxy clusters in the survey: Abell 1413, and MKW 7 . These clusters vary greatly in richness and structure. We show that our surface photometry reliably reaches to a surface brightness of \mu_v = 26.5 mags per arcsec. We find that both clusters show clear excesses over a best-fitting r^{1/4} profile: this was expected for Abell 1413, but not for MKW 7. Both clusters also show evidence of tidal debris in the form of plumes and arc-like structures, but no long tidal arcs were detected. We also find that the central cD galaxy in Abell 1413 is flattened at large radii, with an ellipticity of $\approx 0.8$, the largest measured ellipticity of any cD galaxy to date.

We use numerical simulations to study the evolution of triaxial elliptical galaxies with central black holes. In contrast to earlier numerical studies which used galaxy models with central density ``cores,'' our galaxies have steep central cusps, like those observed in real ellipticals. As a black hole grows in these cuspy triaxial galaxies, the inner regions become rounder owing to chaos induced in the orbit families which populate the model. At larger radii, however, the models maintain their triaxiality, and orbital analyses show that centrophilic orbits there resist stochasticity over many dynamical times. While black hole induced evolution is strong in the inner regions of these galaxies, and reaches out beyond the nominal ``sphere of influence'' of a black hole, our simulations do not show evidence for a rapid {\it global} transformation of the host. The triaxiality of observed elliptical galaxies is therefore not inconsistent with the presence of supermassive black holes at their centers.

We present first results from an HST WFPC2 imaging and STIS spectroscopy program to investigate the structural and star forming properties in the nuclear regions of the Toomre Sequence of merging galaxies. Here we discuss V-band, I-band and H-alpha images of the nuclei. We comment briefly on the connection between the nuclear morphology of the ionized gas and the merger stage.

We present the results of HI line and radio continuum observations of five nearby E+A galaxies. These galaxies have spectra that are dominated by a young stellar component but lack the emission lines characteristic of significant, on-going star formation. They are selected from a unique sample of 21 E+A's identified by Zabludoff et al.(1996) in their spectroscopic search for E+A galaxies using the Las Campanas Redshift Survey, where over 11,000 nearby galaxies were examined. The five E+A galaxies span a range of environments: three are in the field and two are in clusters. Only one system was detected in HI emission, the field E+A galaxy EA1, with a total flux of 0.30 +/- 0.02 Jy km/s and an HI mass of (3.5 +/- 0.2) 10^9 h^(-2) M_sun. The HI morphology and kinematics of EA 1 suggest a galaxy-galaxy interaction, with a dynamical age of about 6 x 10^8 h^(-1) yr inferred from the HI tail lengths and velocities. This age estimate is consistent with the interpretation drawn from optical spectroscopy that starbursts in E+A galaxies began (and subsequently ended) within the last 10^9 yr. Our HI detection limits are such that if the other E+A's in our sample had the HI properties of EA 1, we would have detected (or marginally detected) them. We conclude that E+A galaxies have a range of HI properties. None of the galaxies were detected in radio continuum emission, with upper limits to the radio power of about 10^21 h^(-2) W/Hz. Our limits exclude the possibility that these E+A's are dust-enshrouded massive starburst galaxies, but are insufficient to exclude modest star formation rates of less than a few h^(-2) M_sun per yr.

We construct numerical models of mildly triaxial elliptical galaxies with central density cusps. Using a technique we call ``adiabatic squeezing'', we begin with a spherical gamma=1 Hernquist model and apply a drag to the velocities of the particles along each principle axis. The final models are stable in isolation, preserving their density structure and figure shape over many dynamical timescales. The density profile and axial ratios compare well to the observed properties of elliptical galaxies. The orbital structure of these models show a mixture of tubes, boxes, and boxlets, as expected for triaxial systems, with very few chaotic orbits. These N-body realizations of cuspy triaxial galaxies provide a basis for the study of the dynamical evolution of elliptical galaxies.

Using high resolution N-body simulations with hydrodynamics and star formation, we investigate interactions and the resulting starbursts in galaxies with properties typical of $z\sim 3$. We apply spectral population models to produce mock-HST images, and discuss the observed magnitude, color, and morphological appearance of our simulated galaxies in both the rest-UV and rest-visual bands.

We present an extensive parameter survey to study the influence of halo mass profiles on the development of tidal tails in interacting disk galaxies. We model the galaxies using a fixed exponential disk with a central bulge and vary the halo potential over a range of parameters using both the Hernquist and NFW mass distributions, probing the effect of the halo mass, extent and concentration. We examine the consistency of the results against both observational and theoretical constraints on halo profiles and comment on the failures and weaknesses of different models. Galaxies with rising or flat rotation curves dominated by the halo are inhibited from forming tidal tails unless the halo is abruptly cut off just beyond the disk edge. Conversely, models with declining rotation curves -- resulting either from compact, low mass halos, or from massive disk components in low concentration dark halos -- produce tidal tails very similar to those observed in well-studied interacting systems. As argued by Springel & White (1998), a unifying, quantitative relation for all cases is that the ratio of the escape velocity to the circular velocity at around the solar radius must be v_e/v_c<2.5 for tidal tails to be produced. The galaxy models which appear to fit most of the observational constraints are those which have disk-dominated rotation curves and low concentration halos. We discuss our results in a cosmological context using recent studies which link halo properties to cosmological models.

We perform fully self-consistent stellar dynamical simulations of the accretion of a companion ("satellite") galaxy by a large disk galaxy to investigate the interaction between the disk, halo, and satellite components of the system during a merger. Our fiducial encounter begins with a satellite in a prograde, circular orbit inclined thirty degrees with respect to the disk plane at a galactocentric distance of six disk scalelengths. The satellite's mass is 10% of the disk's mass and its half-mass radius is about 1.3 kpc. The system is modelled with 500 000 particles, sufficient to mitigate numerical relaxation noise over the merging time. The satellite sinks in only ~1 Gyr and a core containing ~45% of its initial mass reaches the centre of the disk. With so much of the satellite's mass remaining intact, the disk sustains significant damage as the satellite passes through. At the solar circle we find that the disk thickens ~60%, the velocity dispersions increase by $\Delta\mbox{\boldmath$\sigma$} \simeq (10,8,8)$ km/s to $(\sigma_R, \sigma_\phi, \sigma_z) \simeq (48, 42, 38)$ km/s, and the asymmetric drift is unchanged at ~18 km/s. Although the disk is not destroyed by these events (hence "minor" mergers), its final state resembles a disk galaxy of earlier Hubble type than its initial state, thicker and hotter, with the satellite's core enhancing the bulge. Thus minor mergers continue to be a promising mechanism for driving galaxy evolution.

We use simulations of merging galaxies to explore the sensitivity of the morphology of tidal tails to variations of the halo mass distributions in the parent galaxies. Our goal is to constrain the mass of dark halos in well-known merging pairs. We concentrate on prograde encounters between equal mass galaxies which represent the best cases for creating tidal tails, but also look at systems with different relative orientations, orbital energies and mass ratios. As the mass and extent of the dark halo increase in the model galaxies, the resulting tidal tails become shorter and less massive, even under the most favorable conditions for producing these features. Our simulations imply that the observed merging galaxies with long tidal tails ($\sim 50-100$ kpc) such as NGC 4038/39 (the Antennae) and NGC 7252 probably have halo:disk+bulge mass ratios less than 10:1. These results conflict with the favored values of the dark halo mass of the Milky Way derived from satellite kinematics and the timing argument which give a halo:disk+bulge mass ratio of $\sim 30:1$. However, the lower bound of the estimated dark halo mass in the Milky Way (mass ratio $\sim 10:1$) is still consistent with the inferred tidal tail galaxy masses. Our results also conflict with the expectations of $\Omega=1$ cosmologies such as CDM which predict much more massive and extended dark halos.