Short gamma-ray bursts are thought to result from the mergers of two neutron stars or a neutron star and stellar mass black hole. The final stages of the merger are generally accompanied by the production of one or more tidal "tails" of ejecta, which fall back onto the remnant-disc system at late times. Using the results of a linear stability analysis, we show that if the material comprising these tails is modeled as adiabatic and the effective adiabatic index satisfies $\gamma \ge 5/3$, then the tails are gravitationally unstable and collapse to form small-scale knots. We analytically estimate the properties of these knots, including their spacing along the tidal tail and the total number produced, and their effect on the mass return rate to the merger remnant. We perform hydrodynamical simulations of the disruption of a polytropic (with the polytropic and adiabatic indices $\gamma$ equal), $\gamma =2$ neutron star by a black hole, and find agreement between the predictions of the linear stability analysis and the distribution of knots that collapse out of the instability. The return of these knots to the black hole induces variability in the fallback rate, which can manifest as variability in the lightcurve of the GRB and -- depending on how rapidly the instability operates -- the prompt emission. The late-time variability induced by the return of these knots is also consistent with the extended emission observed in some GRBs.

A tidal disruption event (TDE) occurs when a star plunges through a supermassive black hole's tidal radius, at which point the star's self-gravity is overwhelmed by the tidal gravity of the black hole. In a partial TDE, where the star does not reach the full disruption radius, only a fraction of the star's mass is tidally stripped while the rest remains intact in the form of a surviving core. Analytical arguments have recently suggested that the temporal scaling of the fallback rate of debris to the black hole asymptotes to $t^{-9/4}$ for partial disruptions, effectively independently of the mass of the intact core. We present hydrodynamical simulations that verify the existence of this predicted, $t^{-9/4}$ scaling. We also define a break timescale -- the time at which the fallback rate transitions from a $t^{-5/3}$ scaling to the characteristic $t^{-9/4}$ scaling -- and measure this break timescale as a function of the impact parameter and the surviving core mass. These results deepen our understanding of the properties and breadth of possible fallback curves expected from TDEs and will therefore facilitate more accurate interpretation of data from wide-field surveys.

White dwarfs are routinely observed to have polluted atmospheres, and sometimes significant infrared excesses, that indicate ongoing accretion of circumstellar dust and rocky debris. Typically this debris is assumed to be in the form of a (circular) disc, and to originate from asteroids that passed close enough to the white dwarf to be pulled apart by tides. However, theoretical considerations suggest that the circularisation of the debris, which initially occupies highly eccentric orbits, is very slow. We therefore hypothesise that the observations may be readily explained by the debris remaining on highly eccentric orbits, and we explore the properties of such debris. For the generic case of an asteroid originating at several au from the white dwarf, we find that all of the tidal debris is always bound to the white dwarf and that the orbital energy distribution of the debris is narrow enough that it executes similar elliptical orbits with only a narrow spread. Assuming that the tidal field of the white dwarf is sufficient to minimise the effects of self-gravity and collisions within the debris, we estimate the time over which the debris spreads into a single elliptical ring, and we generate toy spectra and lightcurves from the initial disruption to late times when the debris distribution is essentially time steady. Finally we speculate on the connection between these simple considerations and the observed properties of these systems, and on additional physical processes that may change this simple picture.

Many accretion discs are thought to be warped. Recent hydrodynamical simulations show that (i) discs can break into distinct planes when the amplitude of an imposed warp is sufficiently high and the viscosity sufficiently low, and that (ii) discs can tear up into discrete rings when an initially planar disc is subject to a forced precession. Previously, we investigated the local stability of isolated, Keplerian, warped discs in order to understand the physics causing an accretion disc to break into distinct planes, finding that anti-diffusion of the warp amplitude is the underlying cause. Here, we explore the behaviour of this instability in disc regions where the rotation profile deviates from Keplerian. We find that at small warp amplitudes non-Keplerian rotation can stabilize the disc by increasing the critical warp amplitude for instability, while at large warp amplitudes non-Keplerian rotation can lead to an increased growth rate for discs that are unstable. Tidal effects on discs in binary systems are typically weak enough such that the disc remains close to Keplerian rotation. However, the inner regions of discs around black holes are strongly affected, with the smallest radius at which the disc can break into discrete planes being a function of the black hole spin. We suggest that interpreting observed frequencies in the power spectra of light curves from accreting compact objects as nodal and apsidal precession of discrete orbits requires an instability that can break the disc into discrete rings such as the one explored here.

In order to provide an explanation for the unexpected radial brightness distribution of the steady accretion discs seen in nova-like variables, Nixon & Pringle (2019) proposed that the accretion energy is redistributed outwards by means of strong, magnetically driven, surface flows. In this paper we note that the "powerful, rotating disc winds" observed in the soft states of black hole X-ray binaries, and also in the disc around a magnetized neutron star in Her X-1, have the properties of the outflows postulated by Nixon & Pringle to exist in the nova-like variable accretion discs around white dwarfs. The relevant properties are that the flows are not winds, but are, instead, bound flows (traveling at less than the escape velocity) and that the mass fluxes in the flows are a substantial fraction of the accretion rate in the disc.

Filamentary structures, or long and narrow streams of material, arise in many areas of astronomy. Here we investigate the stability of such filaments by performing an eigenmode analysis of adiabatic and polytropic fluid cylinders, which are the cylindrical analog of spherical polytropes. We show that these cylinders are gravitationally unstable to perturbations along the axis of the cylinder below a critical wavenumber $k_{\rm crit} \simeq few$, where $k_{\rm crit}$ is measured relative to the radius of the cylinder. Below this critical wavenumber perturbations grow as $\propto e^{\sigma_{\rm u}\tau}$, where $\tau$ is time relative to the sound crossing time across the diameter of the cylinder, and we derive the growth rate $\sigma_{\rm u}$ as a function of wavenumber. We find that there is a maximum growth rate $\sigma_{\rm max} \sim 1$ that occurs at a specific wavenumber $k_{\rm max} \sim 1$, and we derive the growth rate $\sigma_{\rm max}$ and the wavenumbers $k_{\rm max}$ and $k_{\rm crit}$ for a range of adiabatic indices. To the extent that filamentary structures can be approximated as adiabatic and fluid-like, our results imply that these filaments are unstable without the need to appeal to magnetic fields or external media. Further, the objects that condense out of the instability of such filaments are separated by a preferred length scale, form over a preferred timescale, and possess a preferred mass scale.

We present the results of a multi-wavelength follow up campaign for the luminous nuclear transient Gaia16aax, which was first identified in January 2016. The transient is spatially consistent with the nucleus of an active galaxy at z=0.25, hosting a black hole of mass $\rm \sim6\times10^8M_\odot$. The nucleus brightened by more than 1 magnitude in the Gaia G-band over a timescale of less than one year, before fading back to its pre-outburst state over the following three years. The optical spectra of the source show broad Balmer lines similar to the ones present in a pre-outburst spectrum. During the outburst, the $\rm H\alpha$ and $\rm H\beta$ emission lines develop a secondary peak. We also report on the discovery of two transients with similar light curve evolution and spectra: Gaia16aka and Gaia16ajq. We consider possible scenarios to explain the observed outbursts. We exclude that the transient event could be caused by a microlensing event, variable dust absorption or a tidal encounter between a neutron star and a stellar mass black hole in the accretion disk. We consider variability in the accretion flow in the inner part of the disk, or a tidal disruption event of a star $\geq 1 M_{\odot}$ by a rapidly spinning supermassive black hole as the most plausible scenarios. We note that the similarity between the light curves of the three Gaia transients may be a function of the Gaia alerts selection criteria.

Mergers are fundamental to the standard paradigm of galaxy evolution, and provide a natural formation mechanism for supermassive black hole binaries. The formation process of such a binary can have a direct impact on the rate at which stars are tidally disrupted by one or the other black hole, and the luminous signature of the tidal disruption itself can have distinct imprints of a binary companion. In this chapter we review our current understanding of the influence of black hole binarity on the properties of tidal disruption events. We discuss the rates of tidal disruption by supermassive black hole binaries, the impact of a second black hole on the fallback of debris and the formation of an accretion flow, and the prospects for detection of tidal disruption events by supermassive black hole binaries.

The tidal disruption of a star by a supermassive black hole, and the subsequent accretion of the disrupted debris by that black hole, offers a direct means to study the inner regions of otherwise-quiescent galaxies. These tidal disruption events (TDEs) are being discovered at an ever-increasing rate. We present a model for the evolution of the tidally-disrupted debris from a partial TDE, in which a stellar core survives the initial tidal encounter and continues to exert a gravitational influence on the expanding stream of tidally-stripped debris. We use this model to show that the asymptotic fallback rate of material to the black hole in a partial TDE scales as $\propto t^{-2.26\pm0.01}$, and is effectively independent of the mass of the core that survives the encounter; we also estimate the rate at which TDEs approach this asymptotic scaling as a function of the core mass. These findings suggest that the late-time accretion rate onto a black hole from a TDE either declines as $t^{-5/3}$ if the star is completely disrupted or $t^{-9/4}$ if a core is left behind. We emphasize that previous investigations have not recovered this result due to the assumption of a Keplerian energy-period relationship for the debris orbits, which is no longer valid when a surviving core generates a time-dependent, gravitational potential. This dichotomy of fallback rates has important implications for the characteristic signatures of TDEs in the current era of wide-field surveys.

We investigate the impact of time-resolved `gradual' stellar feedback processes in high redshift dwarf spheroidal galaxies. Here `gradual' feedback refers to individual stellar feedback events which deposit energy over a period of time. We conduct high resolution hydrodynamical simulations of dwarf spheroidal galaxies with halo masses of 10$^7$ M$_{\odot}$ - 10$^8$ M$_{\odot}$, based on z = 6 progenitors of the Milky Way's dwarf spheroidal galaxies. We also include a novel feedback prescription for individual massive stars, which includes stellar winds and a HMXB (High Mass X-ray Binary) phase, on top of supernovae. We find the mass of gas unbound across a 1 Gyr starburst is uniformly lowered if gradual feedback mechanisms are included across the range of metallicities, halo concentration parameters and galaxy masses studied here. Furthermore, we find including gradual feedback in the smallest galaxies delays the unbinding of the majority of the gas and facilitates the production of `chimneys' in the dense shell surrounding the feedback generated hot, pressurised `superbubble'. These `chimneys' vent hot gas from the galaxy interior, lowering the temperature of the central 10 kpc of the gaseous halo. Additionally, we find radiative cooling has little effect on the energetics of simulations that include a short, violent starburst compared with those that have a longer, less concentrated starburst. Finally, we investigate the relative impact of HMXB feedback and stellar winds on our results, finding the ubiquity of stellar winds throughout each starburst makes them a defining factor in the final state of the interstellar medium.

The atmosphere of Titan, the largest moon of Saturn, is rich in organic molecules, and it has been suggested that the moon may serve as an analog for the pre-biotic Earth due to its highly reducing chemistry and existence of global hazes. Photochemical models of Titan have predicted the presence of propadiene (historically referred to as allene), CH$_{2}$CCH$_{2}$, an isomer of the well-measured propyne (also called methylacetylene) CH$_{3}$CCH, but its detection has remained elusive due to insufficient spectroscopic knowledge of the molecule - which has recently been remedied with an updated spectral line list. Here we present the first unambiguous detection of the molecule in any astronomical object, observed with the Texas Echelle Cross Echelle Spectrograph (TEXES) on the NASA Infrared Telescope Facility (IRTF) in July 2017. We model its emission line near 12 $\mu$m and measure a volume mixing ratio (VMR) of (6.9 $\pm$ 0.8) $\times$10$^{-10}$ at 175 km, assuming a vertically increasing abundance profile as predicted in photochemical models. Cassini measurements of propyne made during April 2017 indicate that the abundance ratio of propyne to propadiene is 8.2$\pm$1.1 at the same altitude. This initial measurement of the molecule in Titan's stratosphere paves the way towards constraining the amount of atomic hydrogen available on Titan, as well as future mapping of propadiene on Titan from 8 meter and larger ground based observatories, and future detection on other planetary bodies.

