When two black holes merge, the asymmetric emission of gravitational waves provides an impulse to the merged system; this gravitational wave recoil velocity can be up to 4000 km s$^{-1}$, easily fast enough for the black hole to escape its host galaxy. We combine semi-analytic modeling with cosmological zoom-in simulations of a Milky Way-type galaxy to investigate the role of black hole spin and gravitational recoil in the epoch of massive black hole seeding. We sample four different spin distributions (random, aligned, anti-aligned, and zero spin), and compare the resulting merger rates, occupation fractions, and MBH-host relations with what is expected by excluding the effect of recoil. The inclusion of gravitational recoil and MBH spin in the assembly of MBH seeds can reduce the final $z=5$ MBH mass by up to an order of magnitude. The MBH occupation fraction, however, remains effectively unaltered due to episodes of black hole formation following a recoil event. While electromagnetic detections of these events are unlikely, LISA is ideally suited to detect gravitational wave signals from such events.

Gravitational wave measurements provide the most robust constraints on the mass of astrophysical black holes. Using state-of-the-art parameter estimation techniques and gravitational waveform models, we infer the source parameters of the loudest marginal trigger, 170502, found in LIGO from 2015-2017. If this trigger is assumed to be a binary black hole merger, we find it corresponds to a total mass in the source frame of $157^{+55}_{-41}~\rm{M}_\odot$. This is a factor of $\sim3/2$ higher than that of the heaviest confirmed detection, GW170729. The primary and secondary black hole masses are constrained to $94^{+44}_{-28}~\rm{M}_{\odot}$ and $62^{+30}_{-25}~\rm{M}_{\odot}$ respectively with 90% confidence. We find $\chi_\mathrm{eff}= 0.49^{+0.31}_{-0.63}$, the effective spin aligned with the total angular momentum vector, which has most of its support at larger values than the other O1/O2 events. Further, we find the inclusion of higher order harmonics narrows the confidence region for the total binary mass by 10%. Our study illustrates the necessary techniques for inferring the physical parameters of intermediate mass black hole binary candidates $(\gtrsim 100~\mathrm{M}_\odot)$ in the current gravitational wave network.

Harnessing the sheer discovery potential of gravitational wave astronomy will require bold, deliberate, and sustained efforts to train and develop the requisite workforce. The next decade requires a strategic plan to build -- from the ground up -- a robust, open, and well-connected gravitational wave astronomy community with deep participation from traditional astronomers, physicists, data scientists, and instrumentalists. This basic infrastructure is sorely needed as an enabling foundation for research. We outline a set of recommendations for funding agencies, universities, and professional societies to help build a thriving, diverse, and inclusive new field.

Understanding how supermassive black holes (SMBHs) pair and merge helps to inform predictions of off-center, dual, and binary AGN, and provides key insights into how SMBHs grow and co-evolve with their galaxy hosts. As the loudest known gravitational wave source, binary SMBH mergers also hold centerstage for the Laser Interferometer Space Antenna (LISA), a joint ESA/NASA gravitational wave observatory set to launch in 2034. Here, we continue our work to characterize SMBH binary formation and evolution through increasingly more realistic high resolution direct $N$-body simulations, focusing on the effect of SMBH mass ratio, orientation, and eccentricity within a rotating and flattened stellar host. During the dynamical friction phase, we found a prolonged orbital decay for retrograde SMBHs and swift pairing timescales for prograde SMBHs compared to their counterparts in non-rotating models, an effect that becomes more pronounced for smaller mass ratios $M_{\rm sec}/M_{\rm prim} = q$. During this pairing phase, the eccentricity dramatically increases for retrograde configurations, but as the binary forms, the orbital plane flips so that it is almost perfectly prograde, which stifles the rapid eccentricity growth. In prograde configurations, SMBH binaries form and remain at comparatively low eccentricities. As in our prior work, we note that the center of mass of a prograde SMBH binary itself settles into an orbit about the center of the galaxy. Since even the initially retrograde binaries flip their orbital plane, we expect few binaries in rotating systems to reside at rest in the dynamic center of the host galaxy, though this effect is smaller as $q$ decreases.

We study the present-day connection between galaxy morphology and angular momentum using the {\sc Dark Sage} semi-analytic model of galaxy formation. For galaxies between $ 10^{11}-10^{12} \mathrm{M}_{\odot}$ in stellar mass, the model successfully predicts the observed trend whereby galaxies with more prominent disks exhibit higher {\em stellar} disk specific angular momentum ($j_{\rm stellar, disk}$) at fixed stellar mass. However, when we include the gas in the disk, bulge-dominated galaxies have the highest {\em total} disk specific angular momentum ($j_{\rm total, disk}$). We attribute this to a large contribution from an extended disk of cold gas in typical bulge-dominated galaxies. We find the relationship between $j_{\rm dark matter}$ and morphology to be quite complex. Surprisingly, in this stellar mass range, not only do bulge-dominated galaxies tend to live in halos with higher $j_{\rm dark matter}$ than disk-dominated galaxies, but intermediate galaxies (those with roughly equal fractions of bulge and disk mass) have the lowest $j_{\rm dark matter}$ of all. Yet, when controlling for halo mass, rather than stellar mass, the relationship between $j_{\rm dark matter}$ and morphology vanishes. Based on these results, halo mass rather than angular momentum is the main driver of the predicted morphology sequence at high masses. In fact, in our stellar mass range, disk-dominated galaxies live in dark matter halos that are roughly 1/10th the mass of their bulge-dominated counterparts.

The first terrestrial gravitational wave interferometers have dramatically underscored the scientific value of observing the Universe through an entirely different window, and of folding this new channel of information with traditional astronomical data for a multimessenger view. The Laser Interferometer Space Antenna (LISA) will broaden the reach of gravitational wave astronomy by conducting the first survey of the millihertz gravitational wave sky, detecting tens of thousands of individual astrophysical sources ranging from white-dwarf binaries in our own galaxy to mergers of massive black holes at redshifts extending beyond the epoch of reionization. These observations will inform - and transform - our understanding of the end state of stellar evolution, massive black hole birth, and the co-evolution of galaxies and black holes through cosmic time. LISA also has the potential to detect gravitational wave emission from elusive astrophysical sources such as intermediate-mass black holes as well as exotic cosmological sources such as inflationary fields and cosmic string cusps.

The Laser Interferometer Space Antenna (LISA) will open three decades of gravitational wave (GW) spectrum between 0.1 and 100 mHz, the mHz band. This band is expected to be the richest part of the GW spectrum, in types of sources, numbers of sources, signal-to-noise ratios and discovery potential. When LISA opens the low-frequency window of the gravitational wave spectrum, around 2034, the surge of gravitational-wave astronomy will strongly compel a subsequent mission to further explore the frequency bands of the GW spectrum that can only be accessed from space. The 2020s is the time to start developing technology and studying mission concepts for a large-scale mission to be launched in the 2040s. The mission concept would then be proposed to Astro2030. Only space based missions can access the GW spectrum between 10 nHz and 1 Hz because of the Earths seismic noise. This white paper surveys the science in this band and mission concepts that could accomplish that science. The proposed small scale activity is a technology development program that would support a range of concepts and a mission concept study to choose a specific mission concept for Astro2030. In this white paper, we will refer to a generic GW mission beyond LISA as bLISA.

The new field of gravitational wave astrophysics requires a growing pool of students and researchers with unique, interdisciplinary skill sets. It also offers an opportunity to build a diverse, inclusive astronomy community from the ground up. We describe the efforts used by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) NSF Physics Frontiers Center to foster such growth by involving students at all levels in low-frequency gravitational wave astrophysics with pulsar timing arrays (PTAs) and establishing collaboration policies that ensure broad participation by diverse groups. We describe and illustrate the impact of these techniques on our collaboration as a case study for other distributed collaborations.

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.

A space-based interferometer operating in the previously unexplored mHz gravitational band has tremendous discovery potential. If history is any guide, the opening of a new spectral band will lead to the discovery of entirely new sources and phenomena. The mHz band is ideally suited to exploring beyond standard model processes in the early universe, and with the sensitivities that can be reached with current technologies, the discovery space for exotic astrophysical systems is vast.

The discoveries made over the past 20 years by Chandra and XMM-Newton surveys in conjunction with multiwavelength imaging and spectroscopic data available in the same fields have significantly changed the view of the supermassive black hole (SMBH) and galaxy connection. These discoveries have opened up several exciting questions that are beyond the capabilities of current X-ray telescopes and will need to be addressed by observatories in the next two decades. As new observatories peer into the early Universe, we will begin to understand the physics and demographics of SMBH infancy (at $z>6$) and investigate the influence of their accretion on the formation of the first galaxies ($\S$ 2.1). We will also be able to understand the accretion and evolution over the cosmic history (at $z\sim$1-6) of the full population of black holes in galaxies, including low accretion rate, heavily obscured AGNs at luminosities beyond the reach of current X-ray surveys ($\S$2.2 and $\S$2.3), enabling us to resolve the connection between SMBH growth and their environment.

