We present the discovery of the second binary with a Roche lobe-filling hot subdwarf transferring mass to a white dwarf (WD) companion. This 56 minute binary was discovered using data from the Zwicky Transient Facility. Spectroscopic observations reveal an He-sdOB star with an effective temperature of $T_{\rm eff}=33,700\pm1000$ K and a surface gravity of $log(g)=5.54\pm0.11$. The GTC+HiPERCAM light curve is dominated by the ellipsoidal deformation of the He-sdOB star and shows an eclipse of the He-sdOB by an accretion disk as well as a weak eclipse of the WD. We infer a He-sdOB mass of $M_{\rm sdOB}=0.41\pm0.04$ M$_\odot$ and a WD mass of $M_{\rm WD}=0.68\pm0.05$ M$_\odot$. The weak eclipses imply a WD black-body temperature of $63,000\pm10,000$ K and a radius $R_{\rm WD}=0.0148\pm0.0020$ M$_\odot$ as expected for a WD of such high temperature. The He-sdOB star is likely undergoing hydrogen shell burning and will continue transferring mass for $\approx1$ Myrs at a rate of $10^{-9} M_\odot {\rm yr}^{-1}$ which is consistent with the high WD temperature. The hot subdwarf will then turn into a WD and the system will merge in $\approx30$ Myrs. We suggest that Galactic reddening could bias discoveries towards preferentially finding Roche lobe-filling systems during the short-lived shell burning phase. Studies using reddening corrected samples should reveal a large population of helium core-burning hot subdwarfs with $T_{\rm eff}\approx25,000$ K in binaries of 60-90 minutes with WDs. Though not yet in contact, these binaries would eventually come into contact through gravitational wave emission and explode as a sub-luminous thermonuclear supernova or evolve into a massive single WD.

Early observations of Type Ia supernovae (SNe$\,$Ia) provide essential clues for understanding the progenitor system that gave rise to the terminal thermonuclear explosion. We present exquisite observations of SN$\,$2019yvq, the second observed SN$\,$Ia, after iPTF$\,$14atg, to display an early flash of emission in the ultraviolet (UV) and optical. Our analysis finds that SN$\,$2019yvq was unusual, even when ignoring the initial flash, in that it was moderately underluminous for an SN$\,$Ia ($M_g \approx -18.5\,$mag at peak) yet featured very high absorption velocities ($v \approx 15,000\,\mathrm{km\,s}^{-1}$ for Si II $\lambda$6355 at peak). We find that many of the observational features of SN$\,$2019yvq, aside from the flash, can be explained if the explosive yield of radioactive $^{56}\mathrm{Ni}$ is relatively low (we measure $M_{^{56}\mathrm{Ni}} = 0.31 \pm 0.05\,M_\odot$) and it and other iron-group elements are concentrated in the innermost layers of the ejecta. To explain both the UV/optical flash and peak properties of SN$\,$2019yvq we consider four different models: interaction between the SN ejecta and a nondegenerate companion, extended clumps of $^{56}\mathrm{Ni}$ in the outer ejecta, a double-detonation explosion, and the violent merger of two white dwarfs. Each of these models has shortcomings when compared to the observations; it is clear additional tuning is required to better match SN$\,$2019yvq. In closing, we predict that the nebular spectra of SN$\,$2019yvq will feature either H or He emission, if the ejecta collided with a companion, strong [Ca II] emission, if it was a double detonation, or narrow [O I] emission, if it was due to a violent merger.

The Zwicky Transient Facility has begun to discover binary systems with orbital periods that are less than 1\,hr. Combined with dedicated follow-up systems, which allow for high-cadence photometry of these sources, systematic confirmation and characterization of these sources are now possible. Here, we report the discovery of ZTF J190125.42+530929.5, a 40.6\,min orbital period, eclipsing double white-dwarf binary. Both photometric and spectroscopic modeling confirm its nature, yielding an estimated inclination of $i = 86.2^{+0.6}_{-0.2}\,\rm degrees$ and primary and secondary effective temperatures of $\textrm{T}_\textrm{eff} = 28,000^{+500}_{-500}\,K$ and $\textrm{T}_\textrm{eff} = 17,600^{+400}_{-400}\,K$ respectively. This system adds to a growing list of sources for future gravitational-wave detectors and contributes to the demographic analysis of double degenerates.

When fast radio burst (FRB) waves propagate through the local (< 1 pc) environment of the FRB source, electrons in the plasma undergo large-amplitude oscillations. The finite-amplitude effects cause the effective plasma frequency and cyclotron frequency to be dependent on the wave strength. The dispersion measure and rotation measure should therefore vary slightly from burst to burst for a repeating source, depending on the luminosity and frequency of the individual burst. Furthermore, free-free absorption of strong waves is suppressed due to the accelerated electrons' reduced energy exchange in Coulomb collisions. This allows bright low-frequency bursts to propagate through an environment that would be optically thick to low-amplitude waves. Given a large sample of bursts from a repeating source, it would be possible to use the deficit of low-frequency and low-luminosity bursts to infer the emission measure of the local intervening plasma and its distance from the source. Information about the local environment will shed light on the nature of FRB sources.