The Cassini Composite Infrared Spectrometer (CIRS) observed thermal emission in the far- and mid-infrared (from 10 cm$^{-1}$ to 1500 cm$^{-1}$), enabling spatiotemporal studies of ethane on Titan across the span of the Cassini mission from 2004 through 2017. Many previous measurements of ethane on Titan have relied on modeling the molecule's mid-infrared $\nu_{12}$ band, centered on 822 cm$^{-1}$. Other bands of ethane at shorter and longer wavelengths were seen, but have not been modeled to measure ethane abundance. Spectral line lists of the far-infrared $\nu_{4}$ torsional band at 289 cm$^{-1}$ and the mid-infrared $\nu_{8}$ band centered ay 1468 cm$^{-1}$ have recently been studied in the laboratory. We model CIRS observations of each of these bands (along with the $\nu_{12}$ band) separately and compare retrieved mixing ratios from each spectral region. Nadir observations of of the $\nu_{4}$ band probe the low stratosphere below 100 km. Our equatorial measurements at 289 cm$^{-1}$ show an abundance of (1.0$\pm$0.4) $\times$10$^{-5}$ at 88 km, from 2007 to 2017. This mixing ratio is consistent with measurements at higher altitudes, in contrast to the depletion that many photochemical models predict. Measurements from the $\nu_{12}$ and $\nu_{8}$ bands are comparable to each other, with the $\nu_{12}$ band probing an altitude range that extends deeper in the atmosphere. We suggest future studies of planetary atmospheres may observe the $\nu_{8}$ band, enabling shorter wavelength studies of ethane. There may also be an advantage to observing both the ethane $\nu_{8}$ band and nearby methane $\nu_{4}$ band in the same spectral window.

The tidal disruption of stars by supermassive black holes (SMBHs) can be used to probe the SMBH mass function, the properties of individual stars, and stellar dynamics in galactic nuclei. Upcoming missions will detect thousands of TDEs, and accurate theoretical modeling is required to interpret the data with precision. Here we analyze the influence of more realistic stellar structure on the outcome of TDEs; in particular, we compare the fallback rates -- being the rate at which tidally-disrupted debris returns to the black hole -- from progenitors generated with the stellar evolution code {\sc mesa} to $\gamma = 4/3$ and $\gamma = 5/3$ polytropes. We find that {\sc mesa}-generated density profiles yield qualitatively-different fallback rates as compared to polytropic approximations, and that only the fallback curves from low-mass ($1M_{\odot}$ or less), zero-age main-sequence stars are well fit by either a $\gamma = 4/3$ or $5/3$ polytrope. Stellar age has a strong affect on the shape of the fallback curve, and can produce characteristic timescales (e.g., the time to the peak of the fallback rate) that greatly differ from the polytropic values. We use these differences to assess the degree to which the inferred black hole mass from the observed lightcurve can deviate from the true value, and find that the discrepancy can be at the order of magnitude level. Accurate stellar structure also leads to a substantial variation in the critical impact parameter at which the star is fully disrupted, and can increase the susceptibility of the debris stream to fragmentation under its own self-gravity. These results suggest that detailed modeling is required to accurately interpret observed lightcurves of TDEs.

In this white paper, we present a cross-section of important scientific questions that remain partially or completely unanswered, ranging from Titan exosphere to the deep interior, and we detail which instrumentation and mission scenarios should be used to answer them. Our intention is to formulate the science goals for the next generation of planetary missions to Titan in order to prepare the future exploration of the moon. The ESA L-class mission concept that we propose is composed of a Titan orbiter and at least an in situ element (lake lander and/or drone(s)).

From 2004 to 2017, the Cassini spacecraft orbited Saturn, completing 127 close flybys of its largest moon, Titan. Cassini's Composite Infrared Spectrometer (CIRS), one of 12 instruments carried on board, profiled Titan in the thermal infrared (7-1000 microns) throughout the entire 13-year mission. CIRS observed on both targeted encounters (flybys) and more distant opportunities, collecting 8.4 million spectra from 837 individual Titan observations over 3633 hours. Observations of multiple types were made throughout the mission, building up a vast mosaic picture of Titan's atmospheric state across spatial and temporal domains. This paper provides a guide to these observations, describing each type and chronicling its occurrences and global-seasonal coverage. The purpose is to provide a resource for future users of the CIRS data set, as well as those seeking to put existing CIRS publications into the overall context of the mission, and to facilitate future inter-comparison of CIRS results with those of other Cassini instruments, and ground-based observations.

In a standard, steady, thin accretion disc, the radial distribution of the dissipation of the accretion energy is determined simply by energy considerations. Here we draw attention to the fact that while the (quasi-)steady discs in dwarf novae in outburst are in agreement with the expected emission distribution, the steady discs in the nova-like variables are not. We note that essentially the only difference between these two sets of discs is the time for which they have been in the high viscosity, high accretion rate state. In such discs, the major process by which angular momentum is transported outwards is MHD turbulence. We speculate that such turbulence gives rise to corona-like structures (here called magnetically controlled zones, or MCZs) which are also able to provide non-negligible angular momentum transport, the magnitude of which depends on the spatial scale $L$ of the magnetic field structures in such zones. For short-lived, high accretion rate discs (such as those in dwarf novae) we expect $L \sim H$ and the MCZ to have little effect. But, with time (such as in the nova-like variables) an inverse cascade in the MHD turbulence enables $L$, and the net effect of the MCZ, to grow. We present a simple toy model which demonstrates that such ideas can provide an explanation for the difference between the dwarf novae and the nova-like variable discs.

The Cassini mission performed 127 targeted flybys of Titan during its 13-year mission to Saturn, culminating in the Grand Finale between April-September 2017. Here we demonstrate the use of the Atacama Large Millimeter/submillimeter Array (ALMA) to continue Cassini's legacy for chemical and climatological studies of Titan's atmosphere. Whole-hemisphere, interferometric spectral maps of HCN, HNC, HC3N, CH3CN, C2H3CN, C2H5CN and C3H8 were obtained using ALMA in May 2017 at moderate (~0.2'', or 1300 km) spatial resolution, revealing the effects of seasonally-variable chemistry and dynamics on the distribution of each species. The ALMA sub-mm observations of HCN and HC3N are consistent with Cassini infrared data on these species, obtained in the same month. Chemical/dynamical lifetimes of a few years are inferred for C2H3CN and C2H5CN, in reasonably close agreement with the latest chemical models incorporating sticking of C2H5CN to stratospheric aerosol particles. ALMA radial limb flux profiles provide column density information as a function of altitude, revealing maximum abundances in the thermosphere (above 600 km) for HCN, HNC, HC3N and C2H5CN. This constitutes the first detailed measurement of the spatial distribution of HNC, which is found to be confined predominantly to altitudes above 730 $\pm$ 60 km. The HNC emission map shows an east-west hemispheric asymmetry of (13$\pm$3)%. These results are consistent with very rapid production (and loss) of HNC in Titan's uppermost atmosphere, making this molecule an effective probe of short-timescale (diurnal) ionospheric processes.

A star approaching a supermassive black hole (SMBH) can be torn apart in a tidal disruption event (TDE). We examine ultra-deep TDEs, a new regime in which the disrupted debris approaches close to the black hole's Schwarzschild radius, and the leading part intersects the trailing part at the first pericenter passage. We calculate the range of penetration factors $\beta$ vs SMBH masses $M$ that produce these prompt self-intersections using a Newtonian analytic estimate and a general relativistic (GR) geodesic model. We find that significant self-intersection of Solar-type stars requires $\beta \sim 50 - 127$ for $M/M_\odot = 10^4$, down to $\beta \sim 5.6 - 5.9$ for $M/M_\odot = 10^6$. We run smoothed-particle hydrodynamic (SPH) simulations to corroborate our calculations and find close agreement, with a slightly shallower dependence on $M$. We predict that the shock from the collision emits an X-ray flare lasting $t \sim 2$ s with $L \sim 10^{47}$ ergs/s at $E \sim 2$ keV, and the debris has a prompt accretion episode lasting $t \sim$ several min. The events are rare and occur with a rate $\dot{N} \lesssim 10^{-7}$ Mpc$^{-3}$ yr$^{-1}$. Ultra-deep TDEs can probe the strong gravity and demographics of low-mass SMBHs.

We present the first radio/submillimeter detection of monodeuterated methane (CH$_3$D) in Titan's atmosphere, using archival data from of the Atacama Large Millimeter/submillimeter Array (ALMA). The $J_K=2_1-1_1$ and $J_K=2_0-1_0$ transitions at 465.235 and 465.250 GHz ($\sim0.644$ mm) were measured at significance levels of $4.6\sigma$ and $5.7\sigma$, respectively. These two lines were modeled using the Non-linear optimal Estimator for MultivariatE spectral analySIS (NEMESIS) radiative transfer code to determine the disk-averaged CH$_3$D volume mixing ratio = $6.157\times10^{-6}$ in Titan's stratosphere (at altitudes $\gt130$ km). By comparison with the CH$_4$ vertical abundance profile measured by Cassini-Huygens mass spectrometry, the resulting value for D/H in CH$_4$ is $(1.033\pm0.081)\times10^{-4}$. This is consistent with previous ground-based and in-situ measurements from the Cassini-Huygens mission, though slightly lower than the average of the previous values. Additional CH$_3$D observations at higher spatial resolution will be required to determine a value truly comparable with the Cassini-Huygens CH$_4$ measurements, by measuring CH$_3$D with ALMA close to Titan's equator. In the post-Cassini era, spatially resolved observations of CH$_3$D with ALMA will enable the latitudinal distribution of methane to be determined, making this an important molecule for further studies.

Several of NASA missions (TESS, JWST, WFIRST, etc.) and mission concepts (LUVOIR, HabEx, and OST) emphasize the exploration and characterization of exoplanets, and the study of the interstellar medium. We anticipate that a much broader set of chemical environments exists on exoplanets, necessitating data from a correspondingly broader set of chemical reactions. Similarly, the conditions that exist in astrophysical environments are very different from those traditionally probed in laboratory chemical kinetics studies. These are areas where quantum mechanical theory, applied to important reactions via well-validated chemical kinetics models, can fill a critical knowledge gap. Quantum chemical calculations are also introduced to study interior of planets, photochemical escape, and many critical chemical pathways (e.g. prebiotic environments, contaminations, etc.) After years of development of the relevant quantum chemical theories and significant advances in computational power, quantum chemical simulations have currently matured enough to describe real systems with an accuracy that competes with experiments. These approaches, therefore, have become the best possible alternative in many circumstances where performing experiments is too difficult, too expensive, or too dangerous, or simply not possible. In this white paper, several existing quantum chemical studies supporting exoplanetary science, planetary astronomy, and astrophysics are described, and the potential positive impacts of improved models associated with scientific goals of missions are addressed. In the end, a few recommendations from the scientific community to strengthen related research efforts at NASA are provided.