Observational evidence has been mounting for the existence of intermediate mass black holes (IMBHs, 10^2-10^5 Msun), but observing them at all, much less constraining their masses, is very challenging. In one theorized formation channel, IMBHs are the seeds for supermassive black holes in the early universe. As a result, IMBHs are predicted to exist in the local universe in dwarf galaxies, as well as wandering in more massive galaxy halos. However, these environments are not conducive to the accretion events or dynamical signatures that allow us to detect IMBHs. The Laser Interferometer Space Antenna (LISA) will demystify IMBHs by detecting the mergers of these objects out to extremely high redshifts, while measuring their masses with extremely high precision. These observations of merging IMBHs will allow us to constrain the formation mechanism and subsequent evolution of massive black holes, from the 'dark ages' to the present day, and reveal the role that IMBHs play in hierarchical galaxy evolution.

We have compelling evidence for stellar-mass black holes (BHs) of ~5-80 M_sun that form through the death of massive stars. We also have compelling evidence for so-called supermassive BHs (10^5-10^10 M_sun) that are predominantly found in the centers of galaxies. We have very good reason to believe there must be BHs with masses in the gap between these ranges: the first ~10^9 M_sun BHs are observed only hundreds of millions of years after the Big Bang, and all theoretically viable paths to making supermassive BHs require a stage of "intermediate" mass. However, no BHs have yet been reliably detected in the 100-10}^5 M_sun mass range. Uncovering these intermediate-mass BHs of 10^3-10^5 M_sun is within reach in the coming decade. In this white paper we highlight the crucial role that 30-m class telescopes will play in dynamically detecting intermediate-mass black holes, should they exist.

Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumental in identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) observations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with their host galaxies, the dynamical interactions between binaries and their galactic environments, and the fundamental physics of accretion. Multi-messenger observations can also make SMBHBs 'standard sirens' for cosmological distance measurements out to $z\simeq0.5$. LIGO has already ushered in breakthrough insights in our knowledge of black holes. The multi-messenger detection of SMBHBs with PTAs will be a breakthrough in the years $2020-2030$ and beyond, and prepare us for LISA to help complete our views of black hole demographics and evolution at higher redshifts.

Coalescing, massive black-hole (MBH) binaries are the most powerful sources of gravitational waves (GWs) in the Universe, which makes MBH science a prime focus for ongoing and upcoming GW observatories. The Laser Interferometer Space Antenna (LISA) -- a gigameter scale space-based GW observatory -- will grant us access to an immense cosmological volume, revealing MBHs merging when the first cosmic structures assembled in the Dark Ages. LISA will unveil the yet unknown origin of the first quasars, and detect the teeming population of MBHs of $10^4 - 10^7$ solar masses. forming within protogalactic halos. The Pulsar Timing Array, a galactic-scale GW survey, can access the largest MBHs the Universe, detecting the cosmic GW foreground from inspiraling MBH binaries of about 10^9 solar masses. LISA can measure MBH spins and masses with precision far exceeding that from electromagnetic (EM) probes, and together, both GW observatories will provide the first full census of binary MBHs, and their orbital dynamics, across cosmic time. Detecting the loud gravitational signal of these MBH binaries will also trigger alerts for EM counterpart searches, from decades (PTAs) to hours (LISA) prior to the final merger. By witnessing both the GW and EM signals of MBH mergers, precious information will be gathered about the rich and complex environment in the aftermath of a galaxy collision. The unique GW characterization of MBHs will shed light on the deep link between MBHs of $10^4-10^{10}$ solar masses and the grand design of galaxy assembly, as well as on the complex dynamics that drive MBHs to coalescence.

Ultra-compact binaries (UCBs) are systems containing compact or degenerate stars with orbital periods less than one hour. Tens of millions of UCBs are predicted to exist within theGalaxy emitting gravitational waves (GWs) at mHz frequencies. Combining GW searches with electromagnetic (EM) surveys like Gaia and LSST will yield a comprehensive, multimessenger catalog of UCBs in the galaxy. Joint EM and GW observations enable measurements of masses, radii, and orbital dynamics far beyond what can be achieved by independent EM or GW studies. GW+EM surveys of UCBs in the galaxy will yield a trove of unique insights into the nature of white dwarfs, the formation of compact objects, dynamical interactions in binaries, and energetic, accretion-driven phenomena like Type Ia superonovae.

LISA will open the mHz band of gravitational waves (GWs) to the astronomy community. The strong gravity which powers the variety of GW sources in this band is also crucial in a number of important astrophysical processes at the current frontiers of astronomy. These range from the beginning of structure formation in the early universe, through the origin and cosmic evolution of massive black holes in concert with their galactic environments, to the evolution of stellar remnant binaries in the Milky Way and in nearby galaxies. These processes and their associated populations also drive current and future observations across the electromagnetic (EM) spectrum. We review opportunities for science breakthroughs, involving either direct coincident EM+GW observations, or indirect multimessenger studies. We argue that for the US community to fully capitalize on the opportunities from the LISA mission, the US efforts should be accompanied by a coordinated and sustained program of multi-disciplinary science investment, following the GW data through to its impact on broad areas of astrophysics. Support for LISA-related multimessenger observers and theorists should be sized appropriately for a flagship observatory and may be coordinated through a dedicated mHz GW research center.

Electromagnetic observations of the sky have been the basis for our study of the Universe for millennia, cosmic ray studies are now entering their second century, the first neutrinos from an astrophysical source were identified three decades ago, and gravitational waves were directly detected only four years ago. Detections of these messengers are now common. Astrophysics will undergo a revolution in the 2020s as multimessenger detections become routine. The 8th Astro2020 Thematic Area is Multimessenger Astronomy and Astrophysics, which includes the identification of the sources of gravitational waves, astrophysical and cosmogenic neutrinos, cosmic rays, and gamma-rays, and the coordinated multimessenger and multiwavelength follow-ups. Identifying and characterizing multimessenger sources enables science throughout and beyond astrophysics. Success in the multimessenger era requires: (i) sensitive coverage of the non-electromagnetic messengers, (ii) full coverage of the electromagnetic spectrum, with either fast-response observations or broad and deep high-cadence surveys, and (iii) improved collaboration, communication, and notification platforms.

The identification of the electromagnetic counterpart of the gravitational wave event, GW170817, and discovery of neutrinos and gamma-rays from TXS 0506+056 heralded the new era of multi-messenger astrophysics. As the number of multi-messenger events rapidly grow over the next decade, the cyberinfrastructure requirements to handle the increase in data rates, data volume, need for event follow up, and analysis across the different messengers will also explosively grow. The cyberinfrastructure requirements to enhance multi-messenger astrophysics will both be a major challenge and opportunity for astronomers, physicists, computer scientists and cyberinfrastructure specialists. Here we outline some of these requirements and argue for a distributed cyberinfrastructure institute for multi-messenger astrophysics to meet these challenges.

While still in its infancy, the budding field of gravitational-wave astronomy has so far exceeded most expectations, and the achievements that have already been made bode well for the decade to come. While the discoveries made possible by LIGO have captured the imagination of experts and nonexperts alike, it is important when looking ahead to consider those discoveries in the context of the field as a whole. Just as radio, optical, and x-ray radiation probe different physical phenomena occurring on a range of length and energy scales, the future of gravitational-wave astrophysics depends on our ability to open up the entire spectrum. We will describe the scientific prospects for the field of gravitational-wave astronomy as a whole as we enter the coming decade, and we will place the specific contributions from a future space-based gravitational-wave observatory within this context.

There are two big questions cosmologists would like to answer -- How does the Universe work, and what are its origin and destiny? A long wavelength gravitational wave detector -- with million km interferometer arms, achievable only from space -- gives a unique opportunity to address both of these questions. A sensitive, mHz frequency observatory could use the inspiral and merger of massive black hole binaries as standard sirens, extending our ability to characterize the expansion history of the Universe from the onset of dark energy-domination out to a redshift z ~ 10. A low-frequency detector, furthermore, offers the best chance for discovery of exotic gravitational wave sources, including a primordial stochastic background, that could reveal clues to the origin of our Universe.

The LIGO/Virgo gravitational-wave (GW) interferometers have to-date detected ten merging black hole (BH) binaries, some with masses considerably larger than had been anticipated. Stellar-mass BH binaries at the high end of the observed mass range (with "chirp mass" ${\cal M} \gtrsim 25 M_{\odot}$) should be detectable by a space-based GW observatory years before those binaries become visible to ground-based GW detectors. This white paper discusses some of the synergies that result when the same binaries are observed by instruments in space and on the ground. We consider intermediate-mass black hole binaries (with total mass $M \sim 10^2 -10^4 M_{\odot}$) as well as stellar-mass black hole binaries. We illustrate how combining space-based and ground-based data sets can break degeneracies and thereby improve our understanding of the binary's physical parameters. While early work focused on how space-based observatories can forecast precisely when some mergers will be observed on the ground, the reverse is also important: ground-based detections will allow us to "dig deeper" into archived, space-based data to confidently identify black hole inspirals whose signal-to-noise ratios were originally sub-threshold, increasing the number of binaries observed in both bands by a factor of $\sim 4 - 7$.

