In order to investigate the importance of different proposed quenching mechanisms, we use an indirect method to estimate gas masses for ~62,000 SDSS DR7 galaxies. We infer gas surface densities from dust column densities as traced by extinction within the fibre, applying a metallicity correction to account for varying dust-to-gas ratios. We find that both gas fraction and star formation efficiency (SFE) decrease moving away from the star forming main sequence (MS) towards quiescence for all galaxy masses. We further show that both quantities correlate similarly strongly with the departure from the MS, implying the need for any physical model of quenching to invoke a change in $\textit{both}$ gas fraction and SFE. Our results call for a better understanding of the physical processes driving the decrease in star formation efficiency, which has received relatively little attention in the theory of quenching until now.

We present an analysis of star formation and quenching in the SDSS-IV MaNGA-DR15, utilising over 5 million spaxels from $\sim$3500 local galaxies. We estimate star formation rate surface densities ($\Sigma_{\rm SFR}$) via dust corrected $H\alpha$ flux where possible, and via an empirical relationship between specific star formation rate (sSFR) and the strength of the 4000 Angstrom break (D4000) in all other cases. We train a multi-layered artificial neural network (ANN) and a random forest (RF) to classify spaxels into `star forming' and `quenched' categories given various individual (and groups of) parameters. We find that global parameters (pertaining to the galaxy as a whole) perform collectively the best at predicting when spaxels will be quenched, and are substantially superior to local/ spatially resolved and environmental parameters. Central velocity dispersion is the best single parameter for predicting quenching in central galaxies. We interpret this observational fact as a probable consequence of the total integrated energy from AGN feedback being traced by the mass of the black hole, which is well known to correlate strongly with central velocity dispersion. Additionally, we train both an ANN and RF to estimate $\Sigma_{\rm SFR}$ values directly via regression in star forming regions. Local/ spatially resolved parameters are collectively the most predictive at estimating $\Sigma_{\rm SFR}$ in these analyses, with stellar mass surface density at the spaxel location ($\Sigma_*$) being by far the best single parameter. Thus, quenching is fundamentally a global process but star formation is governed locally by processes within each spaxel.

Pseudospectral schemes are a class of numerical methods capable of solving smooth problems with high accuracy thanks to their exponential convergence to the true solution. When applied to discontinuous problems, such as fluid shocks and material interfaces, due to the Gibbs phenomenon, pseudospectral solutions lose their superb convergence and suffer from spurious oscillations across the entire computational domain. Luckily, there exist theoretical remedies for these issues which have been successfully tested in practice for cases of well defined discontinuities. We focus on one piece of this procedure---detecting a discontinuity in spectral data. We show that realistic applications require treatment of discontinuities dynamically developing in time and that it poses challenges associated with shock detection. More precisely, smoothly steepening gradients in the solution spawn spurious oscillations due to insufficient resolution, causing premature shock identification and information loss. We improve existing spectral shock detection techniques to allow us to automatically detect true discontinuities and identify cases for which post-processing is required to suppress spurious oscillations resulting from the loss of resolution. We then apply these techniques to solve an inviscid Burgers' equation in 1D, demonstrating that our method correctly treats genuine shocks caused by wave breaking and removes oscillations caused by numerical constraints.

The shape of low-frequency radio continuum spectra of normal galaxies is not well understood, the key question being the role of physical processes such as thermal absorption in shaping them. In this work we take advantage of the LOFAR Multifrequency Snapshot Sky Survey (MSSS) to investigate such spectra for a large sample of nearby star-forming galaxies. Using the measured 150MHz flux densities from the LOFAR MSSS survey and literature flux densities at various frequencies we have obtained integrated radio spectra for 106 galaxies. The spectra are explained through the use of a three-dimensional model of galaxy radio emission, and radiation transfer dependent on the galaxy viewing angle and absorption processes. Spectra of our galaxies are generally flatter at lower compared to higher frequencies but as there is no tendency for the highly inclined galaxies to have more flattened low-frequency spectra, we argue that the observed flattening is not due to thermal absorption, contradicting the suggestion of Israel & Mahoney (1990). According to our modelled radio maps for M51-like galaxies, the free-free absorption effects can be seen only below 30MHz and in the global spectra just below 20MHz, while in the spectra of starburst galaxies, like M82, the flattening due to absorption is instead visible up to higher frequencies of about 150MHz. Locally, within galactic disks, the absorption effects are distinctly visible in M51-like galaxies as spectral flattening around 100-200MHz in the face-on objects, and as turnovers in the edge-on ones, while in M82-like galaxies there are strong turnovers at frequencies above 700MHz, regardless of viewing angle. Our modelling suggests that the weak spectral flattening observed in the nearby galaxies studied here results principally from synchrotron spectral curvature due to cosmic ray energy losses and propagation effects.

We present a new 300 ks Chandra observation of M87 that limits pileup to only a few per cent of photon events and maps the hot gas properties closer to the nucleus than has previously been possible. Within the supermassive black hole's gravitational sphere of influence, the hot gas is multiphase and spans temperatures from 0.2 to 1 keV. The radiative cooling time of the lowest temperature gas drops to only 0.1-0.5 Myr, which is comparable to its free fall time. Whilst the temperature structure is remarkably symmetric about the nucleus, the density gradient is steep in sectors to the N and S, with $\rho{\propto}r^{-1.5\pm0.1}$, and significantly shallower along the jet axis to the E, where $\rho{\propto}r^{-0.93\pm0.07}$. The density structure within the Bondi radius is therefore consistent with steady inflows perpendicular to the jet axis and an outflow directed E along the jet axis. By putting limits on the radial flow speed, we rule out Bondi accretion on the scale resolved at the Bondi radius. We show that deprojected spectra extracted within the Bondi radius can be equivalently fit with only a single cooling flow model, where gas cools from 1.5 keV down below 0.1 keV at a rate of 0.03 M$_{\odot}$/yr. For the alternative multi-temperature spectral fits, the emission measures for each temperature component are also consistent with a cooling flow model. The lowest temperature and most rapidly cooling gas in M87 is therefore located at the smallest radii at ~100 pc and may form a mini cooling flow. If this cooling gas has some angular momentum, it will feed into the cold gas disk around the nucleus, which has a radius of ~80 pc and therefore lies just inside the observed transition in the hot gas structure.