The reionization of the second electron of helium shapes the physical state of intergalactic gas at redshifts between 2 < z < 5. Because performing full in situ radiative transfer in hydrodynamic simulations is computationally expensive for large volumes, the physics of HeII reionization is often approximated by a uniform UV background model that does not capture the patchiness of reionization. We have devised a model that implements the effects of HeII reionization using semi-analytic calculations of the thermal state of intergalactic gas-- a way to bypass a full radiative transfer simulation while still realizing the physics of HeII reionization that affects observables such as the Lyman alpha forest. Here we present a publicly-available code that flexibly models inhomogeneous HeII reionization in simulations at a negligible computational cost. Because many of the parameters of HeII reionization are uncertain, our model is customizable from a set of free parameters. We show results from this code in MP-Gadget, where this model is implemented. We demonstrate the resulting temperature evolution and temperature-density relation of intergalactic gas-- consistent with recent measurements and previous radiative transfer simulations. We show that the impact of HeII reionization gives rise to subtle signatures in the one-dimensional statistics of the Lyman alpha forest at the level of several per cent, in agreement with previous findings. The flexible nature of these simulations are ideal for studies of HeII reionization and future observations of the HeII Lyman alpha forest.

Spatial fluctuations in ultraviolet backgrounds can subtly modulate the distribution of extragalactic sources, a potential signal and systematic for large-scale structure surveys. While this modulation has been shown to be significant for 3D Ly alpha forest surveys, its relevance for other large-scale structure probes has been hardly explored, despite being the only astrophysical process that likely can affect clustering measurements on the scales of >~Mpc. We estimate that the background fluctuations, modulating the amount of H I, have a fractional effect of (0.03-0.3) x (k/[10^-2 Mpc^-1])^-1 on the power spectrum of 21 cm intensity maps at z = 1-3. We find a smaller effect for H alpha and Ly alpha intensity mapping surveys of (0.001-0.1) x (k/[10^-2 Mpc^-1])-1 and even smaller effect for more traditional surveys that correlate the positions of individual H alpha or Ly alpha emitters. We also estimate the effect of backgrounds on low-redshift galaxy surveys in general based on a simple model in which background fluctuations modulate the rate halo gas cools, modulating star formation: We estimate a maximum fractional effect on the power of ~0.01 (k/[10-2 Mpc-1])-1 at z = 1. We compare sizes of these imprints to cosmological parameter benchmarks for the next generation of redshift surveys: We find that ionizing backgrounds could result in a bias on the squeezed triangle non-Gaussianity parameter fNL that can be larger than unity for power spectrum measurements with a SPHEREx-like galaxy survey, and typical values of intensity bias. Marginalizing over a shape of the form k-1PL, where PL is the linear matter power spectrum, removes much of this bias at the cost of ~ 40% larger statistical errors.

We explore scenarios for reionizing the intergalactic medium with low galaxy ionizing photon escape fractions. We combine simulation-based halo-mass dependent escape fractions with an extrapolation of the observed galaxy rest-ultraviolet luminosity functions to solve for the reionization history from z=20 to z=4. We explore the posterior distributions for key unknown quantities, including the limiting halo mass for star-formation, the ionizing photon production efficiency, and a potential contribution from active galactic nuclei (AGN). We marginalize over the allowable parameter space using a Markov Chain Monte Carlo method, finding a solution which satisfies the most model-independent constraints on reionization. Our fiducial model can match observational constraints with an average escape fraction of <5% throughout the bulk of the epoch of reionization if: i) galaxies form stars down to the atomic cooling limit before reionization and a photosuppression mass of log(M_h/Msol)~9 during/after reionization (-13<M_UV,lim<-11); ii) galaxies become more efficient producers of ionizing photons at higher redshifts and fainter magnitudes, and iii) there is a significant, but sub-dominant, contribution by AGN at z < 7. In this model the faintest galaxies (M_UV>-15) dominate the ionizing emissivity, leading to an earlier start to reionization and a smoother evolution of the ionized volume filling fraction than models which assume a single escape fraction at all redshifts and luminosities. The ionizing emissivity from this model is consistent with observations at z=4-5 (and below, when extrapolated), in contrast to some models which assume a single escape fraction. Our predicted ionized volume filling fraction at z=7 of Q_HII=78% (+\- 8%) is in ~1-2 sigma tension with observations of Lya emitters at z~7 and the damping wing analyses of the two known z>7 quasars, which prefer Q_HII,z=7~40-50%.