Unexpected dynamic phenomena have surprised solar system observers in the past and have led to important discoveries about solar system workings. Observations at the initial stages of these events provide crucial information on the physical processes at work. We advocate for long-term/permanent programs on ground-based and space-based telescopes of all sizes - including Extremely Large Telescopes (ELTs) - to conduct observations of high-priority dynamic phenomena, based on a predefined set of triggering conditions. These programs will ensure that the best initial dataset of the triggering event are taken; separate additional observing programs will be required to study the temporal evolution of these phenomena. While not a comprehensive list, the following are notional examples of phenomena that are rare, that cannot be anticipated, and that provide high-impact advances to our understandings of planetary processes. Examples include: new cryovolcanic eruptions or plumes on ocean worlds; impacts on Jupiter, Saturn, Uranus, or Neptune; extreme eruptions on Io; convective superstorms on Saturn, Uranus, or Neptune; collisions within the asteroid belt or other small-body populations; discovery of an interstellar object passing through our solar system (e.g. 'Oumuamua); and responses of planetary atmospheres to major solar flares or coronal mass ejections.

Imaging and resolved spectroscopy reveal varying environmental conditions in our dynamic solar system. Many key advances have focused on how these conditions change over time. Observatory-level commitments to conduct annual observations of solar system bodies would establish a long-term legacy chronicling the evolution of dynamic planetary atmospheres, surfaces, and rings. Science investigations will use these temporal datasets to address potential biosignatures, circulation and evolution of atmospheres from the edge of the habitable zone to the ice giants, orbital dynamics and planetary seismology with ring systems, exchange between components in the planetary system, and the migration and processing of volatiles on icy bodies, including Ocean Worlds. The common factor among these diverse investigations is the need for a very long campaign duration, and temporal sampling at an annual cadence.

The manner in which warps in accretion disks evolve depends on the magnitude of the viscosity. For small viscosity $(\alpha < H/R)$ the warp evolves in a wave-like manner; for large viscosity $H/R<\alpha \ll 1$ it evolves diffusively. Here $\alpha$ is the viscosity parameter and $H/R$ the disk aspect ratio. Currently there is no simple set of equations which describes the evolution in both regimes. In this paper we describe a possible solution to this problem and introduce a set of one-dimensional equations that describe the evolution of a warped disk that are applicable in both high and low viscosity regimes for arbitrary tilts, but small warps.

The tidal force from a supermassive black hole can rip apart a star that passes close enough in what is known as a Tidal Disruption Event. Typically half of the destroyed star remains bound to the black hole and falls back on highly eccentric orbits, forming an accretion flow which powers a luminous flare. In this paper we use analytical and numerical calculations to explore the effect of stellar rotation on the fallback rate of material. We find that slowly spinning stars ($\Omega_* \lesssim 0.01 \Omega_{\rm{breakup}}$) provide only a small perturbation to fallback rates found in the non-spinning case. However when the star spins faster, there can be significant effects. If the star is spinning retrograde with respect to its orbit the tidal force from the black hole has to spin down the star first before disrupting it, causing delayed and sometimes only partial disruption events. However, if the star is spinning prograde this works with the tidal force and the material falls back sooner and with a higher peak rate. We examine the power-law index of the fallback curves, finding that in all cases the fallback rate overshoots the canonical $t^{-5/3}$ rate briefly after the peak, with the depth of the overshoot dependent on the stellar spin. We also find that in general the late time evolution is slightly flatter than the canonical $t^{-5/3}$ rate. We therefore conclude that considering the spin of the star may be important in modelling observed TDE lightcurves.

We use well-established observational evidence to draw conclusions about the fundamental nature of the viscosity in accretion discs. To do this, we first summarise the observational evidence for the value of the dimensionless accretion disc viscosity parameter $\alpha$, defined by Shakura & Sunyaev (1973, 1976). We find that, for fully ionized discs, the value of $\alpha$ is readily amenable to reliable estimation and that the observations are consistent with the hypothesis that $\alpha \sim 0.2 - 0.3$. In contrast in discs that are not fully ionized, estimates of the value of $\alpha$ are generally less direct and the values obtained are generally $ < 0.01$ and often $ \ll 0.01$. We conclude that this gives us crucial information about the nature of viscosity in accretion discs. First, in fully ionized discs the strength of the turbulence is always limited by being at most trans-sonic. This implies that it is necessary that credible models of the turbulence reflect this fact. Second, the smaller values of $\alpha$ found for less ionized, and therefore less strongly conducting, discs imply that magnetism plays a dominant role. This provides important observational support for the concept of magneto-rotational instability (MRI) driven hydromagnetic turbulence.

Previous investigations have employed more than 100 close observations of Titan by the Cassini orbiter to elucidate connections between the production and distribution of Titan's vast, organic-rich chemical inventory and its atmospheric dynamics. However, as Titan transitions into northern summer, the lack of incoming data from the Cassini orbiter presents a potential barrier to the continued study of seasonal changes in Titan's atmosphere. In our previous work (Thelen et al., 2018), we demonstrated that the Atacama Large Millimeter/submillimeter Array (ALMA) is well suited for measurements of Titan's atmosphere in the stratosphere and lower mesosphere (~100-500 km) through the use of spatially resolved (beam sizes <1'') flux calibration observations of Titan. Here, we derive vertical abundance profiles of four of Titan's trace atmospheric species from the same 3 independent spatial regions across Titan's disk during the same epoch (2012 to 2015): HCN, HC$_3$N, C$_3$H$_4$, and CH$_3$CN. We find that Titan's minor constituents exhibit large latitudinal variations, with enhanced abundances at high latitudes compared to equatorial measurements; this includes CH$_3$CN, which eluded previous detection by Cassini in the stratosphere, and thus spatially resolved abundance measurements were unattainable. Even over the short 3-year period, vertical profiles and integrated emission maps of these molecules allow us to observe temporal changes in Titan's atmospheric circulation during northern spring. Our derived abundance profiles are comparable to contemporary measurements from Cassini infrared observations, and we find additional evidence for subsidence of enriched air onto Titan's south pole during this time period. Continued observations of Titan with ALMA beyond the summer solstice will enable further study of how Titan's atmospheric composition and dynamics respond to seasonal changes.

Submillimeter emission lines of carbon monoxide (CO) in Titan's atmosphere provide excellent probes of atmospheric temperature due to the molecule's long chemical lifetime and stable, well constrained volume mixing ratio. Here we present the analysis of 4 datasets obtained with the Atacama Large Millimeter/Submillimeter Array (ALMA) from 2012 to 2015 that contain strong CO rotational transitions. Utilizing ALMA's high spatial resolution in the 2012, 2014, and 2015 observations, we extract spectra from 3 separate regions on Titan's disk using datasets with beam sizes of ~0.3''. Temperature profiles retrieved by the NEMESIS radiative transfer code are compared to Cassini Composite Infrared Spectrometer (CIRS) and radio occultation science results from similar latitude regions. Small seasonal variations in atmospheric temperature are present from 2012 to 2015 in the stratosphere and mesosphere (~100-500 km) of spatially resolved regions. We measure the stratopause (320 km) to increase in temperature by 5 K in northern latitudes from 2012-2015, while temperatures rise throughout the stratosphere at lower latitudes. While retrieved temperature profiles cover a range of latitudes in these observations, deviations from CIRS nadir maps and radio occultation measurements convolved with the ALMA beam-footprint are not found to be statistically significant, and discrepancies are often found to be less than 5 K throughout the atmosphere. ALMA's excellent sensitivity in the lower stratosphere (60-300 km) provides a highly complementary dataset to contemporary CIRS and radio science observations. The demonstrated utility of CO emission lines in the submillimeter as a tracer of Titan's atmospheric temperature lays the groundwork for future studies of other molecular species, as temperature profiles are found to consistently vary with latitude in all three years by up to 15 K.

Magnetic fields are observed in star forming regions. However simulations of the late stages of star formation that do not include magnetic fields provide a good fit to the properties of young stars including the initial mass function (IMF) and the multiplicity. We argue here that the simulations that do include magnetic fields are unable to capture the correct physics, in particular the high value of the magnetic Prandtl number, and the low value of the magnetic diffusivity. The artificially high (numerical and uncontrolled) magnetic diffusivity leads to a large magnetic flux pervading the star forming region. We argue further that in reality the dynamics of high magnetic Prandtl number turbulence may lead to local regions of magnetic energy dissipation through reconnection, meaning that the regions of molecular clouds which are forming stars might be essentially free of magnetic fields. Thus the simulations that ignore magnetic fields on the scales on which the properties of stellar masses, stellar multiplicities and planet-forming discs are determined, may be closer to reality than those which include magnetic fields, but can only do so in an unrealistic parameter regime.

Of the C$_{3}$H$_{x}$ hydrocarbons, propane (C${_3}$H$_{8}$) and propyne (methylacetylene, CH$_{3}$C$_{2}$H) were first detected in Titan's atmosphere during the Voyager 1 flyby in 1980. Propene (propylene, C$_{3}$H$_{6}$) was first detected in 2013 with data from the Composite InfraRed Spectrometer (CIRS) instrument on Cassini. We present the first measured abundance profiles of propene on Titan from radiative transfer modeling, and compare our measurements to predictions derived from several photochemical models. Near the equator, propene is observed to have a peak abundance of 10 ppbv at a pressure of 0.2 mbar. Several photochemical models predict the amount at this pressure to be in the range 0.3 - 1 ppbv and also show a local minimum near 0.2 mbar which we do not see in our measurements. We also see that propene follows a different latitudinal trend than the other C$_{3}$ molecules. While propane and propyne concentrate near the winter pole, transported via a global convective cell, propene is most abundant above the equator. We retrieve vertical abundances profiles between 125 km and 375 km for these gases for latitude averages between 60$^{\circ}$S to 20$^{\circ}$S, 20$^{\circ}$S to 20$^{\circ}$N, and 20$^{\circ}$N to 60$^{\circ}$N over two time periods, 2004 through 2009 representing Titan's atmosphere before the 2009 equinox, and 2012 through 2015 representing time after the equinox. Additionally, using newly corrected line data, we determined an updated upper limit for allene (propadiene, CH$_{2}$CCH$_{2}$, the isomer of propyne). We claim a 3-$\sigma$ upper limit mixing ratio of 2.5$\times$10$^{-9}$ within 30$^\circ$ of the equator. The measurements we present will further constrain photochemical models by refining reaction rates and the transport of these gases throughout Titan's atmosphere.

Blueshifted absorption lines in the X-ray spectra of AGN show that ultra-fast outflows with typical velocities $v \sim 0.1c$ are a common feature of these luminous objects. Such powerful AGN winds offer an explanation of the observed M-$\sigma$ relation linking the mass of the supermassive black hole and the velocity dispersion in the galaxy's stellar bulge. An extended XMM-Newton study of the luminous Seyfert galaxy PG1211+143 recently revealed a variable multi-velocity wind. Here we report the detection of a short-lived, ultra-fast inflow during the same observation. Previous reports of inflows used single absorption lines with uncertain identifications, but this new result identifies an array of resonance absorption lines of highly ionised Fe, Ca, Ar, S and Si, sharing a common redshift when compared with a grid of realistic photoionization spectra. The redshifted absorption arises in a column of highly ionized matter close to the black hole, with a line-of-sight velocity, $v \sim 0.3c$, inconsistent with the standard picture of a plane circular accretion disc. This may represent the first direct evidence for chaotic accretion in AGN, where accretion discs are generally misaligned to the black hole spin. For sufficient inclinations, the Lense-Thirring effect can break the discs into discrete rings, which then precess, collide and shock, causing near free-fall of gas towards the black hole. The observed accretion rate for the reported infall is comparable to the hard X-ray luminosity in PG1211+143, suggesting that direct infall may be a significant contributor to inner disc accretion.

A circumbinary disc around a pair of merging stellar-mass black holes may be shocked and heated during the recoil of the merged hole, causing a near-simultaneous electromagnetic counterpart to the gravitational wave event. The shocks occur around the recoil radius, where the disc orbital velocity is equal to the recoil velocity. The amount of mass present near this radius at the time of the merger is critical in determining how much radiation is released. We explore the evolution of a circumbinary disc in two limits. First, we consider an accretion disc that feels no torque from the binary. The disc does not survive until the merger unless there is a dead zone, a region of low turbulence. Even with the dead zone, the surface density in this case may be small. Second, we consider a disc that feels a strong binary torque that prevents accretion on to the binary. In this case there is significantly more mass in regions of interest at the time of the merger. A dead zone in this disc increases the mass close to the recoil radius. For typical binary-disc parameters we expect accretion to be significantly slowed by the resonant torque from the binary, and for a dead zone to be present. We conclude that provided significant mass orbits the binary after the formation of the black hole binary and that the radiation produced in recoil shocks can escape the flow efficiently, there is likely to be an observable electromagnetic signal from black hole binary mergers.

Recent observations of several protoplanetary discs have found evidence of departures from flat, circular motion in the inner regions of the disc. One possible explanation for these observations is a disc warp, which could be induced by a planet on a misaligned orbit. We present three-dimensional numerical simulations of the tidal interaction between a protoplanetary disc and a misaligned planet. For low planet masses we show that our simulations accurately model the evolution of inclined planet orbit (up to moderate inclinations). For a planet massive enough to carve a gap, the disc is separated into two components and the gas interior and exterior to the planet orbit evolve separately, forming an inner and outer disc. Due to the inclination of the planet, a warp develops across the planet orbit such that there is a relative tilt and twist between these discs. We show that when other parameters are held constant, the relative inclination that develops between the inner and outer disc depends on the outer radius of the total disc modelled. For a given disc mass, our results suggest that the observational relevance of the warp depends more strongly on the mass of the planet rather than the inclination of the orbit.

We present a community-led assessment of the solar system investigations achievable with NASA's next-generation space telescope, the Wide Field InfraRed Survey Telescope (WFIRST). WFIRST will provide imaging, spectroscopic, and coronagraphic capabilities from 0.43-2.0 $\mu$m and will be a potential contemporary and eventual successor to JWST. Surveys of irregular satellites and minor bodies are where WFIRST will excel with its 0.28 deg$^2$ field of view Wide Field Instrument (WFI). Potential ground-breaking discoveries from WFIRST could include detection of the first minor bodies orbiting in the Inner Oort Cloud, identification of additional Earth Trojan asteroids, and the discovery and characterization of asteroid binary systems similar to Ida/Dactyl. Additional investigations into asteroids, giant planet satellites, Trojan asteroids, Centaurs, Kuiper Belt Objects, and comets are presented. Previous use of astrophysics assets for solar system science and synergies between WFIRST, LSST, JWST, and the proposed NEOCam mission are discussed. We also present the case for implementation of moving target tracking, a feature that will benefit from the heritage of JWST and enable a broader range of solar system observations.

We investigate the effects of gradual heating on the evolution of turbulent molecular clouds of mass $2\times 10^6$ M$_\odot$ and virial parameters ranging between $0.7-1.2$. This gradual heating represents the energy output from processes such as winds from massive stars or feedback from High Mass X-ray binaries (HMXBs), contrasting the impulsive energy injection from supernovae (SNe). For stars with a mass high enough that their lifetime is shorter than the life of the cloud, we include a SN feedback prescription. Including both effects, we investigate the interplay between slow and fast forms of feedback and their effectiveness at triggering/suppressing star formation. We find that SN feedback can carve low density chimneys in the gas, offering a path of least resistance for the energy to escape. Once this occurs the more stable, but less energetic, gradual feedback is able to keep the chimneys open. By funneling the hot destructive gas away from the centre of the cloud, chimneys can have a positive effect on both the efficiency and duration of star formation. Moreover, the critical factor is the number of high mass stars and SNe (and any subsequent HMXBs) active within the free-fall time of each cloud. This can vary from cloud to cloud due to the stochasticity of SN delay times and in HMXB formation. However, the defining factor in our simulations is the efficiency of the cooling, which can alter the Jeans mass required for sink particle formation, along with the number of massive stars in the cloud.

After the tidal disruption of a star by a massive black hole, disrupted stellar debris can fall back to the hole at a rate significantly exceeding its Eddington limit. To understand how black hole mass affects the duration of super-Eddington accretion in tidal disruption events, we first run a suite of simulations of the disruption of a Solar-like star by a supermassive black hole of varying mass to directly measure the fallback rate onto the hole, and we compare these fallback rates to the analytic predictions of the "frozen-in" model. Then, adopting a Zero-Bernoulli Accretion flow as an analytic prescription for the accretion flow around the hole, we investigate how the accretion rate onto the black hole evolves with the more accurate fallback rates calculated from the simulations. We find that numerically-simulated fallback rates yield accretion rates onto the hole that can, depending on the black hole mass, be nearly an order of magnitude larger than those predicted by the frozen-in approximation. Our results place new limits on the maximum black hole mass for which super-Eddington accretion occurs in tidal disruption events.

We investigate the formation of discs within binary systems where at least one component has left the main sequence. In particular we calculate the occurrence rates of systems which can host long-lived, massive discs that may be able to support the formation of planets. We synthesize a population based on Milky Way properties, using both theoretical and observational inputs to constrain key properties such as the shape of the initial mass function, binary fraction, and mass transfer physics. We predict 0.26\% of binary systems will host Second generation discs (where the primary has evolved), and 0.13\% of systems will host Third generation discs (where the secondary also evolves). For the Milky Way, this translates into 130 million and 90 million Second and Third generation systems respectively from an estimated total of 50 billion binary systems. Of these systems that form discs, we estimate approximately 20\% of Second and 3.8\% of Third generation discs have enough mass to form a planetary system. We speculate on how the process of planet formation differs in these systems from conventional planet formation in protostellar discs.

We present Phantom, a fast, parallel, modular and low-memory smoothed particle hydrodynamics and magnetohydrodynamics code developed over the last decade for astrophysical applications in three dimensions. The code has been developed with a focus on stellar, galactic, planetary and high energy astrophysics and has already been used widely for studies of accretion discs and turbulence, from the birth of planets to how black holes accrete. Here we describe and test the core algorithms as well as modules for magnetohydrodynamics, self-gravity, sink particles, H_2 chemistry, dust-gas mixtures, physical viscosity, external forces including numerous galactic potentials as well as implementations of Lense-Thirring precession, Poynting-Robertson drag and stochastic turbulent driving. Phantom is hereby made publicly available.

We present the first maps of cyanoacetylene isotopologues in Titan's atmosphere, including H$^{13}$CCCN and HCCC$^{15}$N, detected in the 0.9 mm band using the Atacama Large Millimeter/submillimeter array (ALMA) around the time of Titan's (southern winter) solstice in May 2017. The first high-resolution map of HC$_3$N in its $v_7=1$ vibrationally excited state is also presented, revealing a unique snapshot of the global HC$_3$N distribution, free from the strong optical depth effects that adversely impact the ground-state ($v=0$) map. The HC$_3$N emission is found to be strongly enhanced over Titan's south pole (by a factor of 5.7 compared to the north pole), consistent with rapid photochemical loss of HC$_3$N from the summer hemisphere combined with production and transport to the winter pole since the April 2015 ALMA observations. The H$^{13}$CCCN/HCCC$^{15}$N flux ratio is derived at the southern HC$_3$N peak, and implies an HC$_3$N/HCCC$^{15}$N ratio of $67\pm14$. This represents a significant enrichment in $^{15}$N compared with Titan's main molecular nitrogen reservoir, which has a $^{14}$N/$^{15}$N ratio of 167, and confirms the importance of photochemistry in determining the nitrogen isotopic ratio in Titan's organic inventory.

Disc accretion on to a black hole is often misaligned from its spin axis. If the disc maintains a significant magnetic field normal to its local plane, we show that dipole radiation from Lense-Thirring precessing disc annuli can extract a significant fraction of the accretion energy, sharply peaked towards small disc radii $R$ (as $R^{-17/2}$ for fields with constant equipartition ratio). This low-frequency emission is immediately absorbed by surrounding matter or refracted towards the regions of lowest density. The resultant mechanical pressure, dipole angular pattern, and much lower matter density towards the rotational poles create a strong tendency to drive jets along the black hole spin axis, similar to the spin-axis jets of radio pulsars, also strong dipole emitters. The coherent primary emission may explain the high brightness temperatures seen in jets. The intrinsic disc emission is modulated at Lense-Thirring frequencies near the inner edge, providing a physical mechanism for low-frequency QPOs. Dipole emission requires nonzero hole spin, but uses only disc accretion energy. No spin energy is extracted, unlike the Blandford-Znajek process. MHD/GRMHD formulations do not directly give radiation fields, but can be checked post-process for dipole emission and so self-consistency, given sufficient resolution. Jets driven by dipole radiation should be more common in AGN than in X-ray binaries, and in low accretion rate states than high, agreeing with observation. In non-black-hole accretion, misaligned disc annuli precess because of the accretor's mass quadrupole moment, similarly producing jets and QPOs.

Current planet formation theories provide successful frameworks with which to interpret the array of new observational data in this field. However, each of the two main theories (core accretion, gravitational instability) is unable to explain some key aspects. In many planet formation calculations, it is usual to treat the initial properties of the planet forming disc (mass, radius, etc.) as free parameters. In this paper, we stress the importance of setting the formation of planet forming discs within the context of the formation of the central stars. By exploring the early stages of disc formation, we introduce the concept of the Maximum Mass Solar Nebula (MMSN), as opposed to the oft-used Minimum Mass Solar Nebula (here mmsn). It is evident that almost all protoplanetary discs start their evolution in a strongly self-gravitating state. In agreement with almost all previous work in this area, we conclude that on the scales relevant to planet formation these discs are not gravitationally unstable to gas fragmentation, but instead form strong, transient spiral arms. These spiral arms can act as efficient dust traps allowing the accumulation and subsequent fragmentation of the dust (but not the gas). This phase is likely to populate the disc with relatively large planetesimals on short timescales while the disc is still veiled by a dusty-gas envelope. Crucially, the early formation of large planetesimals overcomes the main barriers remaining within the core accretion model. A prediction of this picture is that essentially all observable protoplanetary discs are already planet hosting.

This whitepaper discusses the diversity of exoplanets that could be detected by future observations, so that comparative exoplanetology can be performed in the upcoming era of large space-based flagship missions. The primary focus will be on characterizing Earth-like worlds around Sun-like stars. However, we will also be able to characterize companion planets in the system simultaneously. This will not only provide a contextual picture with regards to our Solar system, but also presents a unique opportunity to observe size dependent planetary atmospheres at different orbital distances. We propose a preliminary scheme based on chemical behavior of gases and condensates in a planet's atmosphere that classifies them with respect to planetary radius and incident stellar flux.

The large-scale magnetic fields observed in spiral disc galaxies are often thought to result from dynamo action in the disc plane. However, the increasing importance of Faraday depolarization along any line of sight towards the galactic plane suggests that the strongest polarization signal may come from well above (~0.3-1 kpc) this plane, from the vicinity of the warm interstellar medium (WIM)/halo interface. We propose (see also Henriksen & Irwin 2016) that the observed spiral fields (polarization patterns) result from the action of vertical shear on an initially poloidal field. We show that this simple model accounts for the main observed properties of large-scale fields. We speculate as to how current models of optical spiral structure may generate the observed arm/interarm spiral polarization patterns.

In some galaxies, the stars orbiting the supermassive black hole take the form of an eccentric nuclear disk, in which every star is on a coherent, apsidally-aligned orbit. The most famous example of an eccentric nuclear disk is the double nucleus of Andromeda, and there is strong evidence for many more in the local universe. Despite their apparent ubiquity however, a dynamical explanation for their longevity has remained a mystery: differential precession should wipe out large-scale apsidal-alignment on a short timescale. Here we identify a new dynamical mechanism which stabilizes eccentric nuclear disks, and explain for first time the negative eccentricity gradient seen in the Andromeda nucleus. The stabilizing mechanism drives oscillations of the eccentricity vectors of individual orbits, both in direction (about the mean body of the disk) and in magnitude. Combined with the negative eccentricity gradient, the eccentricity oscillations push some stars near the inner edge of the disk extremely close to the black hole, potentially leading to tidal disruption events. Order of magnitude calculations predict extremely high rates in recently-formed eccentric nuclear disks ($\sim0.1 - 1$ ${\rm yr}^{-1} {\rm gal}^{-1}$). Unless the stellar disks are replenished, these rates should decrease with time as the disk depletes in mass. If eccentric nuclear disks form during gas-rich major mergers, this may explain the preferential occurrence of tidal disruption events in recently-merged and post-merger (E+A/K+A) galaxies.

Accretion discs are generally warped. If a warp in a disc is too large, the disc can `break' apart into two or more distinct planes, with only tenuous connections between them. Further if an initially planar disc is subject to a strong differential precession, then it can be torn apart into discrete annuli that precess effectively independently. In previous investigations, torque-balance formulae have been used to predict where and when the disc breaks into distinct parts. In this work, focusing on discs with Keplerian rotation and where the shearing motions driving the radial communication of the warp are damped locally by turbulence (the `diffusive' regime), we investigate the stability of warped discs to determine the precise criterion for an isolated warped disc to break. We find and solve the dispersion relation, which in general yields three roots. We provide a comprehensive analysis of this viscous-warp instability and the emergent growth rates and their dependence on disc parameters. The physics of the instability can be understood as a combination of (1) a term which would generally encapsulate the classical Lightman-Eardley instability in planar discs (given by $\partial(\nu\Sigma)/\partial\Sigma < 0$) but is here modified by the warp to include $\partial(\nu_1|\psi|)/\partial|\psi| < 0$ and (2) a similar condition acting on the diffusion of the warp amplitude given in simplified form by $\partial(\nu_2|\psi|)/\partial|\psi| < 0$. We discuss our findings in the context of discs with an imposed precession, and comment on the implications for different astrophysical systems.

Vinyl cyanide (C$_2$H$_3$CN) is theorized to form in Titan's atmosphere via high-altitude photochemistry and is of interest regarding the astrobiology of cold planetary surfaces due to its predicted ability to form cell membrane-like structures (azotosomes) in liquid methane. In this work, we follow up on the initial spectroscopic detection of C$_2$H$_3$CN on Titan by Palmer et al. (2017) with the detection of three new C$_2$H$_3$CN rotational emission lines at submillimeter frequencies. These new, high-resolution detections have allowed for the first spatial distribution mapping of C$_2$H$_3$CN on Titan. We present simultaneous observations of C$_2$H$_5$CN, HC$_3$N, and CH$_3$CN emission, and obtain the first (tentative) detection of C$_3$H$_8$ (propane) at radio wavelengths. We present disk-averaged vertical abundance profiles, two-dimensional spatial maps, and latitudinal flux profiles for the observed nitriles. Similarly to HC$_3$N and C$_2$H$_5$CN, which are theorized to be short-lived in Titan's atmosphere, C$_2$H$_3$CN is most abundant over the southern (winter) pole, whereas the longer-lived CH$_3$CN is more concentrated in the north. This abundance pattern is consistent with the combined effects of high-altitude photochemical production, poleward advection, and the subsequent reversal of Titan's atmospheric circulation system following the recent transition from northern to southern winter. We confirm that C$_2$H$_3$CN and C$_2$H$_5$CN are most abundant at altitudes above 200 km. Using a 300 km step model, the average abundance of C$_2$H$_3$CN is found to be $3.03\pm0.29$ ppb, with a C$_2$H$_5$CN/C$_2$H$_3$CN abundance ratio of $2.43\pm0.26$. Our HC$_3$N and CH$_3$CN spectra can be accurately modeled using abundance gradients above the tropopause, with fractional scale-heights of $2.05\pm0.16$ and $1.63\pm0.02$, respectively.

The ice giants Uranus and Neptune are the least understood class of planets in our solar system but the most frequently observed type of exoplanets. Presumed to have a small rocky core, a deep interior comprising ~70% heavy elements surrounded by a more dilute outer envelope of H2 and He, Uranus and Neptune are fundamentally different from the better-explored gas giants Jupiter and Saturn. Because of the lack of dedicated exploration missions, our knowledge of the composition and atmospheric processes of these distant worlds is primarily derived from remote sensing from Earth-based observatories and space telescopes. As a result, Uranus's and Neptune's physical and atmospheric properties remain poorly constrained and their roles in the evolution of the Solar System not well understood. Exploration of an ice giant system is therefore a high-priority science objective as these systems (including the magnetosphere, satellites, rings, atmosphere, and interior) challenge our understanding of planetary formation and evolution. Here we describe the main scientific goals to be addressed by a future in situ exploration of an ice giant. An atmospheric entry probe targeting the 10-bar level, about 5 scale heights beneath the tropopause, would yield insight into two broad themes: i) the formation history of the ice giants and, in a broader extent, that of the Solar System, and ii) the processes at play in planetary atmospheres. The probe would descend under parachute to measure composition, structure, and dynamics, with data returned to Earth using a Carrier Relay Spacecraft as a relay station. In addition, possible mission concepts and partnerships are presented, and a strawman ice-giant probe payload is described. An ice-giant atmospheric probe could represent a significant ESA contribution to a future NASA ice-giant flagship mission.

