A photometric redshift sample of Luminous Red Galaxies (hereafter LRGs) obtained from The DECam Legacy Survey (DECaLS) is analysed to probe cosmic distances by exploiting the wedge approach of the two-point correlation function. Although the cosmological information is highly contaminated by the uncertainties existing in the photometric redshifts from the galaxy map, an angular diameter distance can be probed at the perpendicular configuration in which the measured correlation function is minimally contaminated. An ensemble of wedged correlation functions selected up to a given threshold based on having the least contamination was studied in the previous work (arXiv:1903.09651v2 [astro-ph.CO]) using simulations, and the extracted cosmological information was unbiased within this threshold. We apply the same methodology for analysing the LRG sample from DECaLS which will provide the optical imaging for targeting two-thirds of the DESI footprint and measure the angular diameter distances at $z=0.69$ and $z=0.87$ to be $D_{A}(0.697)=(1499 \pm 77\,\mathrm{Mpc})(r_{d}/r_{d,fid})$ and $D_{A}(0.874)=(1680 \pm 109\,\mathrm{Mpc})(r_{d}/r_{d,fid})$ with a fractional error of 5.14% and 6.48% respectively. We obtain a value of $H_{0}=67.59\pm5.52$ km/s/Mpc which supports the $H_0$ measured by all other BAO results and is consistent with $\Lambda$CDM model.

Context. Solar radio bursts originate mainly from high energy electrons accelerated in solar eruptions like solar flares, jets, and coronal mass ejections. A subcategory of solar radio bursts with short time duration may be used as a proxy to understand the wave generation and propagation within the corona. Aims. Complete case studies of the source size, position and kinematics of short term bursts are very rare due to instrumental limitations. A comprehensive multi-frequency spectroscopic and imaging study was carried out of a clear example of a solar type IIIb-III pair. Methods. In this work, the source of the radio burst was imaged with the interferometric mode, using the remote baselines of the LOw Frequency ARray (LOFAR). A detailed analysis of the fine structures in the spectrum and of the radio source motion with imaging was conducted. Results. The study shows how the fundamental and harmonic components have a significantly different source motion. The apparent source of the fundamental emission at 26.56MHz displaces away from the solar disk center at about 4 times the speed of light, while the apparent source of the harmonic emission at the same frequency shows a speed of < 0.02c. The source size of the harmonic emission, observed in this case, is smaller than that in previous studies, indicating the importance of the use of the remote baselines.

We present LOFAR observations at 150 MHz of the borderline FRI/FRII giant radio galaxy NGC 6251. This paper presents the most sensitive and highest-resolution images of NGC 6251 at these frequencies to date, revealing for the first time a low-surface-brightness extension to the northern lobe, and a possible backflow associated with the southern lobe. The integrated spectra of components of NGC 6251 are consistent with previous measurements at higher frequencies, similar to results from other LOFAR studies of nearby radio galaxies. We find the outer structures of NGC 6251 to be either at equipartition or slightly electron dominated, similar to those of FRII sources rather than FRIs; but this conclusion remains tentative because of uncertainties associated with the geometry and the extrapolation of X-ray measurements to determine the external pressure distribution on the scale of the outer lobes. We place lower limits on the ages of the extension of the northern lobe and the backflow of the southern lobe of $t \gtrsim 250$ Myr and $t \gtrsim 210$ Myr respectively. We present the first detection of polarisation at 150 MHz in NGC 6251. Taking advantage of the high Faraday resolution of LOFAR, we place an upper limit on the magnetic field in the group of $B < 0.2 (\Lambda_B / 10 {\rm kpc})^{-0.5} \mu$G for a coherence scale of $\Lambda_B < 60 {\rm kpc}$ and $B < 13 \mu$G for $\Lambda_B = 240$ kpc.

Low-power radio sources dominate the radio sky. They tend to be small in size and core-dominated, but the origin of their properties and the evolution of their radio plasma are not well constrained. Interestingly, there is mounting evidence that low-power radio sources can significantly impact their surrounding gaseous medium and, therefore, may be more relevant for galaxy evolution than previously thought. In this paper, we present low radio frequency observations obtained with LOFAR at 147 MHz of the radio source hosted by NGC 3998. This is a rare example of a low-power source which is extremely core-dominated, but which has two large-scale lobes of low surface brightness. We combine the new 147 MHz image with available 1400 MHz data to derive the spectral index over the source. Despite the low surface brightness, reminiscent of remnant structures, the lobes show an optically thin synchrotron spectral index (~ 0.6). We interpret this as being due to rapid decollimation of the jets close to the core, to high turbulence of the plasma flow, and entrainment of thermal gas. This could be the result of intermittent activity of the central AGN, or, more likely, temporary disruption of the jet due to the interaction of the jet with the rich circumnuclear ISM. Both would result in sputtering energy injection from the core which would keep the lobes fed, albeit at a low rate. We discuss these results in connection with the properties of low-power radio sources in general. Our findings show that amorphous, low surface brightness lobes should not be interpreted, by default, as remnant structures. Large, deep surveys (in particular the LOFAR 150 MHz LoTSS and the recently started 1400 MHz Apertif survey) will identify a growing number of objects similar to NGC 3998 where these ideas can be further tested.

Spiral structure in disk galaxies could arise from transient modes that create conditions conducive for their regeneration; this is the proposal of Sellwood and Carlberg, based on simulations of stellar disks. The linear response of an axisymmetric stellar disk, to an adiabatic non-axisymmetric transient mode, gives a final distribution function (DF) that is equal to the initial DF everywhere in phase space, except at the Lindblad and corotation resonances where the final DF is singular. We use the nonlinear theory of adiabatic capture into resonance to resolve the singularities and calculate the finite changes in the DF. These take the form of axisymmetric "scars" concentrated around resonances, whose DFs have simple general forms. Global changes in the physical properties are explored for a cool Mestel disk: we calculate the DFs of scars, and estimate the changes in the disk angular momentum, surface density and orbital frequencies leading to shifts of resonances. Resonant torques between disk stars and any new linear non-axisymmetric mode are suppressed within a scar, as is epicyclic heating. Since all the resonances of a linear mode with the same angular wavenumber and pattern speed as its precursor lie inside the scars of the precursor, it suffers less damping. Hence, scars filter the spectrum of noise-generated modes, promoting the renewal of a few select modes. Relic scars sustained by a galaxy disk, due to past tidal interaction with a passing companion, may still be active enablers of non-axisymmetric modes, such as the two-armed "grand design" spiral patterns.

Cosmic rays play a pivotal role in launching galactic winds, particularly in quiescently star-forming galaxies where the hot gas alone is not sufficient to drive a wind. Except for the Milky Way, not much is known about the transport of cosmic rays in galaxies. In this Letter, we present low-frequency observations of the nearby edge-on spiral galaxy NGC 4565 using the LOw-Frequency ARray (LOFAR). With our deep 144-MHz observations, we obtain a clean estimate of the emission originating from old cosmic-ray electrons (CRe), which is almost free from contamination by thermal emission. We measured vertical profiles of the non-thermal radio continuum emission that we fitted with Gaussian and exponential functions. The different profile shapes correspond to 1D cosmic-ray transport models of pure diffusion and advection, respectively. We detect a warp in the radio continuum that is reminiscent of the previously known HI warp. Because the warp is not seen at GHz-frequencies in the radio continuum, its minimum age must be about 100 Myr. The warp also explains the slight flaring of the thick radio disc that can otherwise be well described by a Gaussian profile with an FWHM of 65 arcsec (3.7 kpc). The diffusive radio halo together with the extra-planar X-ray emission may be remnants of enhanced star-forming activity in the past where the galaxy had a galactic wind, as GHz-observations indicate only a weak outflow in the last 40 Myr. NGC 4565 could be in transition from an outflow- to an inflow-dominated phase.

