We explore the effect of a supermassive black hole (SMBH) on the density profile of a fuzzy dark matter (FDM) soliton core at the centre of a dark matter halo. We numerically solve the Schr\"odinger-Poisson equations, treating the black hole as a gravitational point mass, and demonstrate that this additional perturbing term has a `squeezing' effect on the soliton density profile, decreasing the core radius and increasing the central density. In the limit of large black hole mass, the solution approaches one akin to the hydrogen atom, with radius inversely proportional to the black hole mass. By applying our analysis to two specific galaxies (M87 and the Milky Way) and pairing it with known observational limits on the amount of centrally concentrated dark matter, we obtain a constraint on the FDM particle mass, finding that the range $10^{-22.12}$ eV $\lesssim m \lesssim 10^{-22.06}$ eV should be forbidden (taking into account additional factors concerning the life-time of the soliton in the vicinity of a black hole). Improved observational mass measurements of the black hole and total enclosed masses will significantly extend the lower-bound on the excluded FDM mass region, while self-consistent theoretical modeling of the soliton--black hole system can extend the upper-bound.

We present an in-depth exploration of the phenomenon of dynamical friction in a universe where the dark matter is composed entirely of so-called Fuzzy Dark Matter (FDM), ultralight bosons of mass $m\sim\mathcal{O}(10^{-22})\,$eV. We review the classical treatment of dynamical friction before presenting analytic results in the case of FDM for point masses, extended mass distributions, and FDM backgrounds with finite velocity dispersion. We then test these results against a large suite of fully non-linear simulations that allow us to assess the regime of applicability of the analytic results. We apply these results to a variety of astrophysical problems of interest, including infalling satellites in a galactic dark matter background, and determine that \emph{(1)}~for FDM masses $m\gtrsim 10^{-21}\, {\rm eV}\, c^{-2}$, the timing problem of the Fornax dwarf spheroidal's globular clusters is no longer solved and \emph{(2)}~the effects of FDM on the process of dynamical friction for satellites of total mass $M$ and relative velocity $v_{\rm rel}$ should require detailed numerical simulations for $\left(M/10^9~M_{\odot}\right) \left(m/10^{-22}~{\rm eV}\right)\left(100~{\rm km}~{\rm s}^{-1}/v_{\rm rel}\right) \sim 1$, parameters which would lie outside the validated range of applicability of any currently developed analytic theory, due to transient wave structures in the time-dependent regime.

Bose-Einstein Condensate Dark Matter (BECDM; also known as Fuzzy Dark Matter) is motivated by fundamental physics and has recently received significant attention as a serious alternative to the established Cold Dark Matter (CDM) model. We perform cosmological simulations of BECDM gravitationally coupled to baryons and investigate structure formation at high redshifts ($z \gtrsim 5$) for a boson mass $m=2.5\cdot 10^{-22}~{\rm eV}$, exploring the dynamical effects of its wavelike nature on the cosmic web and the formation of first galaxies. Our BECDM simulations are directly compared to CDM as well as to simulations where the dynamical quantum potential is ignored and only the initial suppression of the power spectrum is considered -- a Warm Dark Matter-like ("WDM") model often used as a proxy for BECDM. Our simulations confirm that "WDM" is a good approximation to BECDM on large cosmological scales even in the presence of the baryonic feedback. Similarities also exist on small scales, with primordial star formation happening both in isolated haloes and continuously along cosmic filaments; the latter effect is not present in CDM. Global star formation and metal enrichment in these first galaxies are delayed in BECDM/"WDM" compared to the CDM case: in BECDM/"WDM" first stars form at $z\sim 13$/$13.5$ while in CDM star formation starts at $z\sim 35$. The signature of BECDM interference, not present in "WDM", is seen in the evolved dark matter power spectrum: although the small scale structure is initially suppressed, power on kpc scales is added at lower redshifts. Our simulations lay the groundwork for realistic simulations of galaxy formation in BECDM.