During reionization, the intergalactic medium is heated impulsively by supersonic ionization fronts (I-fronts). The peak gas temperatures behind the I-fronts, $T_\mathrm{reion}$, are a key uncertainty in models of the thermal history after reionization. Here we use high-resolution radiative transfer simulations to study the parameter space of $T_\mathrm{reion}$. We show that $T_\mathrm{reion}$ is only mildly sensitive to the spectrum of incident radiation over most of the parameter space, with temperatures set primarily by I-front speeds. We also explore what current models of reionization predict for $T_\mathrm{reion}$ by measuring I-front speeds in cosmological radiative transfer simulations. We find that the post-I-front temperatures evolve toward hotter values as reionization progresses. Temperatures of $T_\mathrm{reion} = 17,000-22,000$ K are typical during the first half of reionization, but $T_\mathrm{reion} = 25,000 - 30,000$ K may be achieved near the end of this process if I-front speeds reach $\sim10^4$ km/s as found in our simulations. Shorter reionization epochs lead to hotter $T_\mathrm{reion}$. We discuss implications for $z>5$ Ly$\alpha$ forest observations, which potentially include sight lines through hot, recently reionized patches of the Universe. Interpolation tables from our parameter space study are made publicly available, along with a simple fit for the dependence of $T_\mathrm{reion}$ on the I-front speed.

Radiation in the extreme ultraviolet (EUV) and soft X-ray holds clues to the location of the missing baryons, the energetics in stellar feedback processes, and the cosmic enrichment history. Additionally, EUV and soft X-ray photons help determine the ionization state of most intergalactic and circumgalactic metals, shaping the rate at which cosmic gas cools. Unfortunately, this band is extremely difficult to probe observationally due to absorption from the Galaxy. In this paper, we model the contributions of various sources to the cosmic EUV and soft X-ray backgrounds. We bracket the contribution from (1) quasars, (2) X-ray binaries, (3) hot interstellar gas, (4) circumgalactic gas, (5) virialized gas, and (6) supersoft sources, developing models that extrapolate into these bands using both empirical and theoretical inputs. While quasars are traditionally assumed to dominate these backgrounds, we discuss the substantial uncertainty in their contribution. Furthermore, we find that hot intrahalo gases likely emit an O(1) fraction of this radiation at low redshifts, and that interstellar and circumgalactic emission potentially contribute tens of percent to these backgrounds at all redshifts. We estimate that uncertainties in the angular-averaged background intensity impact the ionization corrections for common circumgalactic and intergalactic metal absorption lines by ~0.3-1 dex, and we show that local emissions are comparable to the cosmic background only at r_prox = 10-100 kpc from Milky Way-like galaxies.

Cosmological simulations of the low-density intergalactic medium exhibit a strikingly tight power-law relation between temperature and density that holds over two decades in density. It is found that this relation should roughly apply Delta z ~ 1-2 after a reionization event, and this limiting behavior has motivated the power-law parameterizations used in most analyses of the Ly-alpha forest. This relation has been explained by using equations linearized in the baryonic overdensity (which does not address why a tight power-law relation holds over two decades in density) or by equating the photoheating rate with the cooling rate from cosmological expansion (which we show is incorrect). Previous explanations also did not address why recombination cooling and Compton cooling off of the cosmic microwave background, which are never negligible, do not alter the character of this relation. We provide an understanding for why a tight power-law relation arises for unshocked gas at all densities for which collisional cooling is unimportant. We also use our results to comment on (1) how quickly fluctuations in temperature redshift away after reionization processes, (2) how much shock heating occurs in the low-density intergalactic medium, and (3) how the temperatures of collapsing gas parcels evolve.