Intermediate-mass black holes (IMBHs) have masses of about 100 to 100,000 solar masses. They remain elusive. Observing IMBHs in present-day globular clusters (GCs) would validate a formation channel for seed black holes in the early universe and inform event predictions for gravitational wave facilities. Reaching a large number of GCs per galaxy is key, as models predict that only a few percent will have retained their gravitational-wave fostering IMBHs. Related, many galaxies will need to be examined to establish a robust sample of IMBHs in GCs. These needs can be meet by using a next-generation Very Large Array (ngVLA) to search for IMBHs in the GCs of hundreds of galaxies out to a distance of 25 Mpc. These galaxies hold tens of thousands of GCs in total. We describe how to convert an ngVLA signal from a GC to an IMBH mass according to a semi-empirical accretion model. Simulations of gas flows in GCs would help to improve the robustness of the conversion. Also, self-consistent dynamical models of GCs, with stellar and binary evolution in the presence of IMBHs, would help to improve IMBH retention predictions for present-day GCs.

The inspiral of a stellar-mass compact object into a massive ($\sim 10^{4}$-$10^{7} M_{\odot}$) black hole produces an intricate gravitational-wave signal. Due to the extreme-mass ratios involved, these systems complete $\sim 10^{4}$-$10^{5}$ orbits, most of them in the strong-field region of the massive black hole, emitting in the frequency range $\sim10^{-4}-1~$Hz. This makes them prime sources for the space-based observatory LISA (Laser Interferometer Space Antenna). LISA observations will enable high-precision measurements of the physical characteristics of these extreme-mass-ratio inspirals (EMRIs): redshifted masses, massive black hole spin and orbital eccentricity can be determined with fractional errors $\sim 10^{-4}$-$10^{-6}$, the luminosity distance with better than $\sim 10\%$ precision, and the sky localization to within a few square degrees. EMRIs will provide valuable information about stellar dynamics in galactic nuclei, as well as precise data about massive black hole populations, including the distribution of masses and spins. They will enable percent-level measurements of the multipolar structure of massive black holes, precisely testing the strong-gravity properties of their spacetimes. EMRIs may also provide cosmographical data regarding the expansion of the Universe if inferred source locations can be correlated with galaxy catalogs.

Low-frequency gravitational-wave astronomy can perform precision tests of general relativity and probe fundamental physics in a regime previously inaccessible. A space-based detector will be a formidable tool to explore gravity's role in the cosmos, potentially telling us if and where Einstein's theory fails and providing clues about some of the greatest mysteries in physics and astronomy, such as dark matter and the origin of the Universe.

The existence of ultra-faint dwarf (UFD) galaxies highlights the need to push our theoretical understanding of galaxies to extremely low mass. We examine the formation of UFDs by twice running a fully cosmological simulations of dwarf galaxies, but varying star formation. One run uses a temperature-density threshold for star formation, while the other uses an H$_{2}$-based sub-grid star formation model. The total number of dwarf galaxies that forms is different by a factor of 2 between the two runs, but most of these are satellites, leading to a factor of 5 difference in the number of luminous UFD companions around more massive, isolated dwarfs. The first run yields a 47\% chance of finding a satellite around a M$_{halo}$ $\sim 10^{10}$ M$_{\odot}$ host, while the H$_2$ run predicts only a 16\% chance. Metallicity is the primary physical parameter that creates this difference. As metallicity decreases, the formation of H$_2$ is slowed and relegated to higher-density material. Thus, our H$_2$ run is unable to form many (and often, any) stars before reionization removes gas. These results emphasize that predictions for UFD properties made using hydrodynamic simulations, in particular regarding the frequency of satellites around dwarf galaxies, the slope of the stellar mass function at low masses, as well as the properties of ultra-faint galaxies occupying the smallest halos, are extremely sensitive to the subgrid physics of star formation contained within the simulation. However, upcoming discoveries of ultra-faint dwarfs will provide invaluable constraining power on the physics of the first star formation.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

Cunha et al. (2018) recently reexamined the possibility of detecting gravitational waves from exoplanets, claiming that three ultra-short period systems would be observable by LISA. We revisit their analysis and conclude that the currently known exoplanetary systems are unlikely to be detectable, even assuming a LISA observation time $T_{\rm obs}=4$ yrs. Conclusive statements on the detectability of one of these systems, GP Com b, will require better knowledge of the system's properties, as well as more careful modeling of both LISA's response and the galactic confusion noise. Still, the possibility of exoplanet detection with LISA is interesting enough to warrant further study, as gravitational waves could yield dynamical properties that are difficult to constrain with electromagnetic observations.

The past year has witnessed discovery of the first identified counterparts to a gravitational wave transient (GW 170817A) and a very high-energy neutrino (IceCube-170922A). These source identifications, and ensuing detailed studies, have realized longstanding dreams of astronomers and physicists to routinely carry out observations of cosmic sources by other than electromagnetic means, and inaugurated the era of "multi-messenger" astronomy. While this new era promises extraordinary physical insights into the universe, it brings with it new challenges, including: highly heterogeneous, high-volume, high-velocity datasets; globe-spanning cross-disciplinary teams of researchers, regularly brought together into transient collaborations; an extraordinary breadth and depth of domain-specific knowledge and computing resources required to anticipate, model, and interpret observations; and the routine need for adaptive, distributed, rapid-response observing campaigns to fully exploit the scientific potential of each source. We argue, therefore, that the time is ripe for the community to conceive and propose an Institute for Multi-Messenger Astrophysics that would coordinate its resources in a sustained and strategic fashion to efficiently address these challenges, while simultaneously serving as a center for education and key supporting activities. In this fashion, we can prepare now to realize the bright future that we see, beyond, through these newly opened windows onto the universe.

We study the birth of supermassive black holes from the direct collapse process and characterize the sites where these black hole seeds form. In the pre-reionization epoch, molecular hydrogen (H$_2$) is an efficient coolant, causing gas to fragment and form Population III stars, but Lyman-Werner radiation can suppress H$_2$ formation and allow gas to collapse directly into a massive black hole. The critical flux required to inhibit H$_2$ formation, $J_{\rm crit}$, is hotly debated, largely due to the uncertainties in the source radiation spectrum, H$_2$ self-shielding, and collisional dissociation rates. Here, we test the power of the direct collapse model in a self-consistent, time-dependant, non-uniform Lyman-Werner radiation field -- the first time such has been done in a cosmological volume -- using an updated version of the SPH+N-body tree code Gasoline with H$_2$ non-equilibrium abundance tracking, H$_2$ cooling, and a modern SPH implementation. We vary $J_{\rm crit} $ from $30$ to $10^3$ in units of $J_{21}$ to study how this parameter impacts the number of seed black holes and the type of galaxies which host them. We focus on black hole formation as a function of environment, halo mass, metallicity, and proximity of the Lyman-Werner source. Massive black hole seeds form more abundantly with lower $J_{\rm crit}$ thresholds, but regardless of $J_{\rm crit}$, these seeds typically form in halos that have recently begun star formation. Our results do not confirm the proposed atomic cooling halo pair scenario; rather black hole seeds predominantly form in low-metallicity pockets of halos which already host star formation.

Using a new, high-resolution cosmological hydrodynamic simulation of a Milky Way-type (MW-type) galaxy, we explore how a merger-rich assembly history affects the mass budget of the central supermassive black hole (SMBH). We examine a MW-mass halo at the present epoch whose evolution is characterized by several major mergers to isolate the importance of merger history on black hole accretion. This study is an extension of Bellovary et. al. 2013, which analyzed the accretion of high mass, high redshift galaxies and their central black holes, and found that the gas content of the central black hole reflects what is accreted by the host galaxy halo. In this study, we find that a merger-rich galaxy will have a central SMBH preferentially fed by merger gas. Moreover, we find that nearly 30$\%$ of the accreted mass budget of the SMBH enters the galaxy through the two major mergers in its history, which may account for the increase of merger-gas fueling the SMBH. Through an investigation of the angular momentum of the gas entering the host and its SMBH, we determine that merger gas enters the galaxy with lower angular momentum compared to smooth accretion, partially accounting for the preferential fueling witnessed in the SMBH. In addition, the presence of mergers, particularly major mergers, also helps funnel low angular momentum gas more readily to the center of the galaxy. Our results imply that galaxy mergers play an important role in feeding the SMBH in MW-type galaxies with merger-rich histories. Our results imply that galaxy mergers play an important role in feeding the SMBH in MW-type galaxies with merger-rich histories.