The past two decades have seen a tremendous investment in observational facilities that promise to reveal new and unprecedented discoveries about the universe. In comparison, the investment in theoretical work is completely dwarfed, even though theory plays a crucial role in the interpretation of these observations, predicting new types of phenomena, and informing observing strategies. In this white paper, we argue that in order to reach the promised critical breakthroughs in astrophysics over the next decade and well beyond, the national agencies must take a serious approach to investment in theoretical astrophysics research. We discuss the role of theory in shaping our understanding of the universe, and then we provide a multi-level strategy, from the grassroots to the national, to address the current underinvestment in theory relative to observational work.

We report the discovery of a detached double white dwarf binary with an orbital period of $\approx20.6\,\rm minutes$, PTF J053332.05+020911.6. The visible object in this binary, PTF J0533+0209B, is a $\approx 0.17\,M_\odot$ mass white dwarf with a helium-dominated atmosphere containing traces of hydrogen (DBA). This object exhibits ellipsoidal variations due to tidal deformation, and is the visible component in a single-lined spectroscopic binary with a velocity semi-amplitude of $K_B=618.7\pm6.9 \, \rm km \, s^{-1}$. We have detected significant orbital decay due to the emission of gravitational radiation, and we expect that the Laser Interferometer Space Antenna (LISA) will detect this system with a signal to noise of $8.4^{+4.2}_{-3.0}$ after four years of operation. Because this system already has a well determined orbital period, radial velocity semi-amplitude, temperature, atmospheric composition, surface gravity, and orbital decay rate, a LISA signal will help fully constrain the properties of this system by providing a direct measurement of its inclination. Thus, this binary demonstrates the synergy between electromagnetic and gravitational radiation for constraining the physical properties of an astrophysical object.

We present detailed submillimeter- through centimeter-wave observations of the extraordinary extragalactic transient AT2018cow. The apparent characteristics -- the high radio luminosity, the long-lived emission plateau at millimeter bands, and the sub-relativistic velocity -- have no precedent. A basic interpretation of the data suggests $E_k \gtrsim 10^{48}$ erg coupled to a fast but sub-relativistic ($v \approx 0.13c$) shock in a dense ($n_e \approx 3 \times 10^{5}\,$cm$^{-3}$) medium. We find that the X-ray emission is not naturally explained by an extension of the radio-submm synchrotron spectrum, nor by inverse Compton scattering of the dominant blackbody UVOIR photons by energetic electrons within the forward shock. By $\Delta t \approx20\,$days, the X-ray emission shows spectral softening and erratic inter-day variability. Taken together, we are led to invoke an additional source of X-ray emission: the central engine of the event. Regardless of the nature of this central engine, this source heralds a new class of energetic transients shocking a dense medium, which at early times are most readily observed at millimeter wavelengths.

We present the discovery of a nuclear transient with the Caltech-NRAO Stripe 82 Survey (CNSS), a dedicated radio transient survey carried out with the Karl G. Jansky Very Large Array (VLA). This transient, CNSS J001947.3+003527, exhibited a turn-on over a timescale of $\lesssim$1 yr, increasing in flux density at 3 GHz from $<0.14$ mJy in February 2014 to $4.4\pm0.1$ mJy in March 2015, reaching a peak luminosity of $5\times10^{28} \text{erg s}^{-1} \text{Hz}^{-1}$ around October 2015. The association of CNSS J0019+00 with the nucleus (Gaia and our VLBI positions are consistent to within 1 pc) of a nearby S0 Seyfert galaxy at 77 Mpc, together with the radio spectral evolution, implies that this transient is most likely a tidal disruption event (TDE). Our equipartition analysis indicates the presence of a $\sim$15,000 km s$^{-1}$ outflow, having energy $\sim$10$^{49}$ erg. We derive the radial density profile for the circum-nuclear material in the host galaxy to be proportional to R$^{-2.5}$. All of these properties suggest resemblance with radio-detected thermal TDEs like ASASSN-14li and XMMSL1 J0740-85. No significant X-ray or optical emission is detected from CNSS J0019+00, although this may simply be due to the thermal emission being weak during our late-time follow up observations. From the CNSS survey we have obtained the first unbiased measurement of the rate of radio TDEs, $R(>500 \mu{\rm Jy})$ of about $2\times10^{-3}$ deg$^{-2}$, or equivalently a volumetric rate of about 10 Gpc$^{-3}$ yr$^{-1}$. This rate implies that all-sky radio surveys such as the VLA Sky Survey and those planned with ASKAP, will find many tens of radio TDEs over the next few years.

We present the DSA-2000: a world-leading radio survey telescope and multi-messenger discovery engine for the next decade. The array will be the first true radio camera, outputting science-ready image data over the 0.7 - 2 GHz frequency range with a spatial resolution of 3.5 arcsec. With 2000 x 5 m dishes, the DSA-2000 will have an equivalent point-source sensitivity to SKA1-mid, but with ten times the survey speed. The DSA-2000 is envisaged as an all-sky survey instrument complementary to the ngVLA, and as a counterpart to the LSST (optical), SPHEREx (near-infrared) and SRG/eROSITA (X-ray) all-sky surveys. Over a five-year prime phase, the DSA-2000 will image the entire sky above declination -30 degrees every four months, detecting > 1 unique billion radio sources in a combined full-Stokes sky map with 500 nJy/beam rms noise. This all-sky survey will be complemented by intermediate and deep surveys, as well as spectral and polarization image cubes. The array will be a cornerstone for multi-messenger science, serving as the principal instrument for the US pulsar timing array community, and by searching for radio afterglows of compact object mergers detected by LIGO and Virgo. The array will simultaneously detect and localize ~10,000 fast radio bursts each year, realizing their ultimate use as a cosmological tool. The DSA-2000 will be proposed to the NSF Mid-Scale Research Infrastructure-2 program with a view to first light in 2026