The galaxy distributions along the line-of-sight are significantly contaminated by the uncertainty on redshift measurements obtained through multiband photometry, which makes it difficult to get cosmic distance information measured from baryon acoustic oscillations, or growth functions probed by redshift distortions. We investigate the propagation of the uncertainties into large scale clustering by exploiting all known estimators, and propose the wedge approach as a promising analysis tool to extract cosmic distance information still remaining in the photometric galaxy samples. We test our method using simulated galaxy maps with photometric uncertainties of $\sigma_{0} =\left(0.01, 0.02, 0.03\right)$. The measured anisotropy correlation function $\xi$ is binned into the radial direction of $s$ and the angular direction of $\mu$, and the variations of $\xi(s,\mu)$ with perpendicular and radial cosmic distance measures of $D_A$ and $H^{-1}$ are theoretically estimated by an improved RSD model. Although the radial cosmic distance $H^{-1}$ is unable to be probed from any of the three photometric galaxy samples, the perpendicular component of $D_A$ is verified to be accurately measured even after the full marginalisation of $H^{-1}$. We measure $D_A$ with approximately 6% precision which is nearly equivalent to what we can expect from spectroscopic DR12 CMASS galaxy samples.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

While Type Ia Supernovae (SNe Ia) are one of the most mature cosmological probes, the next era promises to be extremely exciting in the number of different ways SNe Ia are used to measure various cosmological parameters. Here we review the experiments in the 2020s that will yield orders of magnitudes more SNe Ia, and the new understandings and capabilities to constrain systematic uncertainties at a level to match these statistics. We then discuss five different cosmological probes with SNe Ia: the conventional Hubble diagram for measuring dark energy properties, the distance ladder for measuring the Hubble constant, peculiar velocities and weak lensing for measuring sigma8 and strong-lens measurements of H0 and other cosmological parameters. For each of these probes, we discuss the experiments that will provide the best measurements and also the SN Ia-related systematics that affect each one.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

We explore the growth of large-scale magnetic fields in a shear flow, due to helicity fluctuations with a finite correlation time, through a study of the Kraichnan-Moffatt model of zero-mean stochastic fluctuations of the $\alpha$ parameter of dynamo theory. We derive a linear integro-differential equation for the evolution of large-scale magnetic field, using the first-order smoothing approximation and the Galilean invariance of the $\alpha$-statistics. This enables construction of a model that is non-perturbative in the shearing rate $S$ and the $\alpha$-correlation time $\tau_\alpha$. After a brief review of the salient features of the exactly solvable white-noise limit, we consider the case of small but non-zero $\tau_\alpha$. When the large-scale magnetic field varies slowly, the evolution is governed by a partial differential equation. We present modal solutions and conditions for the exponential growth rate of the large-scale magnetic field, whose drivers are the Kraichnan diffusivity, Moffatt drift, Shear and a non-zero correlation time. Of particular interest is dynamo action when the $\alpha$-fluctuations are weak; i.e. when the Kraichnan diffusivity is positive. We show that in the absence of Moffatt drift, shear does not give rise to growing solutions. But shear and Moffatt drift acting together can drive large scale dynamo action with growth rate $\gamma \propto |S|$.

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45$^\circ$00$'$00$''$ to 57$^\circ$00$'$00$''$) were mapped using a fully automated direction-dependent calibration and imaging pipeline that we developed. A total of 325,694 sources are detected with a signal of at least five times the noise, and the source density is a factor of $\sim 10$ higher than the most sensitive existing very wide-area radio-continuum surveys. The median sensitivity is S$_{\rm 144 MHz} = 71\,\mu$Jy beam$^{-1}$ and the point-source completeness is 90% at an integrated flux density of 0.45mJy. The resolution of the images is 6$''$ and the positional accuracy is within 0.2$''$. This data release consists of a catalogue containing location, flux, and shape estimates together with 58 mosaic images that cover the catalogued area. In this paper we provide an overview of the data release with a focus on the processing of the LOFAR data and the characteristics of the resulting images. In two accompanying papers we provide the radio source associations and deblending and, where possible, the optical identifications of the radio sources together with the photometric redshifts and properties of the host galaxies. These data release papers are published together with a further $\sim$20 articles that highlight the scientific potential of LoTSS.

We apply the Tremaine-Weinberg theory of dynamical friction to compute the orbital decay of a globular cluster (GC), on an initially circular orbit inside a cored spherical galaxy with isotropic stellar velocities. The retarding torque on the GC, T(rp) < 0 , is a function of its orbital radius rp . The torque is exerted by stars whose orbits are resonant with the GC's orbit, and given as a sum over the infinitely many possible resonances by the Lynden-Bell Kalnajs (LBK) formula. We calculate the LBK torque T(rp) and determine rp(t), for a GC of mass Mp = 2 x 10^5 M_sun and an Isochrone galaxy of core mass Mc = 4 x 10^8 M_sun and core radius b = 1000pc. (i) When rp > 300 pc many strong resonances are active and, as expected, T = T_C , the classical Chandrasekhar torque. (ii) For rp < 300 pc, T comes mostly from stars nearly co-rotating with the GC, trailing or leading it slightly; Trailing resonances exert stronger torques. (iii) As rp decreases the number and strength of resonances drop, so |T| also decreases, with |T| < 10^{-2} |T_C| at rp = r* = (Mp/Mc)^{1/5} b = 220 pc , a characteristic `filtering' radius. (iv) Many resonances cease to exist inside r* ; this includes all Leading and low-order Trailing ones. (v) The higher-order Trailing resonances inside r* are very weak, with |T| < 10^{-4} |T_C| at rp = 150 pc. (vi) Inspiral times for rp(t) to decay from 300 pc to r* far exceed 10 Gyr.

Rotation measure (RM) synthesis is a widely used polarization processing algorithm for reconstructing polarized structures along the line of sight. Performing RM synthesis on large datasets produced by telescopes like LOFAR can be computationally intensive as the computational cost is proportional to the product of the number of input frequency channels, the number of output Faraday depth values to be evaluated and the number of lines of sight present in the data cube. The required computational cost is likely to get worse due to the planned large area sky surveys with telescopes like the Low Frequency Array (LOFAR), the Murchison Widefield Array (MWA), and eventually the Square Kilometre Array (SKA). The massively parallel General Purpose Graphical Processing Units (GPGPUs) can be used to execute some of the computationally intensive astronomical image processing algorithms including RM synthesis. In this paper, we present a GPU-accelerated code, called cuFFS or CUDA-accelerated Fast Faraday Synthesis, to perform Faraday rotation measure synthesis. Compared to a fast single-threaded and vectorized CPU implementation, depending on the structure and format of the data cubes, our code achieves an increase in speed of up to two orders of magnitude. During testing, we noticed that the disk I/O when using the Flexible Image Transport System (FITS) data format is a major bottleneck and to reduce the time spent on disk I/O, our code supports the faster HDFITS format in addition to the standard FITS format. The code is written in C with GPU-acceleration achieved using Nvidia's CUDA parallel computing platform. The code is available at https://github.com/sarrvesh/cuFFS.

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.

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

The new generation of broad-band radio continuum surveys will provide large data sets with polarization information. New algorithms need to be developed to extract reliable catalogs of linearly polarized sources that can be used to characterize those sources and produce a dense rotation measure (RM) grid to probe magneto-ionized structures along the line of sight via Faraday rotation. The aim of the paper is to develop a computationally efficient and rigorously defined source-finding algorithm for linearly polarized sources. We used a calibrated data set from the LOw Frequency ARray (LOFAR) at 150 MHz centered on the nearby galaxy M51 to search for polarized background sources. With a new imaging software, we re-imaged the field at a resolution of 18''x15'' and cataloged a total of about 3000 continuum sources within 2.5 degrees of the center of M51. We made small Stokes Q and U images centered on each source brighter than 100 mJy in total intensity (201 sources) and used RM synthesis to create corresponding Faraday cubes that were analyzed individually. For each source, the noise distribution function was determined from a subset of measurements at high Faraday depths where no polarization is expected; the peaks in polarized intensity in the Faraday spectrum were identified and the p-value of each source was calculated. Finally, the false discovery rate method was applied to the list of p-values to produce a list of polarized sources and quantify the reliability of the detections. We also analyzed sources fainter than 100 mJy but that were reported as polarized in the literature at at least another radio frequency. Of the 201 sources that were searched for polarization, 6 polarized sources were detected confidently (with a false discovery rate of 5 percent). This corresponds to a number density of 1 polarized source per 3.3 square degrees, or 0.3 source/square degree. [Abridged]