In hierarchical models of structure formation, the first galaxies form in low-mass dark matter potential wells, probing the behavior of dark matter on kiloparsec (kpc) scales. Even though these objects are below the detection threshold of current telescopes, future missions will open an observational window into this emergent world. In this Letter we investigate how the first galaxies are assembled in a `fuzzy' dark matter (FDM) cosmology where dark matter is an ultralight $\sim 10^{-22}$~eV boson and the primordial stars are expected to form along dense dark matter filaments. Using a first-of-its-kind cosmological hydrodynamical simulation, we explore the interplay between baryonic physics and unique wavelike features inherent to FDM. In our simulation, the dark matter filaments show coherent interference patterns on the boson de Broglie scale and develop cylindrical soliton-like cores which are unstable under gravity and collapse into kpc-scale spherical solitons. Features of the dark matter distribution are largely unaffected by the baryonic feedback. On the contrary, the distributions of gas and stars, which do form along the entire filament, exhibit central cores imprinted by dark matter -- a smoking gun signature of FDM.

Compressive isothermal turbulence is known to have a near lognormal density probability distribution function (PDF) with a width that scales with the sonic Mach number and nature of the turbulent driving (solenoidal vs compressive). However, the physical processes that mold the extreme high and low density structures in a turbulent medium can be different, with the densest structures being composed of strong shocks that evolve on shorter timescales than the low density fluid. The density PDF in a turbulent medium exhibits deviations from lognormal due to shocks, that increases with the sonic Mach number, which is often ignored in analytic models for turbulence and star formation. We develop a simple model for turbulence by treating it as a continuous Markov process, which explains both the density PDF and the transient timescales of structures as a function of density, using a framework developed in n Scannapieco & Safarzadeh (2018). Our analytic model depends on only a single parameter, the effective compressive sonic Mach number, and successfully describes the non-lognormal behavior seen in both 1D and 3D simulations of supersonic and subsonic compressive isothermal turbulence. The model quantifies the non-lognormal distribution of density structures in turbulent environments, and has application to star forming molecular clouds and star formation efficiencies.

We investigate the formation, gravitational clustering and interactions of solitons in a self-interacting, non-relativistic scalar field in an expanding universe. Rapid formation of large number of solitons is driven by attractive self-interactions of the field, whereas the slower clustering of solitons is driven by gravitational forces. Driven closer together by gravity, we see a rich plethora of dynamics in the soliton "gas" including mergers, scatterings and formation of soliton binaries. The numerical simulations are complemented by analytic calculations and estimates of (i) the relevant instability length and time scales, (ii) individual soliton profiles and their stability, (iii) number density of produced solitons, and (iv) the two point correlation function of soliton positions as evidence for gravitational clustering.

We analytically calculate the star formation efficiency and dense self-gravitating gas fraction in the presence of magneto-gravo-turbulence using the model of Burkhart (2018), which employs a piecewise lognormal and powerlaw density Probability Distribution Function (PDF). We show that the PDF transition density from lognormal to powerlaw forms is a mathematically motivated critical density for star formation and can be physically related to the density where the Jeans length is comparable to the sonic length, i.e. the post-shock critical density for collapse. When the PDF transition density is taken as the critical density, the instantaneous star formation efficiency ($\epsilon_{\rm inst}$) and depletion time ($\tau_{\rm depl}$) can be calculated from the dense self-gravitating gas fraction represented as the fraction of gas in the PDF powerlaw tail. We minimize the number of free parameters in the analytic expressions for $\epsilon_{\rm inst}$ and $\tau_{\rm depl}$ by using the PDF transition density instead of a parameterized critical density for collapse and thus provide a more direct pathway for comparison with observations. We test the analytic predictions for the transition density and self-gravitating gas fraction against AREPO moving mesh gravoturbulent simulations and find good agreement. We predict that, when gravity dominates the density distribution in the star-forming gas, the star formation efficiency should be weakly anti-correlated with the sonic Mach number while the depletion time should increase with increasing sonic Mach number. The star formation efficiency and depletion time depend primarily on the dense self-gravitating gas fraction, which in turn depends on the interplay of gravity, turbulence and stellar feedback. Our model prediction is in agreement with recent observations, such as the M51 PdBI Arcsecond Whirlpool Survey (PAWS).

Although highly successful on cosmological scales, Cold Dark Matter (CDM) models predict unobserved over-dense `cusps' in dwarf galaxies and overestimate their formation rate. We consider an ultra-light axion-like scalar boson which promises to reduce these observational discrepancies at galactic scales. The model, known as Fuzzy Dark Matter (FDM), avoids cusps, suppresses small-scale power, and delays galaxy formation via macroscopic quantum pressure. We compare the substructure and density fluctuations of galactic dark matter haloes comprised of ultra-light axions to conventional CDM results. Besides self-gravitating subhaloes, FDM includes non-virialized over-dense wavelets formed by quantum interference patterns which are an efficient source of heating to galactic discs. We find that, in the solar neighborhood, wavelet heating is sufficient to give the oldest disc stars a velocity dispersion of $\sim 30 \: \mathrm{km} \: \mathrm{s}^{-1}$ within a Hubble time if energy is not lost from the disc, the velocity dispersion increasing with stellar age as $\sigma_D \propto t^{0.4}$ in agreement with observations. Furthermore, we calculate the radius-dependent velocity dispersion and corresponding scale height caused by the heating of this dynamical substructure in both CDM and FDM with the determination that these effects will produce a flaring that terminates the Milky Way disc at $15 - 20 \: \mathrm{kpc}$. Although the source of thickened discs is not known, the heating due to perturbations caused by dark substructure cannot exceed the total disc velocity dispersion. Therefore, this work provides a lower bound on the FDM particle mass of $m_a > 0.6 \times 10^{-22} \mathrm{eV}$. Furthermore, FDM wavelets with this particle mass should be considered a viable mechanism for producing the observed disc thickening with time.

