Synchrotron emission from astrophysical nonthermal sources usually assumes that the emitting particles are isotropic. By means of large-scale two- and three-dimensional particle-in-cell simulations, we demonstrate that the dissipation of magnetically-dominated ($\sigma_0\gg1$) turbulence in pair plasmas leads to strongly anisotropic particle distributions. At Lorentz factors $\sim \sigma_0 \gamma_{th0}$ (here, $\gamma_{th0}$ is the initial Lorentz factor), the particle velocity is preferentially aligned with the local magnetic field; instead, the highest energy particles are roughly isotropic. This energy-dependent anisotropy leads to a synchrotron spectral flux $\nu F_\nu\propto \nu^s$ that is much harder than for isotropic particles. Remarkably, for $\sigma_0\gg1$ we find that the solid-angle-averaged spectral slope in the slow cooling regime is $s\sim 0.5-0.7$ for a wide range of turbulence fluctuations, $0.25\lesssim \delta B_{\rm rms0}^2/B_0^2\lesssim 4$, despite significant variations in the power-law energy spectrum of nonthermal particles. This is because weaker turbulence levels imprint a stronger degree of anisotropy, thereby counteracting the effect of the steeper particle spectrum. The synchrotron spectral slope may be even harder, $s\gtrsim 0.7$, if the observer is in the plane perpendicular to the mean magnetic field. Our results are independent of domain size and dimensionality. Our findings may help explaining the origin of hard synchrotron spectra of astrophysical nonthermal sources, most notably the radio spectrum of Pulsar Wind Nebulae.

The dissipation of turbulent magnetic fields is an appealing scenario to explain the origin of non-thermal particles in high-energy astrophysical sources. However, it has been suggested that the particle distribution may effectively thermalise when the radiative (synchrotron and/or Inverse Compton) losses are severe. Inspired by recent PIC simulations of relativistic turbulence, which show that electrons are impulsively heated in intermittent current sheets by a strong electric field aligned with the local magnetic field, we instead argue that in plasmas where the particle number density is dominated by the pairs (electron-positron and electron-positron-ion plasmas): (i) as an effect of fast cooling and of different injection times, the electron energy distribution is ${\rm d}n_e/{\rm d}\gamma\propto\gamma^{-2}$ for $\gamma\lesssim\gamma_{\rm heat}$ (the Lorentz factor $\gamma_{\rm heat}$ being close to the equipartition value), while the distribution steepens at higher energies; (ii) since the time scales for the turbulent fields to decay and for the photons to escape are of the same order, the magnetic and the radiation energy densities in the dissipation region are comparable; (iii) if the mass energy of the plasma is dominated by the ion component, the pairs with a Lorentz factor smaller than a critical one (of the order of the proton-to-electron mass ratio) become isotropic, while the pitch angle remains small otherwise. The outlined scenario is consistent with the typical conditions required to reproduce the Spectral Energy Distribution of blazars, and allows one to estimate the magnetisation of the emission site. Finally, we show that turbulence within the Crab Nebula may power the observed gamma-ray flares if the pulsar wind is nearly charge-separated at high latitudes.

We investigate the impact on the properties of high-energy emitting BL Lac objects of a population of electrons with an anisotropic momentum distribution. We adopt a simple phenomenological description of the anisotropy, in which the most energetic electrons have a small pitch angle and the least energetic electrons are isotropic, as proposed by Sobacchi \& Lyubarsky (2019). We explore (i) a simple model that assumes a phenomenological shape for the electron energy distribution, and (ii) a self-consistent scheme in which the electrons follow a distribution which is the result of the balance between injection and radiative losses (we include the effects of the anisotropy on the synchrotron cooling rate). Considering the BL Lac object Mkn 421 as representative of the entire class, we show that in both cases the emission can be satisfactorily reproduced under equipartition between the magnetic field and the relativistic electrons. This is in better agreement with the idea that jets are launched as Poynting dominated flows with respect to the standard isotropic scenario, which requires both a low magnetization and a low radiative efficiency to reproduce the observed emission. The hard spectrum predicted for the inverse Compton continuum at TeV energies could be used as a potential test of the anisotropic model.

