We present observations and models of the kinematics and the distribution of the neutral hydrogen (HI) in the isolated dwarf irregular galaxy, Wolf-Lundmark-Melotte (WLM). We observed WLM with the Green Bank Telescope (GBT) and as part of the MeerKAT Early Science Programme, where 16 dishes were available. The HI disc of WLM extends out to a major axis diameter of 30 arcmin (8.5 kpc), and a minor axis diameter of 20 arcmin (5.6 kpc) as measured by the GBT. We use the MeerKAT data to model WLM using the TiRiFiC software suite, allowing us to fit different tilted-ring models and select the one that best matches the observation. Our final best-fitting model is a flat disc with a vertical thickness, a constant inclination and dispersion, and a radially-varying surface brightness with harmonic distortions. To simulate bar-like motions, we include second-order harmonic distortions in velocity in the tangential and the vertical directions. We present a model with only circular motions included and a model with non-circular motions. The latter describes the data better. Overall, the models reproduce the global distribution and the kinematics of the gas, except for some faint emission at the 2-sigma level. We model the mass distribution of WLM with a pseudo-isothermal (ISO) and a Navarro-Frenk-White (NFW) dark matter halo models. The NFW and the ISO models fit the derived rotation curves within the formal errors, but with the ISO model giving better reduced chi-square values. The mass distribution in WLM is dominated by dark matter at all radii.

We investigate the temperature enhancements and formation heights of impulsive heating phenomena in solar active regions such as Ellerman bombs (EBs), ultraviolet bursts (UVBs), and flaring active-region fibrils (FAFs) using interferometric observations in the millimeter (mm) continuum provided by the Atacama Large Millimeter/submillimeter Array (ALMA). We examined 3 mm signatures of heating events identified in Solar Dynamics Observatory (SDO) observations of an active region and compare the results with synthetic spectra from a 3D radiative magnetohydrodynamic simulation. We estimated the contribution from the corona to the mm brightness using differential emission measure analysis. We report the null detection of EBs in the 3 mm continuum at $\sim1.2$" spatial resolution, which is evidence that they are sub-canopy events that do not significantly contribute to heating the upper chromosphere. In contrast, we find the active region to be populated with multiple compact, bright, flickering mm bursts -- reminiscent of UVBs. The high brightness temperatures of up to $\sim14200$ K in some events may have a significant contribution (up to $\sim$7%) from the corona. We also detect FAF-like events in the 3 mm continuum that show rapid motions of $>10000$ K plasma launched with high plane-of-sky velocities ($37-340\rm\,km\,s^{-1}$) from nanoflare kernels. The mm FAFs are the brightest class of warm canopy fibrils that connect magnetic regions of opposite polarity. The simulation confirms that ALMA should be able to detect the mm counterparts of UVBs and small flares and thus provide a complementary diagnostic for impulsive heating in the solar chromosphere.

MeerKATHI is the current development name for a radio-interferometric data reduction pipeline, assembled by an international collaboration. We create a publicly available end-to-end continuum- and line imaging pipeline for MeerKAT and other radio telescopes. We implement advanced techniques that are suitable for producing high-dynamic-range continuum images and spectroscopic data cubes. Using containerization, our pipeline is platform-independent. Furthermore, we are applying a standardized approach for using a number of different of advanced software suites, partly developed within our group. We aim to use distributed computing approaches throughout our pipeline to enable the user to reduce larger data sets like those provided by radio telescopes such as MeerKAT. The pipeline also delivers a set of imaging quality metrics that give the user the opportunity to efficiently assess the data quality.

We present the full panchromatic afterglow light curve data of GW170817, including new radio data as well as archival optical and X-ray data, between 0.5 and 940 days post-merger. By compiling all archival data, and reprocessing a subset of it, we have ensured that the panchromatic dataset is uniform and therefore immune to the differences in data processing or flux determination methods used by different groups. Simple power-law fits to the uniform afterglow light curve indicate a $t^{0.86\pm0.04}$ rise, a $t^{-1.90\pm0.12}$ decline, and a peak occurring at $155\pm4$ days. The afterglow is optically thin throughout its evolution, consistent with a single spectral index ($-0.569\pm0.002$) across all epochs. This gives a precise and updated estimate of the electron power-law index, $p=2.138\pm0.004$. By studying the diffuse X-ray emission from the host galaxy, we place a conservative upper limit on the hot ionized ISM density, $<$0.01 cm$^{-3}$, consistent with previous afterglow studies. Using the late-time afterglow data we rule out any long-lived neutron star remnant having magnetic field strength between 10$^{10.4}$ G and 10$^{16}$ G. Our fits to the afterglow data using an analytical model that includes VLBI proper motion from Mooley et al (2018), and a structured jet model that ignores the proper motion, indicates that the proper motion measurement needs to be considered while seeking an accurate estimate of the viewing angle.

The Murchison Widefield Array (MWA) has observed the entire southern sky (Declination, $\delta <$ 30 deg) at low radio-frequencies, over the range 72-231 MHz. These observations constitute the GaLactic and Extragalactic All-sky MWA (GLEAM) Survey, and we use the extragalactic catalogue (Galactic latitude, $|b| >$ 10 deg) to define the GLEAM 4-Jy (G4Jy) Sample. This is a complete sample of the 'brightest' radio-sources ($S_{\mathrm{151MHz}} >$ 4 Jy), the majority of which are active galactic nuclei with powerful radio-jets. Crucially, low-frequency observations allow the selection of such sources in an orientation-independent way (i.e. minimising the bias caused by Doppler boosting, inherent in high-frequency surveys). We then use higher-resolution radio images, and information at other wavelengths, to morphologically classify the brightest components in GLEAM. We also conduct cross-checks against the literature, and perform internal matching, in order to improve sample completeness (which is estimated to be $>$ 95.5%). This results in a catalogue of 1,863 sources, making the G4Jy Sample over 10 times larger than that of the revised Third Cambridge Catalogue of Radio Sources (3CRR; $S_{\mathrm{178MHz}} >$ 10.9 Jy). Of these G4Jy sources, 78 are resolved by the MWA (Phase-I) synthesised beam ($\sim$2 arcmin at 200 MHz), and we label 67% of the sample as 'single', 26% as 'double', 4% as 'triple', and 3% as having 'complex' morphology at $\sim$1 GHz (45-arcsec resolution). Alongside this, our value-added catalogue provides mid-infrared source associations (subject to 6-arcsec resolution at 3.4 micron) for the radio emission, as identified through visual inspection and thorough checks against the literature. As such, the G4Jy Sample can be used as a reliable training set for cross-identification via machine-learning algorithms. [Abstract abridged for arXiv submission.]

The entire southern sky (Declination, $\delta <$ 30 deg) has been observed using the Murchison Widefield Array (MWA), which provides radio imaging of $\sim$2-arcmin resolution at low frequencies (72-231 MHz). This is the GaLactic and Extragalactic All-sky MWA (GLEAM) Survey, and we have previously used a combination of visual inspection, cross-checks against the literature, and internal matching to identify the 'brightest' radio-sources ($S_{\mathrm{151MHz}} >$ 4 Jy) in the extragalactic catalogue (Galactic latitude, $|b| >$ 10 deg). We refer to these 1,863 sources as the GLEAM 4-Jy (G4Jy) Sample, and use radio images (of $\leq$ 45-arcsec resolution), and multi-wavelength information, to assess their morphology and identify the galaxy that is hosting the radio emission (where appropriate). Details of how to access all of the overlays used for this work are available at https://github.com/svw26/G4Jy. Alongside this we conduct further checks against the literature, which we document in this paper for individual sources. Whilst the vast majority of the G4Jy Sample are active galactic nuclei with powerful radio-jets, we highlight that it also contains a nebula, two nearby, star-forming galaxies, a cluster relic, and a cluster halo. There are also three extended sources for which we are unable to infer the mechanism that gives rise to the low-frequency emission. In the G4Jy catalogue we provide mid-infrared identifications for 86% of the sources, and flag the remainder as: having an uncertain identification (129 sources), having a faint/uncharacterised mid-infrared host (126 sources), or it being inappropriate to specify a host (2 sources). For the subset of 129 sources, there is ambiguity concerning candidate host-galaxies, and this includes four sources (B0424$-$728, B0703$-$451, 3C 198, and 3C 403.1) where we question the existing identification.

At early times dark matter has a thermal velocity dispersion of unknown amplitude which, for warm dark matter models, can influence the formation of nonlinear structure on observable scales. We propose a new scheme to simulate cosmologies with a small-scale suppression of perturbations that combines two previous methods in a way that avoids the numerical artefacts which have so far prevented either from producing fully reliable results. At low densities and throughout most of the cosmological volume, we represent the dark matter phase-sheet directly using high-accuracy interpolation, thereby avoiding the artificial fragmentation which afflicts particle-based methods in this regime. Such phase-sheet methods are, however, unable to follow the rapidly increasing complexity of the denser regions of dark matter haloes, so for these we switch to an N-body scheme which uses the geodesic deviation equation to track phase-sheet properties local to each particle. In addition, we present a novel high-resolution force calculation scheme based on an oct-tree of cubic force resolution elements which is well suited to approximate the force-field of our combined sheet+particle distribution. Our hybrid simulation scheme enables the first reliable simulations of the internal structure of low-mass haloes in a warm dark matter cosmology.

The low-frequency linearly-polarised radio source population is largely unexplored. However, a renaissance in low-frequency polarimetry has been enabled by pathfinder and precursor instruments for the Square Kilometre Array. In this second paper from the POlarised GaLactic and Extragalactic All-Sky Murchison Widefield Array (MWA) Survey -- the POlarised GLEAM Survey, or POGS -- we present the results from our all-sky MWA Phase I Faraday Rotation Measure survey. Our survey covers nearly the entire Southern sky in the Declination range $-82^{\circ}$ to $+30^{\circ}$ at a resolution between around three and seven arcminutes (depending on Declination) using data in the frequency range 169$-$231 MHz. We have performed two targeted searches: the first covering 25,489 square degrees of sky, searching for extragalactic polarised sources; the second covering the entire sky South of Declination $+30^{\circ}$, searching for known pulsars. We detect a total of 517 sources with 200 MHz linearly-polarised flux densities between 9.9 mJy and 1.7 Jy, of which 33 are known radio pulsars. All sources in our catalogues have Faraday rotation measures in the range $-328.07$ rad m$^{-2}$ to $+279.62$ rad m$^{-2}$. The Faraday rotation measures are broadly consistent with results from higher-frequency surveys, but with typically more than an order of magnitude improvement in the precision, highlighting the power of low-frequency polarisation surveys to accurately study Galactic and extragalactic magnetic fields. We discuss the properties of our extragalactic and known-pulsar source population, how the sky distribution relates to Galactic features, and identify a handful of new pulsar candidates among our nominally extragalactic source population.

We investigate whether the galaxy and star formation model used for the Auriga simulations can produce a realistic globular cluster (GC) population. We compare statistics of GC candidate star particles in the Auriga haloes with catalogues of the Milky Way (MW) and Andromeda (M31) GC populations. We find that the Auriga simulations do produce sufficient stellar mass for GC candidates at radii and metallicities that are typical for the MW GC system (GCS). We also find varying mass-ratios of the simulated GC candidates relative to the observed mass in the MW and M31 GC systems for different bins of galactocentric radius-metallicity (r$_{\text{gal}}$ -[Fe/H]). Overall, the Auriga simulations produce GC candidates with higher metallicities than the MW and M31 GCS and they are found at larger radii than observed. The Auriga simulations would require bound cluster formation efficiencies higher than ten percent for the metal-poor GC candidates, and those within the Solar radius should experience negligible destruction rates to be consistent with observations. GC candidates in the outer halo, on the other hand, should either have low formation efficiencies, or experience high mass loss for the Auriga simulations to produce a GCS that is consistent with that of the MW or M31. Finally, the scatter in the metallicity as well as in the radial distribution between different Auriga runs is considerably smaller than the differences between that of the MW and M31 GCSs. The Auriga model is unlikely to give rise to a GCS that can be consistent with both galaxies.

We present MeerKAT 1.28 GHz total-intensity, polarization, and spectral-index images covering the giant (projected length $l \approx 1.57$~Mpc) X-shaped radio source PKS~2014$-$55 with an unprecedented combination of brightness sensitivity and angular resolution. They show the clear "double boomerang" morphology of hydrodynamical backflows from the straight main jets deflected by the large and oblique hot-gas halo of the host galaxy PGC~064440. The magnetic field orientation in PKS~2014$-$55 follows the flow lines from the jets through the secondary wings. The radio source is embedded in faint ($T_\mathrm{b} \approx 0.5 \mathrm{\,K}$) cocoons having the uniform brightness temperature and sharp outer edges characteristic of subsonic expansion into the ambient intra-group medium. The position angle of the much smaller ($l \sim 25$~kpc) restarted central source is within $5^\circ$ of the main jets, ruling out models that invoke jet re-orientation or two independent jets. Compression and turbulence in the backflows probably produce the irregular and low polarization bright region behind the apex of each boomerang as well as several features in the flow with bright heads and dark tails.

We perform a comparison, object-by-object and statistically, between the Munich semi-analytical model, L-Galaxies, and the IllustrisTNG hydrodynamical simulations. By running L-Galaxies on the IllustrisTNG dark matter-only merger trees, we identify the same galaxies in the two models. This allows us to compare the stellar mass, star formation rate and gas content of galaxies, as well as the baryonic content of subhaloes and haloes in the two models. We find that both the stellar mass functions and the stellar masses of individual galaxies agree to better than $\sim0.2\,$dex. On the other hand, specific star formation rates and gas contents can differ more substantially. At $z=0$ the transition between low-mass star-forming galaxies and high-mass, quenched galaxies occurs at a stellar mass scale $\sim0.5\,$dex lower in IllustrisTNG than in L-Galaxies. IllustrisTNG also produces substantially more quenched galaxies at higher redshifts. Both models predict a halo baryon fraction close to the cosmic value for clusters, but IllustrisTNG predicts lower baryon fractions in group environments. These differences are due primarily to differences in modelling feedback from stars and supermassive black holes. The gas content and star formation rates of galaxies in and around clusters and groups differ substantially, with IllustrisTNG satellites less star-forming and less gas-rich. We show that environmental processes such as ram-pressure stripping are stronger and operate to larger distances and for a broader host mass range in IllustrisTNG. We suggest that the treatment of galaxy evolution in the semi-analytic model needs to be improved by prescriptions which capture local environmental effects more accurately.

We present X-ray imaging spectroscopy of one of the weakest active region (AR) microflares ever studied. The microflare occurred at $\sim$11:04 UT on 2018 September 9 and we studied it using the Nuclear Spectroscopic Telescope ARray (NuSTAR) and the Solar Dynamic Observatory's Atmospheric Imaging Assembly (SDO/AIA). The microflare is observed clearly in 2.5-7 keV with NuSTAR and in Fe XVIII emission derived from the hotter component of the 94 $\unicode{x212B}$ SDO/AIA channel. We estimate the event to be three orders of magnitude lower than a GOES A class microflare with an energy of 1.1$\times$10$^{26}$ erg. It reaches temperatures of 6.7 MK with an emission measure of 8.0$\times$10$^{43}$ cm$^{-3}$. Non-thermal emission is not detected but we instead determine upper limits to such emission. We present the lowest thermal energy estimate for an AR microflare in literature, which is at the lower limits of what is still considered an X-ray microflare.

We explore the chemodynamical properties of a sample of barred galaxies in the Auriga magneto-hydrodynamical cosmological zoom-in simulations, which form boxy/peanut (b/p) bulges, and compare these to the Milky Way (MW). We show that the Auriga galaxies which best reproduce the chemodynamical properties of stellar populations in the MW bulge have quiescent merger histories since redshift $z\sim3.5$: their last major merger occurs at $t_{\rm lookback}>12\,\rm Gyrs$, while subsequent mergers have a stellar mass ratio of $\leq$1:20, suggesting an upper limit of a few percent for the mass ratio of the recently proposed Gaia Sausage/Enceladus merger. These Auriga MW-analogues have a negligible fraction of ex-situ stars in the b/p region ($<1\%$), with flattened, thick disc-like metal-poor stellar populations. The average fraction of ex-situ stars in the central regions of all Auriga galaxies with b/p's is 3% -- significantly lower than in those which do not host a b/p or a bar. While the central regions of these barred galaxies contain the oldest populations, they also have stars younger than 5Gyrs (>30%) and exhibit X-shaped age and abundance distributions. Examining the discs in our sample, we find that in some cases a star-forming ring forms around the bar, which alters the metallicity of the inner regions of the galaxy. Further out in the disc, bar-induced resonances lead to metal-rich ridges in the $V_{\phi}-r$ plane -- the longest of which is due to the Outer Lindblad Resonance. Our results suggest the Milky Way has an uncommonly quiet merger history, which leads to an essentially in-situ bulge, and highlight the significant effects the bar can have on the surrounding disc.

We present MeerKAT 1000 MHz and 1400 MHz observations of a bright radio galaxy in the southern hemisphere, ESO~137-006. The galaxy lies at the centre of the massive and merging Norma galaxy cluster. The MeerKAT continuum images (rms ~0.02 mJy/beam at ~10" resolution) reveal new features that have never been seen in a radio galaxy before: collimated synchrotron threads of yet unknown origin, which link the extended and bent radio lobes of ESO~137-006. The most prominent of these threads stretches in projection for about 80 kpc and is about 1 kpc in width. The radio spectrum of the threads is steep, with a spectral index of up to $\alpha\simeq 2$ between 1000 MHz and 1400 MHz.

We present the data analysis pipeline, commissioning observations and initial results from the GREENBURST fast radio burst (FRB) detection system on the Robert C. Byrd Green Bank Telescope (GBT) previously described by Surnis et al. which uses the 21~cm receiver observing commensally with other projects. The pipeline makes use of a state-of-the-art deep learning classifier to winnow down the very large number of false positive single-pulse candidates that mostly result from radio frequency interference. In our observations totalling 156.5 days so far, we have detected individual pulses from 20 known radio pulsars which provide an excellent verification of the system performance. We also demonstrate, through blind injection analyses, that our pipeline is complete down to a signal-to-noise threshold of 12. Depending on the observing mode, this translates to peak flux sensitivities in the range 0.14--0.89~Jy. Although no FRBs have been detected to date, we have used our results to update the analysis of Lawrence et al. to constrain the FRB all-sky rate to be $1140^{+200}_{-180}$ per day above a peak flux density of 1~Jy. We also constrain the source count index $\alpha=0.83\pm0.06$ which indicates that the source count distribution is substantially flatter than expected from a Euclidean distribution of standard candles (where $\alpha=1.5$). We discuss this result in the context of the FRB redshift and luminosity distributions. Finally, we make predictions for detection rates with GREENBURST, as well as other ongoing and planned FRB experiments.

We report the detection of emission from a non-thermal electron distribution in a small solar microflare (GOES class A5.7) observed by the Nuclear Spectroscopic Telescope Array (NuSTAR), with supporting observation by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). The flaring plasma is well accounted for by a thick-target model of accelerated electrons collisionally thermalizing within the loop, akin to the "coronal thick target" behavior occasionally observed in larger flares. This is the first positive detection of non-thermal hard X-rays from the Sun using a direct imager (as opposed to indirectly imaging instruments). The accelerated electron distribution has a spectral index of 6.3 +/- 0.7, extends down to at least 6.5 keV, and deposits energy at a rate of ~2x1027 erg/s, heating the flare loop to at least 10 MK. The existence of dominant non-thermal emission in X-rays down to <5 keV means that RHESSI emission is almost entirely non-thermal, contrary to what is usually assumed in RHESSI spectroscopy. The ratio of non-thermal to thermal energies is similar to that of large flares, in contrast to what has been found in previous studies of small RHESSI flares. We suggest that a coronal thick target may be a common property of many small microflares based on the average electron energy and collisional mean free path. Future observations of this kind will enable understanding of how flare particle acceleration changes across energy scales, and will aid the push toward the observational regime of nanoflares, which are a possible source of significant coronal heating.

We present anisotropic ``separate universe'' simulations which modify the N-body code Gadget-4 in order to represent a large-scale tidal field through an anisotropic expansion factor. These simulations are used to measure the linear, quasi-linear and nonlinear response of the matter power spectrum to a spatially uniform trace-free tidal field up to wavenumber $k = 10 h \text{Mpc}^{-1}$. Together with the response to a large-scale overdensity measured in previous work, this completely describes the nonlinear matter bispectrum in the squeezed limit. We find that the response amplitude does not approach zero on small scales in physical coordinates, but rather a constant value, $R_K\approx 0.5$ for $k \geq 3 h\text{Mpc}^{-1}$ at $z=0$. This shows that even the inner regions of haloes are affected by the large-scale tidal field. We also measure directly the alignment of halo shapes with the tidal field, finding a clear signal which increases with halo mass.

We have updated the Munich galaxy formation model, L-Galaxies, to follow the radial distributions of stars and atomic and molecular gas in galaxy discs. We include an H2-based star-formation law, as well as a detailed chemical-enrichment model with explicit mass-dependent delay times for SN-II, SN-Ia and AGB stars. Information about the star formation, feedback and chemical-enrichment histories of discs is stored in 12 concentric rings. The new model retains the success of its predecessor in reproducing the observed evolution of the galaxy population, in particular, stellar mass functions and passive fractions over the redshift range 0<=z<=3 and mass range 8<=log(M_*/Msun)<=12, the black hole-bulge mass relation at z=0, galaxy morphology as a function of stellar mass and the mass-metallicity relations of both stellar and gas components. In addition, its detailed modelling of the radial structure of discs allows qualitatively new comparisons with observation, most notably with the relative sizes and masses of the stellar, atomic and molecular components in discs. Good agreement is found with recent data. Comparison of results obtained for simulations differing in mass resolution by more than two orders of magnitude shows that all important distributions are numerically well converged even for this more detailed model. An examination of metallicity and surface-density gradients in the stars and gas indicates that our new model, with star formation, chemical enrichment and feedback calculated self-consistently on local disc scales, reproduces some but not all of the trends seen in recent many-galaxy IFU surveys.

We use the Millennium and Millennium-II simulations to illustrate the Tessellation-Level-Tree (TLT), a hierarchical tree structure linking density peaks in a field constructed by voronoi tessellation of the particles in a cosmological N-body simulation. The TLT uniquely partitions the simulation particles into disjoint subsets, each associated with a local density peak. Each peak is a subpeak of a unique higher peak. The TLT can be persistence filtered to suppress peaks produced by discreteness noise. Thresholding a peak's particle list at $\sim 80\left<\rho\right>$ results in a structure similar to a standard friend-of-friends halo and its subhaloes. For thresholds below $\sim 7\left<\rho\right>$, the largest structure percolates and is much more massive than other objects. It may be considered as defining the cosmic web. For a threshold of $5\left<\rho\right>$, it contains about half of all cosmic mass and occupies $\sim 1\%$ of all cosmic volume; a typical external point is then $\sim 7h^{-1}\mathrm{Mpc}$ from the web. We investigate the internal structure and clustering of TLT peaks. Defining the saddle point density $\rho_{\mathrm{lim}}$ as the density at which a peak joins its parent peak, we show the median value of $\rho_{\mathrm{lim}}$ for FoF-like peaks to be similar to the density threshold at percolation. Assembly bias as a function of $\rho_{\mathrm{lim}}$ is stronger than for any known internal halo property. For peaks of group mass and below, the lowest quintile in $\rho_{\mathrm{lim}}$ has $b\approx 0$, and is thus uncorrelated with the mass distribution.