A supermassive black hole SgrA* with the mass ~4x10^6 Msun resides at the centre of our galaxy. Building up such a massive black hole within the ~10^10 year lifetime of our galaxy would require a mean accretion rate of ~4x10^-4 Msun/yr. At present, X-ray observations constrain the rate of hot gas accretion at the Bondi radius (10^5 R_Sch = 0.04 pc at 8kpc) to \dot{M}_Bondi ~ 3x10^-6 Msun/yr, and polarization measurements constrain it near the event horizon to \dot{M}_horizon ~ 10^{-8} Msun/yr. A range of models was developed to describe the accretion gas onto an underfed black hole. However, the exact physics still remains to be understood. One challenge with the radiation inefficient accretion flows is that even if one understands the dynamics there is no accepted prescription for associating emissivity (and absorption) with the flow. The other issue is the lack of model-independent probes of accretion flow at intermediate radii (between few and ~ 10^5 R_Sch), i.e. the constraints that do not assume a model of accretion flow as an input parameter. We report detection and imaging of the 10^4 K ionized gas disk within 2x10^4 R_Sch in a mm hydrogen recombination line H30alpha: n = 31 -> 30 at 231.9 GHz using the ALMA. The emission was detected with a double-peaked line profile spanning full width of 2,200 km/s with the approaching and the receding components straddling Sgr A*, each offset from it by 0.11arcsec= 0.004pc. The red-shifted side is displaced to the north-east, while the blue-shifted side is displaced to the south-west. The limit on the total mass of ionized gas estimated from the emission is 10^-4 - 10^-5 Sun at a mean hydrogen density 10^5-10^6 cm^-3, depending upon whether or not we assume the presence of a uniform density disk or an ensemble of orbiting clouds, and the amplification factor of the mm radiation due to the strong background source which is Sgr A* continuum.

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.

We consider the effects of radio-wave scattering by cool ionized clumps ($T\sim 10^4\,$K) in circumgalactic media (CGM). The existence of such clumps are inferred from intervening quasar absorption systems, but have long been something of a theoretical mystery. We consider the implications for compact radio sources of the `fog-like' two-phase model of the circumgalactic medium recently proposed by McCourt et al.(2018). In this model, the CGM consists of a diffuse coronal gas ($T\gtrsim 10^6\,$K) in pressure equilibrium with numerous $\lesssim 1\,$pc scale cool clumps or `cloudlets' formed by shattering in a cooling instability. The areal filling factor of the cloudlets is expected to exceed unity in $\gtrsim 10^{11.5} M_\odot$ haloes, and the ensuing radio-wave scattering is akin to that caused by turbulence in the Galactic warm ionized medium (WIM). If $30\,$per-cent of cosmic baryons are in the CGM, we show that for a cool-gas volume fraction of $f_{\rm v}\sim 10^{-3}$, sources at $z_{\rm s}\sim 1$ suffer angular broadening by $\sim 15\,\mu$as and temporal broadening by $\sim 1\,$ms at $\lambda = 30\,$cm, due to scattering by the clumps in intervening CGM. The former prediction will be difficult to test (the angular broadening will suppress Galactic scintillation only for $<10\,\mu$Jy compact synchrotron sources). However the latter prediction, of temporal broadening of localized fast radio bursts, can constrain the size and mass fraction of cool ionized gas clumps as function of halo mass and redshift, and thus provides a test of the model proposed by McCourt et al.(2018).

The dynamics of the high-velocity compact molecular cloud CO-0.40-0.22 have been interpreted as evidence for a $\sim10^{5}M_{\odot}$ black hole within 60 pc of Sgr A*. Recently, Oka et al. have identified a compact millimetre-continuum source, CO-0.40-0.22*, with this candidate black hole. Here we present a collation of radio and infrared data at this location. ATCA constraints on the radio spectrum, and the detection of a mid-infrared counterpart, are in tension with an Sgr A*-like model for CO-0.40-0.22* despite the comparable bolometric to Eddington luminosity ratios under the IMBH interpretation. A protostellar-disk scenario is, however, tenable. CO-0.40-0.22(*) is associated with an arrowhead-shaped infrared-dark cloud (which we call the Freccia Rossa). Radio-continuum observations reveal a candidate HII region associated with the system. If the $V_{\rm LSR}\approx70$ km s$^{-1}$ systemic velocity of CO-0.40-0.22 is common to the entire Freccia Rossa system, we hypothesise that it is the remnant of a high-velocity cloud that has plunged into the Milky Way from the Galactic halo.

