The Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) is a next-generation survey of neutral hydrogen (HI) in the Local Universe. It uses the widefield, high-resolution capability of the Australian Square Kilometer Array Pathfinder (ASKAP), a radio interferometer consisting of 36 x 12-m dishes equipped with Phased-Array Feeds (PAFs), located in an extremely radio-quiet zone in Western Australia. WALLABY aims to survey three-quarters of the sky (-90 degr < Dec < +30 degr) to a redshift of z < 0.26, and generate spectral line image cubes at ~30 arcsec resolution and ~1.6 mJy/beam per 4 km/s channel sensitivity. ASKAP's instantaneous field of view at 1.4 GHz, delivered by the PAF's 36 beams, is about 30 sq deg. At an integrated signal-to-noise ratio of five, WALLABY is expected to detect over half a million galaxies with a mean redshift of z ~ 0.05 (~200 Mpc). The scientific goals of WALLABY include: (a) a census of gas-rich galaxies in the vicinity of the Local Group; (b) a study of the HI properties of galaxies, groups and clusters, in particular the influence of the environment on galaxy evolution; and (c) the refinement of cosmological parameters using the spatial and redshift distribution of low-bias gas-rich galaxies. For context we provide an overview of previous large-scale HI surveys. Combined with existing and new multi-wavelength sky surveys, WALLABY will enable an exciting new generation of panchromatic studies of the Local Universe. - First results from the WALLABY pilot survey are revealed, with initial data products publicly available in the CSIRO ASKAP Science Data Archive (CASDA).

The black hole candidate and X-ray binary MAXI J1535-571 was discovered in September 2017. During the decay of its discovery outburst, and before returning to quiescence, the source underwent at least four reflaring events, with peak luminosities of $\sim$10$^{35-36}$ erg s$^{-1}$ (d/4.1 kpc)$^2$. To investigate the nature of these flares, we analysed a sample of NICER observations taken with almost daily cadence. In this work we present the detailed spectral and timing analysis of the evolution of the four reflares. The higher sensitivity of NICER at lower energies, in comparison with other X-ray detectors, allowed us to constrain the disc component of the spectrum at $\sim$0.5 keV. We found that during each reflare the source appears to trace out a q-shaped track in the hardness-intensity diagram similar to those observed in black hole binaries during full outbursts. MAXI J1535-571 transits between the hard state (valleys) and softer states (peaks) during these flares. Moreover, the Comptonised component is undetected at the peak of the first reflare, while the disc component is undetected during the valleys. Assuming the most likely distance of 4.1 kpc, we find that the hard-to-soft transitions take place at the lowest luminosities ever observed in a black hole transient, while the soft-to-hard transitions occur at some of the lowest luminosities ever reported for such systems.

We present a comparison of galaxy atomic and molecular gas properties in three recent cosmological hydrodynamic simulations, Simba, EAGLE, and Illustris-TNG, versus observations from $z\sim 0-2$. These simulations all rely on similar sub-resolution prescriptions to model cold interstellar gas which they cannot represent directly, and qualitatively reproduce the observed $z\approx 0$ HI and H$_2$ mass functions (HIMF, H2MF), CO(1-0) luminosity functions (COLF), and gas scaling relations versus stellar mass, specific star formation rate, and stellar surface density $\mu_*$, with some quantitative differences. To compare to the COLF, we apply an H$_2$-to-CO conversion factor to the simulated galaxies based on their average molecular surface density and metallicity, yielding substantial variations in $\alpha_{\rm CO}$ and significant differences between models. Using this, predicted $z=0$ COLFs agree better with data than predicted H2MFs. Out to $z\sim 2$, EAGLE's and Simba's HIMF and COLF strongly increase, while TNG's HIMF declines and COLF evolves slowly. EAGLE and Simba reproduce high $L_{\rm CO1-0}$ galaxies at $z\sim 1-2$ as observed, owing partly to a median $\alpha_{\rm CO}(z=2)\sim 1$ versus $\alpha_{\rm CO}(z=0)\sim 3$. Examining \HI, H$_2$, and CO scaling relations, their trends with $M_*$ are broadly reproduced in all models, but EAGLE yields too little HI in green valley galaxies, TNG and Simba overproduce cold gas in massive galaxies, and Simba overproduces molecular gas in small systems. Using Simba variants that exclude individual AGN feedback modules, we find that Simba's AGN jet feedback is primarily responsible by lowering cold gas contents from $z\sim 1\to0$ by suppressing cold gas in $M_*> 10^{10}{\rm M}_\odot$ galaxies, while X-ray feedback suppresses the formation of high-$\mu_*$ systems.

We use SHARK, a semi-analytic galaxy formation model, to investigate the physical processes involved in dictating the shape, scatter and evolution of the HI-halo mass relation at $0\leq z \leq 2$. We compare SHARK with HI clustering and spectral stacking of the HI-halo mass relation derived from observations finding excellent agreement with the former and a deficiency of HI in SHARK at $M_{\rm vir}\approx 10^{12-13} M_{\odot}$ in the latter, but otherwise great agreement below and above that mass threshold. In SHARK, we find that the HI mass increases with the halo mass up to a critical mass of $\approx 10^{11.8} M_{\odot}$; between $\sim 10^{11.8}-10^{13}M_{\odot}$, the scatter in the relation increases by 0.7 dex and the HI mass decreases with the halo mass on average; at $M_{\rm vir} \geq 10^{13} M_{\odot}$, the HI content continues to increase with halo mass. We find that the critical halo mass of $\approx 10^{12} M_{\odot}$ is largely set by feedback from Active Galactic Nuclei (AGN), and the exact shape and scatter of the HI-halo mass relation around that mass is extremely sensitive to how AGN feedback is modelled, with other physical processes playing a less significant role. We determine the main secondary parameters responsible for the scatter of the HI-halo mass relation, namely the halo spin parameter at $M_{\rm vir}\leq 10^{11.8} M_{\odot}$, and the fractional contribution from substructure to the total halo mass for $M_{\rm vir}\geq 10^{13} M_{\odot}$. The scatter at $10^{11.8}<M_{\rm vir}<10^{13} M_{\odot}$ is best described by the black-hole-to-stellar mass ratio of the central galaxy, reflecting the AGN feedback relevance. We present a numerical model to populate dark matter-only simulations with HI at $0\leq z \leq 2$ based solely on halo parameters that are measurable in such simulations.