We present a detailed analysis of the radio galaxy PKS 2250-351, a giant of 1.2 Mpc projected size, its host galaxy, and its environment. We use radio data from the Murchison Widefield Array, the upgraded Giant Metre-wavelength Radio Telescope, the Australian Square Kilometre Array Pathfinder, and the Australia Telescope Compact Array to model the jet power and age. Optical and infra-red data come from the Galaxy And Mass Assembly (GAMA) survey and provide information on the host galaxy and environment. GAMA spectroscopy confirms that PKS 2250-351 lies at z=0.2115 in the irregular, and likely unrelaxed, cluster Abell 3936. We find its host is a massive, `red and dead' elliptical galaxy with negligible star formation but with a highly obscured active galactic nucleus dominating the mid-infrared emission. Assuming it lies on the local M-sigma relation it has an Eddington accretion rate of lambda_EDD~0.014. We find that the lobe-derived jet power (a time-averaged measure) is an order of magnitude greater than the hotspot-derived jet power (an instantaneous measure). We propose that over the lifetime of the observed radio emission (~300 Myr) the accretion has switched from an inefficient advection dominated mode to a thin-disc efficient mode, consistent with the decrease in jet power. We also suggest that the asymmetric radio morphology is due to its environment, with the host of PKS 2250-351 lying to the west of the densest concentration of galaxies in Abell 3936.

We present the radio luminosity function (RLF) of optically-selected quasars below 1~mJy, constructed by applying a Bayesian-fitting stacking technique to objects well below the nominal radio flux-density limit. We test the technique using simulated data, confirming that we can reconstruct the RLF over three orders of magnitude below the typical $5\sigma$ detection threshold. We apply our method to 1.4-GHz flux-densities from the Faint Images of the Radio Sky at Twenty-cm survey (FIRST), extracted at the positions of optical quasars from the Sloan Digital Sky Survey (SDSS) over seven redshift bins up to $z=2.15$ {and measure the RLF down to two orders of magnitude below the FIRST detection threshold}. In the lowest redshift bin ($0.2<z<0.45$), we find that our measured RLF agrees well with deeper data from the literature. The RLF for the radio-loud quasars flattens below $\log_{10}[L_{1.4}/{\rm W\,Hz}^{-1}] \approx 25.5$ and becomes steeper again below $\log_{10}[L_{1.4}/{\rm W\,Hz}^{-1}] \approx 24.8$, where radio-quiet quasars start to emerge. The radio luminosity where radio-quiet quasars emerge coincides with the luminosity where star-forming galaxies are expected to start to dominate the radio source counts. This implies that there could be a significant contribution from star formation in the host galaxies, but additional data is required to investigate this further. The higher-redshift bins show a similar behaviour as for the lowest-$z$ bin, implying that the same physical process may be responsible.

Radio observations allow us to identify a wide range of active galactic nuclei (AGN), which play a significant role in the evolution of galaxies. Amongst AGN at low radio-luminosities is the 'radio-quiet' quasar (RQQ) population, but how they contribute to the total radio emission is under debate, with previous studies arguing that it is predominantly through star formation. In this talk, SVW summarised the results of recent papers on RQQs, including the use of far-infrared data to disentangle the radio emission from the AGN and that from star formation. This provides evidence that black-hole accretion, instead, dominates the radio emission in RQQs. In addition, we find that this accretion-related emission is correlated with the optical luminosity of the quasar, whilst a weaker luminosity-dependence is evident for the radio emission connected with star formation. What remains unclear is the process by which this accretion-related emission is produced. Understanding this for RQQs will then allow us to investigate how this type of AGN influences its surroundings. Such studies have important implications for modelling AGN feedback, and for determining the accretion and star-formation histories of the Universe.

The European Space Agency's Planck satellite, which was dedicated to studying the early Universe and its subsequent evolution, was launched on 14 May 2009. It scanned the microwave and submillimetre sky continuously between 12 August 2009 and 23 October 2013, producing deep, high-resolution, all-sky maps in nine frequency bands from 30 to 857GHz. This paper presents the cosmological legacy of Planck, which currently provides our strongest constraints on the parameters of the standard cosmological model and some of the tightest limits available on deviations from that model. The 6-parameter LCDM model continues to provide an excellent fit to the cosmic microwave background data at high and low redshift, describing the cosmological information in over a billion map pixels with just six parameters. With 18 peaks in the temperature and polarization angular power spectra constrained well, Planck measures five of the six parameters to better than 1% (simultaneously), with the best-determined parameter (theta_*) now known to 0.03%. We describe the multi-component sky as seen by Planck, the success of the LCDM model, and the connection to lower-redshift probes of structure formation. We also give a comprehensive summary of the major changes introduced in this 2018 release. The Planck data, alone and in combination with other probes, provide stringent constraints on our models of the early Universe and the large-scale structure within which all astrophysical objects form and evolve. We discuss some lessons learned from the Planck mission, and highlight areas ripe for further experimental advances.

Dark matter haloes are the basic units of all cosmic structure. They grew by gravitational amplification of weak initial density fluctuations that are still visible on large scales in the cosmic microwave background radiation. Galaxies formed within relatively massive haloes as gas cooled and condensed at their centres, but many hypotheses for the nature of dark matter imply that the halo population should extend to masses many orders of magnitude below those where galaxies can form. Here, we use a novel, multi-zoom technique to create the first consistent simulation of the formation of present-day haloes over the full mass range populated when dark matter is aWeakly Interacting Massive Particle (WIMP) of mass ~100 GeV. The simulation has a dynamic range of 30 orders of magnitude in mass, resolving the internal structure of hundreds of Earth-mass haloes just as well as that of hundreds of rich galaxy clusters. Remarkably, halo density profiles are universal over the entire mass range and are well described by simple two-parameter fitting formulae. Halo mass and concentration are tightly related in a way which depends on cosmology and on the nature of the dark matter. At fixed mass, concentration is independent of local environment for haloes less massive than those of typical galaxies. These results are important for predicting annihilation radiation signals from dark matter, since these should be dominated by contributions from the smallest structures.

Using four years of full-disk-integrated coronal differential emission measures calculated in Schonfeld et al. (2017), we investigate the relative contribution of bremsstrahlung and gyroresonance emission in observations of F$_{10.7}$, the 10.7 cm (2.8 GHz) solar microwave spectral flux density and commonly used activity proxy. We determine that the majority of coronal F$_{10.7}$ is produced by the bremsstrahlung mechanism, but the variability observed over individual solar rotations is often driven by gyroresonance sources rotating across the disk. Our analysis suggests that the chromosphere may contribute significantly to F$_{10.7}$ variability and that coronal bremsstrahlung emission accounts for 14.2 $\pm$ 2.1 sfu ($\sim20\%$) of the observed solar minimum level. The bremsstrahlung emission has a power-law relationship to the total F$_{10.7}$ at high activity levels, and this combined with the observed linearity during low activity yields a continuously differentiable piecewise fit for the bremsstrahlung component as a function of F$_{10.7}$. We find that the bremsstrahlung component fit, along with the Mg II index, correlates better with the observed 5-37 nm spectrum than the common 81 day averaged F$_{10.7}$ proxy. The bremsstrahlung component of F$_{10.7}$ is also well approximated by the moderate-strength photospheric magnetic field parameterization from Henney et al. (2012), suggesting that it could be forecast for use in both atmospheric research and operational models.

We present constraints on the formation and evolution of early-type galaxies (ETGs) with the empirical model EMERGE. The parameters of this model are adjusted so that it reproduces the evolution of stellar mass functions, specific star formation rates, and cosmic star formation rates since $z\approx10$ as well as 'quenched' galaxy fractions and correlation functions. We find that at fixed halo mass present-day ETGs are more massive than late-type galaxies, whereas at fixed stellar mass ETGs populate more massive halos in agreement with lensing results. This effect naturally results from the shape and scatter of the stellar-to-halo mass relation and the galaxy formation histories. The ETG stellar mass assembly is dominated by 'in-situ' star formation below a stellar mass of $3\times10^{11}\mathrm{M}_\odot$ and by merging and accretion of 'ex-situ' formed stars at higher mass. The mass dependence is in tension with current cosmological simulations. Lower mass ETGs show extended star formation towards low redshift in agreement with recent estimates from IFU surveys. All ETGs have main progenitors on the 'main sequence of star formation' with the 'red sequence' appearing at $z \approx 2$. Above this redshift, over 95 per cent of the ETG progenitors are star-forming. More than 90 per cent of $z \approx 2$ 'main sequence' galaxies with $m_* > 10^{10}\mathrm{M}_\odot$ evolve into present-day ETGs. Above redshift 6, more than 80 per cent of the observed stellar mass functions above $10^{9}\mathrm{M}_\odot$ can be accounted for by ETG progenitors with $m_* > 10^{10}\mathrm{M}_\odot$. This implies that current and future high redshift observations mainly probe the birth of present-day ETGs. The source code and documentation of EMERGE are available at github.com/bmoster/emerge.

Using a set of 15 high-resolution magnetohydrodynamic cosmological simulations of Milky Way formation, we investigate the origin of the baryonic material found in stars at redshift zero. We find that roughly half of this material originates from subhalo/satellite systems and half is smoothly accreted from the Inter-Galactic Medium (IGM). About $90 \%$ of all material has been ejected and re-accreted in galactic winds at least once. The vast majority of smoothly accreted gas enters into a galactic fountain that extends to a median galactocentric distance of $\sim 20$ kpc with a median recycling timescale of $\sim 500$ Myr. We demonstrate that, in most cases, galactic fountains acquire angular momentum via mixing of low-angular momentum, wind-recycled gas with high-angular momentum gas in the Circum-Galactic Medium (CGM). Prograde mergers boost this activity by helping to align the disc and CGM rotation axes, whereas retrograde mergers cause the fountain to lose angular momentum. Fountain flows that promote angular momentum growth are conducive to smooth evolution on tracks quasi-parallel to the disc sequence of the stellar mass-specific angular momentum plane, whereas retrograde minor mergers, major mergers and bar-driven secular evolution move galaxies towards the bulge-sequence. Finally, we demonstrate that fountain flows act to flatten and narrow the radial metallicity gradient and metallicity dispersion of disc stars, respectively. Thus, the evolution of galactic fountains depends strongly on the cosmological merger history and is crucial for the chemo-dynamical evolution of Milky Way-sized disc galaxies.

We present cosmological parameter results from the final full-mission Planck measurements of the CMB anisotropies. We find good consistency with the standard spatially-flat 6-parameter $\Lambda$CDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted "base $\Lambda$CDM" in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis gives dark matter density $\Omega_c h^2 = 0.120\pm 0.001$, baryon density $\Omega_b h^2 = 0.0224\pm 0.0001$, scalar spectral index $n_s = 0.965\pm 0.004$, and optical depth $\tau = 0.054\pm 0.007$ (in this abstract we quote $68\,\%$ confidence regions on measured parameters and $95\,\%$ on upper limits). The angular acoustic scale is measured to $0.03\,\%$ precision, with $100\theta_*=1.0411\pm 0.0003$. These results are only weakly dependent on the cosmological model and remain stable, with somewhat increased errors, in many commonly considered extensions. Assuming the base-$\Lambda$CDM cosmology, the inferred late-Universe parameters are: Hubble constant $H_0 = (67.4\pm 0.5)$km/s/Mpc; matter density parameter $\Omega_m = 0.315\pm 0.007$; and matter fluctuation amplitude $\sigma_8 = 0.811\pm 0.006$. We find no compelling evidence for extensions to the base-$\Lambda$CDM model. Combining with BAO we constrain the effective extra relativistic degrees of freedom to be $N_{\rm eff} = 2.99\pm 0.17$, and the neutrino mass is tightly constrained to $\sum m_\nu< 0.12$eV. The CMB spectra continue to prefer higher lensing amplitudes than predicted in base -$\Lambda$CDM at over $2\,\sigma$, which pulls some parameters that affect the lensing amplitude away from the base-$\Lambda$CDM model; however, this is not supported by the lensing reconstruction or (in models that also change the background geometry) BAO data. (Abridged)

We analyse the orbital kinematics of the Milky Way (MW) satellite system utilizing the latest systemic proper motions for 38 satellites based on data from Gaia Data Release 2. Combining these data with distance and line-of-sight velocity measurements from the literature, we use a likelihood method to model the velocity anisotropy, $\beta$, as a function of Galactocentric distance and compare the MW satellite system with those of simulated MW-mass haloes from the APOSTLE and Auriga simulation suites. The anisotropy profile for the MW satellite system increases from $\beta\sim -2$ at $r\sim20$ kpc to $\beta\sim 0.5$ at $r\sim200$ kpc, indicating that satellites closer to the Galactic centre have tangentially-biased motions while those farther out have radially-biased motions. The motions of satellites around APOSTLE host galaxies are nearly isotropic at all radii, while the $\beta(r)$ profiles for satellite systems in the Auriga suite, whose host galaxies are substantially more massive in baryons than those in APOSTLE, are more consistent with that of the MW satellite system. This shape of the $\beta(r)$ profile may be attributed to the central stellar disc preferentially destroying satellites on radial orbits, or intrinsic processes from the formation of the Milky Way system.

We report on the implications for cosmic inflation of the 2018 Release of the Planck CMB anisotropy measurements. The results are fully consistent with the two previous Planck cosmological releases, but have smaller uncertainties thanks to improvements in the characterization of polarization at low and high multipoles. Planck temperature, polarization, and lensing data determine the spectral index of scalar perturbations to be $n_\mathrm{s}=0.9649\pm 0.0042$ at 68% CL and show no evidence for a scale dependence of $n_\mathrm{s}.$ Spatial flatness is confirmed at a precision of 0.4% at 95% CL with the combination with BAO data. The Planck 95% CL upper limit on the tensor-to-scalar ratio, $r_{0.002}<0.10$, is further tightened by combining with the BICEP2/Keck Array BK15 data to obtain $r_{0.002}<0.056$. In the framework of single-field inflationary models with Einstein gravity, these results imply that: (a) slow-roll models with a concave potential, $V" (\phi) < 0,$ are increasingly favoured by the data; and (b) two different methods for reconstructing the inflaton potential find no evidence for dynamics beyond slow roll. Non-parametric reconstructions of the primordial power spectrum consistently confirm a pure power law. A complementary analysis also finds no evidence for theoretically motivated parameterized features in the Planck power spectrum, a result further strengthened for certain oscillatory models by a new combined analysis that includes Planck bispectrum data. The new Planck polarization data provide a stringent test of the adiabaticity of the initial conditions. The polarization data also provide improved constraints on inflationary models that predict a small statistically anisotropic quadrupolar modulation of the primordial fluctuations. However, the polarization data do not confirm physical models for a scale-dependent dipolar modulation.

We present measurements of the cosmic microwave background (CMB) lensing potential using the final $\textit{Planck}$ 2018 temperature and polarization data. We increase the significance of the detection of lensing in the polarization maps from $5\,\sigma$ to $9\,\sigma$. Combined with temperature, lensing is detected at $40\,\sigma$. We present an extensive set of tests of the robustness of the lensing-potential power spectrum, and construct a minimum-variance estimator likelihood over lensing multipoles $8 \le L \le 400$. We find good consistency between lensing constraints and the results from the $\textit{Planck}$ CMB power spectra within the $\rm{\Lambda CDM}$ model. Combined with baryon density and other weak priors, the lensing analysis alone constrains $\sigma_8 \Omega_{\rm m}^{0.25}=0.589\pm 0.020$ ($1\,\sigma$ errors). Also combining with baryon acoustic oscillation (BAO) data, we find tight individual parameter constraints, $\sigma_8=0.811\pm0.019$, $H_0=67.9_{-1.3}^{+1.2}\,\text{km}\,\text{s}^{-1}\,\rm{Mpc}^{-1}$, and $\Omega_{\rm m}=0.303^{+0.016}_{-0.018}$. Combining with $\textit{Planck}$ CMB power spectrum data, we measure $\sigma_8$ to better than $1\,\%$ precision, finding $\sigma_8=0.811\pm 0.006$. We find consistency with the lensing results from the Dark Energy Survey, and give combined lensing-only parameter constraints that are tighter than joint results using galaxy clustering. Using $\textit{Planck}$ cosmic infrared background (CIB) maps we make a combined estimate of the lensing potential over $60\,\%$ of the sky with considerably more small-scale signal. We demonstrate delensing of the $\textit{Planck}$ power spectra, detecting a maximum removal of $40\,\%$ of the lensing-induced power in all spectra. The improvement in the sharpening of the acoustic peaks by including both CIB and the quadratic lensing reconstruction is detected at high significance (abridged).

We present spectropolarimetric imaging observations of the solar corona at low frequencies (80 - 240 MHz) using the Murchison Widefield Array (MWA). These images are the first of their kind, and we introduce an algorithm to mitigate an instrumental artefact by which the total intensity signal contaminates the polarimetric images due to calibration errors. We then survey the range of circular polarization (Stokes V) features detected in over 100 observing runs near solar maximum during quiescent periods. First, we detect around 700 compact polarized sources across our dataset with polarization fractions ranging from less than 0.5% to nearly 100%. These sources exhibit a positive correlation between polarization fraction and total intensity, and we interpret them as a continuum of plasma emission noise storm (Type I burst) continua sources associated with active regions. Second, we report a characteristic "bullseye" structure observed for many low-latitude coronal holes in which a central polarized component is surrounded by a ring of the opposite sense. The central component does not match the sign expected from thermal bremsstrahlung emission, and we speculate that propagation effects or an alternative emission mechanism may be responsible. Third, we show that the large-scale polarimetric structure at our lowest frequencies is reasonably well-correlated with the line-of-sight (LOS) magnetic field component inferred from a global potential field source surface (PFSS) model. The boundaries between opposite circular polarization signs are generally aligned with polarity inversion lines in the model at a height roughly corresponding to that of the radio limb. This is not true at our highest frequencies, however, where the LOS magnetic field direction and polarization sign are often not straightforwardly correlated.

We explore the origin of phase-space substructures revealed by the second Gaia data release in the disc of the Milky Way, such as the ridges in the $V_{\phi}$-$r$ plane, the undulations in the $V_{\phi}$-$r$-$V_r$ space and the streams in the $V_{\phi}$-$V_r$ plane. We use a collisionless N-body simulation with co-spatial thin and thick discs, along with orbit integration, to study the orbital structure close to the Outer Lindblad Resonance (OLR) of the bar. We find that a prominent, long-lived ridge is formed in the $V_{\phi}$-$r$ plane due to the OLR which translates to streams in the $V_{\phi}$-$V_r$ plane and examine which closed periodic and trapped librating orbits are responsible for these features. We find that orbits which carry out small librations around the $x_1(1)$ family are preferentially found at negative $V_r$, giving rise to a `horn'-like feature, while orbits with larger libration amplitudes, trapped around the $x_1(2)$ and $x_1(1)$ families, constitute the positive $V_r$ substructure, i.e. the Hercules-like feature. This changing libration amplitude of orbits will translate to a changing ratio of thin/thick disc stars, which could have implications on the metallicity distribution in this plane. We find that a scenario in which the Sun is placed close to the OLR gives rise to a strong asymmetry in $V_r$ in the $V_{\phi}$-$V_r$ plane (i.e. Hercules vs. `the horn') and subsequently to undulations in the $V_{\phi}$-$r$-$V_r$ space. We also explore a scenario in which the Sun is placed closer to the bar corotation and find that the bar perturbation $alone$ cannot give rise to the these features.

We study the galactic bulges in the Auriga simulations, a suite of thirty cosmological magneto-hydrodynamical zoom-in simulations of late-type galaxies in Milky Way-sized dark matter haloes performed with the moving-mesh code AREPO. We aim to characterize bulge formation mechanisms in this large suite of galaxies simulated at high resolution in a fully cosmological context. The bulges of the Auriga galaxies show a large variety in their shapes,sizes and formation histories. According to observational classification criteria, such as Sersic index and degree of ordered rotation, the majority of the Auriga bulges can be classified as pseudo-bulges, while some of them can be seen as composite bulges with a classical component; however, none can be classified as a classical bulge. Auriga bulges show mostly an in-situ origin, 21 percent of them with a negligible accreted fraction (facc < 0.01). In general,their in-situ component was centrally formed, with 75 percent of the bulges forming most of their stars inside the bulge region at z=0. Part of their in-situ mass growth is rapid and is associated with the effects of mergers, while another part is more secular in origin. In 90 percent of the Auriga bulges, the accreted bulge component originates from less than four satellites.We investigate the relation between the accreted stellar haloes and the bulges of the Auriga simulations. The total bulge mass shows no correlation with the accreted stellar halo mass, as in observations. However, the accreted mass of bulges tends to correlate with their respective accreted stellar halo mass.

We present the first joint observation of a small microflare in X-rays with the Nuclear Spectroscopic Telescope ARray (NuSTAR), UV with the Interface Region Imaging Spectrograph (IRIS) and EUV with the Solar Dynamics Observatory/Atmospheric Imaging Assembly (SDO/AIA). These combined observations allows us to study the microflare's hot coronal and cooler chromospheric/transition region emission. This small microflare peaks from SOL2016-07-26T23:35 to 23:36UT, in both NuSTAR, SDO/AIA and IRIS. Spatially this corresponds to a small loop visible in the SDO/AIA Fe XVIII emission, which matches a similar structure lower in the solar atmosphere seen by IRIS in SJI1330{\AA} and 1400\AA. The NuSTAR emission in both 2.5-4 keV and 4-6 keV, is located in a small source at this loop location. The IRIS slit was over the microflaring loop, and fits show little change in Mg II but do show intensity increases, slight width enhancements and redshifts in Si IV andO IV, indicating that this microflare had most significance in and above the upper chromosphere. The NuSTAR microflare spectrum is well fitted by a thermal component of 5.8MK and $2.3\times10^{44}$ cm$^{-3}$, which corresponds to a thermal energy of $10^{26}$ erg, making it considerably smaller than previously studied X-ray microflares. No non-thermal emission was detected but this could be due to the limited effective exposure time of the observation. This observation shows that even ordinary features seen in UV and EUV, can remarkably have a higher energy component that is clear in X-rays.

Action space synthesizes the orbital information of stars and is well-suited to analyse the rich amount of kinematic disk substructure in the Gaia DR2 radial velocity sample (RVS). In this work, we revisit one of the strongest perturbers in the Milky Way (MW) disk: the m=2 bar. We investigate how its resonances affect the actions of individual test particle stars, i.e., (JR,Lz,Jz) estimated in an axisymmetric MW potential. We confirm that the stars' behaviour is well approximated by scattering and oscillation along a slope Delta JR / Delta Lz = l/m centered on the l:m resonance lines. The Outer Lindblad Resonance (OLR, l=+1,m=2) creates signatures in the stellar action space that can be used to identify the Galactic bar's OLR in the Gaia DR2 RVS data: (a) The JR dependence of the oscillation causes an overdensity ridge (underdensity region) at Lz larger (smaller) than the resonance line in the (Lz,JR) plane. (b) For the first time, we demonstrate that the OLR is expected to cause a gradient in average Jz with Lz across the resonance. (c) We show that the change of predominantly outward to inward motions at the OLR occurs along the resonance line in action space. The latter signature allows us to identify three candidates for the bar's OLR - and therefore its pattern speed Omega_bar - in the Gaia data within 3 kpc from the Sun: 1.85 Omega0, 1.2 Omega0, and 1.6 Omega0 (with ~0.1 Omega0 uncertainty). This demonstrates that (i) the local Gaia action data is consistent with both the short-fast and long-slow bar models in the literature, and that (ii) axisymmetrically estimated actions are a powerful diagnostic even in non-axisymmetric systems.