Low-frequency radio continuum observations of edge-on galaxies are ideal to study cosmic-ray electrons (CREs) in halos via radio synchrotron emission and to measure magnetic field strengths. We obtained new observations of the edge-on spiral galaxy NGC 891 at 129-163 MHz with the LOw Frequency ARray (LOFAR) and at 13-18 GHz with the Arcminute Microkelvin Imager (AMI) and combine them with recent high-resolution Very Large Array (VLA) observations at 1-2 GHz, enabling us to study the radio continuum emission over two orders of magnitude in frequency. The spectrum of the integrated nonthermal flux density can be fitted by a power law with a spectral steepening towards higher frequencies or by a curved polynomial. Spectral flattening at low frequencies due to free-free absorption is detected in star-forming regions of the disk. The mean magnetic field strength in the halo is 7 +- 2 $\mu$G. The scale heights of the nonthermal halo emission at 146 MHz are larger than those at 1.5 GHz everywhere, with a mean ratio of 1.7 +- 0.3, indicating that spectral ageing of CREs is important and that diffusive propagation dominates. The halo scale heights at 146 MHz decrease with increasing magnetic field strengths which is a signature of dominating synchrotron losses of CREs. On the other hand, the spectral index between 146 MHz and 1.5 GHz linearly steepens from the disk to the halo, indicating that advection rather than diffusion is the dominating CRE transport process. This issue calls for refined modelling of CRE propagation.

The nuclear star cluster surrounding the massive black hole at the Galactic Centre consists of young and old stars, with most of the stellar mass in an extended, cuspy distribution of old stars. The compact cluster of young stars was probably born in situ in a massive accretion disc around the black hole. We investigate the effect of the growing gravity of the disc on the orbits of the old stars, using an integrable model of the deformation of a spherical star cluster with anisotropic velocity dispersions. A formula for the perturbed phase space distribution function is derived using linear theory, and new density and surface density profiles are computed. The cusp undergoes a spheroidal deformation with the flattening increasing strongly at smaller distances from the black hole; the intrinsic axis ratio ~ 0.8 at ~ 0.15 pc. Stellar orbits are deformed such that they spend more time near the disc plane and sample the dense inner parts of the disc; this could result in enhanced stripping of the envelopes of red giant stars. Linear theory accounts only for orbits whose apsides circulate. The non-linear theory of adiabatic capture into resonance is needed to understand orbits whose apsides librate. The mechanism is a generic dynamical process, and it may be common in galactic nuclei.

We present idealized models of a razor-thin, axisymmetric, Keplerian stellar disc around a massive black hole, and study non-axisymmetric secular instabilities in the absence of either counter-rotation or loss cones. These discs are prograde mono-energetic waterbags, whose phase space distribution functions are constant for orbits within a range of eccentricities (e) and zero outside this range. The linear normal modes of waterbags are composed of sinusoidal disturbances of the edges of distribution function in phase space. Waterbags which include circular orbits (polarcaps) have one stable linear normal mode for each azimuthal wavenumber m. The m = 1 mode always has positive pattern speed and, for polarcaps consisting of orbits with e < 0.9428, only the m = 1 mode has positive pattern speed. Waterbags excluding circular orbits (bands) have two linear normal modes for each m, which can be stable or unstable. We derive analytical expressions for the instability condition, pattern speeds, growth rates and normal mode structure. Narrow bands are unstable to modes with a wide range in m. Numerical simulations confirm linear theory and follow the non-linear evolution of instabilities. Long-time integration suggests that instabilities of different m grow, interact non-linearly and relax collisionlessly to a coarse-grained equilibrium with a wide range of eccentricities.

We present exact solutions of the incompressible Navier-Stokes equations in a background linear shear flow. The method of construction is based on Kelvin's investigations into linearized disturbances in an unbounded Couette flow. We obtain explicit formulae for all three components of a Kelvin mode in terms of elementary functions. We then prove that Kelvin modes with parallel (though time-dependent) wave vectors can be superposed to construct the most general plane transverse shearing wave. An explicit solution is given, with any specified initial orientation, profile and polarization structure, with either unbounded or shear-periodic boundary conditions.

We present a new analysis of the widely used relation between cavity power and radio luminosity in clusters of galaxies with evidence for strong AGN feedback. We study the correlation at low radio frequencies using two new surveys - the First Alternative Data Release of the TIFR GMRT Sky Survey (TGSS ADR1) at 148 MHz and LOFAR's first all-sky survey, the Multifrequency Snapshot Sky Survey (MSSS) at 140 MHz. We find a scaling relation $P_{\rm cav} \propto L_{148}^{\beta}$, with a logarithmic slope of $\beta = 0.51 \pm 0.14$, which is in good agreement with previous results based on data at 327 MHz. The large scatter present in this correlation confirms the conclusion reached at higher frequencies that the total radio luminosity at a single frequency is a poor predictor of the total jet power. We show that including measurements at 148 MHz alone is insufficient to reliably compute the bolometric radio luminosity and reduce the scatter in the correlation. For a subset of four well-resolved sources, we examine the detected extended structures at low frequencies and compare with the morphology known from higher frequency images and Chandra X-ray maps. In Perseus we discuss details in the structures of the radio mini-halo, while in the 2A 0335+096 cluster we observe new diffuse emission associated with multiple X-ray cavities and likely originating from past activity. For A2199 and MS 0735.6+7421, we confirm that the observed low-frequency radio lobes are confined to the extents known from higher frequencies. This new low-frequency analysis highlights the fact that existing cavity power to radio luminosity relations are based on a relatively narrow range of AGN outburst ages. We discuss how the correlation could be extended using low frequency data from the LOFAR Two-metre Sky Survey (LoTSS) in combination with future, complementary deeper X-ray observations.

We aim to use statistical analysis of a large number of various galaxies to probe, model, and understand relations between different galaxy properties and magnetic fields. We have compiled a sample of 55 galaxies including low-mass dwarf and Magellanic-types, normal spirals and several massive starbursts, and applied principal component analysis (PCA) and regression methods to assess the impact of various galaxy properties on the observed magnetic fields. According to PCA the global galaxy parameters (like HI, H2, and dynamical mass, star formation rate (SFR), near-infrared luminosity, size, and rotational velocity) are all mutually correlated and can be reduced to a single principal component. Further PCA performed for global and intensive (not size related) properties of galaxies (such as gas density, and surface density of the star formation rate, SSFR), indicates that magnetic field strength B is connected mainly to the intensive parameters, while the global parameters have only weak relationships with B. We find that the tightest relationship of B is with SSFR, which is described by a power-law with an index of 0.33+-0.03. The observed weaker associations of B with galaxy dynamical mass and the rotational velocity we interpret as indirect ones, resulting from the observed connection of the global SFR with the available total H2 mass in galaxies. Using our sample we constructed a diagram of B across the Hubble sequence which reveals that high values of B are not restricted by the Hubble type. However, weaker fields appear exclusively in later Hubble types and B as low as about 5muG is not seen among typical spirals. The processes of generation of magnetic field in the dwarf and Magellanic-type galaxies are similar to those in the massive spirals and starbursts and are mainly coupled to local star-formation activity involving the small-scale dynamo mechanism.

We report on the discovery in the LOFAR Multifrequency Snapshot Sky Survey (MSSS) of a giant radio galaxy (GRG) with a projected size of $2.56 \pm 0.07$ Mpc projected on the sky. It is associated with the galaxy triplet UGC 9555, within which one is identified as a broad-line galaxy in the Sloan Digital Sky Survey (SDSS) at a redshift of $0.05453 \pm 1 \times 10^{-5} $, and with a velocity dispersion of $215.86 \pm 6.34$ km/s. From archival radio observations we see that this galaxy hosts a compact flat-spectrum radio source, and we conclude that it is the active galactic nucleus (AGN) responsible for generating the radio lobes. The radio luminosity distribution of the jets, and the broad-line classification of the host AGN, indicate this GRG is orientated well out of the plane of the sky, making its physical size one of the largest known for any GRG. Analysis of the infrared data suggests that the host is a lenticular type galaxy with a large stellar mass ($\log~\mathrm{M}/\mathrm{M}_\odot = 11.56 \pm 0.12$), and a moderate star formation rate ($1.2 \pm 0.3~\mathrm{M}_\odot/\mathrm{year}$). Spatially smoothing the SDSS images shows the system around UGC 9555 to be significantly disturbed, with a prominent extension to the south-east. Overall, the evidence suggests this host galaxy has undergone one or more recent moderate merger events and is also experiencing tidal interactions with surrounding galaxies, which have caused the star formation and provided the supply of gas to trigger and fuel the Mpc-scale radio lobes.