Whether High Velocity Clouds (HVCs) can form by condensation of the hot ($T \sim 10^6 \, {\rm K}$) Galactic corona as a consequence of thermal instabilities has been controversial. Here we re-examine this problem and we suggest that rotation of the corona might be a missing key ingredient. We do this by studying the evolution of the models of rotating galactic coronae presented in Sormani et al. (2018) under the presence of cooling and thermal conduction. We combine a linear stability analysis with the results of local and global hydrodynamical simulations. We find that condensations are likely to occur in regions where the corona has substantial rotational support. Under reasonably general assumptions on the rotation profile of the corona, the locations where condensations are expected are in remarkable agreement with the observed location of the major non-magellanic HVCs complexes in our Galaxy (namely, at distances $< 15 \, \rm kpc$ from the Sun and within $30^\circ$ from the disc plane). We conclude that HVCs can form by thermal instabilities provided that (i) the corona is rotating substantially in the inner ($R < 50 \, \rm kpc$) parts, as suggested by current observational data and predicted by cosmological simulations of galaxy formation; (ii) close to the disc the corona is well-represented by a nearly-equilibrium stratified rotating structure (as opposed to a fast cooling flow). Our results also suggest that a better understanding of the disc-halo interface, including supernova feedback, is critical to understand the origin of HVCs.

Observations show that (i) multiple gas phases can coexist in the atmospheres of galaxies and clusters; (ii) these atmospheres may be significantly rotating in the inner parts, with typical velocities that approach or even exceed the local sound speed. The thermal instability is a natural candidate to explain the formation of cold structures via condensation of a hotter gas phase. Here we systematically study the effect of rotation on the thermal stability of stratified plane-parallel atmospheres, using both analytical arguments and numerical simulations. We find that the formation of cold structures starting from small isobaric perturbations is enhanced in the regions where the rotation of the system is dynamically important (i.e. when the rotational velocity becomes comparable to the sound speed). In particular, the threshold value of the ratio between the cooling and dynamical time $t_{\rm cool}/t_{\rm dyn}$ below which condensations can form is increased by a factor up to $\sim 10$ in the presence of significant rotation. We briefly discuss the implications of our results for galaxies and clusters.

Theoretical modelling and observations of AGN jets suggest that the non-thermal electrons emitting the observed radiation should (i) carry an amount of energy comparable to the magnetic fields ($U_e\sim U_B$), which is likely the case if the magnetic fields play a dynamically important role in the jet's acceleration process; (ii) cool efficiently in a dynamical time ($t_{\rm cool}\lesssim t_{\rm dyn}$), which is suggested by the fact that a large fraction of the jet's kinetic energy is promptly converted into radiation. These expectations are at odds with the results of the simplest one-zone Self-Synchro-Compton (SSC) model for the Spectral Energy Distribution (SED) of BL Lacs. Indeed, the model predicts $U_e\gg U_B$ and $t_{\rm cool}\gg t_{\rm dyn}$ for most of the objects. Here we closely investigate one of the key assumptions of this model, namely that the momentum distribution of the non-thermal electrons is isotropic. We find that this assumption may be an oversimplification. If the magnetic energy is dissipated via a turbulent MHD cascade, the highest energy electrons may instead retain a small pitch angle. Since the synchrotron emissivity is suppressed when the pitch angle is small, this effect may importantly affect the modelling of the SED. As an illustrative example, we present an anisotropic model for the electron momentum distribution such that $U_e\sim U_B$ and $t_{\rm cool}\lesssim t_{\rm dyn}$ at the same time. Our model manages to simultaneously solve the two problems with one only more free parameter with respect to the usual isotropic one-zone SSC model.