For astronomers to make a significant contribution to the reduction of climate change-inducing greenhouse gas emissions, we first must quantify our sources of emissions and review the most effective approaches for reducing them. Here we estimate that Australian astronomers' total greenhouse gas emissions from their regular work activities are $\gtrsim$25 ktCO$_2$-e/yr (equivalent kilotonnes of carbon dioxide per year). This can be broken into $\sim$15 ktCO$_2$-e/yr from supercomputer usage, $\sim$4.2 ktCO$_2$-e/yr from flights (where individuals' flight emissions correlate with seniority), $>$3.3 ktCO$_2$-e/yr from the operation of observatories, and 2.6$\pm$0.4 ktCO$_2$-e/yr from powering office buildings. Split across faculty scientists, postdoctoral researchers, and PhD students, this averages to $\gtrsim$37 tCO$_2$-e/yr per astronomer, over 40% more than what the average Australian non-dependant emits in total, equivalent to $\sim$5$\times$ the global average. To combat these environmentally unsustainable practices, we suggest astronomers should strongly preference use of supercomputers, observatories, and office spaces that are predominantly powered by renewable energy sources. Where facilities that we currently use do not meet this requirement, their funders should be lobbied to invest in renewables, such as solar or wind farms. Air travel should also be reduced wherever possible, replaced primarily by video conferencing, which should also promote inclusivity.

We apply the spectral energy distribution-fitting code ProSpect to multiwavelength imaging for $\sim$7,000 galaxies from the GAMA survey at $z<0.06$, in order to extract their star-formation histories. We combine a parametric description of the star formation history with a closed-box evolution of metallicity where the present-day gas-phase metallicity of the galaxy is a free parameter. We show with this approach that we are able to recover the observationally-determined cosmic star formation history (CSFH), an indication that stars are being formed in the correct epoch of the Universe, on average, for the manner in which we are conducting SED fitting. We also show the contribution to the CSFH of galaxies of different present-day visual morphologies, and stellar masses. This analysis suggests that half of the mass in present-day elliptical galaxies was in place 11 Gyr ago, whereas in other morphological types the stellar mass formed later, up to 6 Gyr ago for present-day irregular galaxies. Similarly, the most massive galaxies in our sample were shown to have formed half their stellar mass by 10.5 Gyr ago, whereas the least massive galaxies formed half their stellar mass as late as 4 Gyr ago (the well-known effect of "galaxy downsizing"). Finally, our metallicity approach allows us to follow the average evolution in gas-phase metallicity for populations of galaxies, and extract the evolution of the cosmic metal mass density in stars and in gas, producing results in broad agreement with observations of metal densities in the Universe.

We investigate the influence of environment on the cold-gas properties of galaxies at z=0 within the TNG100 cosmological, magnetohydrodynamic simulation, part of the IllustrisTNG suite. We extend previous post-processing methods for breaking gas cells into their atomic and molecular phases, and build detailed mocks to comprehensively compare to the latest surveys of atomic hydrogen (HI) in nearby galaxies, namely ALFALFA and xGASS. We use TNG100 to explore the HI content, star formation activity, and angular momentum of satellite galaxies, each as a function of environment, and find that satellites are typically a factor of ~3 poorer in HI than centrals of the same stellar mass, with the exact offset depending sensitively on parent halo mass. Due to the large physical scales on which HI measurements are made (~45--245 kpc), contributions from gas not bound to the galaxy of interest but in the same line of sight crucially lead to larger HI mass measurements in the mocks in many cases, ultimately aligning with observations. This effect is mass-dependent and naturally greater for satellites than centrals, as satellites are never isolated by definition. We also show that HI stripping in TNG100 satellites is closely accompanied by quenching, in tension with observational data that instead favour that HI is preferentially stripped before star formation is reduced.

In November 2019, MAXI detected an X-ray outburst from the known Be X-ray binary system RX J0209.6-7427 located in the outer wing of the Small Magellanic Cloud. We followed the outburst of the system with NICER which led to the discovery of X-ray pulsations with a period of 9.3 s. We analyzed simultaneous X-ray data obtained with NuSTAR and NICER allowing us to characterize the spectrum and provide an accurate estimate of its bolometric luminosity. During the outburst the maximum broadband X-ray luminosity of the system reached $1-2\times10^{39}$ erg/s, thus exceeding by about one order of magnitude the Eddington limit for a typical 1.4 $M_{\odot}$ mass neutron star (NS). Monitoring observations with Fermi/GBM and NICER allowed us to study the spin evolution of the NS and compare it with standard accretion torque models. We found that the NS magnetic field should be of the order of $3\times10^{12}$ G. We conclude that RX J0209.6-7427 exhibited one of the brightest outbursts observed from a Be X-ray binary pulsar in the Magellanic Clouds, reaching similar luminosity level to the 2016 outburst of SMC X-3. Despite the super-Eddington luminosity of RX J0209.6-7427, the NS appears to have only a moderate magnetic field strength.

We use three semi-analytic models (SAMs) of galaxy formation and evolution, run on the same 1$h^{-1}$Gpc MultiDark Planck2 cosmological simulation, to investigate the properties of [OII] emission line galaxies in the redshift range $0.6<z<1.2$. We compare model predictions with different observational data sets, including DEEP2--Firefly galaxies with absolute magnitudes. We estimate the [OII] luminosity, L[OII], using simple relations derived both from the models and observations and also using a public code. This code ideally uses as input instantaneous star formation rates (SFRs), which are only provided by one of the SAMs under consideration. We use this SAM to study the feasibility of inferring galaxies' L[OII] for models that only provide average SFRs. We find that the post-processing computation of L[OII] from average SFRs is accurate for model galaxies with dust attenuated L[OII]$\lesssim10^{42.2}$erg s$^{-1}$ ($<5\%$ discrepancy). We also explore how to derive the [OII] luminosity from simple relations using global properties usually output by SAMs. Besides the SFR, the model L[OII] is best correlated with the observed-frame $u$ and $g$ broad-band magnitudes. These correlations have coefficients (r-values) above 0.64 and a dispersion that varies with L[OII]. We use these correlations and an observational one based on SFR and metallicity to derive L[OII]. These relations result in [OII] luminosity functions and halo occupation distributions with shapes that vary depending on both the model and the method used. Nevertheless, for all the considered models, the amplitude of the clustering at scales above 1$h^{-1}$Mpc remains unchanged independently of the method used to derive L[OII].