We investigate to which accuracy it is possible to recover the real-space two-point correlation function of galaxy clusters from cluster catalogues based on photometric redshifts, and test our ability to measure the redshift and mass evolution of the correlation length and the bias parameter as a function of the redshift uncertainty. We calculate the correlation function for cluster sub-samples covering various mass and redshift bins selected from a light-cone catalogue. To simulate the distribution of clusters in photometric redshift space, we assign to each cluster a redshift randomly extracted from a Gaussian distribution. The dispersion is varied in the range $\sigma_{(z=0)} = 0.001$ to $0.050$. The correlation function in real-space is computed through estimation and deprojection of $w_{p}(r_{p})$. Four mass ranges (from $M_{halo}> 2 \times 10^{13}$ to $M_{halo}> 2 \times 10^{14}$) and six redshift slices covering the redshift range [0,2] are investigated, using cosmological redshifts and photo-z configurations. We find a clear increase of the correlation amplitude as a function of redshift and mass for the $z_{c}$ samples. The evolution of the derived bias parameter is in agreement with theoretical expectations. From our pilot sample limited to $M_{halo}> 5 \times 10^{13} (0.4 < z < 0.7)$, we find that the real-space correlation function can be recovered by deprojection of $w_{p}(r_{p})$ within an accuracy of 5% for $\sigma_{z} = 0.001 \times (1 + z_{c})$ and within 10% for $\sigma_{z} = 0.03 \times (1 + z_{c})$. The evolution of the correlation in redshift and mass is clearly detected for all $\sigma_{z}$ tested. The best-fit parameters $(r_{0}$ and ${\gamma})$ as well as the bias obtained from the deprojection method for all $\sigma_{z}$ are within the $1 \sigma$ uncertainty of the $z_{c}$ sample.

The new generation of low-frequency radio telescopes, such as the Low Frequency Array (LOFAR: a Square Kilometre Array-low pathfinder), provides advancements in our capability of probing Galactic magnetism through low-frequency polarimetry. Maps of diffuse polarized radio emission and Faraday rotation can be used to infer properties of, and trace structure in, the magnetic fields in the ISM. However, to date very little of the sky has been probed at high angular and Faraday depth resolution. We observed a 5x5 degree region centred on the nearby galaxy IC342 using LOFAR in the frequency range 115-178 MHz at 4 arcmin resolution and performed Faraday tomography to detect foreground Galactic polarized synchrotron emission separated by Faraday depth (different amounts of Faraday rotation). Our Faraday depth cube shows rich polarized structure, with up to 30 K of polarized emission at 150 MHz. We detect two overlapping diffuse polarized features that are clearly separated in Faraday depth. Faraday-thick structures at such low frequencies would be too strongly depolarized to explain the observations and are therefore rejected. Only Faraday thin structures will not be strongly depolarized; producing such structures requires localized variations in the ratio of synchrotron emissivity to Faraday depth per unit distance, which can arise from several physical phenomena, such as a transition between regions of ionized and neutral gas. We conclude that the observed polarized emission is Faraday thin, and propose that the emission originates from two neutral clouds in the local ISM. We have modeled the Faraday rotation for this line of sight and estimated that the line of sight component of magnetic field of the local ISM for this direction varies between -0.86 and +0.12 uG. We propose that this may be a useful method for mapping magnetic fields within the local ISM.

The LOFAR Two-metre Sky Survey (LoTSS) is a deep 120-168 MHz imaging survey that will eventually cover the entire Northern sky. Each of the 3170 pointings will be observed for 8 hrs, which, at most declinations, is sufficient to produce ~5arcsec resolution images with a sensitivity of ~0.1mJy/beam and accomplish the main scientific aims of the survey which are to explore the formation and evolution of massive black holes, galaxies, clusters of galaxies and large-scale structure. Due to the compact core and long baselines of LOFAR, the images provide excellent sensitivity to both highly extended and compact emission. For legacy value, the data are archived at high spectral and time resolution to facilitate subarcsecond imaging and spectral line studies. In this paper we provide an overview of the LoTSS. We outline the survey strategy, the observational status, the current calibration techniques, a preliminary data release, and the anticipated scientific impact. The preliminary images that we have released were created using a fully-automated but direction-independent calibration strategy and are significantly more sensitive than those produced by any existing large-area low-frequency survey. In excess of 44,000 sources are detected in the images that have a resolution of 25arcsec, typical noise levels of less than 0.5 mJy/beam, and cover an area of over 350 square degrees in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45d00m00s to 57d00m00s).

We present the first wide area (19 deg$^2$), deep ($\approx120-150$ {\mu}Jy beam$^{-1}$), high resolution ($5.6 \times 7.4$ arcsec) LOFAR High Band Antenna image of the Bo\"otes field made at 130-169 MHz. This image is at least an order of magnitude deeper and 3-5 times higher in angular resolution than previously achieved for this field at low frequencies. The observations and data reduction, which includes full direction-dependent calibration, are described here. We present a radio source catalogue containing 6276 sources detected over an area of $19$\,deg$^2$, with a peak flux density threshold of $5\sigma$. As the first thorough test of the facet calibration strategy, introduced by van Weeren et al., we investigate the flux and positional accuracy of the catalogue. We present differential source counts that reach an order of magnitude deeper in flux density than previously achieved at these low frequencies, and show flattening at 150 MHz flux densities below 10 mJy associated with the rise of the low flux density star-forming galaxies and radio-quiet AGN.

A key ingredient in the evolution of galaxies is the star formation cycle. Recent progress in the study of magnetic fields is revealing the close connection between star formation and its effect on the small-scale structure in the magnetized interstellar medium (ISM). In this contribution we describe how the modern generation of radio telescopes is being used to probe the physics of the ISM through sensitive multiwavelength surveys of gas and magnetic fields, from the inner star forming disk and outward into the galaxy outskirts where large-scale magnetic fields may also play a key role. We highlight unique pioneering efforts towards performing and scientifically exploiting large-scale surveys of the type that the SKA will undertake routinely. Looking to the future, we describe plans for using the Square Kilometre Array (SKA) and its pathfinders to gain important new insights into the cosmic history of galaxy evolution.

With Low-Frequency Array (LOFAR) observations, we have discovered a diverse assembly of steep spectrum emission that is apparently associated with the intra cluster medium (ICM) of the merging galaxy cluster Abell 2034. Such a rich variety of complex emission associated with the ICM has been observed in few other clusters. This not only indicates that Abell 2034 is a more interesting and complex system than previously thought but it also demonstrates the importance of sensitive and high-resolution, low-frequency observations. These observations can reveal emission from relativistic particles which have been accelerated to sufficient energy to produce observable emission or have had their high energy maintained by mechanisms in the ICM. The most prominent feature in our maps is a bright bulb of emission connected to two steep spectrum filamentary structures, the longest of which extends perpendicular to the merger axis for 0.5Mpc across the south of the cluster. The origin of these objects is unclear, with no shock detected in the X-ray images and no obvious connection with cluster galaxies or AGNs. We also find that the X-ray bright region of the cluster coincides with a giant radio halo with an irregular morphology and a very steep spectrum. In addition, the cluster hosts up to three possible radio relics, which are misaligned with the cluster X-ray emission. Finally, we have identified multiple regions of emission with a very steep spectral index that seem to be associated with either tailed radio galaxies or a shock.