Supersonic isothermal turbulence establishes a network of transient dense shocks that sweep up material and have a density profile described by balance between ram pressure of the background fluid versus the magnetic and gas pressure gradient behind the shock. These rare, densest regions of a turbulent environment can become Jeans unstable and collapse to form pre-stellar cores. Using numerical simulations of magneto-gravo-turbulence, we describe the structural properties of dense shocks, which are the seeds of gravitational collapse, as a function of magnetic field strength. In the regime of a weak magnetic field, the collapse is isotropic. Strong magnetic field strengths lead to significant anisotropy in the shocked distribution and collapse occurs preferentially parallel to the field lines. Our work provides insight into analyzing the magnetic field topology and density structures of young protostellar collapse, which the theory presented here predicts are associated with large-scale strong shocks that persist for at least a free-fall time.

We investigate the effects of black hole mergers in star clusters on the black hole mass function. As black holes are not produced in pair-instability supernovae, it is suggested that there is a dearth of high mass stellar black holes. This dearth generates a gap in the upper end of the black hole mass function. Meanwhile, parameter fitting of X-ray binaries suggests the existence of a gap in the mass function under 5 solar masses. We show, through evolving a coagulation equation, that black hole mergers can appreciably fill the upper mass gap, and that the lower mass gap generates potentially observable features at larger mass scales. We also explore the importance of ejections in such systems and whether dynamical clusters can be formation sites of intermediate mass black hole seeds.

The Schr\"odinger-Poisson equations describe the behavior of a superfluid Bose-Einstein condensate under self-gravity with a 3D wave function. As $\hbar/m\to 0$, $m$ being the boson mass, the equations have been postulated to approximate the collisionless Vlasov-Poisson equations also known as the collisionless Boltzmann-Poisson equations. The latter describe collisionless matter with a 6D classical distribution function. We investigate the nature of this correspondence with a suite of numerical test problems in 1D, 2D, and 3D along with analytic treatments when possible. We demonstrate that, while the density field of the superfluid always shows order unity oscillations as $\hbar/m\to 0$ due to interference and the uncertainty principle, the potential field converges to the classical answer as $(\hbar/m)^{2}$. Thus, any dynamics coupled to the superfluid potential is expected to recover the classical collisionless limit as $\hbar/m\to 0$. The quantum superfluid is able to capture rich phenomena such as multiple phase-sheets, shell-crossings, and warm distributions. Additionally, the quantum pressure tensor acts as a regularizer of caustics and singularities in classical solutions. This suggests the exciting prospect of using the Schr\"odinger-Poisson equations as a low-memory method for approximating the high-dimensional evolution of the Vlasov-Poisson equations. As a particular example we consider dark matter composed of ultra-light axions, which in the classical limit ($\hbar/m\to 0$) is expected to manifest itself as collisionless cold dark matter.

In certain astrophysical systems the commonly employed ideal magnetohydrodynamics (MHD) approximation breaks down. Here, we introduce novel explicit and implicit numerical schemes of ohmic resistivity terms in the moving-mesh code AREPO. We include these non-ideal terms for two MHD techniques: the Powell 8-wave formalism and a constrained transport scheme, which evolves the cell-centred magnetic vector potential. We test our implementation against problems of increasing complexity, such as one- and two-dimensional diffusion problems, and the evolution of progressive and stationary Alfv\'en waves. On these test problems, our implementation recovers the analytic solutions to second-order accuracy. As first applications, we investigate the tearing instability in magnetized plasmas and the gravitational collapse of a rotating magnetized gas cloud. In both systems, resistivity plays a key role. In the former case, it allows for the development of the tearing instability through reconnection of the magnetic field lines. In the latter, the adopted (constant) value of ohmic resistivity has an impact on both the gas distribution around the emerging protostar and the mass loading of magnetically driven outflows. Our new non-ideal MHD implementation opens up the possibility to study magneto-hydrodynamical systems on a moving mesh beyond the ideal MHD approximation.