15 years of XMM-Newton observations have established that ultra-fast, highly ionized winds are common in radio-quiet AGN. A simple theory of Eddington-limited accretion correctly predicts the typical velocity (~0.1c) and high ionization of such winds, with observed flow energy capable of ejecting star-forming gas. With a recent extended XMM-Newton observation of the archetypal UFO, PG1211+143, revealing a more complex flow pattern, we suggest that targetted observations over the next decade offer unique potential for probing the inner accretion disc structure and SMBH growth.

We investigate the pre-disruption gravitational dynamics and post-disruption hydrodynamics of the tidal disruption of stars by supermassive black hole (SMBH) binaries. We focus on binaries with relatively low mass primaries ($10^6M_{\odot}$), moderate mass ratios, and separations with reasonably long gravitational wave inspiral times (tens of Myr). First, we generate a large ensemble (between 1 and 10 million) of restricted three-body integrations to quantify the statistical properties of tidal disruptions by circular SMBH binaries of initially-unbound stars. Compared to the reference case of a disruption by a single SMBH, the binary potential induces significant variance into the specific energy and angular momentum of the star at the point of disruption. Second, we use Newtonian numerical hydrodynamics to study the detailed evolution of the fallback debris from 120 disruptions randomly selected from the three-body ensemble (excluding only the most deeply penetrating encounters). We find that the overall morphology of the debris is greatly altered by the presence of the second black hole, and the accretion rate histories display a wide range of behaviors, including order of magnitude dips and excesses relative to control simulations that include only one black hole. Complex evolution persists, in some cases, for many orbital periods of the binary. We find evidence for power in the accretion curves on timescales related to the binary orbital period, though there is no exact periodicity. We discuss our results in the context of future wide-field surveys, and comment on the prospects of identifying and characterizing the subset of events occurring in nuclei with binary SMBHs.