We present a first-principles theory of Resonant Relaxation (RR) of a low mass stellar system orbiting a more massive black hole (MBH). We first extend the kinetic theory of Gilbert (1968) to include the Keplerian field of a black hole of mass $M_\bullet$. Specializing to a Keplerian stellar system of mass $M \ll M_\bullet$, we use the orbit-averaging method of Sridhar & Touma (2015; Paper I) to derive a kinetic equation for RR. This describes the collisional evolution of a system of $N \gg 1$ Gaussian Rings in a reduced 5-dim space, under the combined actions of self-gravity, 1 PN and 1.5 PN relativistic effects of the MBH and an arbitrary external potential. In general geometries RR is driven by both apsidal and nodal resonances, so the distinction between scalar-RR and vector-RR disappears. The system passes through a sequence of quasi-steady secular collisionless equilibria, driven by irreversible 2-Ring correlations that accrue through gravitational interactions, both direct and collective. This correlation function is related to a `wake function', which is the linear response of the system to the perturbation of a chosen Ring. The wake function is easier to appreciate, and satisfies a simpler equation, than the correlation function. We discuss general implications for the interplay of secular dynamics and non-equilibrium statistical mechanics in the evolution of Keplerian stellar systems toward secular thermodynamic equilibria, and set the stage for applications to the RR of axisymmetric discs in Paper III.

We present a theory in three parts, of the secular dynamics of a (Keplerian) stellar system of mass $M$ orbiting a black hole of mass $M_\bullet \gg M$. Here we describe the collisionless dynamics; Papers II and III are on the (collisional) theory of Resonant Relaxation. The mass ratio, $\varepsilon = M/M_\bullet \ll 1$, is a natural small parameter implying a separation of time scales between the short Kepler orbital periods and the longer orbital precessional periods. The collisionless Boltzmann equation (CBE) for the stellar distribution function (DF) is averaged over the fast Kepler orbital phase using the method of multiple scales. The orbit-averaged system is described by a secular DF, $F$, in a reduced phase space. $F$ obeys a secular CBE that includes stellar self-gravity, general relativistic corrections up to 1.5 post-Newtonian order, and external sources varying over secular times. Secular dynamics, even with general time dependence, conserves the semi-major axis of every star. This additional integral of motion promotes extra regularity of the stellar orbits, and enables the construction of equilibria, $F_0$, through a secular Jeans theorem. A linearized secular CBE determines the response and stability of $F_0$. Spherical, non-rotating equilibria may support long-lived, warp-like distortions. We also prove that an axisymmetric, zero-thickness, flat disc is secularly stable to all in-plane perturbations, when its DF, $F_0$, is a monotonic function of the angular momentum at fixed energy.

We study the Resonant Relaxation (RR) of an axisymmetric low mass (or Keplerian) stellar disc orbiting a more massive black hole (MBH). Our recent work on the general kinetic theory of RR is simplified in the standard manner by ignoring the effects of `gravitational polarization', and applied to a zero-thickness, flat, axisymmetric disc. The wake of a stellar orbit is expressed in terms of the angular momenta exchanged with other orbits, and used to derive a kinetic equation for RR under the combined actions of self-gravity, 1 PN and 1.5 PN relativistic effects of the MBH and an arbitrary external axisymmetric potential. This is a Fokker-Planck equation for the stellar distribution function (DF), wherein the diffusion coefficients are given self-consistently in terms of contributions from apsidal resonances between pairs of stellar orbits. The physical kinetics is studied for the two main cases of interest. (1) `Lossless' discs in which the MBH is not a sink of stars, and disc mass, angular momentum and energy are conserved: we prove that general H-functions can increase or decrease during RR, but the Boltzmann entropy is (essentially) unique in being a non-decreasing function of time. Therefore secular thermal equilibria are maximum entropy states, with DFs of the Boltzmann form; the two-Ring correlation function at equilibrium is computed. (2) Discs that lose stars to the MBH through an `empty loss-cone': we derive expressions for the MBH feeding rates of mass, angular momentum and energy in terms of the diffusive flux at the loss-cone boundary.

We present deep LOFAR observations between 120-181 MHz of the "Toothbrush" (RX J0603.3+4214), a cluster that contains one of the brightest radio relic sources known. Our LOFAR observations exploit a new and novel calibration scheme to probe 10 times deeper than any previous study in this relatively unexplored part of the spectrum. The LOFAR observations, when combined with VLA, GMRT, and Chandra X-ray data, provide new information about the nature of cluster merger shocks and their role in re-accelerating relativistic particles. We derive a spectral index of $\alpha = -0.8 \pm 0.1$ at the northern edge of the main radio relic, steepening towards the south to $\alpha \approx - 2$. The spectral index of the radio halo is remarkably uniform ($\alpha = -1.16$, with an intrinsic scatter of $\leq 0.04$). The observed radio relic spectral index gives a Mach number of $\mathcal{M} = 2.8^{+0.5}_{-0.3}$, assuming diffusive shock acceleration (DSA). However, the gas density jump at the northern edge of the large radio relic implies a much weaker shock ($\mathcal{M} \approx 1.2$, with an upper limit of $\mathcal{M} \approx 1.5$). The discrepancy between the Mach numbers calculated from the radio and X-rays can be explained if either (i) the relic traces a complex shock surface along the line of sight, or (ii) if the radio relic emission is produced by a re-accelerated population of fossil particles from a radio galaxy. Our results highlight the need for additional theoretical work and numerical simulations of particle acceleration and re-acceleration at cluster merger shocks.

LOFAR, the Low-Frequency Array, is a powerful new radio telescope operating between 10 and 240 MHz. LOFAR allows detailed sensitive high-resolution studies of the low-frequency radio sky. At the same time LOFAR also provides excellent short baseline coverage to map diffuse extended emission. However, producing high-quality deep images is challenging due to the presence of direction dependent calibration errors, caused by imperfect knowledge of the station beam shapes and the ionosphere. Furthermore, the large data volume and presence of station clock errors present additional difficulties. In this paper we present a new calibration scheme, which we name facet calibration, to obtain deep high-resolution LOFAR High Band Antenna images using the Dutch part of the array. This scheme solves and corrects the direction dependent errors in a number of facets that cover the observed field of view. Facet calibration provides close to thermal noise limited images for a typical 8 hr observing run at $\sim$ 5arcsec resolution, meeting the specifications of the LOFAR Tier-1 northern survey.

Most exoplanetary systems in binary stars are of S--type, and consist of one or more planets orbiting a primary star with a wide binary stellar companion. Gravitational forcing of a single planet by a sufficiently inclined binary orbit can induce large amplitude oscillations of the planet's eccentricity and inclination through the Kozai-Lidov (KL) instability. KL cycling was invoked to explain: the large eccentricities of planetary orbits; the family of close--in hot Jupiters; and the retrograde planetary orbits in eccentric binary systems. However, several kinds of perturbations can quench the KL instability, by inducing fast periapse precessions which stabilize circular orbits of all inclinations: these could be a Jupiter--mass planet, a massive remnant disc or general relativistic precession. Indeed, mutual gravitational perturbations in multiplanet S--type systems can be strong enough to lend a certain dynamical rigidity to their orbital planes. Here we present a new and faster process that is driven by this very agent inhibiting KL cycling. Planetary perturbations enable secular oscillations of planetary eccentricities and inclinations, also called Laplace--Lagrange (LL) eigenmodes. Interactions with a remnant disc of planetesimals can make planets migrate, causing a drift of LL mode periods which can bring one or more LL modes into resonance with binary orbital motion. The results can be dramatic, ranging from excitation of large eccentricities and mutual inclinations to total disruption. Not requiring special physical or initial conditions, binary resonant driving is generic and could have profoundly altered the architecture of many S--type multiplanet systems. It can also weaken the multiplanet occurrence rate in wide binaries, and affect planet formation in close binaries.

We measure the projected 2-point correlation function of galaxies in the 180 deg$^2$ equatorial regions of the GAMA II survey, for four different redshift slices between z = 0.0 and z=0.5. To do this we further develop the Cole (2011) method of producing suitable random catalogues for the calculation of correlation functions. We find that more r-band luminous, more massive and redder galaxies are more clustered. We also find that red galaxies have stronger clustering on scales less than ~3 $h^{-1}$ Mpc. We compare to two different versions of the GALFORM galaxy formation model, Lacey et al (in prep.) and Gonzalez-Perez et al. (2014), and find that the models reproduce the trend of stronger clustering for more massive galaxies. However, the models under predict the clustering of blue galaxies, can incorrectly predict the correlation function on small scales and under predict the clustering in our sample of galaxies with ~3$L_r$ . We suggest possible avenues to explore to improve these cluster- ing predictions. The measurements presented in this paper can be used to test other galaxy formation models, and we make the measurements available online to facilitate this.

We present ultraviolet (UV) observations of six nearby Type~Ia supernovae (SNe~Ia) obtained with the Hubble Space Telescope, three of which were also observed in the near-IR (NIR) with Wide-Field Camera~3. UV observations with the Swift satellite, as well as ground-based optical and near-infrared data provide complementary information. The combined data-set covers the wavelength range $0.2$--$2~\mu$m. By also including archival data of SN 2014J, we analyse a sample spanning observed colour excesses up to $E(B-V)=1.4~$mag. We study the wavelength dependent extinction of each individual SN and find a diversity of reddening laws when characterised by the total-to-selective extinction $R_V$. In particular, we note that for the two SNe with $E(B-V)\gtrsim1~$mag, for which the colour excess is dominated by dust extinction, we find $R_V=1.4\pm0.1$ and $R_V=2.8\pm0.1$. Adding UV photometry reduces the uncertainty of fitted $R_V$ by $\sim50\,$% allowing us to also measure $R_V$ of individual low-extinction objects which point to a similar diversity, currently not accounted for in the analyses when SNe~Ia are used for studying the expansion history of the universe.

We present a model of large-scale dynamo action in a shear flow that has stochastic, zero-mean fluctuations of the $\alpha$ parameter. This is based on a minimal extension of the Kraichnan-Moffatt model, to include a background linear shear and Galilean-invariant $\alpha$-statistics. Using the first order smoothing approximation we derive a linear integro-differential equation for the large-scale magnetic field, which is non perturbative in the shearing rate $S\,$, and the $\alpha$-correlation time $\tau_\alpha\,$. The white-noise case, $\tau_\alpha = 0\,$, is solved exactly, and it is concluded that the necessary condition for dynamo action is identical to the Kraichnan-Moffatt model without shear; this is because white-noise does not allow for memory effects, whereas shear needs time to act. To explore memory effects we reduce the integro-differential equation to a partial differential equation, valid for slowly varying fields when $\tau_\alpha$ is small but non zero. Seeking exponential modal solutions, we solve the modal dispersion relation and obtain an explicit expression for the growth rate as a function of the six independent parameters of the problem. A non zero $\tau_\alpha$ gives rise to new physical scales, and dynamo action is completely different from the white-noise case; e.g. even weak $\alpha$ fluctuations can give rise to a dynamo. We argue that, at any wavenumber, both Moffatt drift and Shear always contribute to increasing the growth rate. Two examples are presented: (a) a Moffatt drift dynamo in the absence of shear; (b) a Shear dynamo in the absence of Moffatt drift.

We propose an algorithm that can be used by amateur astronomers to analyze the images acquired during solar eclipses. The proposed algorithm analyzes the image, detects the eclipse and produces results for parameters like magnitude of eclipse, eclipse obscuration and the approximate distance between the Earth and the Moon.

We study mean field dynamo action in renovating flows with finite and non zero correlation time ($\tau$) in the presence of shear. Previous results obtained when shear was absent are generalized to the case with shear. The question of whether the mean magnetic field can grow in the presence of shear and non helical turbulence, as seen in numerical simulations, is examined. We show in a general manner that, if the motions are strictly non helical, then such mean field dynamo action is not possible. This result is not limited to low (fluid or magnetic) Reynolds numbers nor does it use any closure approximation; it only assumes that the flow renovates itself after each time interval $\tau$. Specifying to a particular form of the renovating flow with helicity, we recover the standard dispersion relation of the $\alpha^2 \Omega$ dynamo, in the small $\tau$ or large wavelength limit. Thus mean fields grow even in the presence of rapidly growing fluctuations, surprisingly, in a manner predicted by the standard quasilinear closure, even though such a closure is not strictly justified. Our work also suggests the possibility of obtaining mean field dynamo growth in the presence of helicity fluctuations, although having a coherent helicity will be more efficient.

We formulate the collisionless Boltzmann equation (CBE) for dense star clusters that lie within the radius of influence of a massive black hole in galactic nuclei. Our approach to these nearly Keplerian systems follows that of Sridhar and Touma (1999): Delaunay canonical variables are used to describe stellar orbits and we average over the fast Keplerian orbital phases. The stellar distribution function (DF) evolves on the longer time scale of precessional motions, whose dynamics is governed by a Hamiltonian, given by the orbit-averaged self-gravitational potential of the cluster. We specialize to razor-thin, planar discs and consider two counter-rotating ("$\pm$") populations of stars. To describe discs of small eccentricities, we expand the $\pm$ Hamiltonian to fourth order in the eccentricities, with coefficients that depend self-consistently on the $\pm$ DFs. We construct approximate $\pm$ dynamical invariants and use Jeans' theorem to construct time-dependent $\pm$ DFs, which are completely described by their centroid coordinates and shape matrices. When the centroid eccentricities are larger than the dispersion in eccentricities, the $\pm$ centroids obey a set of 4 autonomous ordinary differential equations. We show that these can be cast as a two-degree of freedom Hamiltonian system which is nonlinear, yet integrable. We study the linear instability of initially circular discs and derive a criterion for the counter-rotating instability. We then explore the rich nonlinear dynamics of counter-rotating discs, with focus on the variety of steadily precessing eccentric configurations that are allowed. The stability and properties of these configurations are studied as functions of parameters such as the disc mass ratios and angular momentum.

We study the linear $m=1$ counter-rotating instability in a two-component, nearly Keplerian disc. Our goal is to understand these \emph{slow} modes in discs orbiting massive black holes in galactic nuclei. They are of interest not only because they are of large spatial scale--and can hence dominate observations--but also because they can be growing modes that are readily excited by accretion events. Self-gravity being nonlocal, the eigenvalue problem results in a pair of coupled integral equations, which we derive for a two-component softened gravity disc. We solve this integral eigenvalue problem numerically for various values of mass fraction in the counter-rotating component. The eigenvalues are in general complex, being real only in the absence of the counter-rotating component, or imaginary when both components have identical surface density profiles. Our main results are as follows: (i) the pattern speed appears to be non negative, with the growth (or damping) rate being larger for larger values of the pattern speed; (ii) for a given value of the pattern speed, the growth (or damping) rate increases as the mass in the counter-rotating component increases; (iii) the number of nodes of the eigenfunctions decreases with increasing pattern speed and growth rate. Observations of lopsided brightness distributions would then be dominated by modes with the least number of nodes, which also possess the largest pattern speeds and growth rates.

We build on the formulation developed in Sridhar & Singh (JFM, 664, 265, 2010), and present a theory of the \emph{shear dynamo problem} for small magnetic and fluid Reynolds numbers, but for arbitrary values of the shear parameter. Specializing to the case of a mean magnetic field that is slowly varying in time, explicit expressions for the transport coefficients, $\alpha_{il}$ and $\eta_{iml}$, are derived. We prove that, when the velocity field is non helical, the transport coefficient $\alpha_{il}$ vanishes. We then consider forced, stochastic dynamics for the incompressible velocity field at low Reynolds number. An exact, explicit solution for the velocity field is derived, and the velocity spectrum tensor is calculated in terms of the Galilean--invariant forcing statistics. We consider forcing statistics that is non helical, isotropic and delta-correlated-in-time, and specialize to the case when the mean-field is a function only of the spatial coordinate $X_3$ and time $\tau\,$; this reduction is necessary for comparison with the numerical experiments of Brandenburg, R{\"a}dler, Rheinhardt & K\"apyl\"a (ApJ, 676, 740, 2008). Explicit expressions are derived for all four components of the magnetic diffusivity tensor, $\eta_{ij}(\tau)\,$. These are used to prove that the shear-current effect cannot be responsible for dynamo action at small $\re$ and $\rem$, but for all values of the shear parameter.

We study large-scale dynamo action due to turbulence in the presence of a linear shear flow, in the low conductivity limit. Our treatment is nonperturbative in the shear strength and makes systematic use of both the shearing coordinate transformation and the Galilean invariance of the linear shear flow. The velocity fluctuations are assumed to have low magnetic Reynolds number (Rm) but could have arbitrary fluid Reynolds number. The magnetic fluctuations are determined to lowest order in Rm by explicit calculation of the resistive Green's function for the linear shear flow. The mean electromotive force is calculated and an integro-differential equation is derived for the time evolution of the mean magnetic field. In this equation, velocity fluctuations contribute to two different kinds of terms, the C and D terms, in which first and second spatial derivatives of the mean magnetic field, respectively, appear inside the spacetime integrals. The contribution of the D terms is such that the time evolution of the cross-shear components of the mean field do not depend on any other components excepting themselves. Therefore, to lowest order in Rm but to all orders in the shear strength, the D terms cannot give rise to a shear-current assisted dynamo effect. Casting the integro-differential equation in Fourier space, we show that the normal modes of the theory are a set of shearing waves, labelled by their sheared wavevectors. The integral kernels are expressed in terms of the velocity spectrum tensor, which is the fundamental quantity that needs to be specified to complete the integro-differential equation description of the time evolution of the mean magnetic field.

We study large-scale dynamo action due to turbulence in the presence of a linear shear flow. Our treatment is quasilinear and equivalent to the standard `first order smoothing approximation'. However it is non perturbative in the shear strength. We first derive an integro-differential equation for the evolution of the mean magnetic field, by systematic use of the shearing coordinate transformation and the Galilean invariance of the linear shear flow. We show that, for non helical turbulence, the time evolution of the cross-shear components of the mean field do not depend on any other components excepting themselves; this is valid for any Galilean-invariant velocity field, independent of its dynamics. Hence, to all orders in the shear parameter, there is no shear-current type effect for non helical turbulence in a linear shear flow, in quasilinear theory in the limit of zero resistivity. We then develop a systematic approximation of the integro-differential equation for the case when the mean magnetic field varies slowly compared to the turbulence correlation time. For non-helical turbulence, the resulting partial differential equations can again be solved by making a shearing coordinate transformation in Fourier space. The resulting solutions are in the form of shearing waves, labeled by the wavenumber in the sheared coordinates. These shearing waves can grow at early and intermediate times but are expected to decay in the long time limit.