The calibration procedures for data collected for Argus and the 4mm Receiver instruments on the GBT are presented. The measured beam size, aperture efficiency, and main-beam efficiency are derived for the range of observing frequencies (66--116 GHz) within the 3mm atmospheric window. The telescope performs well even at the highest frequencies (>100 GHz). The amount of power in the error pattern of the antenna beam is estimated. Approximately 95% of the total antenna power is contained within 1 degree of the pointed direction. The calibration derived using small sources follows theoretical expectations based on the Ruze equation. In general, bright point source calibrators that are monitored regularly by ALMA can be used to accurately derive the main-beam efficiency of the telescope.

We investigate the impact of dynamical streams and substructure on estimates of the local escape speed and total mass of Milky Way-mass galaxies from modelling the high velocity tail of local halo stars. We use a suite of high-resolution, magneto-hydrodynamical cosmological zoom-in simulations, which resolve phase space substructure in local volumes around solar-like positions. We show that phase space structure varies significantly between positions in individual galaxies and across the suite. Substructure populates the high velocity tail unevenly and leads to discrepancies in the mass estimates. We show that a combination of streams, sample noise and truncation of the high velocity tail below the escape speed leads to a distribution of mass estimates with a median that falls below the true value by $\sim 20 \%$, and a spread of a factor of 2 across the suite. Correcting for these biases, we derive a revised value for the Milky Way mass presented in Deason et al. of $1.29 ^{+0.37}_{-0.47} \times 10^{12}$ $\rm M_{\odot}$.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of a quiet-Sun region at a wavelength of 3 mm, obtained during the first solar ALMA cycle on April 27, 2017, and compare them with available chromospheric observations in the UV and visible as well as with photospheric magnetograms. ALMA images clearly reveal the presence of distinct particularly dark/cool areas in the millimeter maps having temperatures as low as 60% of the normal quiet Sun at 3 mm, which are not seen in the other data. We speculate that ALMA is sensing cool chromospheric gas, whose presence had earlier been inferred from infrared CO spectra.

The Sun is an active star that can have a direct impact on the Earth, its magnetosphere, and the technological infrastructure on which modern society depends. Among the phenomena that drive "space weather" are fast solar wind streams and co-rotating interaction regions, solar flares, coronal mass ejections, the shocks they produce, and the energetic particles they accelerate. Radio emission from these and associated phenomena offer unique diagnostic possibilities that complement those available at other wavelengths. Here, the relevant space weather drivers are briefly described, the potential role of radio observations is outlined, and the requirements of an instrument to provide them are provided: specifically, ultrabroadband imaging spectropolarimetry. The insights provided by radio observations of space weather drivers will not only inform the science of space weather, they will pave the way for new tools for forecasting and "nowcasting" space weather. They will also serve as an important touchstone against which local environment of exoplanets and the impact of "exo-space weather" can be evaluated.

Identification of the mechanisms responsible for heating the solar chromosphere and corona remains an outstanding problem, one of great relevance to late-type stars as well. There has been tremendous progress in the past decade, largely driven by new instruments, new observations, and sophisticated modeling efforts. Despite this progress, gaps remain. We briefly discuss the need for radio coverage of the 3D solar atmosphere and discuss the requirements.

Magnetic reconnection is a fundamental physical process in many laboratory, space, and astrophysical plasma contexts. Solar flares serve as an outstanding laboratory to study the magnetic reconnection and the associated energy release and conversion processes under plasma conditions difficult to reproduce in the laboratory, and with considerable spatiotemporal details not possible elsewhere in astrophysics. Here we emphasize the unique power of remote-sensing observations of solar flares at radio wavelengths. In particular, we discuss the transformative technique of broadband radio dynamic imaging spectroscopy in making significant contributions to addressing several outstanding challenges in magnetic reconnection, including the capability of pinpointing magnetic reconnection sites, measuring the time-evolving reconnecting magnetic fields, and deriving the spatially and temporally resolved distribution function of flare-accelerated electrons.

We describe the design and deployment of GREENBURST, a commensal Fast Radio Burst (FRB) search system at the Green Bank Telescope. GREENBURST uses the dedicated L-band receiver tap to search over the 960$-$1920 MHz frequency range for pulses with dispersion measures out to $10^4$ pc cm$^{-3}$. Due to its unique design, GREENBURST will obtain data even when the L-band receiver is not being used for scheduled observing. This makes it a sensitive single pixel detector capable of reaching deeper in the radio sky. While single pulses from Galactic pulsars and rotating radio transients will be detectable in our observations, and will form part of the database we archive, the primary goal is to detect and study FRBs. Based on recent determinations of the all-sky rate, we predict that the system will detect approximately one FRB for every 2$-$3 months of continuous operation. The high sensitivity of GREENBURST means that it will also be able to probe the slope of the FRB source function, which is currently uncertain in this observing band.

We introduce a Local Background Environment (LBE) estimator that can be measured in and around every galaxy or its dark matter subhalo in high-resolution cosmological simulations. The LBE is designed to capture the influence of environmental effects such as ram-pressure stripping on the formation and evolution of galaxies in semi-analytical models. We define the LBE directly from the particle data within an adaptive spherical shell, and devise a Gaussian mixture estimator to separate background particles from previously unidentified subhalo particles. Analyzing the LBE properties, we find that the LBE of satellite galaxies is not at rest with respect to their host halo, in contrast to typical assumptions. The orientations of the velocities of a subhalo and its LBE are well aligned in the outer infall regions of haloes, but decorrelated near halo center. Significantly, there is no abrupt change in LBE velocity or density at the halo virial radius. This suggests that stripping should also happen beyond this radius. Therefore, we use the time-evolving LBE of galaxies to develop a method to better account for ram-pressure stripping within the Munich semi-analytical model, L-Galaxies. Overall, our new approach results in a significant increase in gas stripping across cosmic time. Central galaxies, as well as satellites beyond the virial radius, can lose a significant fraction of their hot halo gas. As a result, the gas fractions and star formation rates of satellite galaxies are suppressed relative to the fiducial model, although the stellar masses and global stellar mass functions are largely unchanged.

We examine the stellar haloes of the Auriga simulations, a suite of thirty cosmological magneto-hydrodynamical high-resolution simulations of Milky Way-mass galaxies performed with the moving-mesh code AREPO. We study halo global properties and radial profiles out to $\sim 150$ kpc for each individual galaxy. The Auriga haloes are diverse in their masses and density profiles; mean metallicity and metallicity gradients; ages; and shapes, reflecting the stochasticity inherent in their accretion and merger histories. A comparison with observations of nearby late-type galaxies shows very good agreement between most observed and simulated halo properties. However, Auriga haloes are typically too massive. We find a connection between population gradients and mass assembly history: galaxies with few significant progenitors have more massive haloes, possess large negative halo metallicity gradients and steeper density profiles. The number of accreted galaxies, either disrupted or under disruption, that contribute 90% of the accreted halo mass ranges from 1 to 14, with a median of 6.5, and their stellar masses span over three orders of magnitude. The observed halo mass--metallicity relation is well reproduced by Auriga and is set by the stellar mass and metallicity of the dominant satellite contributors. This relationship is found not only for the accreted component but also for the total (accreted + in-situ) stellar halo. Our results highlight the potential of observable halo properties to infer the assembly history of galaxies.

In this work we use solar observations with the ALMA radio telescope at the wavelength of 1.21 mm. The aim of the analysis is to improve understanding of the solar chromosphere, a dynamic layer in the solar atmosphere between the photosphere and corona. The study has an observational and a modeling part. In the observational part full-disc solar images are analyzed. Based on a modified FAL atmospheric model, radiation models for various observed solar structures are developed. Finally, the observational and modeling results are compared and discussed.

We present deep CCS and HC$_7$N observations of the L1495-B218 filaments in the Taurus molecular cloud obtained using the K-band focal plane array on the 100m Green Bank Telescope. We observed the L1495-B218 filaments in CCS $J_N$ = 2$_1$$-$1$_0$ and HC$_7$N $J$ = 21$-$20 with a spectral resolution of 0.038 km s$^{-1}$ and an angular resolution of 31$''$. We observed strong CCS emission in both evolved and young regions and weak emission in two evolved regions. HC$_7$N emission is observed only in L1495A-N and L1521D. We find that CCS and HC$_7$N intensity peaks do not coincide with NH$_3$ or dust continuum intensity peaks. We also find that the fractional abundance of CCS does not show a clear correlation with the dynamical evolutionary stage of dense cores. Our findings and chemical modeling indicate that the fractional abundances of CCS and HC$_7$N are sensitive to the initial gas-phase C/O ratio, and they are good tracers of young condensed gas only when the initial C/O is close to solar value. Kinematic analysis using multiple lines including NH$_3$, HC$_7$N, CCS, CO, HCN, \& HCO$^+$ suggests that there may be three different star formation modes in the L1495-B218 filaments. At the hub of the filaments, L1495A/B7N has formed a stellar cluster with large-scale inward flows (fast mode), while L1521D, a core embedded in a filament, is slowly contracting due to its self-gravity (slow mode). There is also one isolated core that appears to be marginally stable and may undergo quasi-static evolution (isolated mode).

The GaLactic and Extragalactic All-sky MWA survey (GLEAM) is a radio continuum survey at 72-231 MHz of the whole sky south of declination +30 deg, carried out with the Murchison Widefield Array (MWA). In this paper, we derive source counts from the GLEAM data at 200, 154, 118 and 88 MHz, to a flux density limit of 50, 80, 120 and 290 mJy respectively, correcting for ionospheric smearing, incompleteness and source blending. These counts are more accurate than other counts in the literature at similar frequencies as a result of the large area of sky covered and this survey's sensitivity to extended emission missed by other surveys. At S_154MHz > 0.5 Jy, there is no evidence of flattening in the average spectral index (alpha approx. -0.8 where S proportional to nu^alpha) towards the lower frequencies. We demonstrate that the SKA Design Study (SKADS) model by Wilman et al. (2008) significantly underpredicts the observed 154 MHz GLEAM counts, particularly at the bright end. Using deeper LOFAR counts and the SKADS model, we find that sidelobe confusion dominates the thermal noise and classical confusion at nu >~ 100 MHz due to both the limited CLEANing depth and undeconvolved sources outside the field-of-view. We show that we can approach the theoretical noise limit using a more efficient and automated CLEAN algorithm.

Polar brightening of the Sun at radio frequencies has been studied for almost fifty years and yet a disagreement persists between solar atmospheric models and observations. Some observations reported brightening values much smaller than the expected values obtained from the models, with discrepancies being particularly large at millimeter wavelengths. New clues to calibrate the atmospheric models can be obtained with the advent of the Atacama Large Millimeter/submillimeter Array (ALMA) radio interferometer. In this work, we analyzed the lower limit of the polar brightening observed at 100 and 230 GHz by ALMA, during its Science Verification period, 2015 December 16-20. We find that the average polar intensity is higher than the disk intensity at 100 and 230 GHz, with larger brightness intensities at the South pole in eight of the nine maps analyzed. The observational results were compared with calculations of the millimetric limb brightnening emission for two semi-empirical atmospheric models, FAL- C (Fontenla et al. 1993) and SSC (Selhorst et al. 2005a). Both models presented larger limb intensities than the averaged observed values. The intensities obtained with the SSC model were closer to the observations, with polar brightenings of 10.5% and 17.8% at 100 and 230 GHz, respectively. This discrepancy may be due to the presence of chromospheric features (like spicules) at regions close to the limb.

The Longitudinal Study of Astronomy Graduate Students (LSAGS) arose from the 2003 Women in Astronomy Conference, where it was noted that a majority of young members of the American Astronomical Society were women. The astronomy community wishes to make every effort to retain young women in astronomy, so they commissioned a longitudinal study to be conducted that would pinpoint the factors that contribute to retention in general, with a focus on differences between women and men. The LSAGS follows a cohort of people who were graduate students in astronomy or astrophysics during 2006-07. The first survey was conducted during 2007-08, the second during 2012-13, and the third during 2015. The analysis presented in this paper, which is an update to our previous paper on this topic, used a subset of the respondents, all of whom had PhDs in astronomy, astrophysics, or a related field at the time of the third survey. We tested the effects of four major concepts on attrition from physics and astronomy. These concepts included: the imposter syndrome, mentoring and advising during graduate school, the so-called "two-body problem" that occurs when a couple needs to find two jobs in the same geographic area, and gender of the respondent. Having a mentor in grad school did not contribute to working outside of physics or astronomy. Showing characteristics of the imposter syndrome and gender of the respondent had indirect effects on working outside the field. Encouragement of the graduate advisor, the two-body problem, and completing a postdoc, had significant direct effects on working in physics or astronomy. This research identifies specific areas of concern that can be addressed by the scientific community to increase the retention of all people, but especially women, in astronomy and astrophysics.

We use cosmological hydrodynamical simulations of the APOSTLE project along with high-quality rotation curve observations to examine the fraction of baryons in {\Lambda}CDM haloes that collect into galaxies. This 'galaxy formation efficiency' correlates strongly and with little scatter with halo mass, dropping steadily towards dwarf galaxies. The baryonic mass of a galaxy may thus be used to place a lower limit on total halo mass and, consequently, on its asymptotic maximum circular velocity. A number of observed dwarfs seem to violate this constraint, having baryonic masses up to ten times higher than expected from their rotation speeds, or, alternatively, rotating at only half the speed expected for their mass. Taking the data at face value, either these systems have formed galaxies with extraordinary efficiency - highly unlikely given their shallow potential wells - or their dark matter content is much lower than expected from {\Lambda}CDM haloes. This 'missing dark matter' is reminiscent of the inner mass deficit of galaxies with slowly-rising rotation curves, but cannot be explained away by star formation-induced 'cores' in the dark mass profile, since the anomalous deficit applies to regions larger than the luminous galaxies themselves. We argue that explaining the structure of these galaxies would require either substantial modification of the standard Lambda cold dark matter paradigm or else significant revision to the uncertainties in their inferred mass profiles, which should be much larger than reported. Systematic errors in inclination may provide a simple resolution to what would otherwise be a rather intractable problem for the current paradigm.

The upcoming next-generation large area radio continuum surveys can expect tens of millions of radio sources, rendering the traditional method for radio morphology classification through visual inspection unfeasible. We present ClaRAN - Classifying Radio sources Automatically with Neural networks - a proof-of-concept radio source morphology classifier based upon the Faster Region-based Convolutional Neutral Networks (Faster R-CNN) method. Specifically, we train and test ClaRAN on the FIRST and WISE images from the Radio Galaxy Zoo Data Release 1 catalogue. ClaRAN provides end users with automated identification of radio source morphology classifications from a simple input of a radio image and a counterpart infrared image of the same region. ClaRAN is the first open-source, end-to-end radio source morphology classifier that is capable of locating and associating discrete and extended components of radio sources in a fast (< 200 milliseconds per image) and accurate (>= 90 %) fashion. Future work will improve ClaRAN's relatively lower success rates in dealing with multi-source fields and will enable ClaRAN to identify sources on much larger fields without loss in classification accuracy.

The authors of this report met on 28-30 March 2018 at the New Jersey Institute of Technology, Newark, New Jersey, for a 3-day workshop that brought together a group of data providers, expert modelers, and computer and data scientists, in the solar discipline. Their objective was to identify challenges in the path towards building an effective framework to achieve transformative advances in the understanding and forecasting of the Sun-Earth system from the upper convection zone of the Sun to the Earth's magnetosphere. The workshop aimed to develop a research roadmap that targets the scientific challenge of coupling observations and modeling with emerging data-science research to extract knowledge from the large volumes of data (observed and simulated) while stimulating computer science with new research applications. The desire among the attendees was to promote future trans-disciplinary collaborations and identify areas of convergence across disciplines. The workshop combined a set of plenary sessions featuring invited introductory talks and workshop progress reports, interleaved with a set of breakout sessions focused on specific topics of interest. Each breakout group generated short documents, listing the challenges identified during their discussions in addition to possible ways of attacking them collectively. These documents were combined into this report-wherein a list of prioritized activities have been collated, shared and endorsed.

We report on a search for low-frequency radio variability in 944 bright (> 4Jy at 154 MHz) unresolved, extragalactic radio sources monitored monthly for several years with the Murchison Widefield Array. In the majority of sources we find very low levels of variability with typical modulation indices < 5%. We detect 15 candidate low frequency variables that show significant long term variability (>2.8 years) with time-averaged modulation indices M = 3.1 - 7.1%. With 7/15 of these variable sources having peaked spectral energy distributions, and only 5.7% of the overall sample having peaked spectra, we find an increase in the prevalence of variability in this spectral class. We conclude that the variability seen in this survey is most probably a consequence of refractive interstellar scintillation and that these objects must have the majority of their flux density contained within angular diameters less than 50 milli-arcsec (which we support with multi-wavelength data). At 154 MHz we demonstrate that interstellar scintillation time-scales become long (~decades) and have low modulation indices, whilst synchrotron driven variability can only produce dynamic changes on time-scales of hundreds of years, with flux density changes less than one milli-jansky (without relativistic boosting). From this work we infer that the low frequency extra-galactic southern sky, as seen by SKA-Low, will be non-variable on time-scales shorter than one year.

Energy stored in the magnetic field in the solar atmosphere above active regions is a key driver of all solar activity (e.g., solar flares and coronal mass ejections), some of which can affect life on Earth. Radio observations provide a unique diagnostic of the coronal magnetic fields that make them a critical tool for the study of these phenomena, using the technique of broadband radio imaging spectropolarimetry. Observations with the ngVLA will provide unique observations of coronal magnetic fields and their evolution, key inputs and constraints for MHD numerical models of the solar atmosphere and eruptive processes, and a key link between lower layers of the solar atmosphere and the heliosphere. In doing so they will also provide practical "research to operations" guidance for space weather forecasting.

We perform simulations of the capabilities of the next generation Very Large Array to image stellar radio photospheres. For very large (in angle) stars, such as red supergiants within a few hundred parsecs, good imaging fidelity results can be obtained on radio photospheric structures at 38 GHz employing standard techniques, such as disk model fitting and subtraction, with hundreds of resolution elements over the star, even with just the ngVLA-classic baselines to 1000 km. Using the ngVLA Rev B plus long baseline configuration (with baselines out to 9000 km, August 2018), we find for main sequence stars within $\sim$ 10 pc, the photospheres can be easily resolved at 85 GHz, with accurate measures of the mean brightness and size, and possibly imaging large surface structures, as might occur on e.g., active M dwarf stars. For more distant main sequence stars, we find that measurements of sizes and brightnesses can be made using disk model fitting to the u,v-data down to stellar diameters $\sim$ 0.4 mas in a few hours. This size would include M0 V stars to a distance of 15 pc, A0 V stars to 60 pc, and Red Giants to 2.4 kpc. Based on the Hipparcos catalog, we estimate that there are at least 10,000 stars that will be resolved by the ngVLA. While the vast majority of these (95\%) are giants or supergiants, there are still over 500 main sequence stars that can be resolved, with $\sim$ 50 to 150 in each spectral type (besides O stars). Note that these are lower limits, since radio photospheres can be larger than optical, and the Hipparcos catalog might not be complete. Our initial look into the Gaia catalog suggests these numbers might be pessimistic by a factor few.

Powerful radio-galaxies feature heavily in our understanding of galaxy evolution. However, when it comes to studying their properties as a function of redshift and/or environment, the most-detailed studies tend to be limited by small-number statistics. During Focus Meeting 3, on "Radio Galaxies: Resolving the AGN phenomenon", SVW presented a new sample of nearly 2,000 of the brightest radio-sources in the southern sky (Dec. $<$ 30 deg). These were observed at low radio-frequencies as part of the GaLactic and Extragalactic All-sky MWA (GLEAM) Survey, which is a continuum survey conducted using the Murchison Widefield Array (MWA). This instrument is the precursor telescope for the low-frequency component of the Square Kilometre Array, and allows us to select radio galaxies in an orientation-independent way (i.e. minimising the bias caused by Doppler boosting, inherent in high-frequency surveys). Being brighter than 4 Jy at 151 MHz, we refer to these objects as the GLEAM 4-Jy (G4Jy) Sample. The G4Jy catalogue is close to being finalised, with SVW describing how multi-wavelength data have been used to determine the morphology of the radio emission, and identify the host galaxy. In addition, the MWA's excellent spectral-coverage and sensitivity to extended/diffuse emission were highlighted. Both of these aspects are important for understanding the physical mechanisms that take place within active galaxies, and how they interact with their environment.

The low-frequency polarisation properties of radio sources are poorly studied, particularly in statistical samples. However, the new generation of low-frequency telescopes, such as the Murchison Widefield Array (MWA; the precursor for the low-frequency component of the Square Kilometre Array) offers an opportunity to probe the physics of radio sources at very low radio frequencies. In this paper, we present a catalogue of linearly-polarised sources detected at 216 MHz, using data from the Galactic and Extragalactic All-sky MWA (GLEAM) survey. Our catalogue covers the Declination range $-17^{\circ}$ to $-37^{\circ}$ and 24 hours in Right Ascension, at a resolution of around 3 arcminutes. We detect 81 sources (including both a known pulsar and new pulsar candidate) with linearly-polarised flux densities in excess of 18 mJy across a survey area of approximately 6400 square degrees, corresponding to a surface density of 1 source per 79 square degrees. The level of Faraday rotation measured for our sources is broadly consistent with those recovered at higher frequencies, with typically more than an order of magnitude improvement in the uncertainty compared to higher-frequency measurements. However, our catalogue is likely incomplete at low Faraday rotation measures, due to our practice of excluding sources in the region where instrumental leakage appears. The majority of sources exhibit significant depolarisation compared to higher frequencies; however, a small sub-sample repolarise at 216 MHz. We also discuss the polarisation properties of four nearby, large-angular-scale radio galaxies, with a particular focus on the giant radio galaxy ESO 422$-$G028, in order to explain the striking differences in polarised morphology between 216 MHz and 1.4 GHz.

We use the Henriques et al. (2015) version of the Munich galaxy formation model (L-GALAXIES) to investigate why the halo and stellar mass scales above which galaxies are quenched are constant with redshift and coincide with the scale where baryons are most efficiently converted into stars. This model assumes that central galaxies are quenched by AGN feedback when hot halo gas accretes onto a supermassive black hole. Nevertheless, we find that supernova (SN) feedback sets both mass scales. As haloes grow above a threshold mass, SNe can no longer eject material so their hot gas content increases, enhancing the cooling rate onto the central galaxy, its cold gas content, its star formation rate and the growth rate of its central black hole. Strong AGN feedback terminates this short-lived phase by suppressing the fuel supply for star formation. Despite strong evolution of the halo mass - temperature relation, quenching occurs at a redshift-independent halo and stellar mass which coincides with the mass where baryons have been converted into stars with maximal efficiency. These regularities and coincidences are a result of the specific parameters selected by MCMC tuning of the model to fit the observed abundance and passive fraction of galaxies over the redshift range 0<z<3. Thus they are required by the observed evolution of the galaxy population, at least in the context of models of this type.