With the recent discovery of a dozen dusty star-forming galaxies and around 30 quasars at z>5 that are hyper-luminous in the infrared ($\mu$$L_{\rm IR}>10^{13}$ L$_{\odot}$, where $\mu$ is a lensing magnification factor), the possibility has opened up for SPICA, the proposed ESA M5 mid-/far-infrared mission, to extend its spectroscopic studies toward the epoch of reionization and beyond. In this paper, we examine the feasibility and scientific potential of such observations with SPICA's far-infrared spectrometer SAFARI, which will probe a spectral range (35-230 $\mu$m) that will be unexplored by ALMA and JWST. Our simulations show that SAFARI is capable of delivering good-quality spectra for hyper-luminous infrared galaxies (HyLIRGs) at z=5-10, allowing us to sample spectral features in the rest-frame mid-infrared and to investigate a host of key scientific issues, such as the relative importance of star formation versus AGN, the hardness of the radiation field, the level of chemical enrichment, and the properties of the molecular gas. From a broader perspective, SAFARI offers the potential to open up a new frontier in the study of the early Universe, providing access to uniquely powerful spectral features for probing first-generation objects, such as the key cooling lines of low-metallicity or metal-free forming galaxies (fine-structure and H2 lines) and emission features of solid compounds freshly synthesized by Population III supernovae. Ultimately, SAFARI's ability to explore the high-redshift Universe will be determined by the availability of sufficiently bright targets (whether intrinsically luminous or gravitationally lensed). With its launch expected around 2030, SPICA is ideally positioned to take full advantage of upcoming wide-field surveys such as LSST, SKA, Euclid, and WFIRST, which are likely to provide extraordinary targets for SAFARI.

Fitting equilibrium dynamical models to observational data is an essential step in understanding the structure of the gaseous hot haloes that surround our own and other galaxies. However, the two main categories of models that are used in the literature are poorly suited for this task: (i) simple barotropic models are analytic and can therefore be adjusted to match the observations, but are clearly unrealistic because the rotational velocity $v_\phi(R,z)$ does not depend on the distance $z$ from the galactic plane, while (ii) models obtained as a result of cosmological galaxy formation simulations are more realistic, but are impractical to fit to observations due to high computational cost. Here we bridge this gap by presenting a general method to construct axisymmetric baroclinic equilibrium models of rotating galactic coronae in arbitrary external potentials. We consider in particular a family of models whose equipressure surfaces in the $(R,z)$ plane are ellipses of varying axis ratio. These models are defined by two one-dimensional functions, the axial ratio of pressure $q_{\rm axis}(z)$ and the value of the pressure $P_{\rm axis}(z)$ along the galaxy's symmetry axis. These models can have a rotation speed $v_\phi(R,z)$ that realistically decreases as one moves away from the galactic plane, and can reproduce the angular momentum distribution found in cosmological simulations. The models are computationally cheap to construct and can thus be used in fitting algorithms. We provide a python code that given $q_{\rm axis}(z)$, $P_{\rm axis}(z)$ and $\Phi(R,z)$ returns $\rho(R,z)$, $T(R,z)$, $P(R,z)$, $v_\phi(R,z)$. We show a few examples of these models using the Milky Way as a case study.

In the paradigm of magnetic acceleration of relativistic outflows, a crucial point is identifying a viable mechanism to convert the Poynting flux into the kinetic energy of the plasma, and eventually into the observed radiation. Since the plasma is hardly accelerated beyond equipartition, MHD instabilities are often invoked to explain the dissipation of the magnetic energy. Motivated by the fast variability that is shown by the gamma-ray flares of both AGN and PWNe, different authors have proposed the Poynting flux to be dissipated in a region where the flow is converging. Here we perform a linear stability analysis of ultra-relativistic, highly magnetised outflows with such a recollimation nozzle, showing that MHD instabilities are indeed induced by the convergence of the flow. The amplitude of the perturbations increases while recollimation gets stronger, and eventually diverges when the flow is focused to a single point. Hence, depending on the geometry of the outflow, instabilities excited while the flow is converging may play an important role to dissipate the magnetic energy of the plasma.