Thin accretion discs around massive compact objects can support slow pressure modes of oscillations in the linear regime that have azimuthal wavenumber $m=1$. We consider finite, flat discs composed of barotropic fluid for various surface density profiles and demonstrate--through WKB analysis and numerical solution of the eigenvalue problem--that these modes are stable and have spatial scales comparable to the size of the disc. We show that the eigenvalue equation can be mapped to a Schr\"odinger-like equation. Analysis of this equation shows that all eigenmodes have discrete spectra. We find that all the models we have considered support negative frequency eigenmodes; however, the positive eigenfrequency modes are only present in power law discs, albeit for physically uninteresting values of the power law index $\beta$ and barotropic index $\gamma$.

I present a review of incompressible magnetohydrodynamic (MHD) turbulence in a strongly magnetised plasma. The approach is phenomenological even where a more rigorous theory is available, so that a reader armed with paper, pencil and some determination may be able to work through most of the physics. The focus is on the inertial-range spectra for very large (fluid and magnetic) Reynolds numbers. These theories of the inertial-range are built on two important facts: (i) Kraichnan's insight that the turbulent cascades are a result of nonlinear interactions between oppositely directed wavepackets of Elsasser fields; (ii) these oppositely directed wavepackets do not exchange energy, but contribute only to changing each other's spatial structures. I begin with a description and critique of the Iroshnikov-Kraichnan theory, and explore the fundamental departures necessitated by the anisotropic nature of the turbulence. Derivations of the inertial-range spectra of four regimes of MHD turbulence -- the balanced weak, balanced strong, imbalanced weak and the imbalanced strong cascades -- are then presented. The need for studying the spectra of imbalanced turbulence when the waves on the outer scale have a short correlation time is briefly discussed.