We present Neutron Star Interior Composition Explorer X-ray and Arcminute Microkelvin Imager Large Array radio observations of a rapid hard-to-soft state transition in the black hole X-ray transient MAXI J1820+070. During the transition from the hard state to the soft state a switch between two particular types of quasiperiodic oscillations (QPOs) was seen in the X-ray power density spectra, from type-C to type-B, along with a drop in the strength of the broadband X-ray variability and a brief flare in the 7-12 keV band. Soon after this switch (~1.5-2.5 hr) a strong radio flare was observed that corresponded to the launch of superluminal ejecta. Although hints of a connection between QPO transitions and radio flares have been seen in other black hole X-ray transients, our observations constitute the strongest observational evidence to date for a link between the appearance of type-B QPOs and the launch of discrete jet ejections.

We describe the design and assembly of the LUX-ZEPLIN experiment, a direct detection search for cosmic WIMP dark matter particles. The centerpiece of the experiment is a large liquid xenon time projection chamber sensitive to low energy nuclear recoils. Rejection of backgrounds is enhanced by a Xe skin veto detector and by a liquid scintillator Outer Detector loaded with gadolinium for efficient neutron capture and tagging. LZ is located in the Davis Cavern at the 4850' level of the Sanford Underground Research Facility in Lead, South Dakota, USA. We describe the major subsystems of the experiment and its key design features and requirements.

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

We analytically derive the observed size-mass relation of galaxies' atomic hydrogen (HI), including limits on its scatter, based on simple assumptions about the structure of HI discs. We trial three generic profiles for HI surface density as a function of radius. Firstly, we assert that HI surface densities saturate at a variable threshold, and otherwise fall off exponentially with radius or, secondly, radius squared. Our third model assumes the total gas surface density is exponential, with the HI fraction at each radius depending on local pressure. These are tested against a compilation of 110 galaxies from the THINGS, LITTLE THINGS, LVHIS, and Bluedisk surveys, whose HI surface density profiles are well resolved. All models fit the observations well and predict consistent size-mass relations. Using an analytical argument, we explain why processes that cause gas disc truncation - such as ram-pressure stripping - scarcely affect the HI size-mass relation. This is tested with the IllustrisTNG(100) cosmological, hydrodynamic simulation and the Dark Sage semi-analytic model of galaxy formation, both of which capture radially resolved disc structure. For galaxies with m_*>10^9 M_solar and m_HI>10^8 M_solar, both simulations predict HI size-mass relations that align with observations, show no difference between central and satellite galaxies, and show only a minor, second-order dependence on host halo mass for satellites. Ultimately, the universally tight HI size-mass relation is mathematically inevitable and robust. Only by completely disrupting the structure of HI discs, e.g. through overly powerful feedback, could a simulation predict the relation poorly.

We present a spectral study of the ultraluminous Be/X-ray transient pulsar Swift J0243.6+6124 using Neutron Star Interior Composition Explorer (NICER) observations during the system's 2017--2018 giant outburst. The 1.2--10~keV energy spectrum of the source can be approximated with an absorbed cut-off power law model. We detect strong, luminosity-dependent emission lines in the 6--7 keV energy range. A narrow 6.42 keV line, observed in the sub-Eddington regime, is seen to evolve into a broad Fe-line profile in the super-Eddington regime. Other features are found at 6.67 and 6.97 keV in the Fe-line complex. An asymmetric broad line profile, peaking at 6.67 keV, is possibly due to Doppler effects and gravitational redshift. The 1.2--79 keV broadband spectrum from NuSTAR and NICER observations at the outburst peak is well described by an absorbed cut-off power law plus multiple Gaussian lines and a blackbody component. Physical reflection models are also tested to probe the broad iron line feature. Depending on the mass accretion rate, we found emission sites that are evolving from ~5000 km to a range closer to the surface of the neutron star. Our findings are discussed in the framework of the accretion disk and its implication on the magnetic field, the presence of optically thick accretion curtain in the magnetosphere, jet emission, and the massive, ultra-fast outflow expected at super-Eddington accretion rates. We do not detect any signatures of a cyclotron absorption line in the NICER or NuSTAR data.

This paper describes the design and implementation of Stingray, a library in Python built to perform time series analysis and related tasks on astronomical light curves. Its core functionality comprises a range of Fourier analysis techniques commonly used in spectral-timing analysis, as well as extensions for analyzing pulsar data, simulating data sets, and statistical modeling. Its modular build allows for easy extensions and incorporation of its methods into data analysis workflows and pipelines. We aim for the library to be a platform for the implementation of future spectral-timing techniques. Here, we describe the overall vision and framework, core functionality, extensions, and connections to high-level command-line and graphical interfaces. The code is well-tested, with a test coverage of currently 95%, and is accompanied by extensive API documentation and a set of step-by-step tutorials.

We have recently developed a post-processing framework to estimate the abundance of atomic and molecular hydrogen (HI and H2, respectively) in galaxies in large-volume cosmological simulations. Here we compare the HI and H2 content of IllustrisTNG galaxies to observations. We mostly restrict this comparison to $z \approx 0$ and consider six observational metrics: the overall abundance of HI and H2, their mass functions, gas fractions as a function of stellar mass, the correlation between H2 and star formation rate, the spatial distribution of gas, and the correlation between gas content and morphology. We find generally good agreement between simulations and observations, particularly for the gas fractions and the HI mass-size relation. The H2 mass correlates with star formation rate as expected, revealing an almost constant depletion time that evolves up to z = 2 as observed. However, we also discover a number of tensions with varying degrees of significance, including an overestimate of the total neutral gas abundance at z = 0 by about a factor of two and a possible excess of satellites with no or very little neutral gas. These conclusions are robust to the modelling of the HI/H2 transition. In terms of their neutral gas properties, the IllustrisTNG simulations represent an enormous improvement over the original Illustris run. All data used in this paper are publicly available as part of the IllustrisTNG data release.

Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

We use TheThreeHundred project, a suite of 324 resimulated massive galaxy clusters embedded in a broad range of environments, to investigate (i) how the gas content of surrounding haloes correlates with phase-space position at $z=0$, and (ii) to investigate the role that ram pressure plays in this correlation. By stacking all 324 normalised phase-space planes containing 169287 haloes and subhaloes, we show that the halo gas content is tightly correlated with phase-space position. At $\sim\,1.5-2\,\text{R}_{\text{200}}$ of the cluster dark matter halo, we find an extremely steep decline in the halo gas content of infalling haloes and subhaloes irrespective of cluster mass, possibly indicating the presence of an accretion shock. We also find that subhaloes are particularly gas-poor, even in the cluster outskirts, which could indicate active regions of ongoing pre-processing. By modelling the instantaneous ram pressure experienced by each halo and subhalo at $z=0$, we show that the ram pressure intensity is also well correlated with phase-space position, which is again irrespective of cluster mass. In fact, we show that regions in the phase-space plane with high differential velocity between a halo or subhalo and its local gas environment, are almost mutually exclusive with high halo gas content regions. This suggests a causal link between the gas content of objects and the instantaneous ram pressure they experience, where the dominant factor is the differential velocity.