We develop a dynamical theory for the origin of nuclear rings in barred galaxies. In analogy with the standard theory of accretion discs, our theory is based on shear viscous forces among nested annuli of gas. However, the fact that gas follows non circular orbits in an external barred potential has profound consequences: it creates a region of reverse shear in which it is energetically favourable to form a stable ring which does not spread despite dissipation. Our theory allows us to approximately predict the size of the ring given the underlying gravitational potential. The size of the ring is loosely related to the location of the Inner Lindblad Resonance in the epicyclic approximation, but the predicted location is more accurate and is also valid for strongly barred potentials. By comparing analytical predictions with the results of hydrodynamical simulations, we find that our theory provides a viable mechanism for ring formation if the effective sound speed of the gas is low ($\cs\lesssim1\kms$), but that nuclear spirals/shocks created by pressure destroy the ring when the sound speed is high ($\cs\simeq10\kms$). We conclude that whether this mechanism for ring formation is relevant for real galaxies ultimately depends on the effective equation of state of the ISM. Promising confirmation comes from simulations in which the ISM is modelled using state-of-the-art cooling functions coupled to live chemical networks, but more tests are needed regarding the role of turbulence driven by stellar feedback. If the mechanism is relevant in real galaxies, it could provide a powerful tool to constrain the gravitational potential, in particular the bar pattern speed.

IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Large Millimeter Array (ALMA) at longer wavelengths. SPICA, with its 2.5-m telescope actively-cooled to below 8K, will obtain the first spectroscopic determination, in the mid-IR rest-frame, of both the star-formation rate and black hole accretion rate histories of galaxies, reaching lookback times of 12 Gyr, for large statistically significant samples. Densities, temperatures, radiation fields and gas-phase metallicities will be measured in dust-obscured galaxies and active galactic nuclei (AGN), sampling a large range in mass and luminosity, from faint local dwarf galaxies to luminous quasars in the distant Universe. AGN and starburst feedback and feeding mechanisms in distant galaxies will be uncovered through detailed measurements of molecular and atomic line profiles. SPICA's large-area deep spectrophotometric surveys will provide mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, out to redshifts of z~6. Furthermore, SPICA spectroscopy will uncover the most luminous galaxies in the first few hundred million years of the Universe, through their characteristic dust and molecular hydrogen features.

In the paradigm of magnetic launching of astrophysical jets, instabilities in the MHD flow are a good candidate to convert the Poynting flux into the kinetic energy of the plasma. If the magnetised plasma fills the almost entire space, the jet is unstable to helical perturbations of its body. However, the growth rate of these modes is suppressed when the poloidal component of the magnetic field has a vanishing gradient, which may be the actual case for a realistic configuration. Here we show that, if the magnetised plasma is confined into a limited region by the pressure of some external medium, the velocity shear at the contact surface excites unstable modes which can affect a significant fraction of the jet's body. We find that when the Lorentz factor of the jet is $\Gamma\sim10$ ($\Gamma\sim 100$), these perturbations typically develop after propagating along the jet for tens (hundreds) of jet's radii. Surface modes may therefore play an important role in converting the energy of the jet from the Poynting flux to the kinetic energy of the plasma, particularly in AGN. The scaling of the dispersion relation with (i) the angular velocity of the field lines and (ii) the sound speed in the confining gas is discussed.

We study the CO line luminosity ($L_{\rm CO}$), the shape of the CO Spectral Line Energy Distribution (SLED), and the value of the CO-to-$\rm H_2$ conversion factor in galaxies in the Epoch of Reionization (EoR). To this aim, we construct a model that simultaneously takes into account the radiative transfer and the clumpy structure of giant molecular clouds (GMCs) where the CO lines are excited. We then use it to post-process state-of-the-art zoomed, high resolution ($30\, \rm{pc}$), cosmological simulation of a main-sequence ($M_{*}\approx10^{10}\, \rm{M_{\odot}}$, $SFR\approx 100\,\rm{M_{\odot}\, yr^{-1}}$) galaxy, "Alth{\ae}a", at $z\approx6$. We find that the CO emission traces the inner molecular disk ($r\approx 0.5 \,\rm{kpc}$) of Alth{\ae}a with the peak of the CO surface brightness co-located with that of the [CII] 158$\rm \mu m$ emission. Its $L_{\rm CO(1-0)}=10^{4.85}\, \rm{L_{\odot}}$ is comparable to that observed in local galaxies with similar stellar mass. The high ($\Sigma_{gas} \approx 220\, \rm M_{\odot}\, pc^{-2}$) gas surface density in Alth{\ae}a, its large Mach number (\mach$\approx 30$), and the warm kinetic temperature ($T_{k}\approx 45 \, \rm K$) of GMCs yield a CO SLED peaked at the CO(7-6) transition, i.e. at relatively high-$J$, and a CO-to-$\rm H_2$ conversion factor $\alpha_{\rm CO}\approx 1.5 \, \rm M_{\odot} \rm (K\, km\, s^{-1}\, pc^2)^{-1} $ lower than that of the Milky Way. The ALMA observing time required to detect (resolve) at 5$\sigma$ the CO(7-6) line from galaxies similar to Alth{\ae}a is $\approx13$ h ($\approx 38$ h).