We present a theoretical analysis of some unexplored aspects of relaxed Bose-Einstein condensate dark matter (BECDM) haloes. This type of ultralight bosonic scalar field dark matter is a viable alternative to the standard cold dark matter (CDM) paradigm, as it makes the same large-scale predictions as CDM and potentially overcomes CDM's small-scale problems via a galaxy-scale de Broglie wavelength. We simulate BECDM halo formation through mergers, evolved under the Schr\"odinger-Poisson equations. The formed haloes consist of a soliton core supported against gravitational collapse by the quantum pressure tensor and an asymptotic $r^{-3}$ NFW-like profile. We find a fundamental relation of the core=to-halo mass with the dimensionless invariant $\Xi \equiv \lvert E \rvert/M^3/(Gm/\hbar)^2$ or $M_{\rm c}/M \simeq 2.6 \Xi^{1/3}$, linking the soliton to global halo properties. For $r \geq 3.5 \,r_{\rm c}$ core radii, we find equipartition between potential, classical kinetic, and quantum gradient energies. The haloes also exhibit a conspicuous turbulent behavior driven by the continuous reconnection of vortex lines due to wave interference. We analyse the turbulence 1D velocity power spectrum and find a $k^{-1.1}$ power-law. This suggests the vorticity in BECDM haloes is homogeneous, similar to thermally-driven counterflow BEC systems from condensed matter physics, in contrast to a $k^{-5/3}$ Kolmogorov power-law seen in mechanically-driven quantum systems. The mode where the power spectrum peaks is approximately the soliton width, implying the soliton-sized granules carry most of the turbulent energy in BECDM haloes.

We report ALMA observations of polarized dust emission from the protostellar source Ser-emb 8 at a linear resolution of 140 au. Assuming models of dust-grain alignment hold, the observed polarization pattern gives a projected view of the magnetic field structure in this source. Contrary to expectations based on models of strongly magnetized star formation, the magnetic field in Ser-emb 8 does not exhibit an hourglass morphology. Combining the new ALMA data with previous observational studies, we can connect magnetic field structure from protostellar core (~80,000 au) to disk (~100 au) scales. We compare our observations with four magnetohydrodynamic gravo-turbulence simulations made with the AREPO code that have initial conditions ranging from super-Alfv\'enic (weakly magnetized) to sub-Alfv\'enic (strongly magnetized). These simulations achieve the spatial dynamic range necessary to resolve the collapse of protostars from the parsec scale of star-forming clouds down to the ~100 au scale probed by ALMA. Only in the very strongly magnetized simulation do we see both the preservation of the field direction from cloud to disk scales and an hourglass-shaped field at < 1000 au scales. We conduct an analysis of the relative orientation of the magnetic field and the density structure in both the Ser-emb 8 ALMA observations and the synthetic observations of the four AREPO simulations. We conclude that the Ser-emb 8 data are most similar to the weakly magnetized simulations, which exhibit random alignment, in contrast to the strongly magnetized simulation, where the magnetic field plays a role in shaping the density structure in the source. In the weak-field case, it is turbulence -- not the magnetic field -- that shapes the material that forms the protostar, highlighting the dominant role that turbulence can play across many orders of magnitude in spatial scale.

Star formation in our Galaxy occurs in molecular clouds that are self-gravitating, highly turbulent, and magnetized. We study the conditions under which cloud cores inherit large-scale magnetic field morphologies and how the field is governed by cloud turbulence. We present four moving-mesh simulations of supersonic, turbulent, isothermal, self-gravitating gas with a range of magnetic mean-field strengths characterized by the Alfv\'enic Mach number $\mathcal{M}_{{\rm A}, 0}$, resolving pre-stellar core formation from parsec to a few AU scales. In our simulations with the turbulent kinetic energy density dominating over magnetic pressure ($\mathcal{M}_{{\rm A}, 0}>1$), we find that the collapse is approximately isotropic with $B\propto\rho^{2/3}$, core properties are similar regardless of initial mean-field strength, and the field direction on $100$ AU scales is uncorrelated with the mean field. However, in the case of a dominant large-scale magnetic field ($\mathcal{M}_{{\rm A}, 0}=0.35$), the collapse is anisotropic with $B\propto\rho^{1/2}$. This transition at $\mathcal{M}_{{\rm A}, 0}\sim1$ is not expected to be sharp, but clearly signifies two different paths for magnetic field evolution in star formation. Based on observations of different star forming regions, we conclude that star formation in the interstellar medium may occur in both regimes. Magnetic field correlation with the mean-field extends to smaller scales as $\mathcal{M}_{{\rm A}, 0}$ decreases, making future ALMA observations useful for constraining $\mathcal{M}_{{\rm A}, 0}$ of the interstellar medium.