The theoretical and observed populations of pre-cataclysmic variables (pre-CVs) are dominated by systems with low-mass white dwarfs (WDs), while the WD masses in CVs are typically high. In addition, the space density of CVs is found to be significantly lower than theoretical models. We investigate the influence of nova outbursts on the formation and (initial) evolution of CVs. In particular, we calculate the stability of the mass transfer in case all the material accreted on the WD is lost in classical novae, and part of the energy to eject the material comes from a common-envelope like interaction with the companion. In addition, we study the effect of an asymmetry in the mass ejection, that may lead to small eccentricities in the orbit. We find that a common-envelope like ejection significantly decreases the stability of the mass transfer, in particular for low-mass WD. Similarly, the influence of asymmetric mass loss can be important for short-period systems and even more so for low-mass WD, but likely disappears long before the next nova outburst due to orbital circularization. In both cases the mass-transfer rates increase, which may lead to observable (and perhaps already observed) consequences for systems that do survive to become CVs. However, a more detailed investigation of the interaction between nova ejecta and the companion and the evolution of slightly eccentric CVs is needed before definite conclusions can be drawn.

PSR J2129-0429 is a "redback" eclipsing millisecond pulsar binary with an unusually long 15.2 hour orbit. It was discovered by the Green Bank Telescope in a targeted search of unidentified Fermi gamma-ray sources. The pulsar companion is optically bright (mean $m_R = 16.6$ mag), allowing us to construct the longest baseline photometric dataset available for such a system. We present ten years of archival and new photometry of the companion from LINEAR, CRTS, PTF, the Palomar 60-inch, and LCOGT. Radial velocity spectroscopy using the Double-Beam Spectrograph on the Palomar 200-inch indicates that the pulsar is massive: $1.74\pm0.18 M_\odot$. The G-type pulsar companion has mass $0.44\pm0.04 M_\odot$, one of the heaviest known redback companions. It is currently 95\% Roche-lobe filling and only mildly irradiated by the pulsar. We identify a clear 13.1 mmag yr$^{-1}$ secular decline in the mean magnitude of the companion as well as smaller-scale variations in the optical lightcurve shape. This behavior may indicate that the companion is cooling. Binary evolution calculations indicate that PSR J2129-0429 has an orbital period almost exactly at the bifurcation period between systems that converge into tighter orbits as black widows and redbacks and those that diverge into wider pulsar--white dwarf binaries. Its eventual fate may depend on whether it undergoes future episodes of mass transfer and increased irradiation.

The last century has seen enormous progress in our understanding of the Universe. We know the life cycles of stars, the structure of galaxies, the remnants of the big bang, and have a general understanding of how the Universe evolved. We have come remarkably far using electromagnetic radiation as our tool for observing the Universe. However, gravity is the engine behind many of the processes in the Universe, and much of its action is dark. Opening a gravitational window on the Universe will let us go further than any alternative. Gravity has its own messenger: Gravitational waves, ripples in the fabric of spacetime. They travel essentially undisturbed and let us peer deep into the formation of the first seed black holes, exploring redshifts as large as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and unprecedented measurements of black hole masses and spins will make it possible to trace the history of black holes across all stages of galaxy evolution, and at the same time constrain any deviation from the Kerr metric of General Relativity. eLISA will be the first ever mission to study the entire Universe with gravitational waves. eLISA 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 early processes at TeV energies, has guaranteed sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales around singularities and black holes, all the way to cosmological dimensions.

We explore the effect of pulsars, in particular those born with millisecond periods, on their surrounding supernova ejectas. While they spin down, fast-spinning pulsars release their tremendous rotational energy in the form of a relativistic magnetized wind that can affect the dynamics and luminosity of the supernova. We estimate the thermal and non thermal radiations expected from these specific objects, concentrating at times a few years after the onset of the explosion. We find that the bolometric light curves present a high luminosity plateau (that can reach $10^{43-44}\,$erg/s) over a few years. An equally bright TeV gamma-ray emission, and a milder X-ray peak (of order $10^{40-42}\,$erg/s) could also appear a few months to a few years after the explosion, as the pulsar wind nebula emerges, depending on the injection parameters. The observations of these signatures by following the emission of a large number of supernovae could have important implications for the understanding of core-collapse supernovae and reveal the nature of the remnant compact object.

The time-variable electromagnetic sky has been well-explored at a wide range of wavelengths. Numerous high-energy space missions take advantage of the dark Gamma-ray and X-ray sky and utilize very wide field detectors to provide almost continuous monitoring of the entire celestial sphere. In visible light, new wide-field ground-based surveys cover wide patches of sky with ever decreasing cadence, progressing from monthly-weekly time scale surveys to sub-night sampling. In the radio, new powerful instrumentation offers unprecedented sensitivity over wide fields of view, with pathfinder experiments for even more ambitious programs underway. In contrast, the ultra-violet (UV) variable sky is relatively poorly explored, even though it offers exciting scientific prospects. Here, we review the potential scientific impact of a wide-field UV survey on the study of explosive and other transient events, as well as known classes of variable objects, such as active galactic nuclei and variable stars. We quantify our predictions using a fiducial set of observational parameters which are similar to those envisaged for the proposed ULTRASAT mission. We show that such a mission would be able to revolutionize our knowledge about massive star explosions by measuring the early UV emission from hundreds of events, revealing key physical parameters of the exploding progenitor stars. Such a mission would also detect the UV emission from many tens of tidal-disruption events of stars by super massive black holes at galactic nuclei and enable a measurement of the rate of such events. The overlap of such a wide-field UV mission with existing and planned gravitational-wave and high-energy neutrino telescopes makes it especially timely.