The geometry of the accretion flow around stellar-mass black holes can change on timescales of days to months. When a black hole emerges from quiescence (that is, it "turns on" after accreting material from its companion) it has a very hard (high-energy) X-ray spectrum produced by a hot corona positioned above its accretion disk, and then transitions to a soft (lower-energy) spectrum dominated by emission from the geometrically thin accretion disk, which extends to the innermost stable circular orbit. Much debate persists over how this transition occurs and whether it is driven largely by a reduction in the truncation radius of the disk or by a reduction in the spatial extent of the corona. Observations of X-ray reverberation lags in supermassive black-hole systems suggest that the corona is compact and that the disk extends nearly to the central black hole. Observations of stellar-mass black holes, however, reveal equivalent (mass-scaled) reverberation lags that are much larger, leading to the suggestion that the accretion disk in the hard X-ray state of stellar-mass black holes is truncated at a few hundreds of gravitational radii from the black hole. Here we report X-ray observations of the black-hole transient MAXI J1820+070. We find that the reverberation time lags between the continuum-emitting corona and the irradiated accretion disk are 6 to 20 times shorter than previously seen. The timescale of the reverberation lags shortens by an order of magnitude over a period of weeks, whereas the shape of the broadened iron K emission line remains remarkably constant. This suggests a reduction in the spatial extent of the corona, rather than a change in the inner edge of the accretion disk.

Recently a number of studies have found a similarity between the passive fraction of central and satellite galaxies when controlled for both stellar and halo mass. These results suggest that the quenching processes that affect galaxies are largely agnostic to central/satellite status, which contradicts the traditional picture of increased satellite quenching via environmental processes such as stripping, strangulation and starvation. Here we explore this further using the Galaxy And Mass Assembly (GAMA) survey which extends to ~2dex lower in stellar mass than SDSS, is more complete for closely-separated galaxies (>95% compared to >70%), and identifies lower-halo-mass groups outside of the very local Universe (M$_{\mathrm{halo}}\sim10^{12}$M$_{\odot}$ at $0.1<z<0.2$). As far as possible we aim to replicate the selections, completeness corrections and central/satellite division of one of the previous studies but find clear differences between passive fractions of centrals and satellites. We also find that our passive fractions increase with both halo-to-satellite mass ratio and central-to-second rank mass ratio. This suggests that quenching is more efficient in satellites that are low-mass for their halo ($i.e$ at high halo-to-satellite mass ratio in comparison to low halo-to-satellite mass ratio) and are more likely to be passive in older groups - forming a consistent picture of environmental quenching of satellites. We then discuss potential explanations for the previously observed similarity, such as dependence on the group-finding method.

By analysing a sample of galaxies selected from the HI Parkes All Sky Survey (HIPASS) to contain more than 2.5 times their expected HI content based on their optical properties, we investigate what drives these HI eXtreme (HIX) galaxies to be so HI-rich. We model the HI kinematics with the Tilted Ring Fitting Code TiRiFiC and compare the observed HIX galaxies to a control sample of galaxies from HIPASS as well as simulated galaxies built with the semi-analytic model Dark Sage. We find that (1) HI discs in HIX galaxies are more likely to be warped and more likely to host HI arms and tails than in the control galaxies, (2) the average HI and average stellar column density of HIX galaxies is comparable to the control sample, (3) HIX galaxies have higher HI and baryonic specific angular momenta than control galaxies, (4) most HIX galaxies live in higher-spin haloes than most control galaxies. These results suggest that HIX galaxies are HI-rich because they can support more HI against gravitational instability due to their high specific angular momentum. The majority of the HIX galaxies inherits their high specific angular momentum from their halo. The HI content of HIX galaxies might be further increased by gas-rich minor mergers. This paper is based on data obtained with the Australia Telescope Compact Array (ATCA) through the large program C 2705.

In this paper we describe the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to accretion flows in the strong field gravity regime around both stellar-mass and supermassive black-holes. eXTP has the unique capability of using advanced 'spectral-timing-polarimetry' techniques to analyze the rapid variations with three orthogonal diagnostics of the flow and its geometry, yielding unprecedented insight into the inner accreting regions, the effects of strong field gravity on the material within them and the powerful outflows which are driven by the accretion process.

Large-scale cosmological simulations of galaxy formation currently do not resolve the densities at which molecular hydrogen forms, implying that the atomic-to-molecular transition must be modeled either on the fly or in postprocessing. We present an improved postprocessing framework to estimate the abundance of atomic and molecular hydrogen and apply it to the IllustrisTNG simulations. We compare five different models for the atomic-to-molecular transition, including empirical, simulation-based, and theoretical prescriptions. Most of these models rely on the surface density of neutral hydrogen and the ultraviolet (UV) flux in the Lyman-Werner band as input parameters. Computing these quantities on the kiloparsec scales resolved by the simulations emerges as the main challenge. We show that the commonly used Jeans length approximation to the column density of a system can be biased and exhibits large cell-to-cell scatter. Instead, we propose to compute all surface quantities in face-on projections and perform the modeling in two dimensions. In general, the two methods agree on average, but their predictions diverge for individual galaxies and for models based on the observed midplane pressure of galaxies. We model the UV radiation from young stars by assuming a constant escape fraction and optically thin propagation throughout the galaxy. With these improvements, we find that the five models for the atomic-to-molecular transition roughly agree on average but that the details of the modeling matter for individual galaxies and the spatial distribution of molecular hydrogen. We emphasize that the estimated molecular fractions are approximate due to the significant systematic uncertainties.