We investigate systematically whether accreting black hole systems are likely to reach global alignment of the black hole spin and its accretion disc with the binary plane. In low-mass X-ray binaries (LMXBs) there is only a modest tendency to reach such global alignment, and it is difficult to achieve fully: except for special initial conditions we expect misalignment of the spin and orbital planes by ~1 radian for most of the LMXB lifetime. The same is expected in high-mass X-ray binaries (HMXBs). A fairly close approach to global alignment is likely in most stellar-mass ultraluminous X-ray binary systems (ULXs) where the companion star fills its Roche lobe and transfers on a thermal timescale to a black hole of lower mass. These systems are unlikely to show orbital eclipses, as their emission cones are close to the hole's spin axis. This offers a potential observational test, as models for ULXs invoking intermediate-mass black holes do predict eclipses for ensembles of > ~10 systems. Recent observational work shows that eclipses are either absent or extremely rare in ULXs, supporting the picture that most ULXs are stellar-mass binaries with companion stars more massive than the accretor.

We present sub-millimeter spectra of HCN isotopologues on Titan, derived from publicly available ALMA flux calibration observations of Titan taken in early 2014. We report the detection of a new HCN isotopologue on Titan, H$^{13}$C$^{15}$, and confirm an earlier report of detection of DCN. We model high signal-to-noise observations of HCN, H$^{13}$CN, HC$^{15}$N, DCN, and H$^{13}$C$^{15}$ to derive abundances and infer the following isotopic ratios: $^{12}$C/$^{13}$C = 89.8 $\pm$ 2.8, $^{14}$N/$^{15}$N = 72.3 $\pm$ 2.2, D/H = (2.5 $\pm$ 0.2)$\times$10$^{-4}$, and HCN/H$^{13}$C$^{15}$ = 5800 $\pm$ 270 (1$\sigma$ errors). The carbon and nitrogen ratios are consistent with and improve on the precision of previous results, confirming a factor of $\sim$2.3 elevation in $^{14}$N/$^{15}$N in HCN compared to N$_2$ and a lack of fractionation in $^{12}$C/$^{13}$C from the protosolar value. This is the first published measurement of D/H in a nitrile species on Titan, and we find evidence for a factor of $\sim$2 deuterium enrichment in hydrogen cyanide compared to methane. The isotopic ratios we derive may be used as constraints for future models to better understand the fractionation processes occurring in Titan's atmosphere.

With hydrodynamical simulations we determine the conditions under which an initially coplanar planet-disc system that orbits a member of a misaligned binary star evolves to form a planet that undergoes Kozai-Lidov (KL) oscillations once the disc disperses. These oscillations may explain the large orbital eccentricities, as well as the large misalignments with respect to the spin of the central star, observed for some exoplanets. The planet is assumed to be massive enough to open a gap in the disc. The planet's tilt relative to the binary orbital plane is subject to two types of oscillations. The first type, present at even small inclination angles relative to the binary orbital plane, is due to the interaction of the planet with the disc and binary companion and is amplified by a secular resonance. The second type of oscillation is the KL oscillation that operates on both the planet and disc at larger binary inclination angles. We find that for a sufficiently massive disc, even a relatively low inclination planet-disc system can force a planet to an inclination above the critical KL angle, as a consequence of the first type of tilt oscillation, allowing it to undergo the second type of oscillation. We conclude that the hydrodynamical evolution of a sufficiently massive and inclined disc in a binary system broadens the range of systems that form eccentric and misaligned giant planets to include a wide range of initial misalignment angles (20 to 160 degrees).