We have investigated the relationship between the kinematics and mass of young (<3x10^8 years) white dwarfs using proper motions. Our sample is taken from the colour selected catalogues of SDSS (Eisenstein et al. 2006) and the Palomar-Green Survey (Liebert, Bergeron & Holberg 2005), both of which have spectroscopic temperature and gravity determinations. We find that the dispersion decreases with increasing white dwarf mass. This can be explained as a result of less scattering by objects in the Galactic disk during the shorter lifetime of their more massive progenitors. A direct result of this is that white dwarfs with high mass have a reduced scale height, and hence their local density is enhanced over their less massive counterparts. In addition, we have investigated whether the kinematics of the highest mass white dwarfs (>0.95Msun) are consistent with the expected relative contributions of single star evolution and mergers. We find that the kinematics are consistent with the majority of high-mass white dwarfs being formed through single star evolution.

We report on the discovery and observations of the extremely luminous optical transient CSS100217:102913+404220 (CSS100217 hereafter). Spectroscopic observations show this transient was coincident with a galaxy at redshift z=0.147, and reached an apparent magnitude of V ~ 16.3. After correcting for foreground Galactic extinction we determine the absolute magnitude to be M_V =-22.7 approximately 45 days after maximum light. Based on our unfiltered optical photometry the peak optical emission was L = 1.3 x 10^45 erg s^-1, and over a period of 287 rest-frame days had an integrated bolometric luminosity of 1.2 x 10^52 erg. Analysis of the pre-outburst SDSS spectrum of the source shows features consistent with a Narrow-line Seyfert1 (NLS1) galaxy. High-resolution HST and Keck followup observations show the event occurred within 150pc of nucleus of the galaxy, suggesting a possible link to the active nuclear region. However, the rapid outburst along with photometric and spectroscopic evolution are much more consistent with a luminous supernova. Line diagnostics suggest that the host galaxy is undergoing significant star formation. We use extensive follow-up of the event along with archival CSS and SDSS data to investigate the three most likely sources of such an event; 1) an extremely luminous supernova; 2) the tidal disruption of a star by the massive nuclear black hole; 3) variability of the central AGN. We find that CSS100217 was likely an extremely luminous type IIn supernova that occurred within range of the narrow-line region of an AGN. We discuss how similar events may have been missed in past supernova surveys because of confusion with AGN activity.

Black holes of mass M must have a spin angular momentum S below the Kerr limit chi = S/M^2 < 1, but whether astrophysical black holes can attain this limiting spin depends on their accretion history. Gas accretion from a thin disk limits the black-hole spin to chi_gas < 0.9980 +- 0.0002, as electromagnetic radiation from this disk with retrograde angular momentum is preferentially absorbed by the black hole. Extrapolation of numerical-relativity simulations of equal-mass binary black-hole mergers to maximum initial spins suggests these mergers yield a maximum spin chi_eq < 0.95. Here we show that for smaller mass ratios q = m/M << 1, the superradiant extraction of angular momentum from the larger black hole imposes a fundamental limit chi_lim < 0.9979 +- 0.0001 on the final black-hole spin even in the test-particle limit q -> 0 of binary black-hole mergers. The nearly equal values of chi_gas and chi_lim imply that measurement of supermassive black-hole spins cannot distinguish a black hole built by gas accretion from one assembled by the gravitational inspiral of a disk of compact stellar remnants. We also show how superradiant scattering alters the mass and spin predicted by models derived from extrapolating test-particle mergers to finite mass ratios.

Electromagnetic observations over the last 15 years have yielded a growing appreciation for the importance of supermassive black holes (SMBH) to the evolution of galaxies, and for the intricacies of dynamical interactions in our own Galactic center. Here we show that future low-frequency gravitational wave observations, alone or in combination with electromagnetic data, will open up unique windows to these processes. In particular, gravitational wave detections in the 10^{-5}-10^{-1} Hz range will yield SMBH masses and spins to unprecedented precision and will provide clues to the properties of the otherwise undetectable stellar remnants expected to populate the centers of galaxies. Such observations are therefore keys to understanding the interplay between SMBHs and their environments.

The principal goal of this whitepaper is not so much to demonstrate that gravitational wave detectors like LIGO and LISA will help answer many central questions in astronomy and astrophysics, but to make the case that they can help answer a far greater range of questions if we prepare to make the (sometimes substantial) effort to identify electromagnetic counterparts to the gravitational wave sources.

This White Paper to the National Academy of Sciences Astro2010 Decadal Review Committee outlines some of the outstanding questions regarding the assembly history of Massive Black Holes in the nuclei of galaxies and the revolutionary contributions anticipated in this field from low-frequency gravitational wave astronomy.

This white paper briefly describes the astrophysics of ultra-compact binaries, with emphasis of the challenges and opportunities in the next decade.

We study a non-dissipative hydrodynamical mechanism that can stabilize the spin of the accretor in an ultra-compact double white dwarf binary. This novel synchronization mechanism relies on a nonlinear wave interaction spinning down the background star. The essential physics of the synchronization mechanism is summarized as follows. As the compact binary coalesces due to gravitational wave emission, the largest star eventually fills its Roche lobe and accretion starts. The accretor then spins up due to infalling material and eventually reaches a spin frequency where a normal mode of the star is resonantly driven by the gravitational tidal field of the companion. If the resonating mode satisfies a set of specific criteria, which we elucidate in this paper, it exchanges angular momentum with the background star at a rate such that the spin of the accretor locks at this resonant frequency, even though accretion is ongoing. Some of the accreted angular momentum that would otherwise spin up the accretor is fed back into the orbit through this resonant tidal interaction. Two modes capable of stabilizing the accretor's spin are the l=4,m=2 and l=5,m=3 CFS unstable hybrid r-modes, which stabilize the spin of the accretor at frequency 2.6 and 1.5 times the binary's orbital frequency respectively. Since the stabilization mechanism relies on continuously driving a mode at resonance, its lifetime is limited since eventually the mode amplitude saturates due to non-linear mode-mode coupling. Rough estimates of the lifetime of the effect lie from a few orbits to millions of years.