Motivated by the claimed detection of a large population of faint active galactic nuclei (AGN) at high redshift, recent studies have proposed models in which AGN contribute significantly to the z > 4 H I ionizing background. In some models, AGN are even the chief sources of reionization. If correct, these models would make necessary a complete revision to the standard view that galaxies dominated the high-redshift ionizing background. It has been suggested that AGN-dominated models can better account for two recent observations that appear to be in conflict with the standard view: (1) large opacity variations in the z ~ 5.5 H I Lyman-alpha forest, and (2) slow evolution in the mean opacity of the He II Lyman-alpha forest. Large spatial fluctuations in the ionizing background from the brightness and rarity of AGN may account for the former, while the earlier onset of He II reionization in these models may account for the latter. Here we show that models in which AGN emissions source >~ 50 % of the ionizing background generally provide a better fit to the observed H I Lyman-alpha forest opacity variations compared to standard galaxy-dominated models. However, we argue that these AGN-dominated models are in tension with constraints on the thermal history of the intergalactic medium (IGM). Under standard assumptions about the spectra of AGN, we show that the earlier onset of He II reionization heats up the IGM well above recent temperature measurements. We further argue that the slower evolution of the mean opacity of the He II Lyman-alpha forest relative to simulations may reflect deficiencies in current simulations rather than favor AGN-dominated models as has been suggested.

Recent years have brought more precise temperature measurements of the low-density intergalactic medium (IGM). These new measurements constrain the processes that heated the IGM, such as the reionization of H I and of He II. We present a semi-analytical model for the thermal history of the IGM that follows the photoheating history of primordial gas. Our model adopts a multizone approach that, compared to previous models, more accurately captures the inhomogeneous heating and cooling of the IGM during patchy reionization. We compare our model with recent temperature measurements spanning z = 1.6-4.8, finding that these measurements are consistent with scenarios in which the He II was reionized at z = 3-4 by quasars. Significantly longer duration or higher redshift He II reionization scenarios are ruled out by the measurements. For hydrogen reionization, we find that only low-redshift and high-temperature scenarios are excluded. For example, a model in which the IGM was heated to 30 000 K when an ionization front passed, and with hydrogen reionization occurring over 6 < z < 9, is ruled out. Finally, we place constraints on how much heating could owe to TeV blazars, cosmic rays, and other non-standard mechanisms. We find that by z = 2, a maximum of 1 eV of additional heat could be injected per baryon over standard photoheating-only models, with this limit becoming <~0.5 eV at z > 3.

We present results from the Radio Interferometric Planet (RIPL) search for compan- ions to the nearby star GJ 896A. We present 11 observations over 4.9 years. Fitting astrometric parameters to the data reveals a residual with peak-to-peak amplitude of ~ 3 mas in right ascension. This residual is well-fit by an acceleration term of 0.458 \pm 0.032 mas/y^2. The parallax is fit to an accuracy of 0.2 mas and the proper motion terms are fit to accuracies of 0.01 mas/y. After fitting astrometric and acceleration terms residuals are 0.26 mas in each coordinate, demonstrating that stellar jitter does not limit the ability to carry out radio astrometric planet detection and characterization. The acceleration term originates in part from the companion GJ 896B but the amplitude of the acceleration in declination is not accurately predicted by the orbital model. The acceleration sets a mass upper limit of 0.15 MJ at a semi-major axis of 2 AU for a planetary companion to GJ 896A. For semi-major axes between 0.3 and 2 AU upper limits are determined by the maximum angular separation; the upper limits scale from the minimum value in proportion to the inverse of the radius. Upper limits at larger radii are set by the acceleration and scale as the radius squared. An improved solution for the stellar binary system could improve the exoplanet mass sensitivity by an order of magnitude.

Deep Very Large Array imaging of the binary X-ray source SS 433, sometimes classified as a microquasar, has been used to study the intrinsic brightness distribution and evolution of its radio jets. The intrinsic brightness of the jets as a function of age at emission of the jet material tau is recovered by removal of the Doppler boosting and projection effects. We find that intrinsically the two jets are remarkably similar when compared for equal tau, and that they are best described by Doppler boosting of the form D^{2+alpha}, as expected for continuous jets. The intrinsic brightnesses of the jets as functions of age behave in complex ways. In the age range 60 < tau < 150 days, the jet decays are best represented by exponential functions of tau, but linear or power law functions are not statistically excluded. This is followed by a region out to tau ~ 250 days during which the intrinsic brightness is essentially constant. At later times the jet decay can be fit roughly as exponential or power law functions of tau.