Among the potential gravitational wave (GW) sources for the upcoming space-based interferometer LISA are extreme- or intermediate-mass ratio inspirals (EMRI/IMRIs). These events involve the coalescence of supermassive black holes in the mass range $10^5M_{\odot}\lesssim M \lesssim 10^7M_{\odot}$ with companion BHs of much lower masses. A subset of E/IMRIs are expected to occur in the accretion discs of active galactic nuclei (AGN). Previous work has shown that torques exerted by the disc can interfere with the inspiral and cause a phase shift in the GW waveform detectable by LISA. Here we use a suite of two-dimensional hydrodynamical simulations with the moving-mesh code DISCO to present a systematic study of disc torques. We measure torques on an inspiraling BH and compute the corresponding waveform deviations as a function of the binary mass ratio $q\equiv M_2/M_1$, the disc viscosity ($\alpha$), and gas temperature (or equivalently Mach number; $\mathcal{M}$). We find that the absolute value of the gas torques is within an order of magnitude of previously determined planetary migration torques, but their precise value and sign depends non-trivially on the combination of these parameters, the inspiral rate, and the accretion rate onto the satellite BH. The gas imprint is generally detectable by LISA for binaries embedded in AGN discs with surface densities above $\Sigma_0\ge10^{4-6} \rm \, g \, cm^{-2}$, depending on $q$, $\alpha$ and $\mathcal{M}$. We find that the deviations are most pronounced in discs with higher viscosities, and for E/IMRIs detected at frequencies where LISA is most sensitive. Torques in colder discs exhibit a noticeable dependence on the GW-driven inspiral rate as well as strong fluctuations at late stages of the inspiral. Our results further suggest that LISA may be able to place constraints on AGN disc parameters and the physics of disc-satellite interaction.

Using numerical hydrodynamics calculations and a novel method for densely sampling parameter space, we measure the accretion and torque on a binary system from a circumbinary disk. In agreement with previous studies, we find that the net torque on the binary is positive for mass ratios close to unity, and that accretion always drives the binary towards equal mass. Accretion variability depends sensitively on the numerical sink prescription, but the torque and relative accretion onto each component do not depend on the sink timescale. Positive torque and highly variable accretion occurs only for mass ratios greater than around $0.05$. This means that for mass ratios below $0.05$, the binary would migrate inward until the secondary accreted sufficient mass, after which it would execute a U-turn and migrate outward. We explore a range of viscosities, from $\alpha = 0.03$ to $\alpha = 0.15$, and find that this outward torque is proportional to the viscous torque, simply proportional to viscosity in this range. Dependence of accretion and torque on mass ratio is explored in detail, densely sampling mass ratios between $0.01$ and unity. For mass ratio $q > 0.6$, accretion variability is found to exhibit a distinct sawtooth pattern, typically with a five-orbit cycle that provides a "smoking gun" prediction for variable quasars observed over long periods, as a potential means to confirm the presence of a binary.

The coalescence of a compact object with a $10^{4}-10^{7} {\rm M_\odot}$ supermassive black hole (SMBH) produces mHz gravitational waves (GWs) detectable by the future Laser Interferometer Space Antenna (LISA). If such an inspiral occurs in the accretion disc of an active galactic nucleus (AGN), the gas torques imprint a small deviation in the GW waveform. Here we present two-dimensional hydrodynamical simulations with the moving-mesh code DISCO of a BH inspiraling at the GW rate in a binary system with a mass ratio $q\!=\!M_2/M_1\!=\!10^{-3}$, embedded in an accretion disc. We assume a locally isothermal equation of state for the gas (with Mach number $\mathcal{M}=20$) and implement a standard $\alpha$-prescription for its viscosity (with $\alpha = 0.03$). We find disc torques on the binary that are weaker than in previous semi-analytic toy models, and are in the opposite direction: the gas disc slows down, rather than speeds up the inspiral. We compute the resulting deviations in the GW waveform, which scale linearly with the mass of the disc. The SNR of these deviations accumulates mostly at high frequencies, and becomes detectable in a 5-year LISA observation if the total phase shift exceeds a few radians. We find that this occurs if the disc surface density exceeds $\Sigma_0 \gtrsim 10^{2-3}\rm g\,cm^{-2}$, as may be the case in thin discs with near-Eddington accretion rates. Since the characteristic imprint on the GW signal is strongly dependent on disc parameters, a LISA detection of an intermediate mass ratio inspiral would probe the physics of AGN discs and migration.

The discovery of GeV gamma-rays from classical novae has led to a reassessment of these garden-variety explosions, and highlighted their importance for understanding radiative shocks, particle acceleration, and dust formation in more exotic, distant transients. Recent collaboration between observers and theorists has revealed that shocks in novae are energetically important, and can even dominate their bolometric luminosity. Shocks may also explain long-standing mysteries in novae such as dust production, super-Eddington luminosities, and `flares' in optical light curves. Here, we highlight the multi-wavelength facilities of the next decade that will further test our nova shock model and fulfill the promise of novae as powerful astrophysical laboratories.

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.

Classical novae commonly show evidence of rapid dust formation within months of the outburst. However, it is unclear how molecules and grains are able to condense within the ejecta, given the potentially harsh environment created by ionizing radiation from the white dwarf. Motivated by the evidence for powerful radiative shocks within nova outflows, we propose that dust formation occurs within the cool, dense shell behind these shocks. We incorporate a simple molecular chemistry network and classical nucleation theory with a model for the thermodynamic evolution of the post-shock gas, in order to demonstrate the formation of both carbon and forsterite ($\rm Mg_2SiO_4$) grains. The high densities due to radiative shock compression ($n \sim 10^{14}$ cm$^{-3}$) result in CO saturation and rapid dust nucleation. Grains grow efficiently to large sizes $\gtrsim 0.1\mu$m, in agreement with IR observations of dust-producing novae, and with total dust masses sufficient to explain massive extinction events such as V705 Cas. As in dense stellar winds, dust formation is CO-regulated, with carbon-rich flows producing carbon-rich grains and oxygen-rich flows primarily forming silicates. CO is destroyed by non-thermal particles accelerated at the shock, allowing additional grain formation at late times, but the efficiency of this process appears to be low. Given observations showing that individual novae produce both carbonaceous and silicate grains, we concur with previous works attributing this bimodality to chemical heterogeneity of the ejecta. Nova outflows are diverse and inhomogeneous, and the observed variety of dust formation events can be reconciled by different abundances, the range of shock properties, and the observer viewing angle. The latter may govern the magnitude of extinction, with the deepest extinction events occurring for observers within the binary equatorial plane.