It has been well-established that particular centrophilic orbital families in non-spherical galaxies can, in principle, drive a black hole binary to shrink its orbit through three-body scattering until the black holes are close enough to strongly emit gravitational waves. Most of these studies rely on orbital analysis of a static SMBH-embedded galaxy potential to support this view; it is not clear, however, how these orbits transform as the second SMBH enters the center, so our understanding of which orbits actually interact with a SMBH binary is not ironclad. Here, we analyze two flattened galaxy models, one with a single SMBH and one with a binary, to determine which orbits actually do interact with the SMBH binary and how they compare with the set predicted in single SMBH-embedded models. We find close correspondence between the centrophilic orbits predicted to interact with the binary and those that are actually scattered by the binary, in terms of energy and Lz distribution, where Lz is the z component of a stellar particle's angular momentum. Of minor note: because of the larger mass, the binary SMBH has a radius of influence about 4 times larger than in the single SMBH model, which allows the binary to draw from a larger reservoir of orbits to scatter. Of the prediction particles and scattered particles, nearly half have chaotic orbits, 40% have fx:fy=1:1 orbits, 10% have other resonant orbits.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

We present Chandra and XMM-Newton X-ray observations of the Abell 3391/Abell 3395 intercluster filament. It has been suggested that the galaxy clusters Abell 3395, Abell 3391, and the galaxy group ESO-161 located between the two clusters, are in alignment along a large-scale intercluster filament. We find that the filament is aligned close to the plane of the sky, in contrast to previous results. We find a global projected filament temperature kT = $4.45_{-0.55}^{+0.89}$~keV, electron density $n_e=1.08^{+0.06}_{-0.05} \times 10^{-4}$~cm$^{-3}$, and $M_{\rm gas} = 2.7^{+0.2}_{-0.1} \times 10^{13}$~M$_\odot$. The thermodynamic properties of the filament are consistent with that of intracluster medium (ICM) of Abell 3395 and Abell 3391, suggesting that the filament emission is dominated by ICM gas that has been tidally disrupted during an early stage merger between these two clusters. We present temperature, density, entropy, and abundance profiles across the filament. We find that the galaxy group ESO-161 may be undergoing ram pressure stripping in the low density environment at or near the virial radius of both clusters due to its rapid motion through the filament.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

Supermassive black hole (SMBH) binaries residing at the core of merging galaxies are recently found to be strongly affected by the rotation of their host galaxies. The highly eccentric orbits that form when the host is counterrotating emit strong bursts of gravitational waves that propel rapid SMBH binary coalescence. Most prior work, however, focused on planar orbits and a uniform rotation profile, an unlikely interaction configuration. However, the coupling between rotation and SMBH binary evolution appears to be such a strong dynamical process that it warrants further investigation. This study uses direct N-body simulations to isolate the effect of galaxy rotation in more realistic interactions. In particular, we systematically vary the SMBH orbital plane with respect to the galaxy rotation axis, the radial extent of the rotating component, and the initial eccentricity of the SMBH binary orbit. We find that the initial orbital plane orientation and eccentricity alone can change the inspiral time by an order of magnitude. Because SMBH binary inspiral and merger is such a loud gravitational wave source, these studies are critical for the future gravitational wave detector, LISA, an ESA/NASA mission currently set to launch by 2034.

Following the selection of The Gravitational Universe by ESA, and the successful flight of LISA Pathfinder, the LISA Consortium now proposes a 4 year mission in response to ESA's call for missions for L3. The observatory will be based on three arms with six active laser links, between three identical spacecraft in a triangular formation separated by 2.5 million km. LISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using Gravitational Waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the infant Universe at TeV energy scales, has known sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales near the horizons of black holes, all the way to cosmological scales. The LISA mission will scan the entire sky as it follows behind the Earth in its orbit, obtaining both polarisations of the Gravitational Waves simultaneously, and will measure source parameters with astrophysically relevant sensitivity in a band from below $10^{-4}\,$Hz to above $10^{-1}\,$Hz.

We measure the two-point correlation function of G-dwarf stars within 1-3 kpc of the Sun in multiple lines-of-sight using the Schlesinger et al. G-dwarf sample from the SDSS SEGUE survey. The shapes of the correlation functions along individual SEGUE lines-of-sight depend sensitively on both the stellar-density gradients and the survey geometry. We fit smooth disk galaxy models to our SEGUE clustering measurements, and obtain strong constraints on the thin- and thick-disk components of the Milky Way. Specifically, we constrain the values of the thin- and thick-disk scale heights with 3% and 2% precision, respectively, and the values of the thin- and thick-disk scale lengths with 20% and 8% precision, respectively. Moreover, we find that a two-disk model is unable to fully explain our clustering measurements, which exhibit an excess of clustering at small scales (< 50 pc). This suggests the presence of small-scale substructure in the disk system of the Milky Way.

In the second paper of this series, we present results from cosmological simulations on the demographics of flyby interactions to gauge their potential impact on galaxy evolution. In a previous paper, we demonstrated that flybys -- an interaction where two independent halos inter-penetrate but detach at a later time and do not merge -- occur much more frequently than previously believed. In particular, we found that the frequency of flybys increases at low redshift and is comparable to or even greater than the frequency of mergers for halos $\gtrsim 10^{11} M_\odot/h$. In this paper, we classify flybys according to their orbits and the level of perturbation exacted on both the halos involved. We find that the majority of flybys penetrate deeper than $\sim R_{half}$ of the primary and have an initial relative speed $\sim 1.6\times V_{vir}$ of the primary. The typical flyby mass-ratio is $\sim 0.1$ at high $z$ for all halos, while at low $z$, massive primary halos undergo flybys with small secondary halos. We estimate the perturbation from the flyby on both the primary and the secondary and find that a typical flyby is mostly non-perturbative for the primary halo. However, since a massive primary experiences so many flybys at any given time, they are nearly continually a victim of a perturbative event. In particular, we find flybys that cause $\sim 1\%$ change in the binding energy of a primary halo occurs $\gtrsim 1 $ Gyr$^{-1}$ for halos $> 10^{10} M_\odot/h$ for $z \lesssim 4$. Secondary halos, on the other hand, are highly perturbed by the typical encounter, experiencing a change in binding energy of nearly order unity. Our results imply that flybys can drive a significant part of galaxy transformation at moderate to lower redshifts ($z \lesssim 4$). We touch on implications for observational surveys, mass-to-light ratios, and galaxy assembly bias.

During a galaxy merger, the supermassive black hole (SMBH) in each galaxy is thought to sink to the center of the potential and form a supermassive black hole binary; this binary can eject stars via 3-body scattering, bringing the SMBHs ever closer. In a static spherical galaxy model, the binary stalls at a separation of about a parsec after ejecting all the stars in its loss cone -- this is the well-known final parsec problem. However it has been shown that SMBH binaries in non-spherical galactic nuclei harden at a nearly constant rate until reaching the gravitational wave regime. Here we use a suite of direct N-body simulations to follow SMBH binary evolution in both corotating and counterrotating flattened galaxy models. For N larger than 500K, we find that the evolution of the SMBH binary is convergent, and is independent of the particle number. Rotation in general increases the hardening rate of SMBH binaries even more effectively than galaxy geometry alone. SMBH binary hardening rates are similar for co- and counterrotating galaxies. In the corotating case, the center of mass of SMBH binary settles into an orbit that is in a corotation resonance with the background rotating model, and the coalescence time is roughly few hundred Myr faster than a non-rotating flattened model. We find that counterrotation drives SMBHs to coalesce on a nearly radial orbit promptly after forming a hard binary. We discuss the implications for gravitational wave astronomy, hypervelocity star production, and the effect on the structure of the host galaxy.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

We present and test TesseRACt, a non-parametric technique for recovering the concentration of simulated dark matter halos using Voronoi tessellation. TesseRACt is tested on idealized N-body halos that are axisymmetric, triaxial, and contain substructure and compared to traditional least-squares fitting as well as two non-parametric techniques that assume spherical symmetry. TesseRACt recovers halo concentrations within 0.3% of the true value regardless of whether the halo is spherical, axisymmetric, or triaxial. Traditional fitting and non-parametric techniques that assume spherical symmetry can return concentrations that are systematically off by as much as 10% from the true value for non-spherical halos. TesseRACt also performs significantly better when there is substructure present outside $0.5R_{200}$. Given that cosmological halos are rarely spherical and often contain substructure, we discuss implications for studies of halo concentration in cosmological N-body simulations including how choice of technique for measuring concentration might bias scaling relations.

The evolving Laser Interferometer Space Antenna (eLISA) will revolutionize our understanding of the formation and evolution of massive black holes along cosmic history by probing massive black hole binaries in the $10^3-10^7$ solar mass range out to redshift $z\gtrsim 10$. High signal-to-noise ratio detections of $\sim 10-100$ binary coalescences per year will allow accurate measurements of the parameters of individual binaries (such as their masses, spins and luminosity distance), and a deep understanding of the underlying cosmic massive black hole parent population. This wealth of unprecedented information can lead to breakthroughs in many areas of physics, including astrophysics, cosmology and fundamental physics. We review the current status of the field, recent progress and future challenges.