We present the simplest model that permits a largely analytical exploration of the m=1 counter-rotating instability in a "hot" nearly Keplerian disc of collisionless self-gravitating matter. The model consists of a two-component softened gravity disc, whose linear modes are analysed using WKB. The modes are slow in the sense that their (complex) frequency is smaller than the Keplerian orbital frequency by a factor which is of order the ratio of the disc mass to the mass of the central object. Very simple analytical expressions are derived for the precession frequencies and growth rates of local modes; it is shown that a nearly Keplerian disc must be unrealistically hot to avoid an overstability. Global modes are constructed for the case of zero net rotation.

Large-scale dynamo action due to turbulence in the presence of a linear shear flow is studied. Our treatment is quasilinear and kinematic but is non perturbative in the shear strength. We derive the integro-differential equation for the evolution of the mean magnetic field, by systematic use of the shearing coordinate transformation and the Galilean invariance of the linear shear flow. For non helical turbulence the time evolution of the cross-shear components of the mean field do not depend on any other components excepting themselves. This is valid for any Galilean-invariant velocity field, independent of its dynamics. Hence the shear-current assisted dynamo is essentially absent, although large-scale non helical dynamo action is not ruled out.

We present a phenomenological model of imbalanced MHD turbulence in an incompressible magnetofluid. The steady-state cascades, of waves traveling in opposite directions along the mean magnetic field, carry unequal energy fluxes to small length scales, where they decay due to viscous and resistive dissipation. The inertial-range scalings are well-understood when both cascades are weak. We study the case when both cascades are, in a sense, strong. The inertial-range of this imbalanced cascade has the following properties: (i) the ratio of the r.m.s. Elsasser amplitudes is independent of scale, and is equal to the ratio of the corresponding energy fluxes; (ii) in common with the balanced strong cascade, the energy spectra of both Elsasser waves are of the anisotropic Kolmogorov form, with their parallel correlation lengths equal to each other on all scales, and proportional to the two-thirds power of the transverse correlation length; (iii) the equality of cascade time and waveperiod (critical balance) that characterizes the strong balanced cascade does not apply to the Elsasser field with the larger amplitude. Instead, the more general criterion that always applies to both Elsasser fields is that the cascade time is equal to the correlation time of the straining imposed by oppositely-directed waves. Our results are particularly relevant for turbulence in the solar wind. Spacecraft measurements have established that, in the inertial range of solar wind turbulence, waves travelling away from the sun have higher amplitudes than those travelling towards it. Result (i) allows us to infer the turbulent flux ratios from the amplitude ratios, thus providing insight into the origin of the turbulence.

A strong case has been made that radio waves from sources within about half a degree of the Galactic Center undergo extreme diffractive scattering. However, problems arise when standard (``Kolmogorov'') models of electron density fluctuations are employed to interpret the observations of scattering in conjunction with those of free-free radio emission. Specifically, the outer scale of a Kolmogorov spectrum of electron density fluctuations is constrained to be so small that it is difficult to identify an appropriate astronomical setting. Moreover, an unacceptably high turbulent heating rate results if the outer scale of the velocity field coincides with that of the density fluctuations. We propose an alternative model based on folded magnetic field structures that have been reported in numerical simulations of small-scale dynamos. Nearly isothermal density variations across thin current sheets suffice to account for the scattering. There is no problem of excess turbulent heating because the outer scale for the velocity fluctuations is much larger than the widths of the current sheets. We speculate that interstellar magnetic fields could possess geometries that reflect their origins: fields maintained by the galactic dynamo could have large correlation lengths, whereas those stirred by local energetic events might exhibit folded structures.

The transverse velocities of steadily rotating, non--axisymmetric patterns in flat galaxies may be determined by a purely kinematical method, using two dimensional maps of a tracer surface brightness and radial current density. The data--maps could be viewed as the zeroth and first velocity moments of the line--of--sight velocity distribution, which is the natural output of integral--field spectrographs. Our method is closely related to the Tremaine--Weinberg method of estimating pattern speeds of steadily rotating patterns, when the tracer surface brightness satisfies a source--free continuity equation. We prove that, under identical assumptions about the pattern, two dimensional maps may be used to recover not just one number (the pattern speed), but the full vector field of tracer flow in the disc plane. We illustrate the recovery process by applying it to simulated data, and test its robustness by including the effects of noise.