Over a decade ago it was established that the remarkable high energy transients, known as soft gamma-ray repeaters (SGRs), are a Galactic population and originate from neutron stars with intense (<~ 10^15 G) magnetic fields ("magnetars"). On 27 December 2004 a giant flare (fluence >~ 0.3 erg/cm^2) was detected from SGR 1806-20. Here we report the discovery of a fading radio counterpart. We began a monitoring program from 0.2GHz to 250GHz and obtained a high resolution 21-cm radio spectrum which traces the intervening interstellar neutral Hydrogen clouds. Analysis of the spectrum yields the first direct distance measurement of SGR 1806-20. The source is located at a distance greater than 6.4 kpc and we argue that it is nearer than 9.8 kpc. If true, our distance estimate lowers the total energy of the explosion and relaxes the demands on theoretical models. The energetics and the rapid decay of the radio source are not compatible with the afterglow model that is usually invoked for gamma-ray bursts. Instead we suggest that the rapidly decaying radio emission arises from the debris ejected during the explosion.

The final merger of a pair of massive black holes in a galactic nucleus is compelled by gravitational radiation. Gravitational waves from the mergers of black holes of masses (10^5-10^7)(1+z)^{-1} Msun at redshifts of 1-20 will be readily detectable by the Laser Interferometer Space Antenna (LISA), but an electromagnetic afterglow would be helpful in pinpointing the source and its redshift. Long before the merger, the binary "hollows out" any surrounding gas and shrinks slowly compared to the viscous timescale of a circumbinary disk. The inner gas disk is truncated at the radius where gravitational torque from the binary balances the viscous torque, and accretion onto the black holes is diminished. Initially, the inner truncation radius is able to follow the shrinking binary inward. But eventually the gravitational radiation timescale becomes shorter than the viscous timescale in the disk, leading to a merged black hole surrounded by a hollow disk of gas. We show that the subsequent viscous evolution of the hollow, radiation-pressure dominated disk will create a ~10^{43.5}(M/10^6Msun) ergs s^{-1} X-ray source on a timescale \~7(1+z)(M/10^6Msun)^{1.32} yr. This justifies follow-up monitoring of gravitational wave events with next-generation X-ray observatories. Analysis of the detailed light curve of these afterglows will yield new insights into the subtle physics of accretion onto massive black holes.

Gravitational waves and X-ray flares are expected from tidal disruption of stars by a massive black hole. Using a relativistic smoothed particle hydrodynamics code, we investigate the fate of main sequence and Helium stars in plunge orbits passing near Schwarzschild or a Kerr black holes of 10^5-6 Msun. We show that quadrupole gravitational waves emitted during the tidal disruption process are described reasonable well by a point particle approximation even in the strong encounter case. An additional hydrodynamic calculation based on the Godunov method indicates that shocks develop for sufficiently high tidal compressions. The shock-heating results in an X-ray flare, which for solar-type stars disrupted by 10^6Msun black holes is in the keV range, associated with the gravitational wave signal. The hardness and duration of the X-ray flare may serve as a diagnostic of the mass of the central black hole.

One of the most exciting prospects for the LISA gravitational wave observatory is the detection of gravitational radiation from the inspiral of a compact object into a supermassive black hole. The large inspiral parameter space and low amplitude of the signal makes detection of these sources computationally challenging. We outline here a first cut data analysis scheme that assumes realistic computational resources. In the context of this scheme, we estimate the signal-to-noise ratio that a source requires to pass our thresholds and be detected. Combining this with an estimate of the population of sources in the Universe, we estimate the number of inspiral events that LISA could detect. The preliminary results are very encouraging -- with the baseline design, LISA can see inspirals out to a redshift z=1 and should detect over a thousand events during the mission lifetime.

We use a population synthesis approach to characterise, as a function of cosmic time, the extragalactic close binary population descended from stars of low to intermediate initial mass. The unresolved gravitational wave (GW) background due to these systems is calculated for the 0.1-10 mHz frequency band of the planned Laser Interferometer Space Antenna (LISA). This background is found to be dominated by emission from close white dwarf-white dwarf pairs. The spectral shape can be understood in terms of some simple analytic arguments. To quantify the astrophysical uncertainties, we construct a range of evolutionary models which produce populations consistent with Galactic observations of close WD-WD binaries. The models differ in binary evolution prescriptions as well as initial parameter distributions and cosmic star formation histories. We compare the resulting background spectra, whose shapes are found to be insensitive to the model chosen, and different to those found recently by Schneider et al. (2001). From this set of models, we constrain the amplitude of the extragalactic background to be 1E-12 < Omega(1 mHz) < 6E-12, in terms of Omega(f), the fraction of closure density received in gravitational waves in the logarithmic frequency interval around f.