We study the structure and evolution of dark matter halos from z = 300 to z = 6 for two cosmological N-body simulation initialization techniques. While the second order Lagrangian perturbation theory (2LPT) and the Zel'dovich approximation (ZA) both produce accurate present day halo mass functions, earlier collapse of dense regions in 2LPT can result in larger mass halos at high redshift. We explore the differences in dark matter halo mass and concentration due to initialization method through three 2LPT and three ZA initialized cosmological simulations. We find that 2LPT induces more rapid halo growth, resulting in more massive halos compared to ZA. This effect is most pronounced for high mass halos and at high redshift. Halo concentration is, on average, largely similar between 2LPT and ZA, but retains differences when viewed as a function of halo mass. For both mass and concentration, the difference between typical individual halos can be very large, highlighting the shortcomings of ZA-initialized simulations for high-z halo population studies.

It is known that two supermassive black holes (SMBHs) cannot merge in a spherical galaxy within a Hubble time; an emerging picture is that galaxy geometry, rotation, and large potential perturbations may usher the SMBH binary through the critical three-body scattering phase and ultimately drive the SMBH to coalesce. We explore the orbital content within an N-body model of a mildly- flattened, non-rotating, SMBH-embedded elliptical galaxy. When used as the foundation for a study on the SMBH binary coalescence, the black holes bypassed the binary stalling often seen within spherical galaxies and merged on Gyr timescales (Khan et al. 2013). Using both frequency-mapping and angular momentum criteria, we identify a wealth of resonant orbits in the axisymmetric model, including saucers, that are absent from an otherwise identical spherical system and that can potentially interact with the binary. We quantified the set of orbits that could be scattered by the SMBH binary, and found that the axisymmetric model contained nearly seven times the number of these potential loss cone orbits compared to our equivalent spherical model. In this flattened model, the mass of these orbits is roughly 3 times of that of the SMBH, which is consistent with what the SMBH binary needs to scatter to transition into the gravitational wave regime.

The relation of central black hole mass and stellar spheroid velocity dispersion (the M-$\sigma$ relation) is one of the best-known and tightest correlations linking black holes and their host galaxies. There has been much scrutiny concerning the difficulty of obtaining accurate black hole measurements, and rightly so; however, it has been taken for granted that measurements of velocity dispersion are essentially straightforward. We examine five disk galaxies from cosmological SPH simulations and find that line-of-sight effects due to galaxy orientation can affect the measured $\sigma$ by 30%, and consequently black hole mass predictions by up to 1.0 dex. Face-on orientations correspond to systematically lower velocity dispersion measurements, while more edge-on orientations give higher velocity dispersions, due to contamination by disk stars when measuring line of sight quantities. We caution observers that the uncertainty of velocity dispersion measurements is at least 20 km/s, and can be much larger for moderate inclinations. This effect may account for some of the scatter in the locally measured M-$\sigma$ relation, particularly at the low-mass end. We provide a method for correcting observed $\sigma_{\rm los}$ values for inclination effects based on observable quantities.

We present GPU implementations of two fast force calculation methods, based on series expansions of the Poisson equation. One is the Self-Consistent Field (SCF) method, which is a Fourier-like expansion of the density field in some basis set; the other is the Multipole Expansion (MEX) method, which is a Taylor-like expansion of the Green's function. MEX, which has been advocated in the past, has not gained as much popularity as SCF. Both are particle-field method and optimized for collisionless galactic dynamics, but while SCF is a "pure" expansion, MEX is an expansion in just the angular part; it is thus capable of capturing radial structure easily, where SCF needs a large number of radial terms. We show that despite the expansion bias, these methods are more accurate than direct techniques for the same number of particles. The performance of our GPU code, which we call ETICS, is profiled and compared to a CPU implementation. On the tested GPU hardware, a full force calculation for one million particles took ~ 0.1 seconds (depending on expansion cutoff), making simulations with as many as $10^8$ particles fast on a comparatively small number of nodes.

Although supermassive black holes (SMBHs) correlate well with their host galaxies, there is an emerging view that outliers exist. Henize 2-10, NGC 4889, and NGC1277 are examples of SMBHs at least an order of magnitude more massive than their host galaxy suggests. The dynamical effects of such ultramassive central black holes is unclear. Here, we perform direct N-body simulations of mergers of galactic nuclei where one black hole is ultramassive to study the evolution of the remnant and the black hole dynamics in this extreme regime. We find that the merger remnant is axisymmetric near the center, while near the large SMBH influence radius, the galaxy is triaxial. The SMBH separation shrinks rapidly due to dynamical friction, and quickly forms a binary black hole; if we scale our model to the most massive estimate for the NGC1277 black hole, for example, the timescale for the SMBH separation to shrink from nearly a kiloparsec to less than a parsec is roughly 10 Myr. By the time the SMBHs form a hard binary, gravitational wave emission dominates, and the black holes coalesce in a mere few Myr. Curiously, these extremely massive binaries appear to nearly bypass the 3-body scattering evolutionary phase. Our study suggests that in this extreme case, SMBH coalescence is governed by dynamical friction followed nearly directly by gravitational wave emission, resulting in an rapid and efficient SMBH coalescence timescale. We discuss the implications for gravitational wave event rates and hypervelocity star production.

Recently, both simulations and observations have revealed that flybys - fast, one-time interactions between two galaxy halos - are surprisingly common, nearing/comparable to galaxy mergers. Since these are rapid, transient events with the closest approach well outside the galaxy disk, it is unclear if flybys can transform the galaxy in a lasting way. We conduct collisionless N-body simulations of three co-planer flyby interactions between pure-disk galaxies to take a first look at the effects flybys have on disk structure, with particular focus on stellar bar formation. We find that some flybys are capable of inciting a bar with bars forming in both galaxies during our 1:1 interaction and in the secondary during our 10:1 interaction. The bars formed have ellipticities >0.5, sizes on the order of the host disk's scale length, and persist to the end of our simulations, ~5 Gyr after pericenter. The ability of flybys to incite bar formation implies that many processes associated with secular bar evolution may be more closely tied with interactions than previously though.

We present 20 candidate hypervelocity stars from the Sloan Extension for Galactic Understanding and Exploration (SEGUE) G and K dwarf samples. Previous searches for hypervelocity stars have only focused on large radial velocities; in this study we also use proper motions to select the candidates. We determine the hypervelocity likelihood of each candidate by means of Monte Carlo simulations, considering the significant errors often associated with high proper motion stars. We find that nearly half of the candidates exceed their escape velocities with at least 98% probability. Every candidate also has less than a 25% chance of being a high-velocity fluke within the SEGUE sample. Based on orbits calculated using the observed six-dimensional positions and velocities, few, if any, of these candidates originate from the Galactic center. If these candidates are truly hypervelocity stars, they were not ejected by interactions with the Milky Way's supermassive black hole. This calls for a more serious examination of alternative hypervelocity-star ejection scenarios.

Hierarchical structure formation theory is based on the notion that mergers drive galaxy evolution, so a considerable framework of semi-analytic models and N-body simulations has been constructed to calculate how mergers transform a growing galaxy. However, galaxy mergers are only one type of major dynamical interaction between halos -- another class of encounter, a close flyby, has been largely ignored. We analyze a 50 Mpc/h, $1024^3$ collisionless cosmological simulation and find that the number of close flyby interactions is comparable to, or even surpasses, the number of mergers for halo masses $\ga 10^{11}\,{h^{-1} M_\odot}$ at $z \la 2$. Halo flybys occur so frequently to high mass halos that they are continually perturbed, unable to reach a dynamical equilibrium. We also find tentative evidence that at high redshift, $z \ga 14$, flybys are as frequent as mergers. Our results suggest that close halo flybys can play an important role in the evolution of the earliest dark matter halos and their galaxies, and can still influence galaxy evolution at the present epoch.

During a galaxy merger, the supermassive black hole (SMBH) in each galaxy is thought to sink to the center of the potential and form a supermassive black hole binary; this binary can eject stars via 3-body scattering, bringing the SMBHs ever closer. In a static spherical galaxy model, the binary stalls at a separation of about a parsec after ejecting all the stars in its loss cone -- this is the well-known final parsec problem. Earlier work has shown that the centrophilic orbits in triaxial galaxy models are key in refilling the loss cone at a high enough rate to prevent the black holes from stalling. However, the evolution of binary SMBHs has never been explored in axisymmetric galaxies, so it is not clear if the final parsec problem persists in these systems. Here we use a suite of direct N-body simulations to follow SMBH binary evolution in galaxy models with a range of ellipticity. For the first time, we show that mere axisymmetry can solve the final parsec problem; we find the the SMBH evolution is independent of N for an axis ratio of c/a=0.8, and that the SMBH binary separation reaches the gravitational radiation regime for c/a=0.75.