Long-duration, spectrally-soft Gamma-Ray Bursts (GRBs) are associated with Type Ic Core Collapse (CC) Supernovae (SNe), and thus arise from the death of massive stars. In the collapsar model, the jet launched by the central engine must bore its way out of the progenitor star before it can produce a GRB. Most of these jets do not break out, and are instead "choked" inside the star, as the central-engine activity time, $t_{\rm e}$, is not long enough. Modelling the long-soft GRB duration distribution assuming a power-law distribution for their central-engine activity times, $\propto t_{\rm e}^{-\alpha}$ for $t_{\rm e}>t_{\rm b}$, we find a steep distribution ($\alpha\sim4$) and a typical GRB jet breakout time of $t_{\rm b}\sim 60\text{ s}$ in the star's frame. The latter suggests the presence of a low-density, extended envelope surrounding the progenitor star, similar to that previously inferred for low-luminosity GRBs. Extrapolating the range of validity of this power law below what is directly observable, to $t_{\rm e}<t_{\rm b}$, by only a factor of $\sim$4-5 produces enough events to account for all Type Ib/c SNe. Such extrapolation is necessary to avoid fine-tuning the distribution of central engine activity times with the breakout time, which are presumably unrelated. We speculate that central engines launching relativistic jets may operate in all Type Ib/c SNe. In this case, the existence of a common central engine would imply that (i) the jet may significantly contribute to the energy of the SN; (ii) various observational signatures, like the asphericity of the explosion, could be directly related to jet's interaction with the star.

We study the dynamical stability of stationary galactic spiral shocks. The steady-state equilibrium flow contains a shock of the type derived by Roberts in the tightly wound approximation. We find that boundary conditions are critical in determining whether the solutions are stable or not. Shocks are unstable if periodic boundary conditions are imposed. For intermediate strengths of the spiral potential, the instability disappears if boundary conditions are imposed such that the upstream flow is left unperturbed as in the classic analysis of D'yakov and Kontorovich. This reconciles apparently contradictory findings of previous authors regarding the stability of spiral shocks. This also shows that the instability is distinct from the Kelvin-Helmholtz instability, confirming the findings of Kim et al. We suggest that instability is a general characteristics of periodic shocks, regardless of the presence of shear, and provide a physical picture as to why this is the case. For strong spiral potentials, high post-shock shear makes the system unstable also to parasitic Kelvin-Helmholtz instability regardless of the boundary conditions. Our analysis is performed in the context of a simplified problem that, while preserving all the important characteristics of the original problem, strips it from unnecessary complications, and assumes that the gas is isothermal, non self-gravitating, non-magnetised.

Core-collapse supernovae (SNe) of Types Ib and Ic arise from hydrogen-stripped stars, while the latter are also stripped of their helium. Both SN types have a similar temporal evolution, suggesting broadly similar progenitors. However, while some Ic SNe are associated with relativistic jets of gamma-ray bursts (GRBs), no GRB has yet been found in association with Type Ib SNe. Here we find that, even if GRB-like central engines operate in both SNe Ib and Ic, different properties of their envelopes may accommodate this potential tension. In particular, we focus on the case of a low-mass, extended envelope surrounding the progenitor star (as produced, for example, by strong mass losses prior to explosion). If the envelopes of Type Ib SNe are sufficiently massive ($M_{\rm ext}\sim(0.3-1)M_\odot$) and extended ($R_{\rm ext}\sim 10^{13}\text{ cm}$), we show that (i) GRB-like jets cannot break out of the star; (ii) the SN light curve is compatible with current observations. Different envelope properties of Type Ib SNe with respect to (at least some) Type Ic SNe may be connected to the presence of a helium layer surrounding the star.