We formulate a smoothed-particle hydrodynamics numerical method, traditionally used for the Euler equations for fluid dynamics in the context of astrophysical simulations, to solve the non-linear Schrodinger equation in the Madelung formulation. The probability density of the wavefunction is discretized into moving particles, whose properties are smoothed by a kernel function. The traditional fluid pressure is replaced by a quantum pressure tensor, for which a novel, robust discretization is found. We demonstrate our numerical method on a variety of numerical test problems involving the simple harmonic oscillator, Bose-Einstein condensates, collapsing singularities, and dark matter halos governed by the Gross-Pitaevskii-Poisson equation. Our method is conservative, applicable to unbounded domains, and is automatically adaptive in its resolution, making it well suited to study problems with collapsing solutions.

We revisit the integer lattice (IL) method to numerically solve the Vlasov-Poisson equations, and show that a slight variant of the method is a very easy, viable, and efficient numerical approach to study the dynamics of self-gravitating, collisionless systems. The distribution function lives in a discretized lattice phase-space, and each time-step in the simulation corresponds to a simple permutation of the lattice sites. Hence, the method is Lagrangian, conservative, and fully time-reversible. IL complements other existing methods, such as N-body/particle mesh (computationally efficient, but affected by Monte-Carlo sampling noise and two-body relaxation) and finite volume (FV) direct integration schemes (expensive, accurate but diffusive). We also present improvements to the FV scheme, using a moving mesh approach inspired by IL, to reduce numerical diffusion and the time-step criterion. Being a direct integration scheme like FV, IL is memory limited (memory requirement for a full 3D problem scales as N^6, where N is the resolution per linear phase-space dimension). However, we describe a new technique for achieving N^4 scaling. The method offers promise for investigating the full 6D phase-space of collisionless systems of stars and dark matter.

We present a constrained transport (CT) algorithm for solving the 3D ideal magnetohydrodynamic (MHD) equations on a moving mesh, which maintains the divergence-free condition on the magnetic field to machine-precision. Our CT scheme uses an unstructured representation of the magnetic vector potential, making the numerical method simple and computationally efficient. The scheme is implemented in the moving mesh code Arepo. We demonstrate the performance of the approach with simulations of driven MHD turbulence, a magnetized disc galaxy, and a cosmological volume with primordial magnetic field. We compare the outcomes of these experiments to those obtained with a previously implemented Powell divergence-cleaning scheme. While CT and the Powell technique yield similar results in idealized test problems, some differences are seen in situations more representative of astrophysical flows. In the turbulence simulations, the Powell cleaning scheme artificially grows the mean magnetic field, while CT maintains this conserved quantity of ideal MHD. In the disc simulation, CT gives slower magnetic field growth rate and saturates to equipartition between the turbulent kinetic energy and magnetic energy, whereas Powell cleaning produces a dynamically dominant magnetic field. Such difference has been observed in adaptive-mesh refinement codes with CT and smoothed-particle hydrodynamics codes with divergence-cleaning. In the cosmological simulation, both approaches give similar magnetic amplification, but Powell exhibits more cell-level noise. CT methods in general are more accurate than divergence-cleaning techniques, and, when coupled to a moving mesh can exploit the advantages of automatic spatial/temporal adaptivity and reduced advection errors, allowing for improved astrophysical MHD simulations.

Accurate numerical solutions of the equations of hydrodynamics play an ever more important role in many fields of astrophysics. In this work, we reinvestigate the accuracy of the moving-mesh code \textsc{Arepo} and show how its convergence order can be improved for general problems. In particular, we clarify that for certain problems \textsc{Arepo} only reaches first-order convergence for its original formulation. This can be rectified by simple modifications we propose to the time integration scheme and the spatial gradient estimates of the code, both improving the accuracy of the code. We demonstrate that the new implementation is indeed second-order accurate under the $L^1$ norm, and in particular substantially improves conservation of angular momentum. Interestingly, whereas these improvements can significantly change the results of smooth test problems, we also find that cosmological simulations of galaxy formation are unaffected, demonstrating that the numerical errors eliminated by the new formulation do not impact these simulations. In contrast, simulations of binary stars followed over a large number of orbital times are strongly affected, as here it is particularly crucial to avoid a long-term build up of errors in angular momentum conservation.