Observations of the Galactic Center (GC) have accumulated a multitude of "forensic" evidence indicating that several million years ago the center of the Milky Way galaxy was teaming with starforming and accretion-powered activity -- this paints a rather different picture from the GC as we understand it today. We examine a possibility that this epoch of activity could have been triggered by the infall of a satellite galaxy into the Milky Way which began at the redshift of 10 and ended few million years ago with a merger of the Galactic supermassive black hole with an intermediate mass black hole brought in by the inspiralling satellite.

We present a wide (8.5x6.7 degree, 1050x825 kpc), deep (sigma(N_HI)=10^(16.8-17.5) cm^-2) neutral hydrogen (HI) map of the M101 galaxy group. We identify two new HI sources in the group environment, one an extremely low surface brightness (and hitherto unknown) dwarf galaxy, and the other a starless HI cloud, possibly primordial in origin. Our data show that M101's extended HI envelope (Huchtmeier & Witzel 1979) takes the form of a ~100 kpc long tidal loop or plume of HI extending to the southwest of the galaxy. The plume has an HI mass ~ 10^8 Msun and a peak column density of N_HI=5x10^17 cm^-2, and while it rotates with the main body of M101, it shows kinematic peculiarities suggestive of a warp or flaring out of the rotation plane of the galaxy. We also find two new HI clouds near the plume with masses ~ 10^7 Msun, similar to HI clouds seen in the M81/M82 group, and likely also tidal in nature. Comparing to deep optical imaging of the M101 group, neither the plume nor the clouds have any extended optical counterparts down to a limiting surface brightness of mu_B = 29.5. We also trace HI at intermediate velocities between M101 and NGC 5474, strengthening the case for a recent interaction between the two galaxies. The kinematically complex HI structure in the M101 group, coupled with the optical morphology of M101 and its companions, suggests that the group is in a dynamically active state that is likely common for galaxies in group environments.

We explore structure formation in the dark ages ($z\sim 30-6$) using two well-known methods for initializing cosmological $N$-body simulations. Overall, both the Zel'dovich approximation (\za) and second order Lagrangian perturbation theory (\lpt) are known to produce accurate present-day dark matter halo mass functions. However, since the \lpt method drives more rapid evolution of dense regions, it increases the occurrence of rare massive objects -- an effect that is most pronounced at high redshift. We find that \lpt produces more halos that could harbor Population III stars and their black hole remnants, and they produce them earlier. Although the differences between the \lpt and \za mass functions are nearly erased by $z=6$, this small boost to the number and mass of black holes more than doubles the reionized volume of the early Universe. We discuss the implications for reionization and massive black hole growth.

We present a comparative study of recent works on merger-timescales with dynamical friction and find a strong contrast between idealized/isolated mergers (Boylan-Kolchin et al. 2008) and mergers from a cosmological volume (Jiang et al. 2008). Our study measures the duration of mergers in a cosmological N-body simulation of dark matter, with emphasis on higher redshifts (z < 10) and a lower mass range. In our analysis we consider and compare two merger definitions; tidal disruption and coalescence. We find that the merger-time formula proposed by Jiang et al. (2008) describes our results well and conclude that cosmologically motivated merger-time formulae provide a more versatile and statistically robust approximation for practical applications such as semi-analytic/hybrid models.

Hierarchical structure formation theory is based on the notion that mergers drive galaxy evolution, so a considerable framework of semi-analytic models and N-body simulations has been constructed to calculate how mergers transform a growing galaxy. However, galaxy mergers are only one type of major dynamical interaction between halos -- another class of encounter, a close flyby, has been largely ignored. We use cosmological N-body simulations to reconstruct the entire dynamical interaction history of dark matter halos. We present a careful method of identifying and tracking a dark matter halo which resolves the typical classes of anomalies that occur in N-body data. This technique allows us to robustly follow halos and several hierarchical levels of subhalos as they grow, dissolve, merge, and flyby one another -- thereby constructing both a census of the dynamical interactions in a volume and an archive of the dynamical evolution of an individual halo. In addition to a census of mergers, our tool characterizes the frequency of close flyby interactions in the Universe. We find that the number of close flyby interactions is comparable to, or even surpasses, the number of mergers for halo masses $\gtrsim 10^{11}\,\Msun$ at $z \lesssim 2$. Halo flybys occur so frequently to high mass halos that they are continually perturbed, unable to reach a dynamical equilibrium. In particular, we find that Milky Way type halos undergo a similar number of flybys as mergers irrespective of mass-ratio for $z\lesssim 2$. We also find tentative evidence that at high redshift, $z \gtrsim 14$, flybys are as frequent as mergers. Our results suggest that close halo flybys can play an important role in the evolution of the earliest dark matter halos and their galaxies, and can still influence galaxy evolution at the present epoch. [Abridged]