In the paradigm of magnetic acceleration of relativistic jets, one of the key points is identifying a viable mechanism to convert the Poynting flux into the kinetic energy of the plasma beyond equipartition. A promising candidate is the kink instability, which deforms the body of the jet through helical perturbations. Since the detailed structure of real jets is unknown, we explore a large family of cylindrical, force-free equilibria to get robust conclusions. We find that the growth rate of the instability depends primarily on two parameters: (i) the gradient of the poloidal magnetic field; (ii) the Lorentz factor of the perturbation, which is closely related to the velocity of the plasma. We provide a simple fitting formula for the growth rate of the instability. As a tentative application, we use our results to interpret the dynamics of the jet in the nearby active galaxy M87. We show that the kink instability becomes non-linear at a distance from the central black hole comparable to where the jet stops accelerating. Hence (at least for this object), the kink instability of the jet is a good candidate to drive the transition from a Poynting-dominated to a kinetic-energy-dominated flow.

We describe a scenario to explain blazar periodicities with timescales of $\sim$ few years. The scenario is based on a binary super-massive black hole (SMBH) system in which one of the two SMBH carries a jet. We discuss the various mechanisms that can cause the jet to precess and produce corkscrew patterns through space with a scale of $\sim$ few pc. It turns out that the dominant mechanism responsible for the precession is simply the imprint of the jet-carrying SMBH orbital speed on the jet. Gravitational deflection and Lense-Thirring precession (due to the gravitational field of the other SMBH) are second order effects. We complement the scenario with a kinematical jet model which is inspired to the spine-sheath structure observed in M87. One of the main advantages of such a structure is that it allows the peak of the synchrotron emission to scale with frequency according to $\nu F\propto \nu^{\xi}$ as the viewing angle is changed, where $\xi$ is not necessarily 3 or 4 as in the case of jets with uniform velocity, but can be $\xi \sim 1$. Finally, we apply the model to the source PG1553+113, which has been recently claimed to show a $T_{\rm obs}=(2.18\pm 0.08)\text{ yr}$ periodicity. We are able to reproduce the optical and gamma-ray light curves and multiple synchrotron spectra simultaneously. We also give estimates of the source mass and size.

We introduce the Evolution of 21-cm Structure (EOS) project: providing periodic, public releases of the latest cosmological 21-cm simulations. 21-cm interferometry is set to revolutionize studies of the Cosmic Dawn (CD) and epoch of reionization (EoR), eventually resulting in 3D maps of the first billion years of our Universe. Progress will depend on sophisticated data analysis pipelines, which are in turn tested on large-scale mock observations. Here we present the 2016 EOS data release, consisting of the largest (1.6 Gpc on side with a 1024^3 grid), public 21-cm simulations of the CD and EoR. We include calibrated, sub-grid prescriptions for inhomogeneous recombinations and photo-heating suppression of star formation in small mass galaxies. We present two simulation runs that approximately bracket the contribution from faint unseen galaxies. From these two extremes, we predict that the duration of reionization (defined as a change in the mean neutral fraction from 0.9 to 0.1) should be between 2.7 < Delta z < 5.7. The large-scale 21-cm power during the advanced EoR stages can be different by up to a factor of ~10, depending on the model. This difference has a comparable contribution from: (i) the typical bias of sources; and (ii) a more efficient negative feedback in models with an extended EoR driven by faint galaxies. We also make detectability forecasts. With a 1000h integration, HERA and SKA1-low should achieve a signal-to-noise of ~few-hundreds throughout the EoR/CD, while in the maximally optimistic scenario of perfect foreground cleaning, all instruments should make a statistical detection of the cosmic signal. We also caution that our ability to clean foregrounds determines the relative performance of narrow/deep vs. wide/shallow surveys expected with SKA1. Our 21-cm power spectra, simulation outputs and visualizations are publicly available.