We present the discovery of a low-frequency $\approx 5.7$ Hz quasi-periodic oscillation (QPO) feature in observations of the black hole X-ray binary MAXI J1535-571 in its soft-intermediate state, obtained in September-October 2017 by the Neutron Star Interior Composition Explorer (NICER). The feature is relatively broad (compared to other low-frequency QPOs; quality factor $Q\approx 2$) and weak (1.9% rms in 3-10 keV), and is accompanied by a weak harmonic and low-amplitude broadband noise. These characteristics identify it as a weak Type A/B QPO, similar to ones previously identified in the soft-intermediate state of the transient black hole X-ray binary XTE J1550-564. The lag-energy spectrum of the QPO shows increasing soft lags towards lower energies, approaching 50 ms at 1 keV (with respect to a 3-10 keV continuum). This large phase shift has similar amplitude but opposite sign to that seen in Rossi X-ray Timing Explorer data for a Type B QPO from the transient black hole X-ray binary GX 339-4. Previous phase-resolved spectroscopy analysis of the Type B QPO in GX 339-4 pointed towards a precessing jet-like corona illuminating the accretion disk as the origin of the QPO signal. We suggest that this QPO in MAXI J1535-571 may have the same origin, with the different lag sign depending on the scale height of the emitting region and the observer inclination angle.

We introduce the THE THREE HUNDRED project, an endeavour to model 324 large galaxy clusters with full-physics hydrodynamical re-simulations. Here we present the data set and study the differences to observations for fundamental galaxy cluster properties and scaling relations. We find that the modelled galaxy clusters are generally in reasonable agreement with observations with respect to baryonic fractions and gas scaling relations at redshift z = 0. However, there are still some (model-dependent) differences, such as central galaxies being too massive, and galaxy colours (g - r) being bluer (about 0.2 dex lower at the peak position) than in observations. The agreement in gas scaling relations down to 10^{13} h^{-1} M_{\odot} between the simulations indicates that particulars of the sub-grid modelling of the baryonic physics only has a weak influence on these relations. We also include - where appropriate - a comparison to three semi-analytical galaxy formation models as applied to the same underlying dark-matter-only simulation. All simulations and derived data products are publicly available.