We present a single-dish mapping algorithm with a number of advantages over traditional techniques. (1) Our algorithm makes use of weighted modeling, instead of weighted averaging, to interpolate between signal measurements. This smooths the data, but without blurring the data beyond instrumental resolution. Techniques that rely on weighted averaging blur point sources sometimes as much as 40%. (2) Our algorithm makes use of local, instead of global, modeling to separate astronomical signal from instrumental and/or environmental signal drift along the telescope's scans. Other techniques, such as basket weaving, model this drift with simple functional forms (linear, quadratic, etc.) across the entirety of scans, limiting their ability to remove such contaminants. (3) Our algorithm makes use of a similar, local modeling technique to separate astronomical signal from radio-frequency interference (RFI), even if only continuum data are available. (4) Unlike other techniques, our algorithm does not require data to be collected on a rectangular grid or regridded before processing. (5) Data from any number of observations, overlapping or not, may be appended and processed together. (6) Any pixel density may be selected for the final image. We present our algorithm, and evaluate it using both simulated and real data. We are integrating it into the image-processing library of the Skynet Robotic Telescope Network, which includes optical telescopes spanning four continents, and now also Green Bank Observatory's 20-meter diameter radio telescope in West Virginia. Skynet serves hundreds of professional users, and additionally tens of thousands of students, of all ages. Default data products are generated on the fly, but will soon be customizable after the fact.

We report the first science results from the newly completed Expanded Owens Valley Solar Array (EOVSA), which obtained excellent microwave imaging spectroscopy observations of SOL2017-09-10, a classic partially-occulted solar limb flare associated with an erupting flux rope. This event is also well-covered by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) in hard X-rays (HXRs). We present an overview of this event focusing on microwave and HXR data, both associated with high-energy nonthermal electrons, and discuss them within the context of the flare geometry and evolution revealed by extreme ultraviolet (EUV) observations from the Atmospheric Imaging Assembly aboard the Solar Dynamics Observatory (SDO/AIA). The EOVSA and RHESSI data reveal the evolving spatial and energy distribution of high-energy electrons throughout the entire flaring region. The results suggest that the microwave and HXR sources largely arise from a common nonthermal electron population, although the microwave imaging spectroscopy provides information over a much larger volume of the corona.

This is a personal and historical essay based on the Shaw Prize Lecture given in Hong Kong in September 2017 in association with the 2017 Shaw Prize in Astronomy.

The ionosphere contributes time-varying Faraday Rotation (FR) to radio signals passing through it. Correction for the effect of the ionosphere is important for deriving magnetic field information from FR observations of polarized cosmic radio sources, as well as providing valuable diagnostics of the structure of the ionosphere. In this paper, we evaluate the accuracy of models commonly used to correct for its effects using new observations of pulsars at low frequencies, which provide total rotation measures (RM) at better precision than previously available. We evaluate models of the ionosphere derived from modern digital ionosondes that provide electron density information as a function of height, as well as GPS-derived Total Electron Content (TEC) measurements. We combine these density models with reference global magnetic field models to derive ionospheric RM contributions. We find that the models disagree substantially with each other and seek corrections that may explain the differences in RM prediction. Additionally we compare these models to global TEC models and find that local high-cadence TEC measurements are superior to global models for ionospheric RM correction.

The Geostationary Orbital Environmental Satellites (GOES) Soft X-ray (SXR) sensors have provided data relating to, inter alia, the time, intensity and duration of solar flares since the 1970s. The GOES SXR Flare List has become the standard reference catalogue for solar flares and is widely used in solar physics research and space weather. We report here that in the current version of the list there are significant differences between the mean duration of flares which occurred before May 1997 and the mean duration of flares thereafter. Our analysis shows that the reported flare timings for the pre-May 1997 data were not based on the same criteria as is currently the case. This finding has serious implications for all those who used flare duration (or fluence, which depends on the chosen start and end times) as part of their analysis of pre-May 1997 solar events, or statistical analyses of large samples of flares, e.g. as part of the assessment of a Solar Energetic Particle forecasting algorithm.

In this paper we report the experimental implementation of a theoretically-proposed technique for creating a photoionized plasma in the laboratory using X-ray line radiation. Using a Sn laser-plasma to irradiate an Ar gas target, the photoionization parameter, \xi = 4{\pi}F/Ne, reached values of order 50 erg cm/s, where F is the radiation flux in erg/cm2/s. The significance of this is that this technique allows us to mimic effective spectral radiation temperatures in excess of 1 keV. We show that our plasma starts to be collisionally dominated before the peak of the X-ray drive. However, the technique is extendable to higher energy laser systems to create plasmas with parameters relevant to benchmarking codes used to model astrophysical objects.

In this paper we carry out anisotropic "separate universe" simulations by including a large-scale tidal field in the N-body code \textsc{gadget}4 using an anisotropic expansion factor $A_{ij}$. We use the code in a pure \textit{particle-mesh} (PM) mode to simulate the evolution of 16 realizations of an initial density field with and without a large-scale tidal field, which are then used to measure the \textit{response function} describing how the tidal field influences structure formation in the linear and non-linear regimes. Together with the previously measured response to a large scale overdensity, this completely describes the nonlinear matter bispectrum in the squeezed limit. We find that, contrary to the density response, the tidal response never significantly exceeds the large-scale perturbation-theory prediction even on nonlinear scales for the redshift range we discuss. We develop a simple halo model that takes into account the effect of the tidal field and compare it with our direct measurement from the anisotropic N-body simulations.

The explanation of the coronal heating problem potentially lies in the existence of nanoflares, numerous small-scale heating events occuring across the whole solar disk. In this paper, we present the first imaging spectroscopy X-ray observations of three quiet Sun flares during the NuSTAR solar campaigns on 2016 July 26 and 2017 March 21, concurrent with SDO/AIA observations. Two of the three events showed time lags of a few minutes between peak X-ray and extreme ultraviolet (EUV) emissions. Isothermal fits with rather low temperatures in the range $3.2-4.1$ MK and emission measures of $(0.6-15)\times10^{44} \textrm{ cm}^{-3}$ describe their spectra well, resulting in thermal energies in the range $(2-6)\times10^{26}\textrm{ ergs}$. NuSTAR spectra did not show any signs of a nonthermal or higher temperature component. However, since the estimated upper limits of (hidden) nonthermal energy are comparable to the thermal energy estimates, the lack of a nonthermal component in the observed spectra is not a constraining result. The estimated GOES classes from the fitted values of temperature and emission measure fall between $1/1000 \textrm{ and } 1/100$ A class level, making them 8 orders of magnitude fainter in soft X-ray flux than the largest solar flares.

Radio observations are currently the only way to explore accelerated particles in cool stellar environments. We describe how a next generation VLA can contribute to the understanding of the stellar contribution to exo-space weather. This area holds both academic and popular interest, and is expected to grow in the next several decades.

A new 1.4 GHz 19-element, dual-polarization, cryogenic phased array feed (PAF) radio astronomy receiver has been developed for the Robert C. Byrd Green Bank Telescope (GBT) as part of FLAG (Focal L-band Array for the GBT) project. Commissioning observations of calibrator radio sources show that this receiver has the lowest reported beamformed system temperature ($T_{\rm sys}$) normalized by aperture efficiency ($\eta$) of any phased array receiver to date. The measured $T_{\rm sys}/\eta$ is $25.4 \pm 2.5$ K near 1350 MHz for the boresight beam, which is comparable to the performance of the current 1.4 GHz cryogenic single feed receiver on the GBT. The degradation in $T_{\rm sys}/\eta$ at $\sim$ 4 arcmin (required for Nyquist sampling) and $\sim$ 8 arcmin offsets from the boresight is, respectively, $\sim$ 1\% and $\sim$ 20\% of the boresight value. The survey speed of the PAF with seven formed beams is larger by a factor between 2.1 and 7 compared to a single beam system depending on the observing application. The measured performance, both in frequency and offset from boresight, qualitatively agree with predictions from a rigorous electromagnetic model of the PAF. The astronomical utility of the receiver is demonstrated by observations of the pulsar B0329+54 and an extended HII region, the Rosette Nebula. The enhanced survey speed with the new PAF receiver will enable the GBT to carry out exciting new science, such as more efficient observations of diffuse, extended neutral hydrogen emission from galactic in-flows and searches for Fast Radio Bursts.

We present cluster counts and corresponding cosmological constraints from the Planck full mission data set. Our catalogue consists of 439 clusters detected via their Sunyaev-Zeldovich (SZ) signal down to a signal-to-noise ratio of 6, and is more than a factor of 2 larger than the 2013 Planck cluster cosmology sample. The counts are consistent with those from 2013 and yield compatible constraints under the same modelling assumptions. Taking advantage of the larger catalogue, we extend our analysis to the two-dimensional distribution in redshift and signal-to-noise. We use mass estimates from two recent studies of gravitational lensing of background galaxies by Planck clusters to provide priors on the hydrostatic bias parameter, $(1-b)$. In addition, we use lensing of cosmic microwave background (CMB) temperature fluctuations by Planck clusters as an independent constraint on this parameter. These various calibrations imply constraints on the present-day amplitude of matter fluctuations in varying degrees of tension with those from the Planck analysis of primary fluctuations in the CMB; for the lowest estimated values of $(1-b)$ the tension is mild, only a little over one standard deviation, while it remains substantial ($3.7\,\sigma$) for the largest estimated value. We also examine constraints on extensions to the base flat $\Lambda$CDM model by combining the cluster and CMB constraints. The combination appears to favour non-minimal neutrino masses, but this possibility does little to relieve the overall tension because it simultaneously lowers the implied value of the Hubble parameter, thereby exacerbating the discrepancy with most current astrophysical estimates. Improving the precision of cluster mass calibrations from the current 10%-level to 1% would significantly strengthen these combined analyses and provide a stringent test of the base $\Lambda$CDM model.

We investigate how Einstein rings and magnified arcs are affected by small-mass dark-matter haloes placed along the line-of-sight to gravitational lens systems. By comparing the gravitational signature of line-of-sight haloes with that of substructures within the lensing galaxy, we derive a mass-redshift relation that allows us to rescale the detection threshold (i.e. lowest detectable mass) for substructures to a detection threshold for line-of-sight haloes at any redshift. We then quantify the line-of-sight contribution to the total number density of low-mass objects that can be detected through strong gravitational lensing. Finally, we assess the degeneracy between substructures and line-of-sight haloes of different mass and redshift to provide a statistical interpretation of current and future detections, with the aim of distinguishing between CDM and WDM. We find that line-of-sight haloes statistically dominate with respect to substructures, by an amount that strongly depends on the source and lens redshifts, and on the chosen dark matter model. Substructures represent about 30 percent of the total number of perturbers for low lens and source redshifts (as for the SLACS lenses), but less than 10 per cent for high redshift systems. We also find that for data with high enough signal-to-noise ratio and angular resolution, the non-linear effects arising from a double-lens-plane configuration are such that one is able to observationally recover the line-of-sight halo redshift with an absolute error precision of 0.15 at the 68 per cent confidence level.

We compare the transverse velocity dispersions recently measured within the Sculptor dwarf spheroidal galaxy to the predictions of previously published dynamical models. These provide good fits to the observed number count and velocity dispersion profiles of metal-rich and metal-poor stars both in cored and in cusped potentials. At the projected radius where the proper motions were measured, these models predict transverse dispersions in the range 6 to 9.5 km/s, with the tangential dispersion about 1 km/s larger than the (projected) radial dispersion. Both dispersions are predicted to be about 1 km/s larger for metal-poor than for metal-rich stars. At this projected radius, cored and cusped potentials predict almost identical transverse dispersions. The measured tangential dispersion ($8.5 \pm 3.2$ km/s) agrees remarkably well with these predictions, while the measured radial dispersion ($11.5 \pm 4.3$ km/s) differs only at about the $1\sigma$ level. Thus, the proper motion data are in excellent agreement with previous data but do not help to distinguish between cored and cusped potentials. This will require velocity dispersion data (either from proper motions or from radial velocities) with uncertainties well below 1 km/s over a range of projected radii.

We use data from the HST Coma Cluster Treasury program to assess the richness of the Globular Cluster Systems (GCSs) of 54 Coma ultra-diffuse galaxies (UDGs), 18 of which have a half-light radius exceeding 1.5 kpc. We use a hierarchical Bayesian method tested on a large number of mock datasets to account consistently for the high and spatially varying background counts in Coma. These include both background galaxies and intra-cluster GCs (ICGCs), which are disentangled from the population of member GCs in a probabilistic fashion. We find no candidate for a GCS as rich as that of the Milky Way, our sample has GCSs typical of dwarf galaxies. For the standard relation between GCS richness and halo mass, 33 galaxies have a virial mass $M_{vir}\leq10^{11}M_\odot$ at 90% probability. Only three have $M_{vir}>10^{11}M_\odot$ with the same confidence. The mean colour and spread in colour of the UDG GCs are indistinguishable from those of the abundant population of ICGCs. The majority of UDGs in our sample are consistent with the relation between stellar mass and GC richness of 'normal' dwarf galaxies. Nine systems, however, display GCSs that are richer by a factor of 3 or more (at 90% probability). Six of these have sizes $\lesssim1.4$ kpc. Our results imply that the physical mechanisms responsible for the extended size of the UDGs and for the enhanced GC richness of some cluster dwarfs are at most weakly correlated.

We use publicly available data for the Millennium Simulation to explore the implications of the recent detection of assembly bias and splashback signatures in a large sample of galaxy clusters. These were identified in the SDSS/DR8 photometric data by the redMaPPer algorithm and split into high- and low-concentration subsamples based on the projected positions of cluster members. We use simplified versions of these procedures to build cluster samples of similar size from the simulation data. These match the observed samples quite well and show similar assembly bias and splashback signals. Previous theoretical work has found the logarithmic slope of halo density profiles to have a well-defined minimum whose depth decreases and whose radius increases with halo concentration. Projected profiles for the observed and simulated cluster samples show trends with concentration which are opposite to these predictions. In addition, for high-concentration clusters the minimum slope occurs at significantly smaller radius than predicted. We show that these discrepancies all reflect confusion between splashback features and features imposed on the profiles by the cluster identification and concentration estimation procedures. The strong apparent assembly bias is not reflected in the three-dimensional distribution of matter around clusters. Rather it is a consequence of the preferential contamination of low-concentration clusters by foreground or background groups.

Various solar features can be seen on maps of the Sun in the mm and sub-mm wavelength range. The recently installed Atacama Large Millimeter/submillimeter Array (ALMA) is capable of observing the Sun in that wavelength range with an unprecedented spatial, temporal and spectral resolution. To interpret solar observations with ALMA the first important step is to compare ALMA maps with simultaneous images of the Sun recorded in other spectral ranges. First we identify different structures in the solar atmosphere seen in the optical, IR and EUV parts of the spectrum (quiet Sun (QS), active regions (AR), prominences on the disc, magnetic inversion lines (IL), coronal holes (CH) and coronal bright points (CBPs)) in a full disc solar ALMA image. The second aim is to measure the intensities (brightness temperatures) of those structures and compare them with the corresponding QS level. A full disc solar image at 1.21 mm obtained on December 18, 2015 during a CSV-EOC campaign with ALMA is calibrated and compared with full disc solar images from the same day in H\alpha, in He I 1083 nm core, and with SDO images (AIA at 170 nm, 30.4 nm, 21.1 nm, 19.3 nm, and 17.1 nm and HMI magnetogram). The brightness temperatures of various structures are determined by averaging over corresponding regions of interest in the ALMA image. Positions of the QS, ARs, prominences on the disc, ILs, CHs and CBPs are identified in the ALMA image. At 1.21 mm ARs appear as bright areas (but sunspots are dark), while prominences on the disc and CHs are not discernible from the QS background, although having slightly less intensity than surrounding QS regions. ILs appear as large, elongated dark structures and CBPs correspond to ALMA bright points. These results are in general agreement with sparse earlier measurements at similar wavelengths. The identification of CBPs represents the most important new result.

We present the first results of a search for transient hard X-ray (HXR) emission in the quiet solar corona with the \textit{Nuclear Spectroscopic Telescope Array} (\textit{NuSTAR}) satellite. While \textit{NuSTAR} was designed as an astrophysics mission, it can observe the Sun above 2~keV with unprecedented sensitivity due to its pioneering use of focusing optics. \textit{NuSTAR} first observed quiet Sun regions on 2014 November 1, although out-of-view active regions contributed a notable amount of background in the form of single-bounce (unfocused) X-rays. We conducted a search for quiet Sun transient brightenings on time scales of 100 s and set upper limits on emission in two energy bands. We set 2.5--4~keV limits on brightenings with time scales of 100 s, expressed as the temperature T and emission measure EM of a thermal plasma. We also set 10--20~keV limits on brightenings with time scales of 30, 60, and 100 s, expressed as model-independent photon fluxes. The limits in both bands are well below previous HXR microflare detections, though not low enough to detect events of equivalent T and EM as quiet Sun brightenings seen in soft X-ray observations. We expect future observations during solar minimum to increase the \textit{NuSTAR} sensitivity by over two orders of magnitude due to higher instrument livetime and reduced solar background.