A tidal disruption event (TDE) -- when a star is destroyed by the immense gravitational field of a supermassive black hole -- transforms a star into a stream of tidally-shredded debris. The properties of this debris ultimately determine the observable signatures of TDEs. Here we derive a simple, self-similar solution for the velocity profile of the debris streams produced from TDEs, and show that this solution agrees extremely well with numerical results. Using this self-similar solution, we calculate an analytic, approximate expression for the radial density profile of the stream. We show that there is a critical adiabatic index that varies as a function of position along the stream above (below) which the stream is unstable (stable) to gravitational fragmentation. We also calculate the impact of heating and cooling on this stability criterion.

We report interferometric observations of carbon monoxide (CO) and its isotopologues in Titan's atmosphere using the Atacama Large Millimeter/submillimeter Array (ALMA). The following transitions were detected: CO (J = 1-0, 2-1, 3-2, 6-5), $^{13}$CO (J = 2-1, 3-2, 6-5), C$^{18}$O (J = 2-1, 3-2), and C$^{17}$O (J = 3-2). Molecular abundances and the vertical atmospheric temperature profile were derived by modeling the observed emission line profiles using NEMESIS, a line-by-line radiative transfer code. We present the first spectroscopic detection of $^{17}$O in the outer solar system with C$^{17}$O detected at > 8$\sigma$ confidence. The abundance of CO was determined to be 49.6 ${\pm}$ 1.8 ppm, assumed to be constant with altitude, with isotopic ratios $^{12}$C/$^{13}$C = 89.9 ${\pm}$ 3.4, $^{16}$O/$^{18}$O = 486 ${\pm}$ 22, and $^{16}$O/$^{17}$O = 2917${\pm}$359. The measurements of $^{12}$C/$^{13}$C and $^{16}$O/$^{18}$O ratios are the most precise values obtained in Titan's atmospheric CO to date. Our results are in good agreement with previous studies and suggest no significant deviations from standard terrestrial isotopic ratios.

The James Webb Space Telescope (JWST), scheduled for launch in 2018, is the successor to the Hubble Space Telescope (HST) but with a significantly larger aperture (6.5 m) and advanced instrumentation focusing on infrared science (0.6-28.0 $\mu$m ). In this paper we examine the potential for scientific investigation of Titan using JWST, primarily with three of the four instruments: NIRSpec, NIRCam and MIRI, noting that science with NIRISS will be complementary. Five core scientific themes are identified: (i) surface (ii) tropospheric clouds (iii) tropospheric gases (iv) stratospheric composition and (v) stratospheric hazes. We discuss each theme in depth, including the scientific purpose, capabilities and limitations of the instrument suite, and suggested observing schemes. We pay particular attention to saturation, which is a problem for all three instruments, but may be alleviated for NIRCam through use of selecting small sub-arrays of the detectors - sufficient to encompass Titan, but with significantly faster read-out times. We find that JWST has very significant potential for advancing Titan science, with a spectral resolution exceeding the Cassini instrument suite at near-infrared wavelengths, and a spatial resolution exceeding HST at the same wavelengths. In particular, JWST will be valuable for time-domain monitoring of Titan, given a five to ten year expected lifetime for the observatory, for example monitoring the seasonal appearance of clouds. JWST observations in the post-Cassini period will complement those of other large facilities such as HST, ALMA, SOFIA and next-generation ground-based telescopes (TMT, GMT, EELT).

A tidal disruption event, which occurs when a star is destroyed by the gravitational field of a supermassive black hole, produces a stream of debris, the evolution of which ultimately determines the observational properties of the event. Here we show that a post-periapsis caustic -- a location where the locus of gas parcels comprising the stream would collapse into a two-dimensional surface if they evolved solely in the gravitational field of the hole -- occurs when the pericenter distance of the star is on the order of the tidal radius of the hole. It is demonstrated that this "pancake" induces significant density perturbations in the debris stream, and, for stiffer equations of state (adiabatic index $\gamma \gtrsim 5/3$), these fluctuations are sufficient to gravitationally destabilize the stream, resulting in its fragmentation into bound clumps. The results of our findings are discussed in the context of the observational properties of tidal disruption events.

The Hera Saturn entry probe mission is proposed as an M--class mission led by ESA with a contribution from NASA. It consists of one atmospheric probe to be sent into the atmosphere of Saturn, and a Carrier-Relay spacecraft. In this concept, the Hera probe is composed of ESA and NASA elements, and the Carrier-Relay Spacecraft is delivered by ESA. The probe is powered by batteries, and the Carrier-Relay Spacecraft is powered by solar panels and batteries. We anticipate two major subsystems to be supplied by the United States, either by direct procurement by ESA or by contribution from NASA: the solar electric power system (including solar arrays and the power management and distribution system), and the probe entry system (including the thermal protection shield and aeroshell). Hera is designed to perform in situ measurements of the chemical and isotopic compositions as well as the dynamics of Saturn's atmosphere using a single probe, with the goal of improving our understanding of the origin, formation, and evolution of Saturn, the giant planets and their satellite systems, with extrapolation to extrasolar planets. Hera's aim is to probe well into the cloud-forming region of the troposphere, below the region accessible to remote sensing, to the locations where certain cosmogenically abundant species are expected to be well mixed. By leading to an improved understanding of the processes by which giant planets formed, including the composition and properties of the local solar nebula at the time and location of giant planet formation, Hera will extend the legacy of the Galileo and Cassini missions by further addressing the creation, formation, and chemical, dynamical, and thermal evolution of the giant planets, the entire solar system including Earth and the other terrestrial planets, and formation of other planetary systems.

We demonstrate the importance of general relativistic apsidal precession in warped black hole accretion discs by comparing three - dimensional smoothed particle hydrodynamic simulations in which this effect is first neglected, and then included. If apsidal precession is neglected, we confirm the results of an earlier magnetohydrodynamic simulation which made this assumption, showing that at least in this case the $\alpha$ viscosity model produces very similar results to those of simulations where angular momentum transport is due to the magnetorotational instability. Including apsidal precession significantly changes the predicted disc evolution. For moderately inclined discs thick enough that tilt is transported by bending waves, we find a disc tilt which is nonzero at the inner disc edge and oscillates with radius, consistent with published analytic results. For larger inclinations we find disc breaking.

Observational arguments suggest that the growth phases of the supermassive black holes in active galactic nuclei have a characteristic timescale $\sim 10^5$ yr. We show that this is the timescale expected in the chaotic accretion picture of black hole feeding, because of the effect of self-gravity in limiting the mass of any accretion disc feeding event.

We present simulations of the tidal disruption of a solar mass star by a $10^6M_{\odot}$ black hole. These, for the first time, cover the full time evolution of the tidal disruption event, starting well before the initial encounter and continuing until more than 90% of the bound material has returned to the vicinity of the hole. Our results are compared to the analytical prediction for the rate at which tidally-stripped gas falls back. We find that, for our chosen parameters, the overall scaling of the fallback rate, $\dot{M}_{\rm{fb}}$, closely follows the canonical $t^{-5/3}$ power-law. However, our simulations also show that the self-gravity of the tidal stream, which dominates the tidal gravity of the hole at large distances, causes some of the debris to recollapse into bound fragments before returning to the hole. This causes $\dot{M}_{\rm{fb}}$ to vary significantly around the $t^{-5/3}$ average. We discuss the implications of our findings in the context of the event Swift J1644+57.

Astrophysical discs are often warped, that is, their orbital planes change with radius. This occurs whenever there is a non-axisymmetric force acting on the disc, for example the Lense-Thirring precession induced by a misaligned spinning black hole, or the gravitational pull of a misaligned companion. Such misalignments appear to be generic in astrophysics. The wide range of systems that can harbour warped discs - protostars, X-ray binaries, tidal disruption events, quasars and others - allows for a rich variety in the disc's response. Here we review the basic physics of warped discs and its implications.

Observations of warped discs can give insight into the nature of angular momentum transport in accretion discs. Only a few objects are known to show strong periodicity on long timescales, but when such periodicity is present it is often attributed to precession of the accretion disc. The X-ray binary Hercules X-1/HZ Herculis (Her X-1) is one of the best examples of such periodicity and has been linked to disc precession since it was first observed. By using the current best-fitting models to Her X-1, which invoke precession driven by radiation warping, I place a constraint on the effective viscosities that act in a warped disc. These effective viscosities almost certainly arise due to turbulence induced by the magneto-rotational instability. The constraints derived here are in agreement with analytical and numerical investigations into the nature of magneto-hydrodynamic disc turbulence, but at odds with some recent global simulations.

We report the first spectroscopic detection of ethyl cyanide (C$_2$H$_5$CN) in Titan's atmosphere, obtained using spectrally and spatially resolved observations of multiple emission lines with the Atacama Large Millimeter/submillimeter array (ALMA). The presence of C$_2$H$_5$CN in Titan's ionosphere was previously inferred from Cassini ion mass spectrometry measurements of C$_2$H$_5$CNH$^+$. Here we report the detection of 27 rotational lines from C$_2$H$_5$CN (in 19 separate emission features detected at $>3\sigma$ confidence), in the frequency range 222-241 GHz. Simultaneous detections of multiple emission lines from HC$_3$N, CH$_3$CN and CH$_3$CCH were also obtained. In contrast to HC$_3$N, CH$_3$CN and CH$_3$CCH, which peak in Titan's northern (spring) hemisphere, the emission from C$_2$H$_5$CN is found to be concentrated in the southern (autumn) hemisphere, suggesting a distinctly different chemistry for this species, consistent with a relatively short chemical lifetime for C$_2$H$_5$CN. Radiative transfer models show that most of the C$_2$H$_5$CN is concentrated at altitudes 300-600 km, suggesting production predominantly in the mesosphere and above. Vertical column densities are found to be in the range (2-5)$\times10^{14}$ cm$^{-2}$.

Accretion discs are common in binary systems, and they are often found to be misaligned with respect to the binary orbit. The gravitational torque from a companion induces nodal precession in misaligned disc orbits. We calculate whether this precession is strong enough to overcome the internal disc torques communicating angular momentum. For typical parameters precession wins: the disc breaks into distinct planes that precess effectively independently. We run hydrodynamical simulations to check these results, and confirm that disc breaking is widespread and generally enhances accretion on to the central object. This applies in many cases of astrophysical accretion, e.g. supermassive black hole binaries and X--ray binaries.

We analyse the interaction of an eccentric binary with a circular coplanar circumbinary disc that rotates in a retrograde sense with respect to the binary. In the circular binary case, no Lindblad resonances lie within the disc and no Lindblad resonant torques are produced, as was previously known. By analytic means, we show that when the binary orbit is eccentric, there exist components of the gravitational potential of the binary which rotate in a retrograde sense to the binary orbit and so rotate progradely with respect to this disc, allowing a resonant interaction to occur between the binary and the disc. The resulting resonant torques distinctly alter the disc response from the circular binary case. We describe results of three-dimensional hydrodynamic simulations to explore this effect and categorise the response of the disc in terms of modes whose strengths vary as a function of binary mass ratio and eccentricity. These mode strengths are weak compared to the largest mode strengths expected in the prograde case where the binary and disc rotate in the same sense. However, for sufficiently high binary eccentricity, resonant torques open a gap in a retrograde circumbinary disc, while permitting gas inflow on to the binary via gas streams. The inflow results in a time varying accretion rate on to the binary that is modulated over the binary orbital period, as was previously found to occur in the prograde case.