We perform uniformly sampled large-scale cosmological simulations including magnetic fields with the moving mesh code AREPO. We run two sets of MHD simulations: one including adiabatic gas physics only; the other featuring the fiducial feedback model of the Illustris simulation. In the adiabatic case, the magnetic field amplification follows the $B \propto \rho^{2/3}$ scaling derived from `flux-freezing' arguments, with the seed field strength providing an overall normalization factor. At high baryon overdensities the amplification is enhanced by shear flows and turbulence. Feedback physics and the inclusion of radiative cooling change this picture dramatically. In haloes, gas collapses to much larger densities and the magnetic field is amplified strongly and to the same maximum intensity irrespective of the initial seed field of which any memory is lost. At lower densities a dependence on the seed field strength and orientation, which in principle can be used to constrain models of cosmic magnetogenesis, is still present. Inside the most massive haloes magnetic fields reach values of $\sim 10-100\,\,{\rm \mu G}$, in agreement with galaxy cluster observations. The topology of the field is tangled and gives rise to rotation measure signals in reasonable agreement with the observations. However, the rotation measure signal declines too rapidly towards larger radii as compared to observational data.

A method for improving the accuracy of hydrodynamical codes that use a moving Voronoi mesh is described. Our scheme is based on a new regularization scheme that constrains the mesh to be centroidal to high precision while still allowing the cells to move approximately with the local fluid velocity, thereby retaining the quasi-Lagrangian nature of the approach. Our regularization technique significantly reduces mesh noise that is attributed to changes in mesh topology and deviations from mesh regularity. We demonstrate the advantages of our method on various test problems, and note in particular improvements obtained in handling shear instabilities, mixing, and in angular momentum conservation. Calculations of adiabatic jets in which shear excites Kelvin Helmholtz instability show reduction of mesh noise and entropy generation. In contrast, simulations of the collapse and formation of an isolated disc galaxy are nearly unaffected, showing that numerical errors due to the choice of regularization do not impact the outcome in this case.

General relativistic magnetohydrodynamic (GRMHD) simulations of accreting black holes in the radiatively inefficient regime show that systems with sufficient magnetic poloidal flux become magnetically arrested disc (MAD) systems, with a well-defined relationship between the magnetic flux and the mass accretion rate. Recently, Zamaninasab (2014) report that the jet magnetic flux and accretion disc luminosity are tightly correlated over 7 orders of magnitude for a sample of 76 radio-loud active galaxies, concluding that the data are explained by the MAD mode of accretion. Their analysis assumes radiatively efficient accretion, and their sample consists primarily of radiatively efficient sources, while GRMHD simulations of MAD thus far have been carried out in the radiatively inefficient regime. We propose a model to interpret MAD systems in the context of multiple accretion regimes, and apply it to the sample in Zamaninasab (2014), along with additional radiatively inefficient sources from archival data. We show that most of the radiatively inefficient radio-loud galaxies are consistent with being MAD systems. Assuming the MAD relationship found in radiatively inefficient simulations holds at other accretion regimes, a significant fraction of our sample can be candidates for MAD systems. Future GRMHD simulations have yet to verify the validity of this assumption.

Magnetic fields play an important role in many astrophysical systems and a detailed understanding of their impact on the gas dynamics requires robust numerical simulations. Here we present a new method to evolve the ideal magnetohydrodynamic (MHD) equations on unstructured static and moving meshes that preserves the magnetic field divergence-free constraint to machine precision. The method overcomes the major problems of using a cleaning scheme on the magnetic fields instead, which is non-conservative, not fully Galilean invariant, does not eliminate divergence errors completely, and may produce incorrect jumps across shocks. Our new method is a generalization of the constrained transport (CT) algorithm used to enforce the $\nabla\cdot \mathbf{B}=0$ condition on fixed Cartesian grids. Preserving $\nabla\cdot \mathbf{B}=0$ at the discretized level is necessary to maintain the orthogonality between the Lorentz force and $\mathbf{B}$. The possibility of performing CT on a moving mesh provides several advantages over static mesh methods due to the quasi-Lagrangian nature of the former (i.e., the mesh generating points move with the flow), such as making the simulation automatically adaptive and significantly reducing advection errors. Our method preserves magnetic fields and fluid quantities in pure advection exactly.

Double-lobed radio galaxies a few 100s of kpc in extent, like Cygnus A, are common at redshifts of 1 to 2, arising from some 10 per cent of the most powerful Active Galactic Nuclei (AGN). At higher redshifts they are rare, with none larger than a few 10s of kpc known above redshift z~4. Recent studies of the redshift evolution of powerful-jetted objects indicate that they may constitute a larger fraction of the AGN population above redshift 2 than appears from a simple consideration of detected GHz radio sources. The radio band is misleading as the dramatic (1+z)^4 boost in the energy density of the Cosmic Microwave Background (CMB) causes inverse Compton scattering to dominate the energy losses of relativistic electrons in the extended lobes produced by jets, making them strong X-ray, rather than radio, sources. Here we investigate limits to X-ray lobes around two distant quasars, ULASJ112001.48+064124.3 at z = 7.1 and SDSSJ1030+0524 at z=6.3, and find that powerful jets could be operating yet be currently undetectable. Jets may be instrumental in the rapid build-up of billion Msun black hole at a rate that violates the Eddington limit.

A discontinuous Galerkin (DG) method suitable for large-scale astrophysical simulations on Cartesian meshes as well as arbitrary static and moving Voronoi meshes is presented. Most major astrophysical fluid dynamics codes use a finite volume (FV) approach. We demonstrate that the DG technique offers distinct advantages over FV formulations on both static and moving meshes. The DG method is also easily generalized to higher than second-order accuracy without requiring the use of extended stencils to estimate derivatives (thereby making the scheme highly parallelizable). We implement the technique in the AREPO code for solving the fluid and the magnetohydrodynamic (MHD) equations. By examining various test problems, we show that our new formulation provides improved accuracy over FV approaches of the same order, and reduces post-shock oscillations and artificial diffusion of angular momentum. In addition, the DG method makes it possible to represent magnetic fields in a locally divergence-free way, improving the stability of MHD simulations and moderating global divergence errors, and is a viable alternative for solving the MHD equations on meshes where Constrained-Transport (CT) cannot be applied. We find that the DG procedure on a moving mesh is more sensitive to the choice of slope limiter than is its FV method counterpart. Therefore, future work to improve the performance of the DG scheme even further will likely involve the design of optimal slope limiters. As presently constructed, our technique offers the potential of improved accuracy in astrophysical simulations using the moving mesh AREPO code as well as those employing adaptive mesh refinement (AMR).

We develop a simple evolutionary scenario for the growth of supermassive black holes (BHs), assuming growth due to accretion only, to learn about the evolution of the BH mass function from $z=3$ to 0 and from it calculate the energy budgets of different modes of feedback. We tune the parameters of the model by matching the derived X-ray luminosity function (XLF) with the observed XLF of active galactic nuclei. We then calculate the amount of comoving kinetic and bolometric feedback as a function of redshift, derive a kinetic luminosity function and estimate the amount of kinetic feedback and $PdV$ work done by classical double Fanaroff-Riley II (FR II) radio sources. We also derive the radio luminosity function for FR IIs from our synthesized population and set constraints on jet duty cycles. Around 1/6 of the jet power from FR II sources goes into $PdV$ work done in the expanding lobes during the time the jet is on. Anti hierarchical growth of BHs is seen in our model due to addition of an amount of mass being accreted on to all BHs independent of the BH mass. The contribution to the total kinetic feedback by active galaxies in a low accretion, kinetically efficient mode is found to be the most significant at $z<1.5$. FR II feedback is found to be a significant mode of feedback above redshifts $z\sim 1.5$, which has not been highlighted by previous studies.

We construct Tully-Fisher relationships (TFRs) in the $u$, $g$, $r$, $i$ and $z$ bands and stellar mass TFRs (smTFRs) for a sample of $25,698$ late spiral type galaxies (with $0.045<z<0.085$) from the Sloan Digital Sky Survey (SDSS) and study the effects of environment on the relation. We use SDSS-measured Balmer emission line widths, $v_{\rm FWHM}$, as a proxy for disc circular velocity, $v_{\rm circ}$. A priori it is not clear whether we can construct accurate TFRs given the small $3"$ diameter of the fibres used for SDSS spectroscopic measurements. However, we show by modelling the H$\alpha$ emission profile as observed through a $3"$ aperture that for galaxies at appropriate redshifts ($z>0.045$) the fibres sample enough of the disc to obtain a linear relationship between $v_{\rm FWHM}$ and $v_{\rm circ}$, allowing us to obtain a TFR and to investigate dependence on other variables. We also develop a methodology for distinguishing between astrophysical and sample bias in the fibre TFR trends. We observe the well-known steepening of the TFR in redder bands in our sample. We divide the sample of galaxies into four equal groups using projected neighbour density ($\Sigma$) quartiles and find no significant dependence on environment, extending previous work to a wider range of environments and a much larger sample. Having demonstrated that we can construct SDSS-based TFRs is very useful for future applications because of the large sample size available.

We present new GMRT observations of HDF 130, an inverse-Compton (IC) ghost of a giant radio source that is no longer being powered by jets. We compare the properties of HDF 130 with the new and important constraint of the upper limit of the radio flux density at 240 MHz to an analytic model. We learn what values of physical parameters in the model for the dynamics and evolution of the radio luminosity and X-ray luminosity (due to IC scattering of the cosmic microwave background (CMB)) of a Fanaroff-Riley II (FR II) source are able to describe a source with features (lobe length, axial ratio, X-ray luminosity, photon index and upper limit of radio luminosity) similar to the observations. HDF 130 is found to agree with the interpretation that it is an IC ghost of a powerful double-lobed radio source, and we are observing it at least a few Myr after jet activity (which lasted 5--100 Myr) has ceased. The minimum Lorentz factor of injected particles into the lobes from the hotspot is preferred to be $\gamma\sim10^3$ for the model to describe the observed quantities well, assuming that the magnetic energy density, electron energy density, and lobe pressure at time of injection into the lobe are linked by constant factors according to a minimum energy argument, so that the minimum Lorentz factor is constrained by the lobe pressure. We also apply the model to match the features of 6C 0905+3955, a classical double FR II galaxy thought to have a low-energy cutoff of $\gamma\sim10^4$ in the hotspot due to a lack of hotspot inverse-Compton X-ray emission. The models suggest that the low-energy cutoff in the hotspots of 6C 0905+3955 is $\gamma\gtrsim 10^3$, just slightly above the particles required for X-ray emission.

We analyze the \textit{Chandra} High Energy Transmission Grating Spectrometer (HETGS) observation of the Seyfert $2$ AGN IRAS\,18325-5926. We detect a $v=-360^{+41}_{-66}$~km~s$^{-1}$ blueshifted ionized absorber in the X-ray spectrum, with photoionization parameter $\log\xi=2.0^{+0.1}_{-0.1}$ and hydrogen column density $N_{\rm H}=1.55^{+0.75}_{-0.38}\times10^{21}$~cm$^{-2}$. The absorber may be a photoionized wind originating in the obscuring torus/global covering around the black hole or outer edge of the accretion disk. The estimated mass outflow rate suggests that the supermassive black hole in IRAS\,18325-5926 may significantly affect the large-scale environment of the host galaxy, unless the solid angle subtended by the outflow or the gas filling factor is small. A second warm absorber may be needed to explain the absorption features in the vicinity of the iron K edge, although insufficient counts in the data beyond $7.0$~keV make it difficult to assess the nature of the second absorber. Most plausible is a high ionization ($\log\xi\sim 2.3$ to $2.6$), high column density ($N_{\rm H}\sim10^{23}$~cm$^{-2}$) absorber with $v\sim -3000$ to $0$~km~s$^{-1}$, although these parameters are not well constrained. We also examine the broad Fe K emission line in the spectrum, which is likely due to Fe XXV in a highly ionized accretion disk inclined at $25^\circ$, consistent with the \textit{XMM-Newton} EPIC observations of the emission line. Because we are able to view both the obscuring gas and the accretion disk of IRAS\,18325-5926, the surrounding gas of IRAS\,18325-5926 may be patchy or we are viewing the system at an angle just grazing the obscuring torus.

In this study we predict the total distributions of powerful (FR II) active double-lobed radio galaxies and ghost sources, and their observable distribution in the X-ray sky. We develop an analytic model for the evolution of the lobe emission at radio and X-ray energies. During jet activity, a double radio source emits synchrotron radiation in the radio and X-ray emission due to inverse-Compton (IC) upscattering by gamma~10^3 electrons of the cosmic microwave background. After the jets switch off, the radio luminosity (due to higher gamma electrons) falls faster than the X-ray luminosity and for some time the source appears as an IC ghost of a radio galaxy before becoming completely undetectable in the X-ray. With our model, for one set of typical parameters, we predict radio lobes occupy a volume fraction of the universe of 0.01, 0.03, 0.3 at z=2 (during the quasar era) of the filamentary structures in which they are situated, for typical jet lifetimes 5*10^7 yr, 10^8 yr, 5*10^8 yr; however since the inferred abundance of sources depends on how quickly they fall below the radio flux limit the volume filling factor is found to be a strong function of radio galaxy properties such as energy index and minimum gamma factor of injected particles, the latter not well constrained by observations. We test the predicted number density of sources against the Chandra X-ray Deep Field North survey and also find the contribution to the unresolved cosmic X-ray background by the lobes of radio galaxies. 10-30 per cent of observable double-lobed structures in the X-ray are predicted to be IC ghosts. The derived X-ray luminosity function of our synthetic population shows that double-lobed sources have higher space densities than X-ray clusters at redshifts z>2 and X-ray luminosities above 10^44 erg s^-1.