There is an extremely simple relationship between the spectrum of the gravitational wave background produced by a cosmological distribution of discrete gravitational wave sources, the total time-integrated energy spectrum of an individual source, and the present-day comoving number density of remnants. Stated in this way, the background is entirely independent of the cosmology, and only weakly dependent on the evolutionary history of the sources. This relationship allows one easily to compute the amplitude and spectrum of cosmic gravitational wave backgrounds from a broad range of astrophysical sources, and to evaluate the uncertainties therein.

The spatial and temporal coincidence between the gamma-ray burst (GRB) 980425 and supernova (SN) 1998bw has prompted speculation that there exists a class of GRBs produced by SNe (``S-GRBs''). Robust arguments for the existence of a relativistic shock have been presented on the basis of radio observations. A physical model based on the radio observations lead us to propose the following characteristics of supernovae GRBs (S-GRBs): 1) prompt radio emission and implied brightness temperature near or below the inverse Compton limit, 2) high expansion velocity of the optical photosphere as derived from lines widths and energy release larger than usual, 3) no long-lived X-ray afterglow, and 4) a single pulse (SP) GRB profile. Radio studies of previous SNe show that only type Ib and Ic potentially satisfy the first condition. Accordingly we have investigated proposed associations of GRBs and SNe finding no convincing evidence (mainly to paucity of data) to confirm any single connection of a SN with a GRB. If there is a more constraining physical basis for the burst time-history of S-GRBs beyond that of the SP requirement, we suggest the 1% of light curves in the BATSE catalogue similar to that of GRB 980425 may constitute the subclass. Future optical follow-up of bursts with similar profiles should confirm if such GRBs originate from some fraction of SN type Ib/Ic.

Supernova SN 1998bw exploded in the same direction and at the same time as the gamma-ray burst GRB 980425. Here we report radio observations of this type Ic supernova, beginning 4 days after the gamma-ray burst. At its peak the radio source is the most luminous ever seen from a supernova, nu L_nu = 4 x 10^38 erg/s at 5 GHz. More remarkably, the traditional synchrotron interpretation of the radio emission requires that the radio source be expanding at an apparent velocity of at least twice the speed of light, indicating that this supernova was accompanied by a shock wave moving at relativistic speed. The energy U_e associated with the radio-emitting relativistic electrons must lie between 10^49 erg <U_e< 10^52 erg, and thus represents a significant fraction of the total kinetic energy ~10^51 erg associated with supernova explosions. The presence of a relativistic shock in SN 1998bw suggests a physical connection with the gamma-ray burst GRB 980425. We argue that this represents a second class of gamma-ray burst, with much lower emitted energy ~10^48 erg in gamma-rays than the two powerful ~10^53 erg high-redshift gamma-ray bursts previously identified.

RADIO pulsars are thought to born with spin periods of 0.02-0.5 s and space velocities of 100-1000 km/s, and they are inferred to have initial dipole magnetic fields of 10^{11}-10^{13}. The average space velocity of a normal star in the Milky Way is only 30 km/s, which means that pulsars must receive a substantial 'kick' at birth. Here we propose that the birth characteristics of pulsars have a simple physical connection with each other. Magnetic fields maintained by differential rotation between the core and envelope of the progenitor would keep the whole star in a state of approximately uniform rotation until 10 years before the explosion. Such a slowly rotating core has 1000 times less angular momentum than required to explain the rotation of pulsars. Although the specific physical process that 'kicks' the neutron star at birth has not been identified, unless its force is exerted exactly head-on, it will also cause the neutron star to rotate. We identify this process as the origin of the spin of pulsars. Such kicks will cause a correlation between the velocity and spin vectors of pulsars. We predict that many neutron stars are born with periods longer than 2 s, and never become radio pulsars.

Binary stars play a critical role in determining both the dynamical evolution and the observational appearance of globular star clusters. We present a summary of three talks presented at IAU JD15, reviewing the current observational and theoretical status of binaries and related ``exotic'' stars in the Galactic globular cluster system.

We use the grid of models described in paper~I to analyse those millisecond pulsar binaries whose secondaries have been studied optically. In particular, we find cooling ages for these binary systems that range from $< 1 \rm Gyr$ to $\sim \rm 15 Gyr$. Comparison of cooling ages and characteristic spin down ages allows us to constrain the initial spin periods and spin-up histories for individual systems, showing that at least some millisecond pulsars had sub-Eddington accretion rates and long magnetic field decay times.

The presence of low mass, degenerate secondaries in millisecond pulsar binaries offers the opportunity to determine an age for the binary system independent of the rotational properties of the pulsar. To this end, we present here a detailed calculation of the evolution of a grid of low mass ($<0.5 M_{\odot}$) helium core white dwarfs. We investigate the effects of different Hydrogen layer masses and provide results for well-known optical band-passes. We supplement the OPAL opacity calculations with our own calculations for low effective temperatures ($ T_{eff} < 6000 K$) and also provide fitting formulae for the gravity as a function of mass and effective temperature. In paper II we shall apply these results to individual cases.