We investigate the evolution of low mass (Md /Mb = 0.005) misaligned gaseous discs around eccentric supermassive black hole (SMBH) binaries. These are expected to form from randomly oriented accretion events onto a SMBH binary formed in a galaxy merger. When expanding the interaction terms between the binary and a circular ring to quadrupole order and averaging over the binary orbit, we expect four non-precessing disc orientations: aligned or counter-aligned with the binary, or polar orbits around the binary eccentricity vector with either sense of rotation. All other orientations precess around either of these, with the polar precession dominating for high eccentricity. These expectations are borne out by smoothed particle hydrodynamics simulations of initially misaligned viscous circumbinary discs, resulting in the formation of polar rings around highly eccentric binaries in contrast to the co-planar discs around circular binaries. Moreover, we observe disc tearing and violent interactions between differentially precessing rings in the disc significantly disrupting the disc structure and causing gas to fall onto the binary with little angular momentum. While accretion from a polar disc may not promote SMBH binary coalescence (solving the `final-parsec problem'), ejection of this infalling low-angular momentum material via gravitational slingshot is a possible mechanism to reduce the binary separation. Moreover, this process acts on dynamical rather than viscous time scales, and so is much faster.

We investigate the effect of black hole spin on warped or misaligned accretion discs - in particular i) whether or not the inner disc edge aligns with the black hole spin and ii) whether the disc can maintain a smooth transition between an aligned inner disc and a misaligned outer disc, known as the Bardeen-Petterson effect. We employ high resolution 3D smoothed particle hydrodynamics simulations of $\alpha$-discs subject to Lense-Thirring precession, focussing on the bending wave regime where the disc viscosity is smaller than the aspect ratio $\alpha \lesssim H/R$. We first address the controversy in the literature regarding possible steady-state oscillations of the tilt close to the black hole. We successfully recover such oscillations in 3D at both small and moderate inclinations ($\lesssim 15^{\circ}$), provided both Lense-Thirring and Einstein precession are present, sufficient resolution is employed, and provided the disc is not so thick so as to simply accrete misaligned. Second, we find that discs inclined by more than a few degrees in general steepen and break rather than maintain a smooth transition, again in contrast to previous findings, but only once the disc scale height is adequately resolved. Finally, we find that when the disc plane is misaligned to the black hole spin by a large angle, the disc 'tears' into discrete rings which precess effectively independently and cause rapid accretion, consistent with previous findings in the diffusive regime ($\alpha \gtrsim H/R$). Thus misalignment between the disc and the spin axis of the black hole provides a robust mechanism for growing black holes quickly, regardless of whether the disc is thick or thin.

We extend previous studies of the tidal truncation of coplanar disks in binary systems to the more general case of noncoplanar disks. As in the prograde coplanar case, Lindblad resonances play a key role in tidal truncation. We analyze the tidal torque acting on a misaligned nearly circular disk in a circular orbit binary system. We concentrate on the 2:1 inner Lindblad resonance associated with the m=2 tidal forcing (for azimuthal wavenumber m) that plays a major role in the usual coplanar case. We determine the inclination dependence of this torque, which is approximately cos^8(i/2) for misalignment angle i. Compared to the prograde coplanar case (i=0), this torque decreases by a factor of about 2 for i = pi/6 and by a factor of about 20 for i=pi/2. The Lindblad torque decreases to zero for a tilt angle of pi (counter-rotation), consistent with previous investigations. The effects of higher order resonances associated with m>2 tidal forcing may contribute somewhat, but are much more limited than in the i=0 case. These results suggest that misaligned disks in binary systems can be significantly extended compared to their coplanar counterparts. In cases where a disk is sufficiently inclined and viscous, it can overrun all Lindblad resonances and overflow the Roche lobe of the disk central object.

We present spectrally and spatially-resolved maps of HNC and HC$_3$N emission from Titan's atmosphere, obtained using the Atacama Large Millimeter/submillimeter Array (ALMA) on 2013 November 17. These maps show anisotropic spatial distributions for both molecules, with resolved emission peaks in Titan's northern and southern hemispheres. The HC$_3$N maps indicate enhanced concentrations of this molecule over the poles, consistent with previous studies of Titan's photochemistry and atmospheric circulation. Differences between the spectrally-integrated flux distributions of HNC and HC$_3$N show that these species are not co-spatial. The observed spectral line shapes are consistent with HNC being concentrated predominantly in the mesosphere and above (at altitudes $z\gtrsim 400$ km), whereas HC$_3$N is abundant at a broader range of altitudes ($z\approx70$-600 km). From spatial variations in the HC$_3$N line profile, the locations of the HC$_3$N emission peaks are shown to be variable as a function of altitude. The peaks in the integrated emission from HNC and the line core (upper-atmosphere) component of HC$_3$N (at $z\gtrsim300$ km) are found to be asymmetric with respect to Titan's polar axis, indicating that the mesosphere may be more longitudinally-variable than previously thought. The spatially-integrated HNC and HC$_3$N spectra are modeled using the NEMESIS planetary atmosphere code and the resulting best-fitting disk-averaged vertical mixing ratio (VMR) profiles are found to be in reasonable agreement with previous measurements for these species. Vertical column densities of the best-fitting gradient models for HNC and HC$_3$N are $1.9\times10^{13}$ cm$^{-2}$ and $2.3\times10^{14}$ cm$^{-2}$, respectively.

We present the results of high-resolution numerical simulations of gas clouds falling onto binary supermassive black holes to form circumbinary accretion discs, with both prograde and retrograde cloud orbits. We explore a range of clouds masses and cooling rates. We find that for low mass discs that cool fast enough to fragment, prograde discs are significantly shorter-lived than similar discs orbiting retrograde with respect to the binary. For fragmenting discs of all masses, we also find that prograde discs fragment across a narrower radial region. If the cooling is slow enough that the disc enters a self-regulating gravitoturbulent regime, we find that alignment between the disc and binary planes occurs on a timescale primarily dictated by the disc thickness. We estimate realistic cooling times for such discs, and find that in the majority of cases we expect fragmentation to occur. The longer lifetime of low-mass fragmenting retrograde discs allows them to drive significant binary evolution, and may provide a mechanism for solving the 'last parsec problem'.

We use three dimensional hydrodynamical simulations to show that a highly misaligned accretion disk around one component of a binary system can exhibit global Kozai-Lidov cycles, where the inclination and eccentricity of the disk are interchanged periodically. This has important implications for accreting systems on all scales, for example, the formation of planets and satellites in circumstellar and circumplanetary disks, outbursts in X-ray binary systems and accretion on to supermassive black holes.

Be/X-ray binary systems exhibit both periodic (Type I) X-ray outbursts and giant (Type II) outbursts, whose origin has remained elusive. We suggest that Type II X-ray outbursts occur when a highly misaligned decretion disk around the Be star becomes eccentric, allowing the compact object companion to capture a large amount of material at periastron. Using 3D smoothed particle hydrodynamics simulations we model the long term evolution of a representative Be/X-ray binary system. We find that periodic (Type I) X-ray outbursts occur when the neutron star is close to periastron for all disk inclinations. Type II outbursts occur for large misalignment angles and are associated with eccentricity growth that occurs on a timescale of about 10 orbital periods. Mass capture from the eccentric decretion disk results in an accretion disk around the neutron star whose estimated viscous time is long enough to explain the extended duration of Type II outbursts. Previous studies suggested that the outbursts are caused by a warped disk but our results suggest that this is not sufficient, the disk must be both highly misaligned and eccentric to initiate a Type II accretion event.

We present an accretion cycle which can explain state transitions and other observed phenomena in black hole X-ray binaries. This model is based on the process of disc tearing, where individual rings of gas break off the disc and precess effectively independently. This occurs when the Lense-Thirring effect is stronger than the local disc viscosity. We discuss implications of this model for quasi-periodic oscillations and the disc-jet-corona coupling. We also speculate on applying this model to active galactic nuclei and other accreting systems.

The Voyager 1 flyby of Titan in 1980 gave a first glimpse of the chemical complexity of Titan's atmosphere, detecting many new molecules with the infrared spectrometer (IRIS). These included propane (C3H8) and propyne (CH3C2H), while the intermediate-sized C3Hx hydrocarbon (C3H6) was curiously absent. Using spectra from the Composite Infrared Spectrometer (CIRS) on Cassini, we show the first positive detection of propene (C3H6) in Titan's stratosphere (5-sigma significance), finally filling the three-decade gap in the chemical sequence. We retrieve a vertical abundance profile from 100-250 km, that varies slowly with altitude from ~2 ppbv at 100 km, to ~5 ppbv at 200 km. The abundance of C3H6 is less than both C3H8 and CH3C2H, and we remark on an emerging paradigm in Titan's hydrocarbon abundances whereby: alkanes > alkynes > alkenes within the C2Hx and C3Hx chemical families in the lower stratosphere. More generally, there appears to be much greater ubiquity and relative abundance of triple-bonded species than double-bonded, likely due to the greater resistance of triple bonds to photolysis and chemical attack.

Microlensing observations indicate that quasar accretion discs have half-light radii larger than expected from standard theoretical predictions based on quasar fluxes or black hole masses. Blackburne and colleagues have also found a very weak wavelength dependence of these half-light radii. We consider disc temperature profile models that might match these observations. Nixon and colleagues have suggested that misaligned accretion discs around spinning black holes will be disrupted at radii small enough for the Lense-Thirring torque to overcome the disc's viscous torque. Gas in precessing annuli torn off a disc will spread radially and intersect with the remaining disc, heating the disc at potentially large radii. However, if the intersection occurs at an angle of more than a degree or so, highly supersonic collisions will shock-heat the gas to a Compton temperature of T~10^7 K, and the spectral energy distributions (SEDs) of discs with such shock-heated regions are poor fits to observations of quasar SEDs. Torn discs where heating occurs in intermittent weak shocks that occur whenever the intersection angle reaches a tenth of a degree pose less of a conflict with observations, but do not have significantly larger half-light radii than standard discs. We also study two phenomenological disc temperature profile models. We find that discs with a temperature spike at relatively large radii and lowered temperatures at radii inside the spike yield improved and acceptable fits to microlensing sizes in most cases. Such temperature profiles could in principle occur in sub-Keplerian discs partially supported by magnetic pressure. However, such discs overpredict the fluxes from quasars studied with microlensing except in the limit of negligible continuum emission from radii inside the temperature spike.

We review recent progress in studying accretion flows on to supermassive black holes (SMBH). Much of this removes earlier assumptions of symmetry and regularity, such as aligned and prograde disc rotation. This allows a much richer variety of effects, often because cancellation of angular momentum allows rapid infall. Potential applications include lower SMBH spins allowing faster mass growth and suppressing gravitational-wave reaction recoil in mergers, gas-assisted SMBH mergers, and near-dynamical accretion in galaxy centres.

In a recent paper we have shown that the evolution of a misaligned disc around a spinning black hole can result in tearing the disc into many distinct planes. Tearing discs with random orientations produce direct dynamical accretion on to the hole in approximately 70% of all cases. Here we examine the evolution of a misaligned disc around a binary system. We show that these discs are susceptible to tearing for almost all inclinations. We also show that tearing of the disc can result in a significant acceleration of the disc evolution and subsequent accretion on to the binary - by factors up to 10,000 times that of a coplanar prograde disc with otherwise identical parameters. This provides a promising mechanism for driving mergers of supermassive black hole (SMBH) binaries on timescales much shorter than a Hubble time. Disc tearing also suggests new observational signatures of accreting SMBH binaries, and other systems such as protostellar binaries.

In this chapter we describe the remote sensing measurement of nitrogen-bearing species in Titan's atmosphere by the Composite Infrared Spectrometer (CIRS) on the Cassini spacecraft. This instrument, which detects the thermal infrared spectrum from 10-1500 cm-1 (1000-7 microns) is sensitive to vibrational emissions of gases and condensates in Titan's stratosphere and lower mesosphere, permitting the measurement of ambient temperature and the abundances of gases and particulates. Three N-bearing species are firmly detected: HCN, HC3N and C2N2, and their vertical and latitudinal distributions have been mapped. In addition, ices of HC3N and possibly C4N2 are also seen in the far-infrared spectrum at high latitudes during the northern winter. The HC(15)N isotopologue has been measured, permitting the inference of the 14N/15N ratio in this species, which differs markedly (lower) than in the bulk nitrogen reservoir (N2). We also describe the search in the CIRS spectrum, and inferred upper limits, for NH3 and CH3CN. CIRS is now observing seasonal transition on Titan and the gas abundance distributions are changing accordingly, acting as tracers of the changing atmospheric circulation. The prospects for further CIRS science in the remaining five years of the Cassini mission are discussed.

We have searched for the presence of simple P and S-bearing molecules in Titan's atmosphere, by looking for the characteristic signatures of phosphine and hydrogen sulfide in infrared spectra obtained by Cassini CIRS. As a result we have placed the first upper limits on the stratospheric abundances, which are 1 ppb (PH3) and 330 ppb (H2S), at the 2-sigma significance level.

Observations of accreting black holes often provoke suggestions that their jets precess. The precession is usually supposed to result from a combination of the Lense-Thirring effect and accretion disc viscosity. We show that this is unlikely for any type of black hole system, as the disc generally has too little angular momentum compared with a spinning hole to cause any significant movement of the jet direction across the sky on short timescales. Uncorrelated accretion events, as in the chaotic accretion picture of active galactic nuclei, change AGN jet directions only on timescales \gtrsim 10^7 yr. In this picture AGN jet directions are stable on shorter timescales, but uncorrelated with any structure of the host galaxy, as observed. We argue that observations of black-hole jets precessing on timescales short compared to the accretion time would be a strong indication that the accretion disc, and not the standard Blandford-Znajek mechanism, is responsible for driving the jet. This would be particularly convincing in a tidal disruption event. We suggest that additional disc physics is needed to explain any jet precession on timescales short compared with the accretion time. Possibilities include the radiation warping instability, or disc tearing.

We show that in realistic cases of accretion in active galactic nuclei or stellar-mass X-ray binaries, the Lense-Thirring effect breaks the central regions of tilted accretion discs around spinning black holes into a set of distinct planes with only tenuous flows connecting them. If the original misalignment of the outer disc to the spin axis of the hole is $45^{\circ} \lesssim \theta \lesssim 135^{\circ}$, as in $\sim 70$% of randomly oriented accretion events, the continued precession of these discs sets up partially counter-rotating gas flows. This drives rapid infall as angular momentum is cancelled and gas attempts to circularize at smaller radii. Disc breaking close to the black hole leads to direct dynamical accretion, while breaking further out can drive gas down to scales where it can accrete rapidly. For smaller tilt angles breaking can still occur, and may lead to other observable phenomena such as QPOs. For such effects not to appear, the black hole spin must in practice be negligibly small, or be almost precisely aligned with the disc. Qualitatively similar results hold for any accretion disc subject to a forced differential precession, such as an external disc around a misaligned black hole binary.

In general, when gas accretes on to a supermassive black hole binary it is likely to have no prior knowledge of the binary angular momentum. Therefore a circumbinary disc forms with a random inclination angle, theta, to the binary. It is known that for theta < 90 degrees the disc will coalign wrt the binary. If theta > 90 degrees the disc wholly counteraligns if it satisfies cos(theta) < -J_d/2J_b, where J_d and J_b are the magnitudes of the disc and binary angular momentum vectors respectively. If however theta > 90 degrees and this criterion is not satisfied the same disc may counteralign its inner regions and, on longer timescales, coalign its outer regions. I show that for typical disc parameters, describing an accretion event on to a supermassive black hole binary, a misaligned circumbinary disc is likely to wholly co-- or counter--align with the binary plane. This is because the binary angular momentum dominates the disc angular momentum. However with extreme parameters (binary mass ratio M_2/M_1 << 1 or binary eccentricity e ~ 1) the same disc may simultaneously co- and counter-align. It is known that coplanar prograde circumbinary discs are stable. I show that coplanar retrograde circumbinary discs are also stable. A chaotic accretion event on to an SMBH binary will therefore result in a coplanar circumbinary disc that is either prograde or retrograde with respect to the binary plane.

Accretion in the nuclei of active galaxies may occur chaotically. This can produce accretion discs which are counter-rotating or strongly misaligned with respect to the spin of the central supermassive black hole (SMBH), or the axis of a close SMBH binary. Accordingly we consider the cancellation of angular momentum in accretion discs with a significant change of plane (tilt) between inner and outer parts. We estimate analytically the maximum accretion rate through such discs and compare this with the results of Smoothed Particle Hydrodynamics (SPH) simulations. These suggest that accretion rates on to supermassive black holes may be larger by factors $\gtrsim 100$ if the disc is internally tilted in this way rather than planar. This offers a natural way of driving the rapid growth of supermassive black holes, and the coalescence of SMBH binaries.

We simulate the viscous evolution of an accretion disc around a spinning black hole. In general any such disc is misaligned, and warped by the Lense-Thirring effect. Unlike previous studies we use effective viscosities constrained to be consistent with the internal fluid dynamics of the disc. We find that nonlinear fluid effects, which reduce the effective viscosities in warped regions, can promote the breaking of the disc into two distinct planes. This occurs when the Shakura & Sunyaev dimensionless viscosity parameter alpha is <~ 0.3 and the initial angle of misalignment between the disc and hole is >~ 45 degrees. The break can be a long-lived feature, propagating outwards in the disc on the usual alignment timescale, after which the disc is fully co- or counter-aligned with the hole. Such a break in the disc may be significant in systems where we know the inclination of the outer accretion disc to the line of sight, such as some X-ray binaries: the inner disc, and so any jets, may be noticeably misaligned with respect to the orbital plane.

We consider the interaction between a binary system (e.g. two supermassive black holes or two stars) and an external accretion disc with misaligned angular momentum. This situation occurs in galaxy merger events involving supermassive black holes, and in the formation of stellar--mass binaries in star clusters. We work out the gravitational torque between the binary and disc, and show that their angular momenta J_b, J_d stably counteralign if their initial orientation is sufficiently retrograde, specifically if the angle theta between them obeys cos(theta) < -J_d/2J_b, on a time short compared with the mass gain time of the central accretor(s). The magnitude J_b remains unchanged in this process. Counteralignment can promote the rapid merger of supermassive black hole binaries, and possibly the formation of coplanar but retrograde planets around stars in binary systems.

In this paper we describe a first quantitative search for several molecules in Titan's stratosphere in Cassini CIRS infrared spectra. These are: ammonia (NH3), methanol (CH3OH), formaldehyde (H2CO), and acetonitrile (CH3CN), all of which are predicted by photochemical models but only the last of which observed, and not in the infrared. We find non-detections in all cases, but derive upper limits on the abundances from low-noise observations at 25{\deg}S and 75{\deg}N. Comparing these constraints to model predictions, we conclude that CIRS is highly unlikely to see NH3 or CH3OH emissions. However, CH3CN and H2CO are closer to CIRS detectability, and we suggest ways in which the sensitivity threshold may be lowered towards this goal.

We investigate whether a circumbinary gas disc can coalesce a supermassive black hole binary system in the centre of a galaxy. This is known to be problematic for a prograde disc. We show that in contrast, interaction with a retrograde circumbinary disc is considerably more effective in shrinking the binary because there are no orbital resonances. The binary directly absorbs negative angular momentum from the circumbinary disc by capturing gas into a disc around the secondary black hole, or discs around both holes if the binary mass ratio is close to unity. In many cases the binary orbit becomes eccentric, shortening the pericentre distance as the eccentricity grows. In all cases the binary coalesces once it has absorbed the angular momentum of a gas mass comparable to that of the secondary black hole. Importantly, this conclusion is unaffected even if the gas inflow rate through the disc is formally super--Eddington for either hole. The coalescence timescale is therefore always $\sim M_2/\dot M$, where $M_2$ is the secondary black hole mass and $\dot M$ the inflow rate through the circumbinary disc.

The flybys of Jupiter by the Voyager spacecraft in 1979, and over two decades later by Cassini in 2000, have provided us with unique datasets from two different epochs, allowing the investigation of seasonal change in the atmosphere. In this paper we model zonal averages of thermal infrared spectra from the two instruments, Voyager 1 IRIS and Cassini CIRS, to retrieve the vertical and meridional profiles of temperature, and the abundances of the two minor hydrocarbons, acetylene (C2H2) and ethane (C2H6). The spatial variation of these gases is controlled by both chemistry and dynamics, and therefore their observed distribution gives us an insight into both processes. We find that the two gases paint quite different pictures of seasonal change. Whilst the 2-D cross-section of C2H6 abundance is slightly increased and more symmetric in 2000 (northern summer solstice) compared to 1979 (northern fall equinox), the major trend of equator to pole increase remains. For C2H2 on the other hand, the Voyager epoch exhibits almost no latitudinal variation, whilst the Cassini era shows a marked decrease polewards in both hemispheres. At the present time, these experimental findings are in advance of interpretation, as there are no published models of 2-D Jovian seasonal chemical variation available for comparison.

We consider the response of a protostellar disc to a tidally induced warp and the resultant changes in the spectral energy distribution (SED). We argue that for typical protostellar disc parameters the warp is communicated through the disc in a wave-like fashion. We find that the main effects of the warp tend to be at large radii (greater than 30 AU) and, for sufficiently small viscosity, can be quite long-lived. This can result in non-uniform illumination of the disc at these radii and can induce significant changes to the SED at wavelengths greater than 100 microns.

In this paper we select large spectral averages of data from the Cassini Composite Infrared Spectrometer (CIRS) obtained in limb-viewing mode at low latitudes (30S--30N), greatly increasing the path length and hence signal-to-noise ratio for optically thin trace species such as propane. By modeling and subtracting the emissions of other gas species, we demonstrate that at least six infrared bands of propane are detected by CIRS, including two not previously identified in Titan spectra. Using a new line list for the range 1300-1400cm -1, along with an existing GEISA list, we retrieve propane abundances from two bands at 748 and 1376 cm-1. At 748 cm-1 we retrieve 4.2 +/- 0.5 x 10(-7) (1-sigma error) at 2 mbar, in good agreement with previous studies, although lack of hotbands in the present spectral atlas remains a problem. We also determine 5.7 +/- 0.8 x 10(-7) at 2 mbar from the 1376 cm-1 band - a value that is probably affected by systematic errors including continuum gradients due to haze and also an imperfect model of the n6 band of ethane. This study clearly shows for the first time the ubiquity of propane's emission bands across the thermal infrared spectrum of Titan, and points to an urgent need for further laboratory spectroscopy work, both to provide the line positions and intensities needed to model these bands, and also to further characterize haze spectral opacity. The present lack of accurate modeling capability for propane is an impediment not only for the measurement of propane itself, but also for the search for the emissions of new molecules in many spectral regions.