We explore the connection between the atomic gas fraction, f_atm, and `global disc stability' parameter, q, of galaxies within a fully cosmological context by examining galaxies in the Dark Sage semi-analytic model. The q parameter is determined by the ratio of disc specific angular momentum to mass. Dark Sage is well suited to our study, as it includes the numerical evolution of one-dimensional disc structure, making both j_disc and q predicted quantities. We show that Dark Sage produces a clear correlation between gas fraction and j_disc at fixed disc mass, in line with recent results from observations and hydrodynamic simulations. This translates to a tight q--f_atm sequence for star-forming central galaxies, which closely tracks the analytic prediction of Obreschkow et al. The scatter in this sequence is driven by the probability distribution function of mass as a function of j (PDF of j) within discs, specifically where it peaks. We find that halo mass is primarily responsible for the peak location of the PDF of j, at least for low values of q. Two main mechanisms of equal significance are then identified for disconnecting f_atm from q. Mergers in the model can trigger quasar winds, with the potential to blow out most of the gas disc, while leaving the stellar disc relatively unharmed. Ram-pressure stripping of satellite galaxies has a similar effect, where f_atm can drop drastically with only a minimal effect to q. We highlight challenges associated with following these predictions up with observations.

We present a comparison of the observed evolving galaxy stellar mass functions with the predictions of eight semi-analytic models and one halo occupation distribution model. While most models are able to fit the data at low redshift, some of them struggle to simultaneously fit observations at high redshift. We separate the galaxies into 'passive' and 'star-forming' classes and find that several of the models produce too many low-mass star-forming galaxies at high redshift compared to observations, in some cases by nearly a factor of 10 in the redshift range $2.5 < z < 3.0$. We also find important differences in the implied mass of the dark matter haloes the galaxies inhabit, by comparing with halo masses inferred from observations. Galaxies at high redshift in the models are in lower mass haloes than suggested by observations, and the star formation efficiency in low-mass haloes is higher than observed. We conclude that many of the models require a physical prescription that acts to dissociate the growth of low-mass galaxies from the growth of their dark matter haloes at high redshift.

We describe the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept that will provide an unprecedented view of the X-ray sky, performing timing and spectroscopy over both a broad energy band (0.2-30 keV) and a wide range of timescales from microseconds to years. STROBE-X comprises two narrow-field instruments and a wide field monitor. The soft or low-energy band (0.2-12 keV) is covered by an array of lightweight optics (3-m focal length) that concentrate incident photons onto small solid-state detectors with CCD-level (85-175 eV) energy resolution, 100 ns time resolution, and low background rates. This technology has been fully developed for NICER and will be scaled up to take advantage of the longer focal length of STROBE-X. The higher-energy band (2-30 keV) is covered by large-area, collimated silicon drift detectors that were developed for the European LOFT mission concept. Each instrument will provide an order of magnitude improvement in effective area over its predecessor (NICER in the soft band and RXTE in the hard band). Finally, STROBE-X offers a sensitive wide-field monitor (WFM), both to act as a trigger for pointed observations of X-ray transients and also to provide high duty-cycle, high time-resolution, and high spectral-resolution monitoring of the variable X-ray sky. The WFM will boast approximately 20 times the sensitivity of the RXTE All-Sky Monitor, enabling multi-wavelength and multi-messenger investigations with a large instantaneous field of view. This mission concept will be presented to the 2020 Decadal Survey for consideration.

We apply the Jeans Anisotropic MGE (JAM) dynamical modelling method to SAGES Legacy Unifying Globulars and GalaxieS (SLUGGS) survey data of early-type galaxies in the stellar mass range $10^{10}<M_*/{\rm M}_{\odot}<10^{11.6}$ that cover a large radial range of $0.1-4.0$ effective radii. We combine SLUGGS and ATLAS$^{\rm 3D}$ datasets to model the total-mass profiles of a sample of 21 fast-rotator galaxies, utilising a hyperparameter method to combine the two independent datasets. The total-mass density profile slope values derived for these galaxies are consistent with those measured in the inner regions of galaxies by other studies. Furthermore, the total-mass density slopes ($\gamma_{\rm tot}$) appear to be universal over this broad stellar mass range, with an average value of $\gamma_{\rm tot}=-2.12\,\pm\,0.05$, i.e. slightly steeper than isothermal. We compare our results to model galaxies from the Magneticum and EAGLE cosmological hydrodynamic simulations, in order to probe the mechanisms that are responsible for varying total-mass density profile slopes. The simulated-galaxy slopes are shallower than the observed values by $\sim0.1-0.3$, indicating that the physical processes shaping the mass distributions of galaxies in cosmological simulations are still incomplete. For galaxies with $M_*>10^{10.7}{\rm M}_{\odot}$ in the Magneticum simulations, we identify a significant anticorrelation between total-mass density profile slopes and the fraction of stellar mass formed ex situ (i.e. accreted), whereas this anticorrelation is weaker for lower stellar masses, implying that the measured total mass density slopes for low-mass galaxies are less likely to be determined by merger activity.

We present a comparison of nine galaxy formation models, eight semi-analytical and one halo occupation distribution model, run on the same underlying cold dark matter simulation (cosmological box of co-moving width 125$h^{-1}$ Mpc, with a dark-matter particle mass of $1.24\times 10^9 h^{-1}$ Msun) and the same merger trees. While their free parameters have been calibrated to the same observational data sets using two approaches, they nevertheless retain some 'memory' of any previous calibration that served as the starting point (especially for the manually-tuned models). For the first calibration, models reproduce the observed z = 0 galaxy stellar mass function (SMF) within 3-{\sigma}. The second calibration extended the observational data to include the z = 2 SMF alongside the z~0 star formation rate function, cold gas mass and the black hole-bulge mass relation. Encapsulating the observed evolution of the SMF from z = 2 to z = 0 is found to be very hard within the context of the physics currently included in the models. We finally use our calibrated models to study the evolution of the stellar-to-halo mass (SHM) ratio. For all models we find that the peak value of the SHM relation decreases with redshift. However, the trends seen for the evolution of the peak position as well as the mean scatter in the SHM relation are rather weak and strongly model dependent. Both the calibration data sets and model results are publicly available.

We present the public release of the MultiDark-Galaxies: three distinct galaxy catalogues derived from one of the Planck cosmology MultiDark simulations (i.e. MDPL2, with a volume of (1 Gpc/$h$)$^{3}$ and mass resolution of $1.5 \times 10^{9} M_{\odot}/h$) by applying the semi-analytic models GALACTICUS, SAG, and SAGE to it. We compare the three models and their conformity with observational data for a selection of fundamental properties of galaxies like stellar mass function, star formation rate, cold gas fractions, and metallicities - noting that they sometimes perform differently reflecting model designs and calibrations. We have further selected galaxy subsamples of the catalogues by number densities in stellar mass, cold gas mass, and star formation rate in order to study the clustering statistics of galaxies. We show that despite different treatment of orphan galaxies, i.e. galaxies that lost their dark-matter host halo due to the finite mass resolution of the N-body simulation or tidal stripping, the clustering signal is comparable, and reproduces the observations in all three models - in particular when selecting samples based upon stellar mass. Our catalogues provide a powerful tool to study galaxy formation within a volume comparable to those probed by on-going and future photometric and redshift surveys. All model data consisting of a range of galaxy properties - including broad-band SDSS magnitudes - are publicly available.

We use EAGLE to quantify the effect galaxy mergers have on the stellar specific angular momentum of galaxies, $j_{\rm stars}$. We split mergers into: dry (gas-poor)/wet (gas-rich), major/minor, and different spin alignments and orbital parameters. Wet (dry) mergers have an average neutral gas-to-stellar mass ratio of $1.1$ ($0.02$), while major (minor) mergers are those with stellar mass ratios $\ge 0.3$ ($0.1-0.3$). We correlate the positions of galaxies in the $j_{\rm stars}$-stellar mass plane at $z=0$ with their merger history, and find that galaxies of low spins suffered dry mergers, while galaxies of normal/high spins suffered predominantly wet mergers, if any. The radial $j_{\rm stars}$ profiles of galaxies that went through dry mergers are deficient by $\approx 0.3$~dex at $r\lesssim 10\,r_{50}$ (with $r_{50}$ being the half-stellar mass radius), compared to galaxies that went through wet mergers. Studying the merger remnants reveals that dry mergers reduce $j_{\rm stars}$ by $\approx 30$\%, while wet mergers increase it by $\approx 10$\%, on average. The latter is connected to the build-up of the bulge by newly formed stars of high rotational speed. Moving from minor to major mergers accentuates these effects. When the spin vectors of the galaxies prior to the dry merger are misaligned, $j_{\rm stars}$ decreases to a greater magnitude, while in wet mergers co-rotation and high orbital angular momentum efficiently spun-up galaxies. We predict what would be the observational signatures in the $j_{\rm stars}$ profiles driven by dry mergers: (i) shallow radial profiles and (ii) profiles that rise beyond $\approx 10\,r_{50}$, both of which are significantly different from spiral galaxies.

We combine the latest spectrally stacked data of 21-cm emission from the ALFALFA survey with an updated version of the Dark Sage semi-analytic model to investigate the relative contributions of secular and environmental astrophysical processes on shaping the HI fractions and quiescence of galaxies in the local Universe. We calibrate the model to match the observed mean HI fraction of all galaxies as a function of stellar mass. Without consideration of stellar feedback, disc instabilities, and active galactic nuclei, we show how the slope and normalisation of this relation would change significantly. We find Dark Sage can reproduce the relative impact that halo mass is observed to have on satellites' HI fractions and quiescent fraction. However, the model satellites are systematically gas-poor. We discuss how this could be affected by satellite--central cross-contamination from the group-finding algorithm applied to the observed galaxies, but that it is not the full story. From our results, we suggest the anti-correlation between satellites' HI fractions and host halo mass, seen at fixed stellar mass and fixed specific star formation rate, can be attributed almost entirely to ram-pressure stripping of cold gas. Meanwhile, stripping of hot gas from around the satellites drives the correlation of quiescent fraction with halo mass at fixed stellar mass. Further detail in the modelling of galaxy discs' centres is required to solidify this result, however. We contextualise our results with those from other semi-analytic models and hydrodynamic simulations.

We use the hydrodynamic, cosmological EAGLE simulations to investigate how hot gas in haloes condenses to form and grow galaxies. We select haloes from the simulations that are actively cooling and study the temperature, distribution, and metallicity of their hot, cold, and transitioning `cooling' gas, placing these in context of semi-analytic models. Our selection criteria lead us to focus on Milky Way-like haloes. We find the hot-gas density profiles of the haloes form a progressively stronger core over time, the nature of which can be captured by a beta profile that has a simple dependence on redshift. In contrast, the hot gas that will cool over a time-step is broadly consistent with a singular isothermal sphere. We find that cooling gas carries a few times the specific angular momentum of the halo and is offset in spin direction from the rest of the hot gas. The gas loses ~60% of its specific angular momentum during the cooling process, generally remaining greater than that of the halo, and it precesses to become aligned with the cold gas already in the disc. We find tentative evidence that angular-momentum losses are slightly larger when gas cools onto dispersion-supported galaxies. We show that an exponential surface density profile for gas arriving on a disc remains a reasonable approximation, but a cusp containing ~20% of the mass is always present, and disc scale radii are larger than predicted by a vanilla Fall & Efstathiou model. These scale radii are still closely correlated with the halo spin parameter, for which we suggest an updated prescription for galaxy formation models.