Interferometry of the cosmic 21-cm signal is set to revolutionize our understanding of the Epoch of Reionization (EoR), eventually providing 3D maps of the early Universe. Initial detections however will be low signal-to-noise, limited by systematics. To confirm a putative 21-cm detection, and check the accuracy of 21-cm data analysis pipelines, it would be very useful to cross-correlate against a genuine cosmological signal. The most promising cosmological signals are wide-field maps of Lyman alpha emitting galaxies (LAEs), expected from the Subaru Hyper-Suprime Cam (HSC) Ultra-Deep field. Here we present estimates of the correlation between LAE maps at z~7 and the 21-cm signal observed by both the Low Frequency Array (LOFAR) and the planned Square Kilometer Array Phase 1 (SKA1). We adopt a systematic approach, varying both: (i) the prescription of assigning LAEs to host halos; and (ii) the large-scale structure of neutral and ionized regions (i.e. EoR morphology). We find that the LAE-21cm cross-correlation is insensitive to (i), thus making it a robust probe of the EoR. A 1000h observation with LOFAR would be sufficient to discriminate at >1 standard deviation a fully ionized Universe from one with a mean neutral fraction of xHI~0.50, using the LAE-21cm cross-correlation function on scales of R~3-10 Mpc. Unlike LOFAR, whose detection of the LAE-21cm cross-correlation is limited by noise, SKA1 is mostly limited by ignorance of the EoR morphology. However, the planned 100h wide-field SKA1-Low survey will be sufficient to discriminate an ionized Universe from one with xHI~0.25, even with maximally pessimistic assumptions.

The Ly$\alpha$ line of high-redshift galaxies has emerged as a powerful probe of both early galaxy evolution and the epoch of reionization (EoR). Motivated by the upcoming wide-field survey with the Subaru Hyper Suprime-Cam (HSC), we study the angular correlation function (ACF) of narrow-band selected, $z\approx7$ Ly$\alpha$ emitting galaxies (LAEs). The clustering of LAEs is determined by both: (i) their typical host halo masses, $\bar{M}_{\rm h}$; and (ii) the absorption due to a patchy EoR, characterized by an average neutral fraction of the intergalactic medium, $\bar{x}_{\rm HI}$. We bracket the allowed LAE ACFs by exploring extreme scenarios for both the intrinsic Ly$\alpha$ emission and the large-scale pattern (i.e. morphology) of cosmic ionized patches in physical EoR models. Independent of the EoR morphology, current $z\approx 7$ ACF measurements constrain $\bar{x}_{\rm HI}\lesssim 0.5$ (1-$\sigma$). We also find that the low values of the currently-observed ACF imply that LAEs are hosted by relatively small dark matter halos: $\bar{M}_{\rm h} \lesssim 10^{10} M_{\odot}$, with corresponding duty cycles of $\lesssim$ few per cent. These values are over an order of magnitude lower than the analogous ones for color-selected, Lyman break galaxies, suggesting that $z\approx 7$ narrow-band LAEs searches are preferentially selecting young, star-burst galaxies, residing in less massive halos. The upcoming Ultra Deep campaign with the HSC will significantly improve constraints on both the EoR and LAE host halos.

We investigate the impact of sinks of ionizing radiation on the reionization-era 21-cm signal, focusing on 1-point statistics. We consider sinks in both the intergalactic medium and inside galaxies. At a fixed filling factor of HII regions, sinks will have two main effects on the 21-cm morphology: (i) as inhomogeneous absorbers of ionizing photons they result in smaller and more widespread cosmic HII patches; and (ii) as reservoirs of neutral gas they contribute a non-zero 21-cm signal in otherwise ionized regions. Both effects damp the contrast between neutral and ionized patches during reionization, making detection of the epoch of reionization with 21-cm interferometry more challenging. Here we systematically investigate these effects using the latest semi-numerical simulations. We find that sinks dramatically suppress the peak in the redshift evolution of the variance, corresponding to the midpoint of reionization. As previously predicted, skewness changes sign at midpoint, but the fluctuations in the residual HI suppress a late-time rise. Furthermore, large levels of residual HI dramatically alter the evolution of the variance, skewness and power spectrum from that seen at lower levels. In general, the evolution of the large-scale modes provides a better, cleaner, higher signal-to-noise probe of reionization.