Tidal stripping and three-body interactions with the central supermassive black hole may eject stars from the Milky Way. These stars would comprise a set of `intragroup' stars that trace the past history of interactions in our galactic neighborhood. Using the Sloan Digital Sky Survey DR7, we identify candidate solar metallicity red giant intragroup stars using color cuts that are designed to exclude nearby M and L dwarfs. We present 677 intragroup candidates that are selected between 300 kpc and 2 Mpc, and are either the reddest intragroup candidates (M7-M10) or are L dwarfs at larger distances than previously detected.

We combine high resolution N-body simulations with deep observations of neutral hydrogen (HI) in nearby galaxy groups in order to explore two well-known theories of HI cloud formation: HI stripping by galaxy interactions and dark matter minihalos with embedded HI gas. This paper presents new data from three galaxy groups, Canes Venatici I, NGC 672, and NGC 45, and assembles data from our previous galaxy group campaign to generate a rich HI cloud archive to compare to our simulated data. We find no HI clouds in the Canes Venatici I, NGC 672, or NGC 45 galaxy groups. We conclude that HI clouds in our detection space are most likely to be generated through recent, strong galaxy interactions. We find no evidence of HI clouds associated with dark matter halos above M_HI = 10^6 M_Sun, within +/- 700 km/s of galaxies, and within 50 kpc projected distance of galaxies.

We explore the growth of < 10^7 Msun black holes that reside at the centers of spiral and field dwarf galaxies in a Local Group type of environment. We use merger trees from a cosmological N-body simulation known as Via Lactea II (VL-2) as a framework to test two merger-driven semi-analytic recipes for black hole growth that include dynamical friction, tidal stripping, and gravitational wave recoil in over 20,000 merger tree realizations. First, we apply a Fundamental Plane limited (FPL) model to the growth of Sgr A*, which drives the central black hole to a maximum mass limited by the Black Hole Fundamental Plane after every merger. Next, we present a new model that allows for low-level Prolonged Gas Accretion (PGA) during the merger. We find that both models can generate a Sgr A* mass black hole. We predict a population of massive black holes in local field dwarf galaxies - if the VL-2 simulation is representative of the growth of the Local Group, we predict up to 35 massive black holes (< 10^6 Msun) in Local Group field dwarfs. We also predict that hundreds of < 10^5 Msun black holes fail to merge, and instead populate the Milky Way halo, with the most massive of them at roughly the virial radius. In addition, we find that there may be hundreds of massive black holes ejected from their hosts into the nearby intergalactic medium due to gravitational wave recoil. We discuss how the black hole population in the Local Group field dwarfs may help to constrain the growth mechanism for Sgr A*.

We comment on the recent observation of a 115-day modulation in the X-ray flux of the ultraluminous X-ray source (ULX) NGC 5408 X-1, and in particular, the interpretation of this modulation as the orbital period. We suggest that this modulation may instead be due to a precessing jet, and is thus superorbital in nature. Comparing the properties of this ULX with those of the prototypical microquasar SS 433, we argue that NGC 5408 X-1 is very similar to SS 433: a hyper-accreting stellar mass black hole in a shorter-period binary. If the analogy holds, the 115-day modulation is best explained by the still poorly-understood physics of inner-disc/jet precession and a longer observing baseline would be able to reveal an intrinsic phase jitter that is associated with such a precession.

Cosmological galaxy formation models predict the existence of dark matter minihalos surrounding galaxies and in filaments connecting groups of galaxies. The more massive of these minihalos are predicted to host HI gas that should be detectable by current radio telescopes such as the GBT. We observed the region including the M81/M82 and NGC 2403 galaxy groups, searching for observational evidence of an HI component associated with dark matter halos within the "M81 Filament", using the Robert C. Byrd Green Bank Telescope (GBT). The map covers an 8.7 degree x 21.3 degree (480 kpc x 1.2 Mpc) region centered between the M81/M82 and NGC 2403 galaxy groups. Our observations cover a wide velocity range, from -890 to 1320 km/s, which spans much of the range predicted by cosmological N-body simulations for dark matter minihalo velocities. Our search is not complete in the velocity range -210 to 85 km/s, containing Galactic emission and the HVC Complex A. For an HI cloud at the distance of M81, with a size < 10 kpc, our average 5-sigma mass detection limit is 3.2 x 10^6 M_Sun, for a linewidth of 20 km/s. We compare our observations to two large cosmological N-body simulations and find that the simulation predicts a significantly greater number of detectable minihalos than are found in our observations, and that the simulated minihalos do not match the phase space of observed HI clouds. These results place strong constraints on the HI gas that can be associated with dark-matter halos. Our observations indicate that the majority of extragalactic HI clouds with a mass greater than 10^6 M_Sun are likely to be generated through tidal stripping caused by galaxy interactions.

Observations indicate that roughly 60% of the baryons may exist in a Warm-Hot Intergalactic Medium (WHIM) at low redshifts. Following up on previous results showing that gas is released through galaxy mergers, we use a semi-analytic technique to estimate the fraction of gas mass lost from haloes solely due to mergers. We find that up to ~25% of the gas in a halo can unbind over the course of galaxy assembly. This process does not act preferentially on smaller mass haloes; bigger haloes \emph{always} release larger amounts of gas in a given volume of the Universe. However, if we include multi-phase gas accretion onto haloes, we find that only a few percent is unbound. We conclude that either non-gravitational processes may be in play to heat up the gas in the galaxies prior to unbinding by mergers or most of the baryons in the WHIM have never fallen into virialised dark matter haloes. We present a budget for stocking the WHIM compiled from recent work.

We calculate the gravitational wave signal from the growth of 10 million solar mass supermassive black holes (SMBH) from the remnants of Population III stars. The assembly of these lower mass black holes is particularly important because observing SMBHs in this mass range is one of the primary science goals for the Laser Interferometer Space Antenna (LISA), a planned NASA/ESA mission to detect gravitational waves. We use high resolution cosmological N-body simulations to track the merger history of the host dark matter halos, and model the growth of the SMBHs with a semi-analytic approach that combines dynamical friction, gas accretion, and feedback. We find that the most common source in the LISA band from our volume consists of mergers between intermediate mass black holes and SMBHs at redshifts less than 2. This type of high mass ratio merger has not been widely considered in the gravitational wave community; detection and characterization of this signal will likely require a different technique than is used for SMBH mergers or extreme mass ratio inspirals. We find that the event rate of this new LISA source depends on prescriptions for gas accretion onto the black hole as well as an accurate model of the dynamics on a galaxy scale; our best estimate yields about 40 sources with a signal-to-noise ratio greater than 30 occur within a volume like the Local Group during SMBH assembly -- extrapolated over the volume of the universe yields roughly 500 observed events over 10 years, although the accuracy of this rate is affected by cosmic variance.

We have observed the NGC 2403 group of galaxies using the Robert C. Byrd Green Bank Telescope (GBT) in a search for faint, extended neutral hydrogen clouds similar to the clouds found around the M81/M82 group, which is located approximately 250 kpc from the NGC 2403 group along the same filament of galaxies. For an HI cloud with a size < 10 kpc within 50 kpc of a group galaxy, our 7-sigma mass detection limit is 2.2 x 10^6 M_sun for a cloud with a linewidth of 20 km/s, over the velocity range from -890 to 1750 km/s. At this sensitivity level we detect 3 new HI clouds in the direction of the group, as well as the known galaxies. The mean velocity of the new clouds differs from that of the group galaxies by more than 250 km/s, but are in the range of Milky Way High Velocity Clouds (HVCs) in that direction. It is most likely that the clouds are part of the Milky Way HVC population. If HI clouds exist in the NGC 2403 group, their masses are less than 2.2 x 10^6 M_sun. We also compared our results to structures that are expected based on recent cosmological models, and found none of the predicted clouds. If NGC 2403 is surrounded by a population of dark matter halos similar to those proposed for the Milky Way in recent models, our observations imply that their HI content is less than 1% of their total mass.

Galaxy merger simulations have explored the behaviour of gas within the galactic disk, yet the dynamics of hot gas within the galaxy halo has been neglected. We report on the results of high-resolution hydrodynamic simulations of colliding galaxies with hot halo gas. We explore a range of mass ratios, gas fractions and orbital configurations to constrain the shocks and gas dynamics within the progenitor haloes. We find that : (i) A strong shock is produced in the galaxy haloes before the first passage, increasing the temperature of the gas by almost an order of magnitude to $T\sim 10^{6.3}$ K. (ii) The X-ray luminosity of the shock is strongly dependent on the gas fraction; it is $\gtrsim 10^{39}$ erg/s for halo gas fractions larger than 10%. (iii) The hot diffuse gas in the simulation produces X-ray luminosities as large as $10^{42}$ erg/s. This contributes to the total X-ray background in the Universe. (iv) We find an analytic fit to the maximum X-ray luminosity of the shock as a function of merger parameters. This fit can be used in semi-analytic recipes of galaxy formation to estimate the total X-ray emission from shocks in merging galaxies. (v) $\sim$ 10-20% of the initial gas mass is unbound from the galaxies for equal-mass mergers, while $3-5%$ of the gas mass is released for the 3:1 and 10:1 mergers. This unbound gas ends up far from the galaxy and can be a feasible mechanism to enrich the IGM with metals.

Galactic double white dwarfs were postulated as a source of confusion limited noise for LISA, the future space-based gravitational wave observatory. Until very recently, the Galactic population consisted of a relatively well studied disk population, a somewhat studied smaller bulge population and a mostly unknown, but potentially large halo population. It has been argued that the halo population may produce a signal that is much stronger (factor of ~5 in spectral amplitude) than the disk population. However, this surprising result was not based on an actual calculation of a halo white dwarf population but was derived on (i) the assumption that one can extrapolate the halo population properties from those of the disk population and (ii) the postulated (unrealistically) high number of white dwarfs in the halo. We perform the first calculation of a halo white dwarf population using population synthesis models. Our comparison with the signal arising from double white dwarfs in the Galactic disk+bulge clearly shows that it is impossible for the double white dwarf halo signal to exceed that of the rest of the Galaxy. Using microlensing results to give an upper limit on the content of white dwarfs in the halo (~30 % baryonic mass in white dwarfs), our predicted halo signal is a factor of 10 lower than the disk+bulge signal. Even in the implausible case where all of the baryonic halo mass is found in white dwarfs, the halo signal does not become comparable to that of the disk+bulge, and thus would still have a negligible effect on the detection of other LISA sources.

(Abridged) We use high resolution cosmological N-body simulations to study the growth of intermediate to supermassive black holes from redshift 49 to zero. We track the growth of black holes from the seeds of population III stars to black holes in the range of 10^3 < M < 10^7 Msun -- not quasars, but rather IMBH to low-mass SMBHs. These lower mass black holes are the primary observable for the Laser Interferometer Space Antenna (LISA). The large-scale dynamics of the black holes are followed accurately within the simulation down to scales of 1 kpc; thereafter, we follow the merger analytically from the last dynamical friction phase to black hole coalescence. We find that the merger rate of these black holes is R~25 per year between 8 < z < 11 and R = 10 per year at z=3. Before the merger occurs the incoming IMBH may be observed with a next generation of X-ray telescopes as a ULX source with a rate of about ~ 3 - 7 per year for 1 < z < 5. We develop an analytic prescription that captures the most important black hole growth mechanisms: galaxy merger-driven gas accretion and black hole coalescence. Using this, we find that SMBH at the center of Milky Way type galaxy was in place with most of its mass by z = 4.7, and most of the growth was driven by gas accretion excited by major mergers. Hundreds of black holes have failed to coalesce with the SMBH by z=0, some with masses of 10000 Msun, orbiting within the dark matter halo with luminosities up to ~ 30000 Lsun. These X-ray sources can easily be observed with Chandra at ~ 100 kpc.

Extreme mass ratio bursts (EMRBs) have been proposed as a possible source for future space-borne gravitational wave detectors, such as the Laser Interferometer Space Antenna (LISA). These events are characterized by long-period, nearly-radial orbits of compact objects around a central massive black hole. The gravitational radiation emitted during such events consists of a short burst, corresponding to periapse passage, followed by a longer, silent interval. In this paper we investigate the impact of including relativistic corrections to the description of the compact object's trajectory via a geodesic treatment, as well as including higher-order multipole corrections in the waveform calculation. The degree to which the relativistic corrections are important depends on the EMRB's orbital parameters. We find that relativistic EMRBs (v_{max}}/c > 0.25) are not rare and actually account for approximately half of the events in our astrophysical model. The relativistic corrections tend to significantly change the waveform amplitude and phase relative to a Newtonian description, although some of this dephasing could be mimicked by parameter errors. The dephasing over several bursts could be of particular importance not only to gravitational wave detection, but also to parameter estimation, since it is highly correlated to the spin of the massive black hole. Consequently, we postulate that if a relativistic EMRB is detected, such dephasing might be used to probe the relativistic character of the massive black hole and obtain information about its spin.

During the inspiral and merger of a binary black hole, gravitational radiation is emitted anisotropically due to asymmetries in the merger configuration. This anisotropic radiation leads to a gravitational wave kick, or recoil velocity, as large as ~ 4000 km/sec. We investigate the effect gravitational recoil has on the retention of intermediate mass black holes (IMBH) within Galactic globular clusters. Assuming that our current understanding of IMBH-formation is correct and yields an IMBH-seed in every globular cluster, we find a significant problem retaining low mass IMBHs (1000 $\Msun$) in the typical merger-rich globular cluster environment. Given a uniform black hole spin distribution and orientation and a Kroupa IMF, we find that at most 3% of the globular clusters can retain an IMBH larger than 1000 $\Msun$ today. For a population of black holes that better approximates mass loss from winds and supernovae, we find that 16% of globulars can retain an IMBH larger than 1000 $\Msun$. Our calculations show that if there are black holes of mass $M > 60 \Msun$ in a cluster, repeated IMBH-BH encounters will eventually eject a 1000 $\Msun$ IMBH with greater than 30% probability. As a consequence, a large population of rogue black holes may exist in our Milky Way halo. We discuss the dynamical implications of this subpopulation, and its possible connection to ultraluminous X-ray sources (ULXs).

Understanding how seed black holes grow into intermediate and supermassive black holes (IMBHs and SMBHs, respectively) has important implications for the duty-cycle of active galactic nuclei (AGN), galaxy evolution, and gravitational wave astronomy. Most studies of the cosmological growth and merger history of black holes have used semianalytic models and have concentrated on SMBH growth in luminous galaxies. Using high resolution cosmological N-body simulations, we track the assembly of black holes over a large range of final masses -- from seed black holes to SMBHs -- over widely varying dynamical histories. We used the dynamics of dark matter halos to track the evolution of seed black holes in three different gas accretion scenarios. We have found that growth of Sagittarius A* - size SMBH reaches its maximum mass M_{SMBH}~10^6Msun at z~6 through early gaseous accretion episodes, after which it stays at near constant mass. At the same redshift, the duty-cycle of the host AGN ends, hence redshift z=6 marks the transition from an AGN to a starburst galaxy which eventually becomes the Milky Way. By tracking black hole growth as a function of time and mass, we estimate that the IMBH merger rate reaches a maximum of R_{max}=55 yr^-1 at z=11. From IMBH merger rates we calculate N_{ULX}=7 per Milky Way type galaxy per redshift in redshift range 2<z<6.

Stellar mass compact objects in short period orbits about a $10^{4.5}$--$10^{7.5}$ solar mass massive black hole (MBH) are thought to be a significant continuous-wave source of gravitational radiation for the ESA/NASA Laser Interferometer Space Antenna (LISA) gravitational wave detector. However, these extreme mass-ratio inspiral sources began in long-period, nearly parabolic orbits that have multiple close encounters with the MBH. The gravitational radiation emitted during these close encounters may be detectable by LISA as a gravitational wave burst if the characteristic passage timescale is less than $10^5$ seconds. Scaling a static, spherical model to the size and mass of the Milky Way bulge we estimate an event rate of ~ 15 per year for such burst signals, detectable by LISA with signal-to-noise greater than five, originating in our galaxy. When extended to include Virgo cluster galaxies our estimate increases to a gravitational wave burst rate of ~ 18. We conclude that these extreme mass-ratio burst sources may be a steady and significant source of gravitational radiation in the LISA data streams.

Based on a high resolution cosmological n-body simulation, we track the hierarchical growth of black holes in galaxy clusters from z=20 to z=0. We present a census of black holes as function of redshift and will determine their mass assembly history under a variety of assumptions regarding the importance of gas accretion in black hole growth, from early supercritical Eddington accretion to gas-poor hierarchical assembly. Following a galaxy merger, black holes are expected to form, inspiral and merge after strongly radiating energy via gravitational waves. For each binary black hole inspiral and merger, we determine the expected gravitational wave signal for the Laser Interferometer Space Antenna (LISA), and calculate the LISA event rate as a function of time. We will calculate the black hole mass assembly history for several black hole growth scenerios, so that we can explore tests to characterize each model observationally. In particular, we will study how well LISA observations will be able to distinguish between these very different assembly scenarios.

Gravitational radiation from the Galactic population of close white dwarf binaries (CWDBs) is expected to produce a confusion-limited signal at the lower end of the sensitivity band of the Laser Interferometer Space Antenna (LISA). The canonical scale height of the disk population has been taken to be 90 pc for most studies of the nature of this confusion-limited signal. This estimate is probably too low, and the consequences of a more realistic scale height are investigated with a model of the LISA signal due to populations of close white dwarf binaries with different scale heights. If the local space density of CWDBs is held constant, increasing the scale height results in both an increase in the overall strength of the confusion-limited signal as well as in increase in the frequency at which the signals become individually resolvable. If the total number of binaries is held constant, increasing the scale height results in a reduction of the number of expected bright signals above the confusion-limited signal at low frequencies. We introduce an estimator for comparing this transition frequency that takes into account the signal spreading at higher frequencies.

Stars and compact objects that plunge toward a black hole are either 1) captured, emitting gravitational waves as the orbit decays, 2) tidally disrupted, leaving a disc of baryonic material, 3) scattered to a large radius, where they may thereafter avoid encounters with the black hole or 4) swallowed whole, contributing to black hole growth. These processes occur on a dynamical time, which implies that for a static spherically symmetric stellar system, the loss cone is quickly emptied. However, most elliptical galaxies and spiral bulges are thought to be triaxial in shape. The centrophilic orbits comprising the backbone of a triaxial galaxy have been suggested as one way to keep the loss cone around a supermassive black hole filled with stars, stellar remnants, and intermediate mass black holes. We investigate the evolution of the loss cone population in a triaxial galaxy model with high resolution N-body simulations. We find that enough regular orbits flow through angular momentum space to maintain a full loss cone for a Hubble time. This increases the astrophysical capture rate by several orders of magnitude. In the Milky Way, for example, we find that the white dwarf capture rate can be as high as $10^{-5}$ per year, 100 times larger than previous estimates based on spherical models for the bulge.

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.

The evolution of a stellar bar transforms not only the galactic disk, but also the host dark matter halo. We present high resolution, fully self-consistent N-body simulations that clearly demonstrate that dark matter halo central density cusps flatten as the bar torques the halo. This effect is independent of the bar formation mode and occurs even for rather short bars. The halo and bar evolution is mediated by resonant interactions between orbits in the halo and the bar pattern speed, as predicted by linear Hamiltonian perturbation theory. The bar lengthens and slows as it loses angular momentum, a process that occurs even in rather warm disks. We demonstrate that the bar and halo response can be critically underestimated for experiments that are unable to resolve the relevant resonant dynamics; this occurs when the phase space in the resonant region is under sampled or plagued by noise.

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.

Supermassive black holes can be seen as an agent of galaxy transformation. In particular, a supermassive black hole can cause a triaxial galaxy to evolve toward axisymmetry by inducing chaos in centrophilic orbit families. This is one way in which a single supermassive black hole can induce large-scale changes in the structure of its host galaxy -- changes on scales far larger than the Schwarzschild radius ($O(10^{-5}) \rm{pc}$) and the radius of influence of the black hole ($O(1)-O(100) \rm{pc}$). We will discuss the transformative power of supermassive black holes in light of recent high resolution N-body realizations of cuspy triaxial galaxies.

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.

In order for the core Fundamental Plane (cFP) to endure, small ellipticals must not survive mergers with giant ellipticals, despite the fact the small secondary galaxy can be as much as a million times more dense than the primary. However, our previous set of experiments has shown that, for purely stellar galaxies, the secondary does in fact survive mergers with primaries up to 100 times more massive. In this paper, we demonstrate the effect that a massive central black hole has on mergers of cFP galaxies. Our results indicate that the addition of a massive central singularity inside the primary galaxy provides strong enough tidal forces to destroy dense cFP companions when the secondary's orbital decay is sufficiently elongated. The destruction of the secondary acts to preserve the original low central density profile of the primary in the giant merger remnant, which allows the remnant to remain on the cFP. On more circular orbits, though, the secondary is only disrupted near the end of the merger, and the degree to which the secondary particles disperse depends on the amount of orbital energy left in the merger. Hence, there are some mergers for which the cFP is not preserved in our experiments. We find that if the secondary is not dispersed, it forms a spinning stellar disk with a central density that forces the merger remnant off the cFP.

We have developed an effective simulation method and applied it to the problem of the accretion of very dense secondary companions by tenuous primaries. We have studied the accretion of objects of varying luminosity ratios, with sizes and densities drawn from the core Fundamental Plane for elliptical galaxies under the assumption that the mass distribution in the central parts of the galaxies follows the light. The results indicate that in mergers with mass ratios greater than 10, the smaller object is only stripped down to the highest density encountered in the primary during the accretion process. Thus, the form of the core Fundamental Plane suggests that the mass distribution in these galaxy centers is different than the light distribution, or that mergers on the core Fundamental Plane require more than the dynamics of visible stars.