The equation of state of cold supra-nuclear-density matter, such as in neutron stars, is an open question in astrophysics. A promising method for constraining the neutron star equation of state is modelling pulse profiles of thermonuclear X-ray burst oscillations from hotspots on accreting neutron stars. The pulse profiles, constructed using spherical and oblate neutron star models, are comparable to what would be observed by a next-generation X-ray timing instrument like ASTROSAT, NICER, or LOFT. In this paper, we showcase the use of an evolutionary optimization algorithm to fit pulse profiles to determine the best-fit masses and radii. By fitting synthetic data, we assess how well the optimization algorithm can recover the input parameters. Multiple Poisson realizations of the synthetic pulse profiles, constructed with 1.6 million counts and no background, were fitted with the Ferret algorithm to analyze both statistical and degeneracy-related uncertainty, and to explore how the goodness-of-fit depends on the input parameters. For the regions of parameter space sampled by our tests, the best-determined parameter is the projected velocity of the spot along the observer's line-of-sight, with an accuracy of $\le3$% compared to the true value and with $\le5$\% statistical uncertainty. The next best-determined are the mass and radius; for a neutron star with a spin frequency of 600\,Hz, the best-fit mass and radius are accurate to $\le5$%, with respective uncertainties of $\le7$% and $\le10$%. The accuracy and precision depend on the observer inclination and spot co-latitude, with values of $\sim1$ % achievable in mass and radius if both the inclination and co-latitude are > 60 degrees.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

We use the EAGLE cosmological hydrodynamic simulation suite to study the specific angular momentum of galaxies, $j$, with the aims of (i) investigating the physical causes behind the wide range of $j$ at fixed mass and (ii) examining whether simple, theoretical models can explain the seemingly complex and non-linear nature of the evolution of $j$. We find that $j$ of the stars, $j_{\rm stars}$, and baryons, $j_{\rm bar}$, are strongly correlated with stellar and baryon mass, respectively, with the scatter being highly correlated with morphological proxies such as gas fraction, stellar concentration, (u-r) intrinsic colour, stellar age and the ratio of circular velocity to velocity dispersion. We compare with available observations at $z=0$ and find excellent agreement. We find that $j_{\rm bar}$ follows the theoretical expectation of an isothermal collapsing halo under conservation of specific angular momentum to within $\approx 50$%, while the subsample of rotation-supported galaxies are equally well described by a simple model in which the disk angular momentum is just enough to maintain marginally stable disks. We extracted evolutionary tracks of the stellar spin parameter of EAGLE galaxies and found that the fate of their $j_{\rm stars}$ at $z=0$ depends sensitively on their star formation and merger histories. From these tracks, we identified two distinct physical channels behind low $j_{\rm stars}$ galaxies at $z=0$: (i) galaxy mergers, and (ii) early star formation quenching. The latter can produce galaxies with low $j_{\rm stars}$ and early-type morphologies even in the absence of mergers.

We present a new spectral-timing technique for phase-resolved spectroscopy and apply it to the low-frequency Type B quasi-periodic oscillation (QPO) from the black hole X-ray binary GX 339-4. We show that on the QPO time-scale the spectrum changes not only in normalisation, but also in spectral shape. Using several different spectral models which parameterise the blackbody and power-law components seen in the time-averaged spectrum, we find that both components are required to vary, although the fractional rms amplitude of blackbody emission is small, ~1.4 per cent compared to ~25 per cent for the power-law emission. However the blackbody variation leads the power-law variation by ~0.3 in relative phase (~110 degrees), giving a significant break in the Fourier lag-energy spectrum that our phase-resolved spectral models are able to reproduce. Our results support a geometric interpretation for the QPO variations where the blackbody variation and its phase relation to the power-law are explained by quasi-periodic heating of the approaching and receding sides of the disc by a precessing Comptonising region. The small amplitude of blackbody variations suggests that the Comptonising region producing the QPO has a relatively large scale-height, and may be linked to the base of the jet, as has previously been suggested to explain the binary orbit inclination-dependence of Type B QPO amplitudes.

We present the new semi-analytic model of galaxy evolution, DARK SAGE, a heavily modified version of the publicly available SAGE code. The model is designed for detailed evolution of galactic discs. We evolve discs in a series of annuli with fixed specific angular momentum, which allows us to make predictions for the radial and angular-momentum structure of galaxies. Most physical processes, including all channels of star formation and associated feedback, are performed in these annuli. We present the surface density profiles of our model spiral galaxies, both as a function of radius and specific angular momentum, and find the discs naturally build a pseduobulge-like component. Our main results are focussed on predictions relating to the integrated mass--specific angular momentum relation of stellar discs. The model produces a distinct sequence between these properties in remarkable agreement with recent observational literature. We investigate the impact Toomre disc instabilities have on shaping this sequence and find they are crucial for regulating both the mass and spin of discs. Without instabilities, high-mass discs would be systematically deficient in specific angular momentum by a factor of ~2.5, with increased scatter. Instabilities also appear to drive the direction in which the mass--spin sequence of spiral galaxy discs evolves. With them, we find galaxies of fixed mass have higher specific angular momentum at later epochs.

Fitting the thermal X-ray spectra of neutron stars (NSs) in quiescent X-ray binaries can constrain the masses and radii of NSs. The effect of undetected hot spots on the spectrum, and thus on the inferred NS mass and radius, has not yet been explored for appropriate atmospheres and spectra. A hot spot would harden the observed spectrum, so that spectral modeling tends to infer radii that are too small. However, a hot spot may also produce detectable pulsations. We simulated the effects of a hot spot on the pulsed fraction and spectrum of the quiescent NSs X5 and X7 in the globular cluster 47 Tucanae, using appropriate spectra and beaming for hydrogen atmosphere models, incorporating special and general relativistic effects, and sampling a range of system angles. We searched for pulsations in archival Chandra HRC-S observations of X5 and X7, placing 90\% confidence upper limits on their pulsed fractions below 16\%. We use these pulsation limits to constrain the temperature differential of any hot spots, and to then constrain the effects of possible hot spots on the X-ray spectrum and the inferred radius from spectral fitting. We find that hot spots below our pulsation limit could bias the spectroscopically inferred radius downward by up to 28\%. For Cen X-4 (which has deeper published pulsation searches), an undetected hot spot could bias its inferred radius downward by up to 10\%. Improving constraints on pulsations from quiescent LMXBs may be essential for progress in constraining their radii.