We present stellar dynamical models of the lopsided, double-peaked nucleus of M31, derived from Hubble Space Telescope (HST) photometry. A Schwarzscild-type method, in conjunction with Richardson-Lucy deconvolution, was employed to construct steadily rotating, hot, stellar disks. The stars orbit a massive dark object, on prograde and retrograde quasi-periodic loop orbits. Our results support Tremaine's eccentric disk model, extended to include a more massive disk, non zero pattern speed ($\Omega$), and different viewing angle. Most of the disk mass populated prograde orbits, with $\simeq 3.4%$ on retrograde orbits. The best fits to photometric and kinematic maps were disks with $\Omega\approx 16\kmspc $. We speculate on the origins of the lopsidedness, invoking recent work on the linear overstability of nearly Keplerian disks, that possess even a small amount of a counter-rotating component. Accretion of material-no more massive than a globular cluster-onto a preexisting stellar disk, will account for the mass in our retrograde orbits, and could have stimulated the lopsidedness seen in the nucleus of M31.

We investigate the orbital structure of a model triaxial star cluster, centered around a supermassive black hole (BH), appropriate to galactic nuclei. Sridhar and Touma (1999) proved that the presence of the BH enforces some regularity in the dynamics within the radius of influence of the BH. We employ their averaging method to reduce the degrees of freedom from three to two. Numerical orbit integrations, together with Poincar\'e surfaces of section allow us to draw a global portrait of the orbital structure; in our calculations we employ a model cluster potential that is triaxial and harmonic. The averaged dynamics of the axisymmetric case is integrable, and we present a detailed comparison of orbits in oblate and prolate axisymmetric potentials. Both cases support resonant orbits with fixed values of eccentricity, inclination, and periapse, whose lines of nodes rotates steadily. We then systematically explore significantly triaxial potentials, possessing small oblateness, or prolateness. Resonant orbits and their families are studied both numerically, and through secular perturbation theory. Chaos appears to be suppressed for all the cases we studied, and we obtain effective third integrals. Some of the orbits appear to reinforce the shape of the potential; we provide phase space, as well as real space portraits of these orbits. A particularly promising resonant orbit exists in highly prolate, triaxial potentials.

The twin peaks in the nucleus of M31 have been interpreted by Tremaine as a thick, eccentric, disk of stars orbiting a massive dark object; the required alignment of the apoapsides of the stellar orbits could be maintained by self-gravity, and the whole structure might be a discrete, nonlinear eigenmode. The pattern speed of this mode could, in principle, be determined by the Tremaine-Weinberg (TW) method, which requires measurements of the surface brightness, and radial velocity along a strip parallel to the line of nodes. However, spectroscopic observations along the line of nodes are not available. We propose a variant of the TW method, which exploits a basic feature of the eccentric disk model, to extract estimates of the pattern speed from {\it Hubble Space Telescope} spectroscopic data, taken along the line joining the two peaks. Within limitations imposed by the data, we estimate that the pattern rotates in a prograde manner and, for an assumed disk inclination of $77^{\circ}$, the pattern speed $|\Omega_p| < 30 {\rm km s^{-1} pc^{-1}}$, or period more than $200,000$ years.

Violent relaxation during the collapse of a galaxy halo is known to be incomplete in realistic cases such as cosmological infall or mergers. We adopt a physical picture of strong but short lived interactions between potential fluctuations and particle orbits, using the broad framework outlined by Tremaine (1987) for incorporating incompleteness of the relaxation. We are guided by results from plasma physics, viz. the quasilinear theory of Landau damping, but allow for significant differences in our case. Crucially, wave particle scattering does not drive the system to an equilibrium distribution function of the exponential type, even in regions of phase space allowed by the constraints. The physical process is mixing without friction in ``action'' space, for which the simplest possible model is a constant phase space density modulated by the constraints. Our distribution function does not use the exponential functions of the energy prevalent in other work, which we regard as inappropriate to collisionless systems. The halo of the self-consistent, parameter-free solutions show an r^(-4) behavior in density at large r, an r^(1/4) surface brightness profile in the region 0.1-8 r_e, and a radially anisotropic velocity dispersion profile outside an isotropic core. The energy distribution seen in simulations N(E) singles out the pericenter cutoff model as the most realistic among the variants we have explored.

We classify orbits of stars that are bound to central black holes in galactic nuclei. The stars move under the combined gravitational influences of the black hole and the central star cluster. Within the sphere of influence of the black hole, the orbital periods of the stars are much shorter than the periods of precession. We average over the orbital motion and end up with a simpler problem and an extra integral of motion: the product of the black hole mass and the semimajor axis of the orbit. Thus the black hole enforces some degree of regularity in its neighborhood. Well within the sphere of influence, (i) planar, as well as three dimensional, axisymmetric configurations-both of which could be lopsided-are integrable, (ii) fully three dimensional clusters with no spatial symmetry whatsover must have semi-regular dynamics with two integrals of motion. Similar considerations apply to stellar orbits when the black hole grows adiabatically. We introduce a family of planar, non-axisymmetric potential perturbations, and study the orbital structure for the harmonic case in some detail. In the centered potentials there are essentially two main families of orbits: the familiar loops and lenses, which were discussed in Sridhar and Touma (1997, MNRAS, 287, L1-L4). We study the effect of lopsidedness, and identify a family of loop orbits, whose orientation reinforces the lopsidedness, an encouraging sign for the construction of self-consistent models of eccentric, discs around black holes, such as in M31 and NGC 4486B.

We construct three dimensional axisymmetric, cuspy density distributions, whose potentials are of St\"ackel form in parabolic coordinates. As in Sridhar and Touma (1997), a black hole of arbitrary mass may be added at the centre, without destroying the St\"ackel form of the potentials. The construction uses a classic method, originally due to Kuzmin (1956), which is here extended to parabolic coordinates. The models are highly oblate, and the cusps are "weak", with the density, $\rho \propto 1/r^k$, where $0<k<1$.

Kraichnan (1965) proposed that MHD turbulence occurs as a result of collisions between oppositely directed Alfv\'en wave packets. Recent work has generated some controversy over the nature of non linear couplings between colliding Alfv\'en waves. We find that the resolution to much of the confusion lies in the existence of a new type of turbulence, intermediate turbulence, in which the cascade of energy in the inertial range exhibits properties intermediate between those of weak and strong turbulent cascades. Some properties of intermediate MHD turbulence are: (i) in common with weak turbulent cascades, wave packets belonging to the inertial range are long lived; (ii) however, components of the strain tensor are so large that, similar to the situation in strong turbulence, perturbation theory is not applicable; (iii) the breakdown of perturbation theory results from the divergence of neighboring field lines due to wave packets whose perturbations in velocity and magnetic fields are localized, but whose perturbations in displacement are not; (iv) 3--wave interactions dominate individual collisions between wave packets, but interactions of all orders $n\geq 3$ make comparable contributions to the intermediate turbulent energy cascade; (v) successive collisions are correlated since wave packets are distorted as they follow diverging field lines; (vi) in common with the weak MHD cascade, there is no parallel cascade of energy, and the cascade to small perpendicular scales strengthens as it reaches higher wave numbers; (vii) For an appropriate weak excitation, there is a natural progression from a weak, through an intermediate, to a strong cascade.

We present cuspy, non-axisymmetric, scale-free mass models of discs, whose gravitational potentials are of St\"ackel form in parabolic coordinates. A black hole may be added at the centre, without in any way affecting the St\"ackel form; the dynamics in these potentials is, of course, fully integrable. The surface density, $\Sigma_{disc}\propto 1/r^{\gamma}$, where $0 <\gamma < 1$ corresponds to steep cusps for which the central force diverges. Thus cusps, black holes, and non-axisymmetry are not a sure recipe for chaos, as is generally assumed. A new family of orbits, lens orbits, emerges to replace the box orbits of models of elliptical galaxies that have constant-density cores. Loop orbits are conspicuous by their absence. Both lenses and boxlets (the other family of orbits), can be elongated in the direction of the density distribution, a property that is favourable for the construction of non-axisymmetric, self-consistent equilibrium models of elliptical galaxies.