Dwarf, irregular and infrared-luminous starburst galaxies are all known to have "steep" luminosity functions with faint-end behavior roughly phi(L)-propto- L^(-1.8). This form is exactly what is expected if the luminosities of these objects fade with time as L-propto-t^(-1.3), because the objects spend more time at low luminosities than high, even if they form with a wide range of initial masses. Models of young stellar populations show this fading behavior when the star formation has occured in a single, short, recent burst. Steep luminosity functions therefore do not require steep mass functions if the galaxies are powered by fading bursts. The local-galaxy H-alpha luminosity function---which is less steep than L^(-1.8)---is also well-fit by this mechanism, because ionizing photon flux fades much more quickly than broad-band optical luminosity. An age-luminosity relation and a wavelength-dependence of the luminosity function are both predicted. In the context of this mechanism, the slope of the luminosity function provides a constraint on the stellar initial mass function in the bursts.

We analyse the sample of pulsar proper motions, taking detailed account of the selection effects of the original surveys. We treat censored data using survival statistics. From a comparison of our results with Monte Carlo simulations, we find that the mean birth speed of a pulsar is 250-300 km/s, rather than the 450 km/s foundby Lyne & Lorimer (1994). The resultant distribution is consistent with a maxwellian with dispersion $ \sigma_v = 190 km/s$. Despite the large birth velocities, we find that the pulsars with long characteristic ages show the asymmetric drift, indicating that they are dynamically old. These pulsars may result from the low velocity tail of the younger population, although modified by their origin in binaries and by evolution in the galactic potential.

We have used WFPC2 to construct B, V color-magnitude diagrams of four metal-rich globular clusters, NGC 104 (47 Tuc), NGC 5927, NGC 6388, and NGC 6441. All four clusters have well populated red horizontal branches (RHB), as expected for their metallicity. However, NGC 6388 and 6441 also exhibit a prominent blue HB (BHB) extension, including stars reaching as faint in V as the turnoff luminosity. This discovery demonstrates directly for the first time that a major population of hot HB stars can exist in old, metal-rich systems. This may have important implications for the interpretation of the integrated spectra of elliptical galaxies. The cause of the phenomenon remains uncertain. We examine the possibility that NGC 6388 and 6441 are older than the other clusters, but a simple difference in age may not be sufficient to produce the observed distributions along the HB. The high central densities in NGC 6388 and 6441 suggest that the existence of the blue HB (BHB) tails might be caused by stellar interactions in the dense cores of these clusters, which we calculate to have two of the highest collision rates among globular clusters in the Galaxy. Tidal collisions might act in various ways to enhance loss of envelope mass, and therefore populate the blue side of the HB. However, the relative frequency of tidal collisions does not seem large enough (compared to that of the clusters with pure RHBs) to account for such a drastic difference in HB morphology. While a combination of an age difference and dynamical interactions may help, prima facie the lack of a radial gradient in the BHB/RHB star ratio seems to argue against dynamical effects playing a role.

The evolution of a binary star system by various analytic approximations of mass transfer is discussed, with particular attention payed to the stability of these processes against runaway on the thermal and dynamical timescales of the mass donating star. Mass transfer in red giant - neutron star binary systems is used as a specific example of such mass transfer, and is investigated. Hjellming and Webbink's (1987) results on the dynamic timescale response of a convective star with a core to mass loss are applied, with new results. It is found that mass transfer is usually stable, so long as the the wind's specific angular momentum does not exceed the angular momentum per reduced mass of the system. This holds for both dynamical and thermal timescales. Those systems which are not stable will usually transfer mass on the thermal timescale. Included are graphs illustrating the variation of the Roche radius exponent with mass ratio in the binary, for various parameters in the non-conservative mass transfer, as well as evolutionary paths of interacting red giant neutron star binaries.

We present distributions of colors of stars along the horizontal branch of the globular cluster NGC 2808, from Hubble Space Telescope WFPC2 imaging in B, V, and an ultraviolet filter (F218W). This cluster's HB is already known to be strongly bimodal, with approximately equal-sized HB populations widely separated in the color-magnitude diagram. Our images reveal a long blue tail with two gaps, for a total of four nearly distinct HB groups. These gaps are very narrow, corresponding to envelope-mass differences of only \sim 0.01 Msun. This remarkable multimodality may be a signature of mass-loss processes, subtle composition variations, or dynamical effects; we briefly summarize the possibilities. The existence of narrow gaps between distinct clumps on the HB presents a challenge for models that attempt to explain HB bimodality or other peculiar HB structures.

We model the dynamics of test binaries in isotropic, multi-mass models of galactic globular clusters. The evolution of binary orbits through the cluster potentials is modeled, including second order diffusion terms, and probabilities for close encounters with field stars are calculated. We carry out Monte Carlo simulations of the effects of the binary--single star encounters on the binary population and distribution in the cluster, and estimate the collision rate for different stellar populations in globular clusters with different structural parameters. Assuming a Salpeter IMF, for low concentration clusters the core encounter rate is dominated by turnoff mass main--sequence stars and medium mass white dwarfs. For high concentration, high density clusters the encounter probabilities are increasingly dominated by neutron stars and heavy white dwarfs. Hence we predict a smaller ratio of blue stragglers and cataclysmic variables to pulsars in high concentration clusters. The total number of millisecond pulsars, and the ratio of single to binary pulsars, is broadly consistent with the observed population, suggesting the binary--single star encounters contribute significantly to the pulsar formation rate in globular clusters, for the whole range of globular cluster types. The number of millisecond pulsars and the ratio of pulsars in different globular clusters is best explained by a total binary fraction comparable to that of the galaxy, and a modest number of primordial neutron stars in the globular clusters.