This paper describes a new publicly available codebase for modelling galaxy formation in a cosmological context, the "Semi-Analytic Galaxy Evolution" model, or SAGE for short. SAGE is a significant update to that used in Croton et al. (2006) and has been rebuilt to be modular and customisable. The model will run on any N-body simulation whose trees are organised in a supported format and contain a minimum set of basic halo properties. In this work we present the baryonic prescriptions implemented in SAGE to describe the formation and evolution of galaxies, and their calibration for three N-body simulations: Millennium, Bolshoi, and GiggleZ. Updated physics include: gas accretion, ejection due to feedback, and reincorporation via the galactic fountain; a new gas cooling--radio mode active galactic nucleus (AGN) heating cycle; AGN feedback in the quasar mode; a new treatment of gas in satellite galaxies; and galaxy mergers, disruption, and the build-up of intra-cluster stars. Throughout, we show the results of a common default parameterization on each simulation, with a focus on the local galaxy population.

We address the possibility that intelligent civilisations that destroy themselves could present signatures observable by humanity. Placing limits on the number of self-destroyed civilisations in the Milky Way has strong implications for the final three terms in Drake's Equation, and would allow us to identify which classes of solution to Fermi's Paradox fit with the evidence (or lack thereof). Using the Earth as an example, we consider a variety of scenarios in which humans could extinguish their own technological civilisation. Each scenario presents some form of observable signature that could be probed by astronomical campaigns to detect and characterise extrasolar planetary systems. Some observables are unlikely to be detected at interstellar distances, but some scenarios are likely to produce significant changes in atmospheric composition that could be detected serendipitously with next-generation telescopes. In some cases, the timing of the observation would prove crucial to detection, as the decay of signatures is rapid compared to humanity's communication lifetime. In others, the signatures persist on far longer timescales.

A key component of explaining the array of galaxies observed in the Universe is the feedback of active galactic nuclei, each powered by a massive black hole's accretion disc. For accretion to occur, angular momentum must be lost by that which is accreted. Electromagnetic radiation must offer some respite in this regard, the contribution for which is quantified in this paper, using solely general relativity, under the thin-disc regime. Herein, I calculate extremised situations where photons are entirely responsible for energy removal in the disc and then extend and relate this to the standard relativistic accretion disc outlined by Novikov & Thorne, which includes internal angular-momentum transport. While there is potential for the contribution of angular-momentum removal from photons to be >~1% out to ~10^4 Schwarzschild radii if the disc is irradiated and maximally liberated of angular momentum through inverse Compton scattering, it is more likely of order 10^2 Schwarzschild radii if thermal emission from the disc itself is stronger. The effect of radiation/scattering is stronger near the horizons of fast-spinning black holes, but, ultimately, other mechanisms must drive angular-momentum liberation/transport in accretion discs.

Using the suite of high-resolution zoom re-simulations of individual haloes by Martig et al., and the large-scale simulation \emph{MassiveBlack-II}, we examine the differences in measured galaxy properties from techniques with various aperture definitions of where galaxies end. We perform techniques popular in the literature and present a new technique of our own, where the aperture radius is based on the baryonic mass profiles of simulated (sub)haloes. For the average Milky-Way-mass system, we find the two most popular techniques in the literature return differences of order 30 per cent for stellar mass, a factor of 3 for gas mass, 40 per cent for star formation rate, and factors of several for gas accretion and ejection rates. Individual cases can show variations greater than this, with the severity dependent on the concentration of a given system. The average difference in integrated properties for a more general galaxy population are not as striking, but are still significant for stellar and gas mass, especially for optical-limit apertures. The large differences that can occur are problematic for comparing results from various publications. We stress the importance of both defining and justifying a technique choice and discourage using popular apertures that use an exact fraction of the virial radius, due to the unignorable variation in galaxy-to-(sub)halo size. Finally, we note that technique choice does not greatly affect simulated galaxies from lying within the scatter of observed scaling relations, but it can alter the derived best-fit slope for the Kennicutt-Schmidt relation.

Finite-source effects of gravitationally microlensed stars have been well discussed in the literature, but the role that stellar rotation plays has been neglected. A differential magnification map applied to a differentially Doppler-shifted surface alters the profiles of absorption lines, compromising their ordinarily symmetric nature. Herein, we assess the degree to which this finite-source effect of differential limb magnification (DLM), in combination with stellar rotation, alters spectroscopically derived stellar properties. To achieve this, we simulated a grid of high-magnification microlensing events using synthetic spectra. Our analysis shows that rotation of the source generates differences in the measured equivalent widths of absorption lines supplementary to DLM alone, but only of the order of a few percent. Using the wings of H alpha from the same simulated data, we confirmed the result of Johnson et al. (2010) that DLM alters measurements of effective temperature by < 100 K for dwarf stars, while showing rotation to bear no additional effect.

QSO emission-line spectra are compared to predictions based on theoretical ionizing continua of accretion disks. Observed line intensities do not show the expected trend of higher ionization with higher accretion disk temperature as derived from the black hole mass and accretion rate. This suggests that, at least for accretion rates close to the Eddington limit, the inner disk does not reach temperatures as high as expected from standard disk theory. Modified radial temperature profiles, taking account of winds or advection in the inner disk, achieve better agreement with observation. This conclusion agrees with an earlier study of QSO continuum colors as a function of disk temperature. The emission lines of radio-detected and radio-undetected sources show different trends as a function of disk temperature.

AGN with double-peaked narrow lines (DPAGN) may be caused by kiloparsec scale binary AGN, bipolar outflows, or rotating gaseous disks. We examine the class of DPAGN in which the two narrow line components have closely similar intensity as being especially likely to involve disks or jets. Two spectroscopic indicators support this likelihood. For DPAGN from Smith et al. (2010), the "equal-peaked" objects (EPAGN) have [Ne V]/[O III] ratios lower than for a control sample of non-double peaked AGN. This is unexpected for a pair of normal AGN in a galactic merger, but may be consistent with [O III] emission from a rotating ring with relatively little gas at small radii. Also, [O III]/H-beta ratios of the redshifted and blueshifted systems in the EPAGN are more similar to each other than in a control sample, suggestive of a single ionizing source and inconsistent with the binary interpretation.