In this paper, we describe the development of a new L-band (1.4 GHz) Cryogenic Phased Array Feed (PAF) system, referred to as the GBT2 array. Results from initial measurements made with the GBT2 array are also presented. The PAF was developed for the Green Bank Telescope (GBT) as part of the Focal L-band Array for the GBT (FLAG) project. During the first stage of the development work (Phase I), a prototype cryogenic 19 element dual-polarized array with "Kite" dipole elements was developed and tested on the GBT. The measured system temperature over efficiency ($T_{sys}/\eta$) ratio for the bore sight beam of the Kite array was 45.5 K at 1.55 GHz. The off-boresight $T_{sys}/\eta$ shows an increase by 13 K at an offset equal to the half power beam width (7$^{'}$.2 at 1.7 GHz). Our measurements indicate that the off-boresight degradation and field-of-view (FoV) limitation of the Kite array is simply due to the fixed array size. To increase the FoV, a new 19-element GBT2 array with larger array spacing was developed during FLAG Phase II. The frequency response of the array was optimized from 1.2 to 1.6 GHz. A system with larger cryostat, new low noise amplifiers (LNAs), down-conversion and digitization close to the front end, unformatted digital transmission over fiber, ROACH II based polyphase filter banks (PFBs) with bandwidth 150 MHz and a data acquisition system that records voltage samples from one of the PFB channels were all developed. The data presented here is processed off-line. The receiver temperature measured with the new system is 17 K at 1.4 GHz, an improvement $>$ 8 K over the previous Kite array. Measurements with the GBT2 array on the telescope are in progress. A real time 150 MHz beamformer is also being developed as part of an NSF-funded collaboration between NRAO/GBO/BYU \& West Virginia University (Beamformer Project) to support science observations.

We characterize and provide a catalog of thirty >100 MeV sustained gamma-ray emission (SGRE) events observed by Fermi LAT. These events are temporally and spectrally distinct from the associated solar flares. Their spectra are consistent with decay of pions produced by >300 MeV protons and are not consistent with electron bremsstrahlung. SGRE start times range from CME onset to two hours later. Their durations range from about four minutes to twenty hours and appear to be correlated with durations of >100 MeV SEP proton events. The >300 MeV protons producing SGRE have spectra that can be fit with power laws with a mean index of ~4 and RMS spread of 1.8. Gamma-ray line measurements indicate that SGRE proton spectra are steeper above 300 MeV than they are below 300 MeV. The number of SGRE protons >500 MeV is on average about ten times more than then the number in the associated flare and about fifty to one hundred times less than the number in the accompanying SEP. SGRE can extend tens of degrees from the are site. Sustained bremsstrahlung from MeV electrons was observed in one SGRE event. Flare >100 keV X-ray emission appears to be associated with SGRE and with intense SEPs. From this observation, we provide arguments that lead us to propose that sub-MeV to MeV protons escaping from the flare contribute to the seed population that is accelerated by shocks onto open field lines to produce SEPs and onto field lines returning to the Sun to produce SGRE.

Despite the fact that the mean matter density of the universe has been measured to an accuracy of a few percent within the standard $\Lambda$CDM paradigm, its median density is not known even to order of magnitude. Typical points lie in low-density regions and are not part of a collapsed structure of any scale. Locally, the dark matter distribution is then simply a stretched version of that in the early universe. In this single-stream regime, the distribution of unsmoothed density is sensitive to the initial power spectrum on all scales, in particular on very small scales, and hence to the nature of the dark matter. It cannot be estimated reliably using conventional cosmological simulations because of the enormous dynamic range involved, but a suitable excursion set procedure can be used instead. For the Planck cosmological parameters, a 100 GeV WIMP, corresponding to a free-streaming mass $\sim 10^{-6}M_\odot$, results in a median density of $\sim 4\times 10^{-3}$ in units of the mean density, whereas a 10 $\mu$eV axion with free-streaming mass $\sim 10^{-12}M_\odot$ gives $\sim 3\times 10^{-3}$, and Warm Dark Matter with a (thermal relic) mass of 1 keV gives $\sim 8\times 10^{-2}$. In CDM (but not in WDM) universes, single-stream regions are predicted to be topologically isolated by the excursion set formalism. A test by direct N-Body simulations seems to confirm this prediction, although it is still subject to finite size and resolution effects. Unfortunately, it is unlikely that any of these properties is observable and so suitable for constraining the properties of dark matter.

We present the first high-resolution Atacama Large Millimeter/Submillimeter Array (ALMA) observations of a sunspot at wavelengths of 1.3 mm and 3 mm, obtained during the solar ALMA Science Verification campaign in 2015, and compare them with the predictions of semi-empirical sunspot umbral/penumbral atmosphere models. For the first time millimeter observations of sunspots have resolved umbral/penumbral brightness structure at the chromospheric heights, where the emission at these wavelengths is formed. We find that the sunspot umbra exhibits a radically different appearance at 1.3 mm and 3 mm, whereas the penumbral brightness structure is similar at the two wavelengths. The inner part of the umbra is ~600 K brighter than the surrounding quiet Sun (QS) at 3 mm and is ~700 K cooler than the QS at 1.3 mm, being the coolest part of sunspot at this wavelength. On average, the brightness of the penumbra at 3 mm is comparable to the QS brightness, while at 1.3 mm it is ~1000 K brighter than the QS. Penumbral brightness increases towards the outer boundary in both ALMA bands. Among the tested umbral models, that of Severino et al. (1994) provides the best fit to the observational data, including both the ALMA data analyzed in this study and data from earlier works. No penumbral model amongst those considered here gives a satisfactory fit to the currently available measurements. ALMA observations at multiple mm wavelengths can be used for testing existing sunspot models, and serve as an important input to constrain new empirical models.

We report a NuSTAR observation of a solar microflare, SOL2015-09-01T04. Although it was too faint to be observed by the GOES X-ray Sensor, we estimate the event to be an A0.1 class flare in brightness. This microflare, with only 5 counts per second per detector observed by RHESSI, is fainter than any hard X-ray (HXR) flare in the existing literature. The microflare occurred during a solar pointing by the highly sensitive NuSTAR astrophysical observatory, which used its direct focusing optics to produce detailed HXR microflare spectra and images. The microflare exhibits HXR properties commonly observed in larger flares, including a fast rise and more gradual decay, earlier peak time with higher energy, spatial dimensions similar to the RHESSI microflares, and a high-energy excess beyond an isothermal spectral component during the impulsive phase. The microflare is small in emission measure, temperature, and energy, though not in physical size; observations are consistent with an origin via the interaction of at least two magnetic loops. We estimate the increase in thermal energy at the time of the microflare to be 2.4x10^27 ergs. The observation suggests that flares do indeed scale down to extremely small energies and retain what we customarily think of as "flarelike" properties.

NuSTAR is a highly sensitive focusing hard X-ray (HXR) telescope and has observed several small microflares in its initial solar pointings. In this paper, we present the first joint observation of a microflare with NuSTAR and Hinode/XRT on 2015 April 29 at ~11:29 UT. This microflare shows heating of material to several million Kelvin, observed in Soft X-rays (SXRs) with Hinode/XRT, and was faintly visible in Extreme Ultraviolet (EUV) with SDO/AIA. For three of the four NuSTAR observations of this region (pre-, decay, and post phases) the spectrum is well fitted by a single thermal model of 3.2-3.5 MK, but the spectrum during the impulsive phase shows additional emission up to 10 MK, emission equivalent to A0.1 GOES class. We recover the differential emission measure (DEM) using SDO/AIA, Hinode/XRT, and NuSTAR, giving unprecedented coverage in temperature. We find the pre-flare DEM peaks at ~3 MK and falls off sharply by 5 MK; but during the microflare's impulsive phase the emission above 3 MK is brighter and extends to 10 MK, giving a heating rate of about $2.5 \times 10^{25}$ erg s$^{-1}$. As the NuSTAR spectrum is purely thermal we determined upper-limits on the possible non-thermal bremsstrahlung emission. We find that for the accelerated electrons to be the source of the heating requires a power-law spectrum of $\delta \ge 7$ with a low energy cut-off $E_{c} \lesssim 7$ keV. In summary, this first NuSTAR microflare strongly resembles much more powerful flares.

We introduce the Cluster-EAGLE (C-EAGLE) simulation project, a set of cosmological hydrodynamical zoom simulations of the formation of $30$ galaxy clusters in the mass range $10^{14}<M_{200}/\mathrm{M}_{\odot}<10^{15.4}$ that incorporates the Hydrangea sample of Bah\'e et al. (2017). The simulations adopt the state-of-the-art EAGLE galaxy formation model, with a gas particle mass of $1.8\times10^{6}\,\mathrm{M}_{\odot}$ and physical softening length of $0.7\,\mathrm{kpc}$. In this paper, we introduce the sample and present the low-redshift global properties of the clusters. We calculate the X-ray properties in a manner consistent with observational techniques, demonstrating the bias and scatter introduced by using estimated masses. We find the total stellar content and black hole masses of the clusters to be in good agreement with the observed relations. However, the clusters are too gas rich, suggesting that the AGN feedback model is not efficient enough at expelling gas from the high-redshift progenitors of the clusters. The X-ray properties, such as the spectroscopic temperature and the soft-band luminosity, and the Sunyaev-Zel'dovich properties are in reasonable agreement with the observed relations. However, the clusters have too high central temperatures and larger-than-observed entropy cores, which is likely driven by the AGN feedback after the cluster core has formed. The total metal content and its distribution throughout the ICM are a good match to the observations.

Daily differential emission measure (DEM) distributions of the solar corona are derived from spectra obtained by the Extreme-ultraviolet Variability Experiment (EVE) over a 4-year period starting in 2010 near solar minimum and continuing through the maximum of solar cycle 24. The DEMs are calculated using six strong emission features dominated by Fe lines of charge states VIII, IX, XI, XII, XIV, and XVI that sample the non-flaring coronal temperature range 0.3--5 MK. A proxy for the non-XVIII emission in the wavelength band around the 93.9 \AA\ line is demonstrated. There is little variability in the cool component of the corona (T $<$ 1.3 MK) over the four years, suggesting that the quiet-Sun corona does not respond strongly to the solar cycle, whereas the hotter component (T $>$ 2.0 MK) varies by more than an order of magnitude. A discontinuity in the behavior of coronal diagnostics in 2011 February--March, around the time of the first X-class flare of cycle 24, suggests fundamentally different behavior in the corona under solar minimum and maximum conditions. This global state transition occurs over a period of several months. The DEMs are used to estimate the thermal energy of the visible solar corona (of order $10^{31}$ erg), its radiative energy loss rate (2.5--8 $\times 10^{27}$ erg s$^{-1}$), and the corresponding energy turnover timescale (about an hour). The uncertainties associated with the DEMs and these derived values are mostly due to the coronal Fe abundance and density and the CHIANTI atomic line database.

We present EMERGE, an Empirical ModEl for the foRmation of GalaxiEs, describing the evolution of individual galaxies in large volumes from $z\sim10$ to the present day. We assign a star formation rate to each dark matter halo based on its growth rate, which specifies how much baryonic material becomes available, and the instantaneous baryon conversion efficiency, which determines how efficiently this material is converted to stars, thereby capturing the baryonic physics. Satellites are quenched following the delayed-then-rapid model, and they are tidally disrupted once their subhalo has lost a significant fraction of its mass. The model is constrained with observed data extending out to high redshift. The empirical relations are very flexible, and the model complexity is increased only if required by the data, assessed by several model selection statistics. We find that for the same final halo mass galaxies can have very different star formation histories. Nevertheless, the average star formation and accretion rates are in good agreement with models following an abundance matching strategy. Galaxies that are quenched at $z=0$ typically have a higher peak star formation rate compared to their star-forming counterparts. The accretion of stars can dominate the total mass of massive galaxies, but is insignificant for low-mass systems, independent of star-formation activity. EMERGE predicts stellar-to-halo mass ratios for individual galaxies and introduces scatter self-consistently. We find that at fixed halo mass, passive galaxies have a higher stellar mass on average. The intra-cluster-mass in massive haloes can be up to 8 times larger than the mass of the central galaxy. Clustering for star-forming and quenched galaxies is in good agreement with observational constraints, indicating a realistic assignment of galaxies to haloes.

The Atacama Large Millimeter-submillimeter Array (ALMA) radio telescope has commenced science observations of the Sun starting in late 2016. Since the Sun is much larger than the field of view of individual ALMA dishes, the ALMA interferometer is unable to measure the background level of solar emission when observing the solar disk. The absolute temperature scale is a critical measurement for much of ALMA solar science, including the understanding of energy transfer through the solar atmosphere, the properties of prominences, and the study of shock heating in the chromosphere. In order to provide an absolute temperature scale, ALMA solar observing will take advantage of the remarkable fast-scanning capabilities of the ALMA 12m dishes to make single-dish maps of the full Sun. This article reports on the results of an extensive commissioning effort to optimize the mapping procedure, and it describes the nature of the resulting data. Amplitude calibration is discussed in detail: a path that utilizes the two loads in the ALMA calibration system as well as sky measurements is described and applied to commissioning data. Inspection of a large number of single-dish datasets shows significant variation in the resulting temperatures, and based on the temperature distributions we derive quiet-Sun values at disk center of 7300 K at lambda=3 mm and 5900 K at lambda=1.3 mm. These values have statistical uncertainties of order 100 K, but systematic uncertainties in the temperature scale that may be significantly larger. Example images are presented from two periods with very different levels of solar activity. At a resolution of order 25 arcsec, the 1.3 mm wavelength images show temperatures on the disk that vary over about a 2000 K range.

We report the discovery of a brightness enhancement in the center of a large sunspot umbra at a wavelength of 3 mm using the Atacama Large Millimeter/sub-millimeter Array (ALMA). Sunspots are amongst the most prominent features on the solar surface, but many of their aspects are surprisingly poorly understood. We analyzed a {\lambda}=3 mm (100 GHz) mosaic image obtained by ALMA, which includes a large sunspot within the active region AR12470 on December 16, 2015. The 3 mm map has a field-of-view and spatial resolution, which is the highest spatial-resolution map of an entire sunspot in this frequency range. We find a gradient of 3 mm brightness from a high value in the outer penumbra to a low value in the inner penumbra/outer umbra. Within the inner umbra, there is a marked increase in 3mm brightness temperature, which we call an umbral brightness enhancement. This enhanced emission corresponds to a temperature excess of 800 K relative to the surrounding inner penumbral region and coincides with excess brightness in the 1330 and 1400 {\AA} slitjaw images of the Interface Region Imaging Spectrograph (IRIS), adjacent to a partial lightbridge. This {\lambda}=3 mm brightness enhancement may be an intrinsic feature of the sunspot umbra at chromospheric heights, such as a manifestation of umbral flashes, or it could be related to a coronal plume since the brightness enhancement was coincident with the footpoint of a coronal loop observed at 171 {\AA}.

Eruptive phenomena such as plasmoid ejections or jets are an important feature of solar activity with the potential for improving our understanding of the dynamics of the solar atmosphere. Such ejections are often thought to be signatures of the outflows expected in regions of fast magnetic reconnection. The 304 A EUV line of Helium, formed at around 10^5 K, is found to be a reliable tracer of such phenomena, but the determination of physical parameters from such observations is not straightforward. We have observed a plasmoid ejection from an X-ray bright point simultaneously at millimeter wavelengths with ALMA, at EUV wavelengths with AIA, in soft X-rays with Hinode/XRT. This paper reports the physical parameters of the plasmoid obtained by combining the radio, EUV and X-ray data. As a result, we conclude that the plasmoid can consist either of (approximately) isothermal 10^5 K plasma that is optically thin at 100 GHz, or else a 10^4 K core with a hot envelope. The analysis demonstrates the value of the additional temperature and density constraints that ALMA provides, and future science observations with ALMA will be able to match the spatial resolution of space-borne and other high-resolution telescopes.

The magnetic fields observed in the Milky~Way and nearby galaxies appear to be in equipartition with the turbulent, thermal, and cosmic ray energy densities, and hence are expected to be dynamically important. However, the origin of these strong magnetic fields is still unclear, and most previous attempts to simulate galaxy formation from cosmological initial conditions have ignored them altogether. Here, we analyse the magnetic fields predicted by the simulations of the Auriga Project, a set of 30 high-resolution cosmological zoom simulations of Milky~Way-like galaxies, carried out with a moving-mesh magneto-hydrodynamics code and a detailed galaxy formation physics model. We find that the magnetic fields grow exponentially at early times owing to a small-scale dynamo with an e-folding time of roughly $100\,\rm{Myr}$ in the center of halos until saturation occurs around $z=2-3$, when the magnetic energy density reaches about $10\%$ of the turbulent energy density with a typical strength of $10-50\,\rm{\mu G}$. In the galactic centers the ratio between magnetic and turbulent energy remains nearly constant until $z=0$. At larger radii, differential rotation in the disks leads to linear amplification that typically saturates around $z=0.5$ to $z=0$. The final radial and vertical variations of the magnetic field strength can be well described by two joint exponential profiles, and are in good agreement with observational constraints. Overall, the magnetic fields have only little effect on the global evolution of the galaxies as it takes too long to reach equipartition. We also demonstrate that our results are well converged with numerical resolution.

Observations of the Sun at millimeter and submillimeter wavelengths offer a unique probe into the structure, dynamics, and heating of the chromosphere; the structure of sunspots; the formation and eruption of prominences and filaments; and energetic phenomena such as jets and flares. High-resolution observations of the Sun at millimeter and submillimeter wavelengths are challenging due to the intense, extended, low- contrast, and dynamic nature of emission from the quiet Sun, and the extremely intense and variable nature of emissions associated with energetic phenomena. The Atacama Large Millimeter/submillimeter Array (ALMA) was designed with solar observations in mind. The requirements for solar observations are significantly different from observations of sidereal sources and special measures are necessary to successfully carry out this type of observations. We describe the commissioning efforts that enable the use of two frequency bands, the 3 mm band (Band 3) and the 1.25 mm band (Band 6), for continuum interferometric-imaging observations of the Sun with ALMA. Examples of high-resolution synthesized images obtained using the newly commissioned modes during the solar commissioning campaign held in December 2015 are presented. Although only 30 of the eventual 66 ALMA antennas were used for the campaign, the solar images synthesized from the ALMA commissioning data reveal new features of the solar atmosphere that demonstrate the potential power of ALMA solar observations. The ongoing expansion of ALMA and solar-commissioning efforts will continue to enable new and unique solar observing capabilities.

We characterize the contribution from accreted material to the galactic discs of the Auriga Project, a set of high resolution magnetohydrodynamic cosmological simulations of late-type galaxies performed with the moving-mesh code AREPO. Our goal is to explore whether a significant accreted (or ex-situ) stellar component in the Milky Way disc could be hidden within the near-circular orbit population, which is strongly dominated by stars born in-situ. One third of our models shows a significant ex-situ disc but this fraction would be larger if constraints on orbital circularity were relaxed. Most of the ex-situ material ($\gtrsim 50\%$) comes from single massive satellites ($> 6 \times 10^{10}~M_{\odot}$). These satellites are accreted with a wide range of infall times and inclination angles (up to $85^{\circ}$). Ex-situ discs are thicker, older and more metal-poor than their in-situ counterparts. They show a flat median age profile, which differs from the negative gradient observed in the in-situ component. As a result, the likelihood of identifying an ex-situ disc in samples of old stars on near-circular orbits increases towards the outskirts of the disc. We show three examples that, in addition to ex-situ discs, have a strongly rotating dark matter component. Interestingly, two of these ex-situ stellar discs show an orbital circularity distribution that is consistent with that of the in-situ disc. Thus, they would not be detected in typical kinematic studies.

Based on five years of observations with the 40-foot telescope at Green Bank Observatory (GBO), Reichart & Stephens (2000) found that the radio source Cassiopeia A had either faded more slowly between the mid-1970s and late 1990s than Baars et al. (1977) had found it to be fading between the late 1940s and mid-1970s, or that it had rebrightened and then resumed fading sometime between the mid-1970s and mid-1990s, in L band (1.4 GHz). Here, we present 15 additional years of observations of Cas A and Cyg A with the 40-foot in L band, and three and a half additional years of observations of Cas A, Cyg A, Tau A, and Vir A with GBO's recently refurbished 20-meter telescope in L and X (9 GHz) bands. We also present a more sophisticated analysis of the 40-foot data, and a reanalysis of the Baars et al. (1977) data, which reveals small, but non-negligible differences. We find that overall, between the late 1950s and late 2010s, Cas A faded at an average rate of $0.670 \pm 0.019$ %/yr in L band, consistent with Reichart & Stephens (2000). However, we also find, at the 6.3$\sigma$ credible level, that it did not fade at a constant rate. Rather, Cas A faded at a faster rate through at least the late 1960s, rebrightened (or at least faded at a much slower rate), and then resumed fading at a similarly fast rate by, at most, the late 1990s. Given these differences from the original Baars et al. (1977) analysis, and given the importance of their fitted spectral and temporal models for flux-density calibration in radio astronomy, we update and improve on these models for all four of these radio sources. In doing so, we additionally find that Tau A is fading at a rate of $0.102^{+0.042}_{-0.043}$ %/yr in L band.

We introduce the Hydrangea simulations, a suite of 24 cosmological hydrodynamic zoom-in simulations of massive galaxy clusters (M_200c = 10^14-10^15 M_Sun) with baryon particle masses of ~10^6 M_Sun. Designed to study the impact of the cluster environment on galaxy formation, they are a key part of the `Cluster-EAGLE' project (Barnes et al. 2017). They use a galaxy formation model developed for the EAGLE project, which has been shown to yield both realistic field galaxies and hot gas fractions of galaxy groups consistent with observations. The total stellar mass content of the simulated clusters agrees with observations, but central cluster galaxies are too massive, by up to 0.6 dex. Passive satellite fractions are higher than in the field, and at stellar masses Mstar > 10^10 M_Sun this environmental effect is quantitatively consistent with observations. The predicted satellite stellar mass function matches data from local cluster surveys. Normalized to total mass, there are fewer low-mass (Mstar < 10^10 M_Sun) galaxies within the virial radius of clusters than in the field, primarily due to star formation quenching. Conversely, the simulations predict an overabundance of massive galaxies in clusters compared to the field that persists to their far outskirts (> 5r_200c). This is caused by a significantly increased stellar mass fraction of (sub-)haloes in the cluster environment, by up to ~0.3 dex even well beyond r_200c. Haloes near clusters are also more concentrated than equally massive field haloes, but these two effects are largely uncorrelated.

We use state-of-the-art, three-dimensional non-local thermodynamic equilibrium (non-LTE) radiative magnetohydrodynamic simulations of the quiet solar atmosphere to carry out detailed tests of chromospheric magnetic field diagnostics from free-free radiation at millimeter and submillimeter wavelengths (mm/submm). The vertical component of the magnetic field was deduced from the mm/submm brightness spectra and the degree of circular polarization synthesized at millimeter frequencies. We used the frequency bands observed by the Atacama Large Millimeter/Submillimeter Array (ALMA) as a convenient reference. The magnetic field maps obtained describe the longitudinal magnetic field at the effective formation heights of the relevant wavelengths in the solar chromosphere. The comparison of the deduced and model chromospheric magnetic fields at the spatial resolution of both the model and current observations demonstrates a good correlation, but has a tendency to underestimate the model field. The systematic discrepancy of about 10 percent is probably due to averaging of the restored field over the heights contributing to the radiation, weighted by the strength of the contribution. On the whole, the method of probing the longitudinal component of the magnetic field with free-free emission at mm/submm wavelengths is found to be applicable to measurements of the weak quiet-Sun magnetic fields. However, successful exploitation of this technique requires very accurate measurements of the polarization properties (primary beam and receiver polarization response) of the antennas, which will be the principal factor that determines the level to which chromospheric magnetic fields can be measured. Consequently, high-resolution and high-precision observations of circularly polarized radiation at millimeter wavelengths can be a powerful tool for producing chromospheric longitudinal magnetograms.

In order to understand the role of radio-quiet quasars (RQQs) in galaxy evolution, we must determine the relative levels of accretion and star-formation activity within these objects. Previous work at low radio flux-densities has shown that accretion makes a significant contribution to the total radio emission, in contrast with other quasar studies that suggest star formation dominates. To investigate, we use 70 RQQs from the Spitzer-Herschel Active Galaxy Survey. These quasars are all at $z$ ~ 1, thereby minimising evolutionary effects, and have been selected to span a factor of ~100 in optical luminosity, so that the luminosity dependence of their properties can be studied. We have imaged the sample using the Karl G. Jansky Very Large Array (JVLA), whose high sensitivity results in 35 RQQs being detected above 2 $\sigma$. This radio dataset is combined with far-infrared luminosities derived from grey-body fitting to Herschel photometry. By exploiting the far-infrared--radio correlation observed for star-forming galaxies, and comparing two independent estimates of the star-formation rate, we show that star formation alone is not sufficient to explain the total radio emission. Considering RQQs above a 2-$\sigma$ detection level in both the radio and the far-infrared, 92 per cent are accretion-dominated, and the accretion process accounts for 80 per cent of the radio luminosity when summed across the objects. The radio emission connected with accretion appears to be correlated with the optical luminosity of the RQQ, whilst a weaker luminosity-dependence is evident for the radio emission connected with star formation.

We present observations of the occulted active region AR12222 during the third {\em NuSTAR} solar campaign on 2014 December 11, with concurrent {\em SDO/}AIA and {\em FOXSI-2} sounding rocket observations. The active region produced a medium size solar flare one day before the observations, at $\sim18$UT on 2014 December 10, with the post-flare loops still visible at the time of {\em NuSTAR} observations. The time evolution of the source emission in the {\em SDO/}AIA $335\textrm{\AA}$ channel reveals the characteristics of an extreme-ultraviolet late phase event, caused by the continuous formation of new post-flare loops that arch higher and higher in the solar corona. The spectral fitting of {\em NuSTAR} observations yields an isothermal source, with temperature $3.8-4.6$ MK, emission measure $0.3-1.8 \times 10^{46}\textrm{ cm}^{-3}$, and density estimated at $2.5-6.0 \times 10^8 \textrm{ cm}^{-3}$. The observed AIA fluxes are consistent with the derived {\em NuSTAR} temperature range, favoring temperature values in the range $4.0-4.3$ MK. By examining the post-flare loops' cooling times and energy content, we estimate that at least 12 sets of post-flare loops were formed and subsequently cooled between the onset of the flare and {\em NuSTAR} observations, with their total thermal energy content an order of magnitude larger than the energy content at flare peak time. This indicates that the standard approach of using only the flare peak time to derive the total thermal energy content of a flare can lead to a large underestimation of its value.

We analyse the properties of the HI gas distribution in the Auriga project, a set of magnetohydrodynamic cosmological simulations performed with the moving-mesh code AREPO and a physics model for galaxy formation that succeeds in forming realistic late-type galaxies in the 30 Milky Way-sized haloes simulated in this project. We use a simple approach to estimate the neutral hydrogen fraction in our simulation set, which treats low-density and star-forming gas separately, and we explore two different prescriptions to subtract the contribution of molecular hydrogen from the total HI content. The HI gas in the vast majority of the systems forms extended discs although more disturbed morphologies are present. Notwithstanding the general good agreement with observed HI properties -- such as radial profiles and the mass-diameter relation -- the Auriga galaxies are systematically larger and more gas-rich than typical nearby galaxies. Interestingly, the amount of HI gas outside the disc plane correlates with the star formation rate, consistent with a picture where most of this extra-planar HI gas originates from a fountain-like flow. Our findings are robust with respect to the different assumptions adopted for computing the molecular hydrogen fraction and do not vary significantly over a wide range of numerical resolution. The HI modelling introduced in this paper can be used in future work to build artificial interferometric HI data cubes, allowing an even closer comparison of the gas dynamics in simulated galaxies with observations.

We introduce a suite of thirty cosmological magneto-hydrodynamical zoom simulations of the formation of galaxies in isolated Milky Way mass dark haloes. These were carried out with the moving mesh code \textlcsc{AREPO}, together with a comprehensive model for galaxy formation physics, including AGN feedback and magnetic fields, which produces realistic galaxy populations in large cosmological simulations. We demonstrate that our simulations reproduce a wide range of present-day observables, in particular, two component disc dominated galaxies with appropriate stellar masses, sizes, rotation curves, star formation rates and metallicities. We investigate the driving mechanisms that set present-day disc sizes/scale lengths, and find that they are related to the angular momentum of halo material. We show that the largest discs are produced by quiescent mergers that inspiral into the galaxy and deposit high angular momentum material into the pre-existing disc, simultaneously increasing the spin of dark matter and gas in the halo. More violent mergers and strong AGN feedback play roles in limiting disc size by destroying pre-existing discs and by suppressing gas accretion onto the outer disc, respectively. The most important factor that leads to compact discs, however, is simply a low angular momentum for the halo. In these cases, AGN feedback plays an important role in limiting central star formation and the formation of a massive bulge.

How impulsive magnetic energy release leads to solar eruptions and how those eruptions are energized and evolve are vital unsolved problems in Heliophysics. The standard model for solar eruptions summarizes our current understanding of these events. Magnetic energy in the corona is released through drastic restructuring of the magnetic field via reconnection. Electrons and ions are then accelerated by poorly understood processes. Theories include contracting loops, merging magnetic islands, stochastic acceleration, and turbulence at shocks, among others. Although this basic model is well established, the fundamental physics is poorly understood. HXR observations using grazing-incidence focusing optics can now probe all of the key regions of the standard model. These include two above-the-looptop (ALT) sources which bookend the reconnection region and are likely the sites of particle acceleration and direct heating. The science achievable by a direct HXR imaging instrument can be summarized by the following science questions and objectives which are some of the most outstanding issues in solar physics (1) How are particles accelerated at the Sun? (1a) Where are electrons accelerated and on what time scales? (1b) What fraction of electrons is accelerated out of the ambient medium? (2) How does magnetic energy release on the Sun lead to flares and eruptions? A Focusing Optics X-ray Solar Imager (FOXSI) instrument, which can be built now using proven technology and at modest cost, would enable revolutionary advancements in our understanding of impulsive magnetic energy release and particle acceleration, a process which is known to occur at the Sun but also throughout the Universe.

The absolute brightness temperature of the Sun at millimeter wavelengths is an important diagnostic of the solar chromosphere. Because the Sun is so bright, measurement of this property usually involves the operation of telescopes under extreme conditions and requires a rigorous performance assessment of the telescope. In this study, we establish solar observation and calibration techniques at 2.6-mm wavelength for the Nobeyama 45-m telescope and derive the absolute solar brightness temperature accurately. We tune the superconductor-insulator-superconductor (SIS) receiver by inducing different bias voltages onto the SIS mixer to prevent saturation. Then, we examine the linearity of the receiver system by comparing outputs derived from different tuning conditions. Further, we measure the lunar filled beam efficiency of the telescope using the New Moon, and then derive the absolute brightness temperature of the Sun. The derived solar brightness temperature is 7700+-310 K at 115 GHz. The telescope beam pattern is modeled as a summation of three Gaussian functions and derived using the solar limb. The real shape of the Sun is determined via deconvolution of the beam pattern from the observed map. Such well-calibrated single-dish observations are important for high-resolution chromospheric studies because they provide the absolute temperature scale missing from interferometer observations.

The Sculptor dwarf spheroidal galaxy contains two distinct stellar populations: one metal-rich and the other metal-poor. Several authors have argued that in order for these two populations to reside in the same gravitational potential, the dark matter halo must have a core similar to that observed in the stellar count profile. This would rule out the cuspy Navarro-Frenk-White (NFW) density profiles predicted for halos and subhalos by dark matter only simulations of the Lambda Cold Dark Matter (Lambda-CDM) cosmological model. We present a new theoretical framework to analyse stellar count and velocity observations in a self-consistent manner based on separable models, f(E,J)=g(J)h(E), for the distribution function of an equilibrium spherical system. We use this machinery to analyse available photometric and kinematic data for the two stellar populations in Sculptor. We find, contrary to some previous claims, that the data are consistent with populations in equilibrium within an NFW dark matter potential with structural parameters in the range expected in Lambda-CDM. Our solution gives a maximum circular velocity for Sculptor between 20 and 35 km/s. We discuss why some previous authors came to a different conclusion.

Recent studies have revealed an oscillating asymmetry in the vertical structure of the Milky Way's disc. Here we analyze 16 high-resolution, fully cosmological simulations of the evolution of individual Milky Way-sized galaxies, carried out with the MHD code AREPO. At redshift zero, about $70\%$ of our galactic discs show strong vertical patterns, with amplitudes that can exceed 2 kpc. Half of these are typical `integral sign' warps. The rest are oscillations similar to those observed in the Milky Way. Such structures are thus expected to be common. The associated mean vertical motions can be as large as 30~km/s. Cold disc gas typically follows the vertical patterns seen in the stars. These perturbations have a variety of causes: close encounters with satellites, distant flybys of massive objects, accretion of misaligned cold gas from halo infall or from mergers. Tidally induced vertical patterns can be identified in both young and old stellar populations, whereas those originating from cold gas accretion are seen mainly in the younger populations. Galaxies with regular or at most weakly perturbed discs are usually, but not always, free from recent interactions with massive companions, although we have one case where an equilibrium compact disc reforms after a merger.

We study the quenching of star formation as a function of redshift, environment and stellar mass in the galaxy formation simulations of Henriques et al. (2015), which implement an updated version of the Munich semi-analytic model (L-GALAXIES) on the two Millennium Simulations after scaling to a Planck cosmology. In this model massive galaxies are quenched by AGN feedback depending on both black hole and hot gas mass, and hence indirectly on stellar mass. In addition, satellite galaxies of any mass can be quenched by ram-pressure or tidal stripping of gas and through the suppression of gaseous infall. This combination of processes produces quenching efficiencies which depend on stellar mass, host halo mass, environment density, distance to group centre and group central galaxy properties in ways which agree qualitatively with observation. Some discrepancies remain in dense regions and close to group centres, where quenching still seems too efficient. In addition, although the mean stellar age of massive galaxies agrees with observation, the assumed AGN feedback model allows too much ongoing star formation at late times. The fact that both AGN feedback and environmental effects are stronger in higher density environments leads to a correlation between the quenching of central and satellite galaxies which roughly reproduces observed conformity trends inside haloes.

We study the evolution of spectral early-type galaxies in clusters, groups and the field up to redshift 0.9 using the EDisCS dataset. We measure Re, Ie, and sigma for 154 cluster and 68 field galaxies. We study the evolution of the zero point of the fundamental plane (FP) and confirm results in the literature, but now also for the low cluster velocity dispersion regime. The mass-to-light ratio varies as Delta log M/L_B=(-0.54+-0.01)z=(-1.61+-0.01)log(1+z) in clusters, independent of their velocity dispersion. The evolution is stronger (Delta log M/L_B=(-0.76+-0.01)z=(-2.27+-0.03)log(1+z)) for field galaxies. The FP residuals correlate with galaxy mass and become progressively negative at low masses. The effect is visible at z>=0.7 for cluster galaxies and at z>=0.5 for field galaxies. We investigate the size evolution of our galaxy sample. We find that the half-luminosity radius for a galaxy with a dynamical or stellar mass of 2x10^11 Msol varies as (1+z)^{-1.0+-0.3} for both cluster and field galaxies. At the same time, stellar velocity dispersions grow with redshift, as (1+z)^{0.59+-0.10} at constant dynamical mass, and as (1+z)^{0.34+- 0.14} at constant stellar mass. The measured size evolution reduces to Re (1+z)^{-0.5+- 0.2} and sigma (1+z)^{0.41+-0.08}, at fixed dynamical masses, and Re (1+z)^{-0.68+-0.4} and sigma (1+z)^{0.19+-0.10}, at fixed stellar masses, when the progenitor bias (galaxies that locally are of spectroscopic early-type, but not very old, disappear from the EDisCS high-redshift sample; these galaxies tend to be large in size) is taken into account. Taken together, the variations in size and velocity dispersion imply that the luminosity evolution with redshift derived from the zero point of the FP is somewhat milder than that derived without taking these variations into account.

We analyse the mass assembly of central galaxies in the EAGLE hydrodynamical simulations. We build merger trees to connect galaxies to their progenitors at different redshifts and characterize their assembly histories by focusing on the time when half of the galaxy stellar mass was assembled into the main progenitor. We show that galaxies with stellar mass $M_*<10^{10.5}M_{\odot}$ assemble most of their stellar mass through star formation in the main progenitor (`in-situ' star formation). This can be understood as a consequence of the steep rise in star formation efficiency with halo mass for these galaxies. For more massive galaxies, however, an increasing fraction of their stellar mass is formed outside the main progenitor and subsequently accreted. Consequently, while for low-mass galaxies the assembly time is close to the stellar formation time, the stars in high-mass galaxies typically formed long before half of the present-day stellar mass was assembled into a single object, giving rise to the observed anti-hierarchical downsizing trend. In a typical present-day $M_*\geq10^{11}M_{\odot}$ galaxy, around $20\%$ of the stellar mass has an external origin. This fraction decreases with increasing redshift. Bearing in mind that mergers only make an important contribution to the stellar mass growth of massive galaxies, we find that the dominant contribution comes from mergers with galaxies of mass greater than one tenth of the main progenitor's mass. The galaxy merger fraction derived from our simulations agrees with recent observational estimates.

Roughly half of all stars reside in galaxies without significant ongoing star formation. However, galaxy formation models indicate that it is energetically challenging to suppress the cooling of gas and the formation of stars in galaxies that lie at the centers of their dark matter halos. In this Letter, we show that the dependence of quiescence on black hole and stellar mass is a powerful discriminant between differing models for the mechanisms that suppress star formation. Using observations of 91 star-forming and quiescent central galaxies with directly-measured black hole masses, we find that quiescent galaxies host more massive black holes than star-forming galaxies with similar stellar masses. This observational result is in qualitative agreement with models that assume that effective, more-or-less continuous AGN feedback suppresses star formation, strongly suggesting the importance of the black hole in producing quiescence in central galaxies.

The predicted abundance and properties of the low-mass substructures embedded inside larger dark matter haloes differ sharply among alternative dark matter models. Too small to host galaxies themselves, these subhaloes may still be detected via gravitational lensing, or via perturbations of the Milky Way's globular cluster streams and its stellar disk. Here we use the Apostle cosmological simulations to predict the abundance and the spatial and velocity distributions of subhaloes in the range 10^6.5-10^8.5 solar masses inside haloes of mass ~ 10^12 solar masses in LCDM. Although these subhaloes are themselves devoid of baryons, we find that baryonic effects are important. Compared to corresponding dark matter only simulations, the loss of baryons from subhaloes and stronger tidal disruption due to the presence of baryons near the centre of the main halo, reduce the number of subhaloes by ~ 1/4 to 1/2, independently of subhalo mass, but increasingly towards the host halo centre. We also find that subhaloes have non-Maxwellian orbital velocity distributions, with centrally rising velocity anisotropy and positive velocity bias which reduces the number of low-velocity subhaloes, particularly near the halo centre. We parameterise the predicted population of subhaloes in terms of mass, galactocentric distance, and velocities. We discuss implications of our results for the prospects of detecting dark matter substructures and for possible inferences about the nature of dark matter.

We present the public data release of halo and galaxy catalogues extracted from the EAGLE suite of cosmological hydrodynamical simulations of galaxy formation. These simulations were performed with an enhanced version of the GADGET code that includes a modified hydrodynamics solver, time-step limiter and subgrid treatments of baryonic physics, such as stellar mass loss, element-by-element radiative cooling, star formation and feedback from star formation and black hole accretion. The simulation suite includes runs performed in volumes ranging from 25 to 100 comoving megaparsecs per side, with numerical resolution chosen to marginally resolve the Jeans mass of the gas at the star formation threshold. The free parameters of the subgrid models for feedback are calibrated to the redshift z=0 galaxy stellar mass function, galaxy sizes and black hole mass - stellar mass relation. The simulations have been shown to match a wide range of observations for present-day and higher-redshift galaxies. The raw particle data have been used to link galaxies across redshifts by creating merger trees. The indexing of the tree produces a simple way to connect a galaxy at one redshift to its progenitors at higher redshift and to identify its descendants at lower redshift. In this paper we present a relational database which we are making available for general use. A large number of properties of haloes and galaxies and their merger trees are stored in the database, including stellar masses, star formation rates, metallicities, photometric measurements and mock gri images. Complex queries can be created to explore the evolution of more than 10^5 galaxies, examples of which are provided in appendix. (abridged)

Radio emission observations from the Learmonth and Bruny Island radio spectrographs are analysed to determine the nature of a train of discrete, periodic radio \lq sparks\rq (finite-bandwidth, short-duration isolated radio features) which precede a type II burst. We analyse extreme ultraviolet (EUV) imaging from SDO/AIA at multiple wavelengths and identify a series of quasi-periodic rapidly-propagating enhancements, which we interpret as a fast wave train, and link these to the detected radio features. The speeds and positions of the periodic rapidly propagating fast waves and the coronal mass ejection (CME) were recorded using running-difference images and time-distance analysis. From the frequency of the radio sparks the local electron density at the emission location was estimated for each. Using an empirical model for the scaling of density in the corona, the calculated electron density was used to obtain the height above the surface at which the emission occurs, and the propagation velocity of the emission location. The period of the radio sparks, $\delta t_\mathrm{r}$=1.78$\pm$0.04 min, matches the period of the fast wave train observed at 171 $\AA$, $\delta t_\mathrm{EUV}$=1.7 $\pm$ 0.2 min. The inferred speed of the emission location of the radio sparks, 630 km s$^{-1}$, is comparable to the measured speed of the CME leading edge, 500 km s$^{-1}$, and the speeds derived from the drifting of the type II lanes. The calculated height of the radio emission (obtained from the density) matches the observed location of the CME leading edge. From the above evidence we propose that the radio sparks are caused by the quasi-periodic fast waves, and the emission is generated as they catch up and interact with the leading edge of the CME.

We study radial perturbations of a wormhole in $R^2$ gravity to determine regions of stability. We also investigate massive and massless particle orbits and tidal forces in this space-time for a radially infalling observer.

As a result of their internal dynamical coherence, thin stellar streams formed by disrupting globular clusters (GCs) can act as detectors of dark matter (DM) substructure in the Galactic halo. Perturbations induced by close flybys amplify into detectable density gaps, providing a probe both of the abundance and of the masses of DM subhaloes. Here, we use N-body simulations to show that the Galactic population of giant molecular clouds (GMCs) can also produce gaps (and clumps) in GC streams, and so may confuse the detection of DM subhaloes. We explore the cases of streams analogous to the observed Palomar 5 and GD1 systems, quantifying the expected incidence of structure caused by GMC perturbations. Deep observations should detect such disturbances regardless of the substructure content of the Milky Way's halo. Detailed modelling will be needed to demonstrate that any detected gaps or clumps were produced by DM subhaloes rather than by molecular clouds.

We study the vertical structure of a stellar disk obtained from a fully cosmological high-resolution hydrodynamical simulation of the formation of a Milky Way-like galaxy. At the present day, the disk's mean vertical height shows a well-defined and strong pattern, with amplitudes as large as 3 kpc in its outer regions. This pattern is the result of a satellite - host halo - disk interaction and reproduces, qualitatively, many of the observable properties of the Monoceros Ring. In particular we find disk material at the distance of Monoceros ($R \sim$ 12-16 kpc, galactocentric) extending far above the mid plane ( 30$^{\circ}$, $\langle Z \rangle \sim$ 1-2 kpc) in both hemispheres, as well as well-defined arcs of disk material at heliocentric distances $\gtrsim 5$ kpc. The pattern was first excited $\approx 3$ Gyr ago as an $m=1$ mode that later winds up into a leading spiral pattern. Interestingly, the main driver behind this perturbation is a low-mass low-velocity fly-by encounter. The satellite has total mass, pericentre distance and pericentric velocity of $\sim 5\%$ of the host, $\sim 80$ kpc, and 215 km/s, respectively. The satellite is not massive enough to directly perturb the galactic disk but we show that the density field of the host dark matter halo responds to this interaction resulting in a strong amplification of the perturbative effects. This subsequently causes the onset and development of the Monoceros-like feature.

We present results based on full-mission Planck observations of temperature and polarization anisotropies of the CMB. These data are consistent with the six-parameter inflationary LCDM cosmology. From the Planck temperature and lensing data, for this cosmology we find a Hubble constant, H0= (67.8 +/- 0.9) km/s/Mpc, a matter density parameter Omega_m = 0.308 +/- 0.012 and a scalar spectral index with n_s = 0.968 +/- 0.006. (We quote 68% errors on measured parameters and 95% limits on other parameters.) Combined with Planck temperature and lensing data, Planck LFI polarization measurements lead to a reionization optical depth of tau = 0.066 +/- 0.016. Combining Planck with other astrophysical data we find N_ eff = 3.15 +/- 0.23 for the effective number of relativistic degrees of freedom and the sum of neutrino masses is constrained to < 0.23 eV. Spatial curvature is found to be |Omega_K| < 0.005. For LCDM we find a limit on the tensor-to-scalar ratio of r <0.11 consistent with the B-mode constraints from an analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP data leads to a tighter constraint of r < 0.09. We find no evidence for isocurvature perturbations or cosmic defects. The equation of state of dark energy is constrained to w = -1.006 +/- 0.045. Standard big bang nucleosynthesis predictions for the Planck LCDM cosmology are in excellent agreement with observations. We investigate annihilating dark matter and deviations from standard recombination, finding no evidence for new physics. The Planck results for base LCDM are in agreement with BAO data and with the JLA SNe sample. However the amplitude of the fluctuations is found to be higher than inferred from rich cluster counts and weak gravitational lensing. Apart from these tensions, the base LCDM cosmology provides an excellent description of the Planck CMB observations and many other astrophysical data sets.

We present results from the the first campaign of dedicated solar observations undertaken by the \textit{Nuclear Spectroscopic Telescope ARray} ({\em NuSTAR}) hard X-ray telescope. Designed as an astrophysics mission, {\em NuSTAR} nonetheless has the capability of directly imaging the Sun at hard X-ray energies ($>$3~keV) with an increase in sensitivity of at least two magnitude compared to current non-focusing telescopes. In this paper we describe the scientific areas where \textit{NuSTAR} will make major improvements on existing solar measurements. We report on the techniques used to observe the Sun with \textit{NuSTAR}, their limitations and complications, and the procedures developed to optimize solar data quality derived from our experience with the initial solar observations. These first observations are briefly described, including the measurement of the Fe K-shell lines in a decaying X-class flare, hard X-ray emission from high in the solar corona, and full-disk hard X-ray images of the Sun.

We have used the Karl G. Jansky Very Large Array to image ~100 sq. deg. of SDSS Stripe 82 at 1-2 GHz. The survey consists of 1,026 snapshot observations of 2.5 minutes duration, using the hybrid CnB configuration. The survey has good sensitivity to diffuse, low surface brightness structures and extended radio emission, making it highly synergistic with existing 1.4 GHz radio observations of the region. The principal data products are continuum images, with 16 x 10 arcsecond resolution, and a catalogue containing 11,782 point and Gaussian components resulting from fits to the thresholded Stokes-I brightness distribution, forming approximately 8,948 unique radio sources. The typical effective 1{\sigma} noise level is 88 {\mu}Jy / beam. Spectral index estimates are included, as derived from the 1 GHz of instantaneous bandwidth. Astrometric and photometric accuracy are in excellent agreement with existing narrowband observations. A large-scale simulation is used to investigate clean bias, which we extend into the spectral domain. Clean bias remains an issue for snapshot surveys with the VLA, affecting our total intensity measurements at the ~1{\sigma} level. Statistical spectral index measurements are in good agreement with existing measurements derived from matching separate surveys at two frequencies. At flux densities below ~35{\sigma} the median in-band spectral index measurements begin to exhibit a bias towards flatness that is dependent on both flux density and the intrinsic spectral index. In-band spectral curvature measurements are likely to be unreliable for all but the very brightest components. Image products and catalogues are publicly available via an FTP server.

The GaLactic and Extragalactic All-sky MWA survey (GLEAM) is a new relatively low resolution, contiguous 72-231 MHz survey of the entire sky south of declination +25 deg. In this paper, we outline one approach to determine the relative contribution of system noise, classical confusion and sidelobe confusion in GLEAM images. An understanding of the noise and confusion properties of GLEAM is essential if we are to fully exploit GLEAM data and improve the design of future low-frequency surveys. Our early results indicate that sidelobe confusion dominates over the entire frequency range, implying that enhancements in data processing have the potential to further reduce the noise.

This paper outlines how the new GaLactic and Extragalactic All-sky MWA Survey (GLEAM, Wayth et al. 2015), observed by the Murchison Widefield Array covering the frequency range 72 - 231 MHz, allows identification of a new large, complete, sample of more than 2000 bright extragalactic radio sources selected at 151 MHz. With a flux density limit of 4 Jy this sample is significantly larger than the canonical fully-complete sample, 3CRR (Laing, Riley & Longair 1983). In analysing this small bright subset of the GLEAM survey we are also providing a first user check of the GLEAM catalogue ahead of its public release (Hurley-Walker et al. in prep). Whilst significant work remains to fully characterise our new bright source sample, in time it will provide important constraints to evolutionary behaviour, across a wide redshift and intrinsic radio power range, as well as being highly complementary to results from targeted, small area surveys.

This document was created by the Solar Simulations for the Atacama Large Millimeter Observatory Network (SSALMON) in preparation of the first regular observations of the Sun with the Atacama Large Millimeter/submillimeter Array (ALMA), which are anticipated to start in ALMA Cycle 4 in October 2016. The science cases presented here demonstrate that a large number of scientifically highly interesting observations could be made already with the still limited solar observing modes foreseen for Cycle 4 and that ALMA has the potential to make important contributions to answering long-standing scientific questions in solar physics. With the proposal deadline for ALMA Cycle 4 in April 2016 and the Commissioning and Science Verification campaign in December 2015 in sight, several of the SSALMON Expert Teams composed strategic documents in which they outlined potential solar observations that could be feasible given the anticipated technical capabilities in Cycle 4. These documents have been combined and supplemented with an analysis, resulting in recommendations for solar observing with ALMA in Cycle 4. In addition, the detailed science cases also demonstrate the scientific priorities of the solar physics community and which capabilities are wanted for the next observing cycles. The work on this White Paper effort was coordinated in close cooperation with the two international solar ALMA development studies led by T. Bastian (NRAO, USA) and R. Brajsa, (ESO). This document will be further updated until the beginning of Cycle 4 in October 2016. In particular, we plan to adjust the technical capabilities of the solar observing modes once finally decided and to further demonstrate the feasibility and scientific potential of the included science cases by means of numerical simulations of the solar atmosphere and corresponding simulated ALMA observations.

We use the semi-analytic model developed by Henriques et al. (2015) to explore the origin of star formation history diversity for galaxies that lie at the centre of their dark matter haloes and have present-day stellar masses in the range 5-8 $\times$ 10$^{10}$ M$_{\odot}$, similar to that of the Milky Way. In this model, quenching is the dominant physical mechanism for introducing scatter in the growth histories of these galaxies. We find that present-day quiescent galaxies have a larger variety of growth histories than star-formers since they underwent 'staggered quenching' - a term describing the correlation between the time of quenching and present-day halo mass. While halo mass correlates broadly with quiescence, we find that quiescence is primarily a function of black hole mass, where galaxies quench when heating from their active galactic nuclei becomes sufficient to offset the redshift-dependent cooling rate. In this model, the emergence of a prominent quiescent population is the main process that flattens the stellar mass-halo mass relation at mass scales at or above that of the Milky Way.

We present the first observations of quiescent active regions (ARs) using NuSTAR, a focusing hard X-ray telescope capable of studying faint solar emission from high temperature and non-thermal sources. We analyze the first directly imaged and spectrally resolved X-rays above 2~keV from non-flaring ARs, observed near the west limb on 2014 November 1. The NuSTAR X-ray images match bright features seen in extreme ultraviolet and soft X-rays. The NuSTAR imaging spectroscopy is consistent with isothermal emission of temperatures $3.1-4.4$~MK and emission measures $1-8\times 10^{46}$~cm$^{-3}$. We do not observe emission above 5~MK but our short effective exposure times restrict the spectral dynamic range. With few counts above 6~keV, we can place constraints on the presence of an additional hotter component between 5 and 12~MK of $\sim 10^{46}$cm$^{-3}$ and $\sim 10^{43}$ cm$^{-3}$, respectively, at least an order of magnitude stricter than previous limits. With longer duration observations and a weakening solar cycle (resulting in an increased livetime), future NuSTAR observations will have sensitivity to a wider range of temperatures as well as possible non-thermal emission.

We introduce methods which allow observed galaxy clustering to be used together with observed luminosity or stellar mass functions to constrain the physics of galaxy formation. We show how the projected two-point correlation function of galaxies in a large semi-analytic simulation can be estimated to better than ~10% using only a very small subsample of the subhalo merger trees. This allows measured correlations to be used as constraints in a Monte Carlo Markov Chain exploration of the astrophysical and cosmological parameter space. An important part of our scheme is an analytic profile which captures the simulated satellite distribution extremely well out to several halo virial radii. This is essential to reproduce the correlation properties of the full simulation at intermediate separations. As a first application, we use low-redshift clustering and abundance measurements to constrain a recent version of the Munich semi-analytic model. The preferred values of most parameters are consistent with those found previously, with significantly improved constraints and somewhat shifted "best" values for parameters that primarily affect spatial distributions. Our methods allow multi-epoch data on galaxy clustering and abundance to be used as joint constraints on galaxy formation. This may lead to significant constraints on cosmological parameters even after marginalising over galaxy formation physics.

The scaling of disk galaxy rotation velocity with baryonic mass (the "Baryonic Tully-Fisher" relation, BTF) has long confounded galaxy formation models. It is steeper than the M ~ V^3 scaling relating halo virial masses and circular velocities and its zero point implies that galaxies comprise a very small fraction of available baryons. Such low galaxy formation efficiencies may in principle be explained by winds driven by evolving stars, but the tightness of the BTF relation argues against the substantial scatter expected from such vigorous feedback mechanism. We use the APOSTLE/EAGLE simulations to show that the BTF relation is well reproduced in LCDM simulations that match the size and number of galaxies as a function of stellar mass. In such models, galaxy rotation velocities are proportional to halo virial velocity and the steep velocity-mass dependence results from the decline in galaxy formation efficiency with decreasing halo mass needed to reconcile the CDM halo mass function with the galaxy luminosity function. Despite the strong feedback, the scatter in the simulated BTF is smaller than observed, even when considering all simulated galaxies and not just rotationally-supported ones. The simulations predict that the BTF should become increasingly steep at the faint end, although the velocity scatter at fixed mass should remain small. Observed galaxies with rotation speeds below ~40 km/s seem to deviate from this prediction. We discuss observational biases and modeling uncertainties that may help to explain this disagreement in the context of LCDM models of dwarf galaxy formation.

We use galaxy-galaxy lensing to study the dark matter halos surrounding a sample of Locally Brightest Galaxies (LBGs) selected from the Sloan Digital Sky Survey. We measure mean halo mass as a function of the stellar mass and colour of the central galaxy. Mock catalogues constructed from semi-analytic galaxy formation simulations demonstrate that most LBGs are the central objects of their halos, greatly reducing interpretation uncertainties due to satellite contributions to the lensing signal. Over the full stellar mass range, $10.3 < \log [M_*/M_\odot] < 11.6$, we find that passive central galaxies have halos that are at least twice as massive as those of star-forming objects of the same stellar mass. The significance of this effect exceeds $3\sigma$ for $\log [M_*/M_\odot] > 10.7$. Tests using the mock catalogues and on the data themselves clarify the effects of LBG selection and show that it cannot artificially induce a systematic dependence of halo mass on LBG colour. The bimodality in halo mass at fixed stellar mass is reproduced by the astrophysical model underlying our mock catalogue, but the sign of the effect is inconsistent with recent, nearly parameter-free age-matching models. The sign and magnitude of the effect can, however, be reproduced by halo occupation distribution models with a simple (few-parameter) prescription for type-dependence.

The Atacama Large Millimeter/submillimeter Array (ALMA) is a new powerful tool for observing the Sun at high spatial, temporal, and spectral resolution. These capabilities can address a broad range of fundamental scientific questions in solar physics. The radiation observed by ALMA originates mostly from the chromosphere - a complex and dynamic region between the photosphere and corona, which plays a crucial role in the transport of energy and matter and, ultimately, the heating of the outer layers of the solar atmosphere. Based on first solar test observations, strategies for regular solar campaigns are currently being developed. State-of-the-art numerical simulations of the solar atmosphere and modeling of instrumental effects can help constrain and optimize future observing modes for ALMA. Here we present a short technical description of ALMA and an overview of past efforts and future possibilities for solar observations at submillimeter and millimeter wavelengths. In addition, selected numerical simulations and observations at other wavelengths demonstrate ALMA's scientific potential for studying the Sun for a large range of science cases.

The recent discovery of impulsive solar burst emission in the 30 THz band is raising new interpretation challenges. One event associated with a GOES M2 class flare has been observed simultaneously in microwaves, H-alpha, EUV, and soft X-ray bands. Although these new observations confirm some features found in the two prior known events, they exhibit time profile structure discrepancies between 30 THz, microwaves, and hard X-rays (as inferred from the Neupert effect). These results suggest a more complex relationship between 30 THz emission and radiation produced at other wavelength ranges. The multiple frequency emissions in the impulsive phase are likely to be produced at a common flaring site lower in the chromosphere. The 30 THz burst emission may be either part of a nonthermal radiation mechanism or due to the rapid thermal response to a beam of high-energy particles bombarding the dense solar atmosphere.

We use weak gravitational lensing to measure mean mass profiles around Locally Brightest Galaxies (LBGs). These are selected from the SDSS/DR7 spectroscopic and photometric catalogues to be brighter than any neighbour projected within 1.0 Mpc and differing in redshift by $<1000$ km/s. Most ($> 83\%$) are expected to be the central galaxies of their dark matter halos. Previous stacking analyses have used this LBG sample to measure mean Sunyaev-Zeldovich flux and mean X-ray luminosity as a function of LBG stellar mass. In both cases, a simulation of the formation of the galaxy population was used to estimate effective halo mass for LBGs of given stellar mass, allowing the derivation of scaling relations between the gas properties of halos and their mass. By comparing results from a variety of simulations to our lensing data, we show that this procedure has significant model dependence reflecting: (i) the failure of any given simulation to reproduce observed galaxy abundances exactly; (ii) a dependence on the cosmology underlying the simulation; and (iii) a dependence on the details of how galaxies populate halos. We use our lensing results to recalibrate the scaling relations, eliminating most of this model dependence and explicitly accounting both for residual modelling uncertainties and for observational uncertainties in the lensing results. The resulting scaling relations link the mean gas properties of dark halos to their mass over an unprecedentedly wide range, $10^{12.5}<M_{500}/ \mathrm{M_\odot}<10^{14.5}$, and should fairly and robustly represent the full halo population.

We have updated the Munich galaxy formation model to the Planck first-year cosmology, while modifying the treatment of baryonic processes to reproduce recent data on the abundance and passive fractions of galaxies from z= 3 down to z=0. Matching these more extensive and more precise observational results requires us to delay the reincorporation of wind ejecta, to lower the surface density threshold for turning cold gas into stars, to eliminate ram-pressure stripping in haloes less massive than ~10^14 Msun, and to modify our model for radio mode feedback. These changes cure the most obvious failings of our previous models, namely the overly early formation of low-mass galaxies and the overly large fraction of them that are passive at late times. The new model is calibrated to reproduce the observed evolution both of the stellar mass function and of the distribution of star formation rate at each stellar mass. Massive galaxies (M>10^11 [Msun]) assemble most of their mass before z=1 and are predominantly old and passive at z=0, while lower mass galaxies assemble later and, for M<10^9.5 (Msun), are still predominantly blue and star forming at z=0. This phenomenological but physically based model allows the observations to be interpreted in terms of the efficiency of the various processes that control the formation and evolution of galaxies as a function of their stellar mass, gas content, environment and time.

The Local Group of galaxies offer some of the most discriminating tests of models of cosmic structure formation. For example, observations of the Milky Way (MW) and Andromeda satellite populations appear to be in disagreement with N-body simulations of the "Lambda Cold Dark Matter" ({\Lambda}CDM) model: there are far fewer satellite galaxies than substructures in cold dark matter halos (the "missing satellites" problem); dwarf galaxies seem to avoid the most massive substructures (the "too-big-to-fail" problem); and the brightest satellites appear to orbit their host galaxies on a thin plane (the "planes of satellites" problem). Here we present results from APOSTLE (A Project Of Simulating The Local Environment), a suite of cosmological hydrodynamic simulations of twelve volumes selected to match the kinematics of the Local Group (LG) members. Applying the Eagle code to the LG environment, we find that our simulations match the observed abundance of LG galaxies, including the satellite galaxies of the MW and Andromeda. Due to changes to the structure of halos and the evolution in the LG environment, the simulations reproduce the observed relation between stellar mass and velocity dispersion of individual dwarf spheroidal galaxies without necessitating the formation of cores in their dark matter profiles. Satellite systems form with a range of spatial anisotropies, including one similar to that of the MW, confirming that such a configuration is not unexpected in {\Lambda}CDM. Finally, based on the observed velocity dispersion, size, and stellar mass, we provide new estimates of the maximum circular velocity for the halos of nine MW dwarf spheroidals.

ALMA will open a new chapter in the study of the Sun by providing a leap in spatial resolution and sensitivity compared to currently available mm wave- length observations. In preparation of ALMA, we have carried out a large number of observational tests and state-of-the-art radiation MHD simulations. Here we review the best available observations of the Sun at millimeter wavelengths.Using state of the art radiation MHD simulations of the solar atmosphere we demonstrate the huge potential of ALMA observations for uncovering the nature of the solar chromosphere. We show that ALMA will not only provide a reliable probe of the thermal structure and dynamics of the chromosphere, it will also open up a powerful new diagnostic of magnetic field at chromospheric heights, a fundamentally important, but so far poorly known parameter.

We present meter wave solar radio spectra of the highest spectrotemporal resolution achieved to date. The observations, obtained with the first station of the Long Wavelength Array (LWA1), show unprecedented detail of solar emissions across a wide bandwidth during a Type III/IIIb storm. Our flux calibration demonstrates that the LWA1 can detect Type III bursts much weaker than 1 SFU, much lower than previous observations, and that the distribution of fluxes in these bursts varies with frequency. The high sensitivity and low noise in the data provide strong constraints to models of this type of plasma emission. The continuous generation of electron beams in the corona revealed by the high density Type III storm is evidence for ubiquitous magnetic reconnection in the lower corona. Such an abundance of reconnection events not only contributes to the total coronal energy budget, but also provides an engine by which to form the populations of seed particles responsible for proton-rich solar energetic particle events. An active region with such levels of reconnection and the accompanying type III/IIIb storms is here proposed to be associated with an increase of SEP production if a CME erupts. The data's constraints on existing theories of type IIIIb production are used to make an association of the observed type IIIb storm to specific electron beam paths with increased inhomogeneities in density, temperature, and or turbulence. This scenario ties in the observed timing of III and IIIb storms, constrained theories of type III and IIIb emission, and the ability of the emitting AR to produce a strong SEP event. The result requires but a single observable to cement these ideas, the statistical correlation of type III/IIIb activity with SEP-productive AR.

The Planck design and scanning strategy provide many levels of redundancy that can be exploited to provide tests of internal consistency. One of the most important is the comparison of the 70GHz and 100GHz channels. Based on different instrument technologies, with feeds located differently in the focal plane, analysed independently by different teams using different software, and near the minimum of diffuse foreground emission, these channels are in effect two different experiments. The 143GHz channel has the lowest noise level on Planck, and is near the minimum of unresolved foreground emission. In this paper, we analyse the level of consistency achieved in the 2013 Planck data. We concentrate on comparisons between the 70/100/143GHz channel maps and power spectra, particularly over the angular scales of the first and second acoustic peaks, on maps masked for diffuse Galactic emission and for strong unresolved sources. Difference maps covering angular scales from 8deg-15arcmin are consistent with noise, and show no evidence of cosmic microwave background structure. Including small but important corrections for unresolved-source residuals, we demonstrate agreement between 70 and 100GHz power spectra averaged over 70<l<390 at the 0.8% level, and agreement between 143 and 100GHz power spectra of 0.4% over the same l range. These values are within and consistent with the overall uncertainties in calibration given in the Planck 2013 results. We also present results based on the 2013 likelihood analysis showing consistency at the 0.35% between the 100/143/217GHz power spectra. We analyse calibration procedures and beams to determine what fraction of these differences can be accounted for by known approximations or systematic errors that could be controlled even better in the future, reducing uncertainties still further. Several possible small improvements are described...(abridged)

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity.

We present results from the first solar full-disk F$_{10.7}$ (the radio flux at $10.7$ cm, $2.8$ GHz) image taken with the S-band receivers on the recently upgraded Karl G. Jansky Very Large Array (VLA) in order to assess the relationship between the F$_{10.7}$ index and solar extreme ultra-violet (EUV) emission. To identify the sources of the observed $2.8$ GHz emission, we calculate differential emission measures (DEMs) from EUV images collected by the Atmospheric Imaging Assembly (AIA) and use them to predict the bremsstrahlung component of the radio emission. By comparing the bremsstrahlung prediction and radio observation we find that $8.1\pm 0.5\%$ of the variable component of the F$_{10.7}$ flux is associated with the gyroresonance emission mechanism. Additionally, we identify optical depth effects on the radio limb which may complicate the use of F$_{10.7}$ time series as an EUV proxy. Our analysis is consistent with a coronal iron abundance that is $4$ times the photospheric level.

We present a sub-100 pc-scale analysis of the CO molecular gas emission and kinematics of the gravitational lens system SDP.81 at redshift 3.042 using Atacama Large Millimetre/submillimetre Array (ALMA) science verification data and a visibility-plane lens reconstruction technique. We find clear evidence for an excitation dependent structure in the unlensed molecular gas distribution, with emission in CO (5-4) being significantly more diffuse and structured than in CO (8-7). The intrinsic line luminosity ratio is r_8-7/5-4 = 0.30 +/- 0.04, which is consistent with other low-excitation starbursts at z ~ 3. An analysis of the velocity fields shows evidence for a star-forming disk with multiple velocity components that is consistent with a merger/post-coalescence merger scenario, and a dynamical mass of M(< 1.56 kpc) = 1.6 +/- 0.6 x 10^10 M_sol . Source reconstructions from ALMA and the Hubble Space Telescope show that the stellar component is offset from the molecular gas and dust components. Together with Karl G. Jansky Very Large Array CO (1-0) data, they provide corroborative evidence for a complex ~2 kpc-scale starburst that is embedded within a larger ~15 kpc structure.

We examine the circular velocity profiles of galaxies in {\Lambda}CDM cosmological hydrodynamical simulations from the EAGLE and LOCAL GROUPS projects and compare them with a compilation of observed rotation curves of galaxies spanning a wide range in mass. The shape of the circular velocity profiles of simulated galaxies varies systematically as a function of galaxy mass, but shows remarkably little variation at fixed maximum circular velocity. This is especially true for low-mass dark matter-dominated systems, reflecting the expected similarity of the underlying cold dark matter haloes. This is at odds with observed dwarf galaxies, which show a large diversity of rotation curve shapes, even at fixed maximum rotation speed. Some dwarfs have rotation curves that agree well with simulations, others do not. The latter are systems where the inferred mass enclosed in the inner regions is much lower than expected for cold dark matter haloes and include many galaxies where previous work claims the presence of a constant density "core". The "cusp vs core" issue is thus better characterized as an "inner mass deficit" problem than as a density slope mismatch. For several galaxies the magnitude of this inner mass deficit is well in excess of that reported in recent simulations where cores result from baryon-induced fluctuations in the gravitational potential. We conclude that one or more of the following statements must be true: (i) the dark matter is more complex than envisaged by any current model; (ii) current simulations fail to reproduce the effects of baryons on the inner regions of dwarf galaxies; and/or (iii) the mass profiles of "inner mass deficit" galaxies inferred from kinematic data are incorrect.

We simulate the phase-space distribution of stellar mass in 9 massive Lambda-CDM galaxy clusters by applying the semi-analytic particle tagging method of Cooper et al. to the Phoenix suite of high-resolution N-body simulations (M200 = 7.5 to 33 x 10^14 Msol). The resulting surface brightness (SB) profiles of brightest cluster galaxies (BCGs) match well to observations. On average, stars formed in galaxies accreted by the BCG account for ~90 per cent of its total mass (the remainder is formed in situ). In circular BCG-centred apertures, the superposition of multiple debris clouds (each ~10 per cent of the total BCG mass) from different progenitors can result in an extensive outer diffuse component, qualitatively similar to a 'cD envelope'. These clouds typically originate from tidal stripping at z < 1 and comprise both streams and the extended envelopes of other massive galaxies in the cluster. Stars at very low SB contribute a significant fraction of the total cluster stellar mass budget: in the central 1 Mpc^2 of a z ~ 0.15 cluster imaged at SDSS-like resolution, our fiducial model predicts 80-95 per cent of stellar mass below a SB of mu_V = 26.5 mag arcsec^2 is associated with accreted stars in the envelope of the BCG. The ratio of BCG stellar mass (including this diffuse component) to total cluster stellar mass is ~30 per cent.

Impulsive 30 THz continuum bursts have been recently observed in solar flares, utilizing small telescopes with a unique and relatively simple optical setup concept. The most intense burst was observed together with a GOES X2 class event on October 27, 2014, also detected at two sub-THz frequencies, RHESSI X-rays and SDO/HMI and EUV. It exhibits strikingly good correlation in time and in space with white light flare emission. It is likely that this association may prove to be very common. All three 30 THz events recently observed exhibited intense fluxes in the range of 104 solar flux units, considerably larger than those measured for the same events at microwave and sub-mm wavelengths. The 30 THz burst emission might be part of the same spectral burst component found at sub-THz frequencies. The 30 THz solar bursts open a promising new window for the study of flares at their origin

We present a sub-50 pc-scale analysis of the gravitational lens system SDP.81 at redshift 3.042 using Atacama Large Millimetre/submillimetre Array (ALMA) science verification data. We model both the mass distribution of the gravitational lensing galaxy and the pixelated surface brightness distribution of the background source using a novel Bayesian technique that fits the data directly in visibility space. We find the 1 and 1.3 mm dust emission to be magnified by a factor of u_tot = 17.6+/-0.4, giving an intrinsic total star-formation rate of 315+/-60 M_sol/yr and a dust mass of 6.4+/-1.5*10^8 M_sol. The reconstructed dust emission is found to be non-uniform, but composed of multiple regions that are heated by both diffuse and strongly clumped star-formation. The highest surface brightness region is a ~1.9*0.7 kpc disk-like structure, whose small extent is consistent with a potential size-bias in gravitationally lensed starbursts. Although surrounded by extended star formation, with a density of 20-30+/-10 M_sol/yr/kpc^2, the disk contains three compact regions with densities that peak between 120-190+/-20 M_sol/yr/kpc^2. Such star-formation rate densities are below what is expected for Eddington-limited star-formation by a radiation pressure supported starburst. There is also a tentative variation in the spectral slope of the different star-forming regions, which is likely due to a change in the dust temperature and/or opacity across the source.

We report the results of a joint analysis of data from BICEP2/Keck Array and Planck. BICEP2 and Keck Array have observed the same approximately 400 deg$^2$ patch of sky centered on RA 0h, Dec. $-57.5\deg$. The combined maps reach a depth of 57 nK deg in Stokes $Q$ and $U$ in a band centered at 150 GHz. Planck has observed the full sky in polarization at seven frequencies from 30 to 353 GHz, but much less deeply in any given region (1.2 $\mu$K deg in $Q$ and $U$ at 143 GHz). We detect 150$\times$353 cross-correlation in $B$-modes at high significance. We fit the single- and cross-frequency power spectra at frequencies $\geq 150$ GHz to a lensed-$\Lambda$CDM model that includes dust and a possible contribution from inflationary gravitational waves (as parameterized by the tensor-to-scalar ratio $r$), using a prior on the frequency spectral behavior of polarized dust emission from previous \planck\ analysis of other regions of the sky. We find strong evidence for dust and no statistically significant evidence for tensor modes. We probe various model variations and extensions, including adding a synchrotron component in combination with lower frequency data, and find that these make little difference to the $r$ constraint. Finally we present an alternative analysis which is similar to a map-based cleaning of the dust contribution, and show that this gives similar constraints. The final result is expressed as a likelihood curve for $r$, and yields an upper limit $r_{0.05}<0.12$ at 95% confidence. Marginalizing over dust and $r$, lensing $B$-modes are detected at $7.0\,\sigma$ significance.

We present results from thirteen cosmological simulations that explore the parameter space of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) simulation project. Four of the simulations follow the evolution of a periodic cube L = 50 cMpc on a side, and each employs a different subgrid model of the energetic feedback associated with star formation. The relevant parameters were adjusted so that the simulations each reproduce the observed galaxy stellar mass function at z = 0.1. Three of the simulations fail to form disc galaxies as extended as observed, and we show analytically that this is a consequence of numerical radiative losses that reduce the efficiency of stellar feedback in high-density gas. Such losses are greatly reduced in the fourth simulation - the EAGLE reference model - by injecting more energy in higher density gas. This model produces galaxies with the observed size distribution, and also reproduces many galaxy scaling relations. In the remaining nine simulations, a single parameter or process of the reference model was varied at a time. We find that the properties of galaxies with stellar mass <~ M* (the "knee" of the galaxy stellar mass function) are largely governed by feedback associated with star formation, while those of more massive galaxies are also controlled by feedback from accretion onto their central black holes. Both processes must be efficient in order to reproduce the observed galaxy population. In general, simulations that have been calibrated to reproduce the low-redshift galaxy stellar mass function will still not form realistic galaxies, but the additional requirement that galaxy sizes be acceptable leads to agreement with a large range of observables.

We have simulated the Expanded Owens Valley Solar Array (EOVSA) radio images generated at multiple frequencies from a model solar active region, embedded in a realistic solar disk model, and explored the resulting datacube for different spectral analysis schemes to evaluate the potential for realizing one of EOVSA's most important scientific goals--coronal magnetography. In this paper, we focus on modeling the gyroresonance and free-free emission from an on-disk solar active region model with realistic complexities in electron density, temperature and magnetic field distribution. We compare the magnetic field parameters extrapolated from the image datacube along each line of sight after folding through the EOVSA instrumental profile with the original (unfolded) parameters used in the model. We find that even the most easily automated, image-based analysis approach (Level 0) provides reasonable quantitative results, although they are affected by systematic effects due to finite sampling in the Fourier (uv) plane. Finally, we note the potential for errors due to misidentified harmonics of the gyrofrequency, and discuss the prospects for applying a more sophisticated spectrally-based analysis scheme (Level 1) to resolve the issue in cases where improved uv coverage and spatial resolution are available.

We examine a sample of $\sim 250 000$ "locally brightest galaxies" selected from the Sloan Digital Sky Survey to be central galaxies within their dark matter halos. These were previously stacked by the Planck Collaboration to measure the Sunyaev-Zel'dovich signal as a function of central galaxy stellar mass. Here, we stack the X-ray emission from these halos using data from the ROSAT All-Sky Survey. We detect emission across almost our entire sample, including emission which we attribute to hot gas around galaxies spanning a range of 1.2 dex in stellar mass (corresponding to nearly two orders of magnitude in halo mass) down to $M* = 10^{10.8} M_{\odot}$ ($M_{500} \approx 10^{12.6} M_{\odot}$). Over this range, the X-ray luminosity can be fit by a power-law, either of stellar mass or of halo mass. A single unified scaling relation between mass and $L_X$ applies for galaxies, groups, and clusters. This relation has a steeper slope than expected for self-similarity, in contrast to the $Y_{SZ}$-$M_{500}$ relation, showing the importance of non-gravitational heating. If this heating is predominantly due to AGN feedback, the lack of a break in our relation suggests that AGN feedback is tightly self-regulated and fairly gentle, in agreement with recent high-resolution simulations. Our results are consistent with earlier measurements of the $L_X$-$L_K$ relation for elliptical galaxies and of the $L_X$-$M_{500}$ relation for optically-selected galaxy clusters. However, our $L_X$-$M_{500}$ relation lies more than a factor of two below most previous relations based on X-ray-selected cluster samples. We argue that optical selection offers a less biased view of the $L_X$-$M_{500}$ relation.

We present the all-sky Planck catalogue of Sunyaev-Zeldovich (SZ) sources detected from the 29 month full-mission data. The catalogue (PSZ2) is the largest SZ-selected sample of galaxy clusters yet produced and the deepest all-sky catalogue of galaxy clusters. It contains 1653 detections, of which 1203 are confirmed clusters with identified counterparts in external data-sets, and is the first SZ-selected cluster survey containing > $10^3$ confirmed clusters. We present a detailed analysis of the survey selection function in terms of its completeness and statistical reliability, placing a lower limit of 83% on the purity. Using simulations, we find that the Y5R500 estimates are robust to pressure-profile variation and beam systematics, but accurate conversion to Y500 requires. the use of prior information on the cluster extent. We describe the multi-wavelength search for counterparts in ancillary data, which makes use of radio, microwave, infra-red, optical and X-ray data-sets, and which places emphasis on the robustness of the counterpart match. We discuss the physical properties of the new sample and identify a population of low-redshift X-ray under- luminous clusters revealed by SZ selection. These objects appear in optical and SZ surveys with consistent properties for their mass, but are almost absent from ROSAT X-ray selected samples.

We present a 67--93.6 GHz spectral line survey of Orion-KL with the new 4 mm Receiver on the Green Bank Telescope (GBT). The survey reaches unprecedented depths and covers the low-frequency end of the 3 mm atmospheric window which has been relatively unexplored previously. The entire spectral-line survey is published electronically for general use by the astronomical community. The calibration and performance of 4 mm Receiver on the GBT is also summarized.

We update the all-sky Planck catalogue of 1227 clusters and cluster candidates (PSZ1) published in March 2013, derived from Sunyaev-Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. Addendum. We deliver an updated version of the PSZ1 catalogue, reporting the further confirmation of 86 Planck-discovered clusters. In total, the PSZ1 now contains 947 confirmed clusters, of which 214 were confirmed as newly discovered clusters through follow-up observations undertaken by the Planck Collaboration. The updated PSZ1 contains redshifts for 913 systems, of which 736 (~80.6%) are spectroscopic, and associated mass estimates derived from the Y_z mass proxy. We also provide a new SZ quality flag, derived from a novel artificial neural network classification of the SZ signal, for the remaining 280 candidates. Based on this assessment, the purity of the updated PSZ1 catalogue is estimated to be 94%. In this release, we provide the full updated catalogue and an additional readme file with further information on the Planck SZ detections.

Understanding the interplay between black-hole accretion and star formation, and how to disentangle the two, is crucial to our understanding of galaxy formation and evolution. To investigate, we use a combination of optical and near-infrared photometry to select a sample of 74 quasars from the VISTA Deep Extragalactic Observations (VIDEO) Survey, over 1 deg$^2$. The depth of VIDEO allows us to study very low accretion rates and/or lower-mass black holes, and 26 per cent of the candidate quasar sample has been spectroscopically confirmed. We use a radio-stacking technique to sample below the nominal flux-density threshold using data from the Very Large Array at 1.4 GHz and find, in agreement with other work, that a power-law fit to the quasar-related radio source counts is inadequate at low flux density. By comparing with a control sample of galaxies (where we match in terms of stellar mass), and by estimating the star formation rate, we suggest that this radio emission is predominantly caused by accretion activity rather than star-formation activity.

Aims. We use advanced 3D NLTE radiative magnetohydrodynamic simulations of the solar atmosphere to carry out detailed tests of chromospheric diagnostics at millimeter and submillimeter wavelengths. Methods. We focused on the diagnostics of the thermal structure of the chromosphere in the wavelength bands from 0.4 mm up to 9.6 mm that can be accessed with the Atacama Large Millimeter/Submillimeter Array (ALMA) and investigated how these diagnostics are affected by the instrumental resolution. Results. We find that the formation height range of the millimeter radiation depends on the location in the simulation domain and is related to the underlying magnetic structure. Nonetheless, the brightness temperature is a reasonable measure of the gas temperature at the effective formation height at a given location on the solar surface. There is considerable scatter in this relationship, but this is significantly reduced when very weak magnetic fields are avoided. Our results indicate that although instrumental smearing reduces the correlation between brightness and temperature, millimeter brightness can still be used to reliably diagnose electron temperature up to a resolution of 1". If the resolution is more degraded, then the value of the diagnostic diminishes rapidly. Conclusions. We conclude that millimeter brightness can image the chromospheric thermal structure at the height at which the radiation is formed. Thus multiwavelength observations with ALMA with a narrow step in wavelength should provide sufficient information for a tomographic imaging of the chromosphere.

Radio wavelengths offer the unique possibility of tracing the total star-formation rate in galaxies, both obscured and unobscured. As such, they may provide the most robust measurement of the star-formation history of the Universe. In this chapter we highlight the constraints that the SKA can place on the evolution of the star-formation history of the Universe, the survey area required to overcome sample variance, the spatial resolution requirements, along with the multi-wavelength ancillary data that will play a major role in maximising the scientific promise of the SKA. The required combination of depth and resolution means that a survey to trace the star formation in the Universe should be carried out with a facility that has a resolution of at least ~0.5arcsec, with high sensitivity at < 1 GHz. We also suggest a strategy that will enable new parameter space to be explored as the SKA expands over the coming decade.

It has become apparent that active galactic nuclei (AGN) may have a significant impact on the growth and evolution of their host galaxies and vice versa but a detailed understanding of the interplay between these processes remains elusive. Deep radio surveys provide a powerful, obscuration-independent tool for measuring both star formation and AGN activity in high-redshift galaxies. Multiwavelength studies of deep radio fields show a composite population of star-forming galaxies and AGN, with the former dominating at the lowest flux densities (S$_{1.4\mathrm{GHz}}<$100~$\mu$Jy). The sensitivity and resolution of the SKA will allow us to identify, and separately trace, the total star formation in the bulges of individual high-redshift galaxies, the related nuclear activity and any star formation occurring on larger scales within a disc. We will therefore gain a detailed picture of the apparently simultaneous development of stellar populations and black holes in the redshift range where both star-formation and AGN activity peak (1$\leq$z$\leq$4). In this chapter we discuss the role of the SKA in studying the connection between AGN activity and galaxy evolution, and the most critical technical requirements for such of studies

The "Lambda Cold Dark Matter" (LCDM) model of cosmic structure formation is eminently falsifiable: once its parameters are fixed on large scales, it becomes testable in the nearby Universe. Observations within our Local Group of galaxies, including the satellite populations of the Milky Way and Andromeda, appear to contradict LCDM predictions: there are far fewer satellite galaxies than dark matter halos (the "missing satellites" problem), galaxies seem to avoid the largest substructures (the "too big to fail" problem), and the brightest satellites appear to orbit their host galaxies on a thin plane (the "planes of satellites" problem). We present results from the first hydrodynamic simulations of the Local Group that match the observed abundance of galaxies. We find that when baryonic and dark matter are followed simultaneously in the context of a realistic galaxy formation model, all three "problems" are resolved within the LCDM paradigm.

The polarized thermal emission from Galactic dust is the main foreground present in measurements of the polarization of the cosmic microwave background (CMB) at frequencies above 100GHz. We exploit the Planck HFI polarization data from 100 to 353GHz to measure the dust angular power spectra $C_\ell^{EE,BB}$ over the range $40<\ell<600$ well away from the Galactic plane. These will bring new insights into interstellar dust physics and a precise determination of the level of contamination for CMB polarization experiments. We show that statistical properties of the emission can be characterized over large fractions of the sky using $C_\ell$. For the dust, they are well described by power laws in $\ell$ with exponents $\alpha^{EE,BB}=-2.42\pm0.02$. The amplitudes of the polarization $C_\ell$ vary with the average brightness in a way similar to the intensity ones. The dust polarization frequency dependence is consistent with modified blackbody emission with $\beta_d=1.59$ and $T_d=19.6$K. We find a systematic ratio between the amplitudes of the Galactic $B$- and $E$-modes of 0.5. We show that even in the faintest dust-emitting regions there are no "clean" windows where primordial CMB $B$-mode polarization could be measured without subtraction of dust emission. Finally, we investigate the level of dust polarization in the BICEP2 experiment field. Extrapolation of the Planck 353GHz data to 150GHz gives a dust power $\ell(\ell+1)C_\ell^{BB}/(2\pi)$ of $1.32\times10^{-2}\mu$K$_{CMB}^2$ over the $40<\ell<120$ range; the statistical uncertainty is $\pm0.29$ and there is an additional uncertainty (+0.28,-0.24) from the extrapolation, both in the same units. This is the same magnitude as reported by BICEP2 over this $\ell$ range, which highlights the need for assessment of the polarized dust signal even in the cleanest windows of the sky.

We have constructed the first all-sky map of the thermal Sunyaev-Zeldovich (tSZ) effect by applying specifically tailored component separation algorithms to the 100 to 857 GHz frequency channel maps from the Planck survey. These maps show an obvious galaxy cluster tSZ signal that is well matched with blindly detected clusters in the Planck SZ catalogue. To characterize the signal in the tSZ map we have computed its angular power spectrum. At large angular scales ($\ell < 60$), the major foreground contaminant is the diffuse thermal dust emission. At small angular scales ($\ell > 500$) the clustered Cosmic Infrared Background (CIB) and residual point sources are the major contaminants. These foregrounds are carefully modelled and subtracted. We measure the tSZ power spectrum in angular scales, $0.17^{\circ} \lesssim \theta \lesssim 3.0^{\circ}$, that were previously unexplored. The measured tSZ power spectrum is consistent with that expected from the Planck catalogue of SZ sources, with additional clear evidence of signal from unresolved clusters and, potentially, diffuse warm baryons. We use the tSZ power spectrum to obtain the following cosmological constraints: $\sigma_8(\Omega_{\mathrm{m}}/0.28)^{3.2/8.1}=0.784 \pm 0.016 (68% C.L.). Marginalized band-powers of the Planck tSZ power spectrum and the best-fit model are given. The non-Gaussianity of the Compton parameter map is further characterized by computing its 1D probability distribution function and its bispectrum. These are used to place additional independent constraints on $\sigma_{8}$.

Aims. An efficient means of locating calibrator sources for International LOFAR is developed and used to determine the average density of usable calibrator sources on the sky for subarcsecond observations at 140 MHz. Methods. We used the multi-beaming capability of LOFAR to conduct a fast and computationally inexpensive survey with the full International LOFAR array. Sources were pre-selected on the basis of 325 MHz arcminute-scale flux density using existing catalogues. By observing 30 different sources in each of the 12 sets of pointings per hour, we were able to inspect 630 sources in two hours to determine if they possess a sufficiently bright compact component to be usable as LOFAR delay calibrators. Results. Over 40% of the observed sources are detected on multiple baselines between international stations and 86 are classified as satisfactory calibrators. We show that a flat low-frequency spectrum (from 74 to 325 MHz) is the best predictor of compactness at 140 MHz. We extrapolate from our sample to show that the density of calibrators on the sky that are sufficiently bright to calibrate dispersive and non-dispersive delays for the International LOFAR using existing methods is 1.0 per square degree. Conclusions. The observed density of satisfactory delay calibrator sources means that observations with International LOFAR should be possible at virtually any point in the sky, provided that a fast and efficient search using the methodology described here is conducted prior to the observation to identify the best calibrator.

We develop and test a new statistical method to measure the kinematic Sunyaev-Zel'dovich (kSZ) effect. A sample of independently detected clusters is combined with the cosmic flow field predicted from a galaxy redshift survey in order to derive a matched filter that optimally weights the kSZ signal for the sample as a whole given the noise involved in the problem. We apply this formalism to realistic mock microwave skies based on cosmological $N$-body simulations, and demonstrate its robustness and performance. In particular, we carefully assess the various sources of uncertainty, cosmic microwave background primary fluctuations, instrumental noise, uncertainties in the determination of the velocity field, and effects introduced by miscentring of clusters and by uncertainties of the mass-observable relation (normalization and scatter). We show that available data (\plk\ maps and the MaxBCG catalogue) should deliver a $7.7\sigma$ detection of the kSZ. A similar cluster catalogue with broader sky coverage should increase the detection significance to $\sim 13\sigma$. We point out that such measurements could be binned in order to study the properties of the cosmic gas and velocity fields, or combined into a single measurement to constrain cosmological parameters or deviations of the law of gravity from General Relativity.

We introduce the Virgo Consortium's EAGLE project, a suite of hydrodynamical simulations that follow the formation of galaxies and black holes in representative volumes. We discuss the limitations of such simulations in light of their finite resolution and poorly constrained subgrid physics, and how these affect their predictive power. One major improvement is our treatment of feedback from massive stars and AGN in which thermal energy is injected into the gas without the need to turn off cooling or hydrodynamical forces, allowing winds to develop without predetermined speed or mass loading factors. Because the feedback efficiencies cannot be predicted from first principles, we calibrate them to the z~0 galaxy stellar mass function and the amplitude of the galaxy-central black hole mass relation, also taking galaxy sizes into account. The observed galaxy mass function is reproduced to $\lesssim 0.2$ dex over the full mass range, $10^8 < M_*/M_\odot \lesssim 10^{11}$, a level of agreement close to that attained by semi-analytic models, and unprecedented for hydrodynamical simulations. We compare our results to a representative set of low-redshift observables not considered in the calibration, and find good agreement with the observed galaxy specific star formation rates, passive fractions, Tully-Fisher relation, total stellar luminosities of galaxy clusters, and column density distributions of intergalactic CIV and OVI. While the mass-metallicity relations for gas and stars are consistent with observations for $M_* \gtrsim 10^9 M_\odot$, they are insufficiently steep at lower masses. The gas fractions and temperatures are too high for clusters of galaxies, but for groups these discrepancies can be resolved by adopting a higher heating temperature in the subgrid prescription for AGN feedback. EAGLE constitutes a valuable new resource for studies of galaxy formation.

We develop a new method to account for the finite lifetimes of stars and trace individual abundances within a semi-analytic model of galaxy formation. At variance with previous methods, based on the storage of the (binned) past star formation history of model galaxies, our method projects the information about the metals produced by each simple stellar population (SSP) in the future. Using this approach, an accurate accounting of the timings and properties of the individual SSPs composing model galaxies is possible. We analyse the dependence of our chemical model on various ingredients, and apply it to six simulated haloes of roughly Milky Way mass and with no massive close neighbour at z=0. For all models considered, the [Fe/H] distributions of the stars in the disc component are in good agreement with Milky Way data, while for the spheroid component (whose formation we model only through mergers) these are offset low with respect to observational measurements for the Milky Way bulge. This is a consequence of narrow star formation histories, with relatively low rates of star formation. The slow recycling of gas and energy from supernovae in our chemical model has important consequences on the predicted star formation rates, which are systematically lower than the corresponding rates in the same physical model but with an instantaneous recycling approximation. The halo that resembles most our Galaxy in terms of its global properties also reproduces the observed relation between the average metallicity and luminosity of the Milky Way satellites, albeit with a slightly steeper slope.

We present cosmological N-body resimulations of the assembly of the Brightest Cluster Galaxies (BCGs) in rich clusters. At $z=2$ we populate dark matter subhalos with self-gravitating stellar systems whose abundance and structure match observed high-redshift galaxies. By $z=0$, mergers have built much larger galaxies at cluster centre. Their dark matter density profiles are shallower than in corresponding dark-matter-only simulations, but their total mass density profiles (stars + dark matter) are quite similar. Differences are found only at radii where the effects of central black holes may be significant. Dark matter density slopes shallower than $\gamma=1.0$ occur for $r/r_{200} < 0.015$, close to the half-light radii of the BCGs. Our experiments support earlier suggestions that NFW-like profiles are an attractor for the hierarchical growth of structure in collisionless systems -- total mass density profiles asymptote to the solution found in dark-matter-only simulations over the radial range where mergers produce significant mixing between stars and dark matter. Simulated dark matter fractions are substantially higher in BCGs than in field ellipticals, reaching 80\% within the half-light radius. We also estimate that supermassive black hole mergers should create BCG cores as large as $r_{c}\sim 3 \,\mathrm{kpc}$. The good agreement of all these properties with recent observational studies of BCG structure suggests that dissipational processes have not played a dominant role in the assembly of the observed systems.

Measurements of the total logarithmic central slope of the mass profile in galaxy clusters constrain their evolution and assembly history and that of their brightest cluster galaxies. We report the first full surface brightness distribution modelling of the inner region of the galaxy cluster MACS J1149.5+2223. We compare these results with a position-based modelling approach for which we employ more than twice the previously known positional constraints. This is the first time that the detailed lensed image configuration of two non-central cluster galaxies with Einstein rings has been mapped. Due to the extended radial coverage provided by the multiple images in this system, we are able to determine the slope $\partial \log{ \kappa }/\partial \log{R} = -0.33$ of the total projected mass distribution from $8$ to $80~\mathrm{kpc}$. This is within the cluster-to-cluster scatter estimates from previous cluster measurements. Our reconstruction of the image surface brightness distribution of the large central spiral galaxy has a root mean square residual for all image pixels of $1.14~\sigma$, where $\sigma$ is the observational background noise. This corresponds to a reconstruction of the positions of bright clumps in the central galaxy with an rms of $0.063~\mathrm{arcsec}$.

On the arcminute angular scales probed by Planck, the CMB anisotropies are gently perturbed by gravitational lensing. Here we present a detailed study of this effect, detecting lensing independently in the 100, 143, and 217GHz frequency bands with an overall significance of greater than 25sigma. We use the temperature-gradient correlations induced by lensing to reconstruct a (noisy) map of the CMB lensing potential, which provides an integrated measure of the mass distribution back to the CMB last-scattering surface. Our lensing potential map is significantly correlated with other tracers of mass, a fact which we demonstrate using several representative tracers of large-scale structure. We estimate the power spectrum of the lensing potential, finding generally good agreement with expectations from the best-fitting LCDM model for the Planck temperature power spectrum, showing that this measurement at z=1100 correctly predicts the properties of the lower-redshift, later-time structures which source the lensing potential. When combined with the temperature power spectrum, our measurement provides degeneracy-breaking power for parameter constraints; it improves CMB-alone constraints on curvature by a factor of two and also partly breaks the degeneracy between the amplitude of the primordial perturbation power spectrum and the optical depth to reionization, allowing a measurement of the optical depth to reionization which is independent of large-scale polarization data. Discarding scale information, our measurement corresponds to a 4% constraint on the amplitude of the lensing potential power spectrum, or a 2% constraint on the RMS amplitude of matter fluctuations at z~2.

We describe the processing of the 531 billion raw data samples from the High Frequency Instrument (hereafter HFI), which we performed to produce six temperature maps from the first 473 days of Planck-HFI survey data. These maps provide an accurate rendition of the sky emission at 100, 143, 217, 353, 545, and 857 GHz with an angular resolution ranging from 9.7 to 4.6 arcmin. The detector noise per (effective) beam solid angle is respectively, 10, 6, 12 and 39 microKelvin in HFI four lowest frequency channel (100--353 GHz) and 13 and 14 kJy/sr for the 545 and 857 GHz channels. Using the 143 GHz channel as a reference, these two high frequency channels are intercalibrated within 5% and the 353 GHz relative calibration is at the percent level. The 100 and 217 GHz channels, which together with the 143 GHz channel determine the high-multipole part of the CMB power spectrum (50 < l <2500), are intercalibrated at better than 0.2 %.

We study how feedback influences baryon infall onto galaxies using cosmological, zoom-in simulations of haloes with present mass $M_{vir}=6.9\times10^{11} M_{\odot}$ to $1.7\times10^{12} M_{\odot}$. Starting at z=4 from identical initial conditions, implementations of weak and strong stellar feedback produce bulge- and disc-dominated galaxies, respectively. Strong feedback favours disc formation: (1) because conversion of gas into stars is suppressed at early times, as required by abundance matching arguments, resulting in flat star formation histories and higher gas fractions; (2) because 50% of the stars form in situ from recycled disc gas with angular momentum only weakly related to that of the z=0 dark halo; (3) because late-time gas accretion is typically an order of magnitude stronger and has higher specific angular momentum, with recycled gas dominating over primordial infall; (4) because 25-30% of the total accreted gas is ejected entirely before z~1, removing primarily low angular momentum material which enriches the nearby inter-galactic medium. Most recycled gas roughly conserves its angular momentum, but material ejected for long times and to large radii can gain significant angular momentum before re-accretion. These processes lower galaxy formation efficiency in addition to promoting disc formation.

The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the SZ effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter LCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25 sigma. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations derived from CMB data and that derived from SZ data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak.

We use the SDSS/DR8 galaxy sample to study the radial distribution of satellite galaxies around isolated primaries, comparing to semi-analytic models of galaxy formation based on the Millennium and Millennium-II simulations. SDSS satellites behave differently around high- and low-mass primaries: those orbiting objects with $M_*>10^{11}M_\odot$ are mostly red and are less concentrated towards their host than the inferred dark matter halo, an effect that is very pronounced for the few blue satellites. On the other hand, less massive primaries have steeper satellite profiles that agree quite well with the expected dark matter distribution and are dominated by blue satellites, even in the inner regions where strong environmental effects are expected. In fact, such effects appear to be strong only for primaries with $M_* > 10^{11}M_\odot$. This behaviour is not reproduced by current semi-analytic simulations, where satellite profiles always parallel those of the dark matter and satellite populations are predominantly red for primaries of all masses. The disagreement with SDSS suggests that environmental effects are too efficient in the models. Modifying the treatment of environmental and star formation processes can substantially increase the fraction of blue satellites, but their radial distribution remains significantly shallower than observed. It seems that most satellites of low-mass primaries can continue to form stars even after orbiting within their joint halo for 5 Gyr or more.

We analyze the formation histories of 19 galaxies from cosmological smoothed particle hydrodynamics zoom-in resimulations. We construct mock three-colour images and show that the models reproduce observed trends in the evolution of galaxy colours and morphologies. However, only a small fraction of galaxies contains bars. Many galaxies go through phases of central mass growth by in-situ star formation driven by gas-rich mergers or misaligned gas infall. These events lead to accretion of low-angular momentum gas to the centres and leave imprints on the distributions of z=0 stellar circularities, radii and metallicities as functions of age. Observations of the evolution of structural properties of samples of disc galaxies at z=2.5-0.0 infer continuous mass assembly at all radii. Our simulations can only explain this if there is a significant contribution from mergers or misaligned infall, as expected in a LambdaCDM universe. Quiescent merger histories lead to high kinematic disc fractions and inside-out growth, but show little central growth after the last `destructive' merger at z>1.5. For sufficiently strong feedback, as assumed in our models, a moderate amount of merging does not seem to be a problem for the z=0 disc galaxy population, but may rather be a requirement. The average profiles of simulated disc galaxies agree with observations at z>=1.5. At z<=1, there is too much growth in size and too little growth in central mass, possibly due to the under-abundance of bars. The discrepancies may partly be caused by differences between the star formation histories of the simulations and those assumed for observations.

We describe the all-sky Planck catalogue of clusters and cluster candidates derived from Sunyaev--Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. The catalogue contains 1227 entries, making it over six times the size of the Planck Early SZ (ESZ) sample and the largest SZ-selected catalogue to date. It contains 861 confirmed clusters, of which 178 have been confirmed as clusters, mostly through follow-up observations, and a further 683 are previously-known clusters. The remaining 366 have the status of cluster candidates, and we divide them into three classes according to the quality of evidence that they are likely to be true clusters. The Planck SZ catalogue is the deepest all-sky cluster catalogue, with redshifts up to about one, and spans the broadest cluster mass range from (0.1 to 1.6) 10^{15}Msun. Confirmation of cluster candidates through comparison with existing surveys or cluster catalogues is extensively described, as is the statistical characterization of the catalogue in terms of completeness and statistical reliability. The outputs of the validation process are provided as additional information. This gives, in particular, an ensemble of 813 cluster redshifts, and for all these Planck clusters we also include a mass estimated from a newly-proposed SZ-mass proxy. A refined measure of the SZ Compton parameter for the clusters with X-ray counter-parts is provided, as is an X-ray flux for all the Planck clusters not previously detected in X-ray surveys.