We consider the electrodynamics of in-spiraling binary pulsars, showing that there are two distinct ways in which they may emit radiation. On the one hand, even if the pulsars do not rotate, we show that in vacuo orbital rotation generates magnetic quadrupole emission, which, in the late stages of the binary evolution becomes nearly as effective as magnetic dipole emission by a millisecond pulsar. On the other hand, we show that interactions of the two magnetic fields generate powerful induction electric fields, which cannot be screened by a suitable distribution of charges and currents like they are in isolated pulsars. We compute approximate electromotive forces for this case.

By depleting the ionizing photon budget available to expand cosmic HII regions, recombining systems (or Lyman limit systems) can have a large impact during (and following) cosmic reionization. Unfortunately, directly resolving such structures in large-scale reionization simulations is computationally impractical. Instead, here we implement a sub-grid prescription for tracking inhomogeneous recombinations in the intergalactic medium. Building on previous work parameterizing photo-heating feedback on star-formation, we present large-scale, semi-numeric reionization simulations which self-consistently track the local (sub-grid) evolution of both sources and sinks of ionizing photons. Our simple, single-parameter model naturally results in both an extended reionization and a modest, slowly-evolving emissivity, consistent with observations. Recombinations are instrumental in slowing the growth of large HII regions, and damping the rapid rise of the ionizing background in the late stages of (and following) reionization. As a result, typical HII regions are smaller by factors of $\sim 2-3$ throughout reionization. The large-scale ($k\lesssim 0.2\text{ Mpc$^{-1}$}$) ionization power spectrum is suppressed by factors of $\gtrsim 2-3$ in the second half of reionization. Therefore properly modeling recombinations is important in interpreting virtually all reionization observables, including upcoming interferometry with the redshifted 21 cm line. Consistent with previous works, we find the clumping factor of ionized gas to be $C_{\rm HII}\sim 4$ at the end of reionization.

An ionizing UV background (UVB) inhibits gas accretion and photo-evaporates gas from the shallow potential wells of small, dwarf galaxies. During cosmological reionization, this effect can result in negative feedback: suppressing star-formation inside HII regions, thus impeding their continued growth. It is difficult to model this process, given the enormous range of scales involved. We tackle this problem using a tiered approach: combining parameterized results from single-halo collapse simulations with large-scale models of reionization. In the resulting reionization models, the ionizing emissivity of galaxies depends on the local values of the reionization redshift and the UVB intensity. We present a physically-motivated analytic expression for the average minimum mass of star-forming galaxies, which can be readily used in modeling galaxy formation. We find that UVB feedback: (i) delays the end stages of reionization by less than 0.5 in redshift; (ii) results in a more uniform distribution of HII regions, peaked on smaller-scales (with large-scale ionization power suppressed by tens of percent); and (iii) suppresses the global photoionization rate per baryon by a factor of < 2 towards the end of reionization. However, the impact is modest, since the hydrodynamic response of the gas to the UVB occurs on a time-scale comparable to reionization. In particular, the popular approach of modeling UVB feedback with an instantaneous transition in the minimum mass of star-forming galaxies, dramatically overestimates its importance. UVB feedback does not significantly affect reionization unless: (i) molecularly-cooled galaxies contribute significantly to reionization; or (ii) internal feedback processes strongly couple with UVB feedback in the early Universe. Since both are considered unlikely, we conclude that there is no significant self-regulation of reionization by UVB feedback.

The reionization of the intergalactic medium (IGM) was likely inhomogeneous and extended. By heating the IGM and photo-evaporating gas from the outskirts of galaxies, this process can have a dramatic impact on the growth of structures. Using a suite of spherically-symmetric collapse simulations spanning a large parameter space, we study the impact of an ionizing ultraviolet background (UVB) on the condensation of baryons onto dark matter halos. We present an expression for the halo baryon fraction, which is an explicit function of: (i) halo mass; (ii) UVB intensity; (iii) redshift; (iv) redshift at which the halo was exposed to a UVB. We also present a corresponding expression for the characteristic or critical mass, defined as the halo mass which retains half of its baryons compared to the global value. Since our results are general and physically-motivated, they can be broadly applied to inhomogeneous reionization models.