We study the highly magnified arc SGAS J122651.3+215220 caused by a star-forming galaxy at $z_s=2.93$ crossing the lensing caustic cast by the galaxy cluster SDSS J1226+2152 ($z_l=0.43$), using Hubble Space Telescope observations. We report in the arc several asymmetric surface brightness features whose angular separations are a fraction of an arcsecond from the lensing critical curve and appear to be highly but unequally magnified image pairs of underlying compact sources, with one brightest pair having clear asymmetry consistently across four filters. One explanation of unequal magnification is microlensing by intracluster stars, which induces independent flux variations in the images of individual or groups of source stars in the lensed galaxy. For a second possibility, intracluster dark matter subhalos invisible to telescopes effectively perturb lensing magnifications near the critical curve and give rise to persistently unequal image pairs. Our modeling suggests, at least for the most prominent identified image pair, that the microlensing hypothesis is in tension with the absence of notable asymmetry variation over a six-year baseline, while subhalos of $\sim 10^6$--$10^8\,M_\odot$ anticipated from structure formation with Cold Dark Matter typically produce stationary and sizable asymmetries. We judge that observations at additional times and more precise lens models are necessary to stringently constrain temporal variability and robustly distinguish between the two explanations. The arc under this study is a scheduled target of a Director's Discretionary Early Release Science program of the James Webb Space Telescope, which will provide deep images and a high-resolution view with integral field spectroscopy.

We report the detection of new binary black hole merger events in the publicly available data from the second observing run of advanced LIGO and advanced Virgo (O2). The mergers were discovered using the new search pipeline described in Venumadhav et al. [Phys. Rev. D 100, 023011 (2019)], and are above the detection thresholds as defined in Abbott et al. (LIGO Scientific and Virgo Collaborations) [Phys. Rev. X 9, 031040 (2019).] Three of the mergers (GW170121, GW170304, GW170727) have inferred probabilities of being of astrophysical origin $p_{\rm astro} > 0.98$. The remaining three (GW170425, GW170202, GW170403) are less certain, with $p_{\rm astro}$ ranging from 0.5 to 0.8. The newly found mergers largely share the statistical properties of previously reported events, with the exception of GW170403, the least secure event, which has a highly negative effective spin parameter $\chi_{\rm eff}$ . The most secure new event, GW170121 ($p_{\rm astro} > 0.99$), is also notable due to its inferred negative value of $\chi_{\rm eff}$, which is inconsistent with being positive at the ~95.8% confidence level. The new mergers nearly double the sample of gravitational wave events reported from O2, and present a substantial opportunity to explore the statistics of the binary black hole population in the Universe. The number of detected events is not surprising since we estimate that the detection volume of our pipeline may be larger than that of other pipelines by as much as a factor of two (with significant uncertainties in the estimate). The increase in volume is larger when the constituent detectors of the network have very different sensitivities, as is likely to be the case in current and future runs.

We perform a detailed parameter estimation study of binary black hole merger events reported in Zackay et al. and Venumadhav et al.. These are some of the faintest signals reported so far, and hence, relative to the loud events in the GWTC-1 catalog, the data should have lesser constraining power on their intrinsic parameters. Hence we examine the robustness of parameter inference to choices made in the analysis, as well as any potential systematics. We check the impact of different methods of estimating the noise power spectral density, different waveform models, and different priors for the compact object spins. For most of the events, the resulting differences in the inferred values of the parameters are much smaller than their statistical uncertainties. The estimation of the effective spin parameter $\chi_{\mathrm{eff}}$, i.e. the projection of the mass-weighted total spin along the angular momentum, can be sensitive to analysis choices for two of the sources with the largest effective spin magnitudes, GW151216 and GW170403. The primary differences arise from using a 3D isotropic spin prior: the tails of the posterior distributions should be interpreted with care and due consideration of the other data analysis choices.

The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA as a facility for studying physics. We review four areas in which the SKA is expected to make major contributions to our understanding of fundamental physics: cosmic dawn and reionisation; gravity and gravitational radiation; cosmology and dark energy; and dark matter and astroparticle physics. These discussions demonstrate that the SKA will be a spectacular physics machine, which will provide many new breakthroughs and novel insights on matter, energy and spacetime.

We introduce a new technique to search for gravitational wave events from compact binary mergers that produce a clear signal only in a single gravitational wave detector, and marginal signals in other detectors. Such a situation can arise when the detectors in a network have different sensitivities, or when sources have unfavorable sky locations or orientations. We start with a short list of loud single-detector triggers from regions of parameter space that are empirically unaffected by glitches (after applying signal-quality vetoes). For each of these triggers, we compute evidence for astrophysical origin from the rest of the detector network by coherently combining the likelihoods from all detectors and marginalizing over extrinsic geometric parameters. We report the discovery of two new binary black hole (BBH) mergers in the second observing run of Advanced LIGO and Virgo (O2), in addition to the ones that were reported in Abbott et al. (2018) and Venumadhav et al. (2019). We estimate that the two events have false alarm rates of one in 19 years (60 O2) and one in 11 years (36 O2). One of the events, GW170817A, has primary and secondary masses $m_1^{\rm src} = 56_{-10}^{+16} \, M_\odot$ and $m_2^{\rm src} = 40_{-11}^{+10} \, M_\odot$ in the source frame. The existence of GW170817A should be very informative about the theoretically predicted upper mass gap for stellar mass black holes. Its effective spin parameter is measured to be $\chi_{\rm eff} = 0.5 \pm 0.2$, which is consistent with the tendency of the heavier detected BBH systems to have large and positive effective spin parameters. The other event, GWC170402, will be discussed thoroughly in future work.

Searches for gravitational waves crucially depend on exact signal processing of noisy strain data from gravitational wave detectors, which are known to exhibit significant non-Gaussian behavior. In this paper, we study two distinct non-Gaussian effects in the LIGO/Virgo data which reduce the sensitivity of searches: first, variations in the noise power spectral density (PSD) on timescales of more than a few seconds; and second, loud and abrupt transient `glitches' of terrestrial or instrumental origin. We derive a simple procedure to correct, at first order, the effect of the variation in the PSD on the search background. Given the knowledge of the existence of localized glitches in particular segments of data, we also develop a method to insulate statistical inference from these glitches, so as to cleanly excise them without affecting the search background in neighboring seconds. We show the importance of applying these methods on the publicly available LIGO data, and measure an increase in the detection volume of at least $15\%$ from the PSD-drift correction alone, due to the improved background distribution.

We introduce an algorithm for placing template waveforms for the search of compact binary mergers in gravitational wave interferometer data. We exploit the smooth dependence of the amplitude and unwrapped phase of the frequency-domain waveform on the parameters of the binary. We group waveforms with similar amplitude profiles and perform a singular value decomposition of the phase profiles to obtain an orthonormal basis for the phase functions. The leading basis functions span a lower-dimensional linear space in which the unwrapped phase of any physical waveform is well approximated. The optimal template placement is given by a regular grid in the space of linear coefficients. The algorithm is applicable to any frequency-domain waveform model and detector sensitivity curve. It is computationally efficient and requires little tuning. Applying this method, we construct a set of template banks suitable for the search of aligned-spin binary neutron star, neutron-star--black-hole and binary black hole mergers in LIGO--Virgo data.

We report a new binary black hole merger in the publicly available LIGO First Observing Run (O1) data release. The event has an inverse false alarm rate of one per six years in the detector-frame chirp-mass range $\mathcal{M}^{\rm det} \in [20,40]M_\odot$ in a new independent analysis pipeline that we developed. Our best estimate of the probability that the event is of astrophysical origin is $P_{\rm astro} \sim 0.71\, .$ The estimated physical parameters of the event indicate that it is the merger of two massive black holes, $\mathcal{M}^{\rm det} = 31^{+2}_{-3}\,M_\odot$ with an effective spin parameter, $\chi_{\rm eff} = 0.81^{+0.15}_{-0.21}$, making this the most highly spinning merger reported to date. It is also among the two highest redshift mergers observed so far. The high aligned spin of the merger supports the hypothesis that merging binary black holes can be created by binary stellar evolution.

In this paper, we report on the construction of a new and independent pipeline for analyzing the public data from the first observing run of advanced LIGO for mergers of compact binary systems. The pipeline incorporates different techniques and makes independent implementation choices in all its stages including the search design, the method to construct template banks, the automatic routines to detect bad data segments ("glitches") and to insulate good data from them, the procedure to account for the non-stationary nature of the detector noise, the signal-quality vetoes at the single-detector level and the methods to combine results from multiple detectors. Our pipeline enabled us to identify a new binary black-hole merger GW151216 in the public LIGO data. This paper serves as a bird's eye view of the pipeline's important stages. Full details and derivations underlying the various stages will appear in accompanying papers.

We show that recent observations of the compact binary, AM CVn type system, ES Ceti are fully consistent with theoretical predictions of stable mass transfer moderated by angular momentum loss due to gravitational-wave radiation. One of the main predictions of this model (for degenerate donors) is a widening of the binary. The mass transfer rate inferred from the observed rate of change in the orbital frequency is consistent with that inferred from the observed flux using the recent Gaia DR2 parallax

We examine a caustic-straddling arc at $z=0.9397$ in the field of the galaxy cluster MACS J0416.1-2403 ($z=0.397$) using archival multiband HST images and show that its surface brightness exhibits anomalies that can be explained by a single highly magnified star undergoing microlensing. First, we show that the surface brightness pattern is not perfectly symmetric across the cluster critical curve, which is inconsistent with a locally smooth lens model; the location of the candidate star exhibits the most significant asymmetry. Second, our analysis indicates that the asymmetric feature has $\sim 30\%$ higher flux in the 2012 visits compared to the Frontier Fields program visits in 2014. Moreover, the variable asymmetric feature shows an anomalous color between the F814W and F105W filters in 2014. These anomalies are naturally explained by microlensing induced variability of a caustic-transiting blue supergiant in a star-forming region, with a mean magnification factor around $\mu \sim 200$. We extend this study to a statistical analysis of the whole arc image and find tentative evidence of the increased mismatch of the two images in the proximity of the critical line. Robust detection of one or multiple caustic-transiting stars in this arc will enable detailed follow-up studies that can shed light on the small-scale structure of the dark matter inside the cluster halo.

The dynamical assembly of binary black holes (BBHs) in dense star clusters (SCs) is one of the most promising pathways for producing observable gravitational wave (GW) sources, however several other formation scenarios likely operate as well. One of the current outstanding questions is how these different pathways may be distinguished apart. In this paper we suggest a new multi-messenger observable that can be used to constrain the formation of BBH mergers originating from SCs: the electromagnetic signal from tidal disruptions (TDs) of stars by BBHs. Such TDs will show variability in their light curve from the orbital motion of the disruptive BBHs, and can therefore be used to map the BBH orbital period distribution, and thereby also the dynamical mechanisms that eventually drive the BBHs to merger. Using an analytical approach including General Relativistic effects, we find that the orbital period distribution of BBHs within globular clusters peaks on timescales of days, which we argue is unique to this assembly pathway. We propose that the search for variable TDs in current and future EM transient surveys might be used to constrain the merger history of BBHs in SCs.

Luminous stars in background galaxies straddling the lensing caustic of a foreground galaxy cluster can be individually detected due to extreme magnification factors of $\sim 10^2$--$10^3$, as recently observed in deep HST images. We propose a direct method to probe the presence of dark matter subhalos in galaxy clusters by measuring the astrometric perturbation they induce on the image positions of magnified stars or bright clumps: lensing by subhalos breaks the symmetry of a smooth critical curve, traced by the midpoints of close image pairs. For the giant arc at $z = 0.725$ behind the lensing cluster Abell 370 at $z = 0.375$, a promising target for detecting image pairs of stars, we find that subhalos of masses in the range $10^6$--$10^8\,M_\odot$ with the abundance predicted in the cold dark matter theory should typically imprint astrometric distortions at the level of $20$--$80\,{\rm mas}$. We estimate that $\sim 10\,$hr integrations with JWST at $\sim 1$--$3\,\mu{\rm m}$ may uncover several magnified stars whose image doublets will reveal the subhalo-induced structures of the critical curve. This method can probe a dynamic range in the subhalo to cluster halo mass ratio $m/M \sim 10^{-7}$--$10^{-9}$, thereby placing new constraints on the nature of dark matter.

We present a method to accelerate the evaluation of the likelihood in gravitational wave parameter estimation. Parameter estimation codes compute likelihoods of similar waveforms, whose phases and amplitudes differ smoothly with frequency. We exploit this by precomputing frequency-binned overlaps of the best-fit waveform with the data. We show how these summary data can be used to approximate the likelihood of any waveform that is sufficiently probable within the required accuracy. We demonstrate that $\simeq 60$ bins suffice to accurately compute likelihoods for strain data at a sampling rate of $4096\,$Hz and duration of $T=2048\,$s around the binary neutron star merger GW170817. Relative binning speeds up parameter estimation for frequency domain waveform models by a factor of $\sim 10^4$ compared to naive matched filtering and $\sim 10$ compared to reduced order quadrature.

Relative binning is a new method for fast and accurate evaluation of the likelihood of gravitational wave strain data. This technique can be used to produce reliable posterior distributions for compact object mergers with very moderate computational resources. We use a fast likelihood evaluation code based on this technique to estimate the parameters of the double neutron-star merger event GW170817 using publicly available LIGO data. We obtain statistically similar posteriors using either Markov-chain Monte-Carlo or nested sampling. The results do not favor non-zero aligned spins at a statistically significant level. There is no significant sign of non-zero tidal deformability (as quantified by the Bayesian evidence), whether or not high-spin or low-spin priors are adopted. Our posterior samples are publicly available, and we also provide a tutorial Python code to implement fast likelihood evaluation using the relative binning method.

We revisit the 21 cm power spectrum from the epoch of cosmic dawn in light of the recent EDGES detection of the 21 cm global signal at frequencies corresponding to $z\sim20$. The shape of the signal suggests that the spin temperature of neutral hydrogen was coupled to the kinetic temperature of the gas relatively rapidly ($19\lesssim z \lesssim 21$). We therefore consider models in which the UV photons were dominantly produced in the rarest and most massive halos ($M\gtrsim 10^9M_\odot$), since their abundance grows fast enough at those redshifts to account for this feature of the signal. We show that these models predict large power spectrum amplitudes during the inhomogeneous coupling, and then inhomogeneous heating by CMB and Lyman-$\alpha$ photons due to the large shot noise associated with the rare sources. The power spectrum is enhanced by more than an order of magnitude compared to previous models which did not include the shot noise contribution, making it a promising target for upcoming radio interferometers that aim to detect high-redshift 21 cm fluctuations.

We propose a new method to measure the tensor-to-scalar ratio $r$ using the circular polarization of the 21 cm radiation from the pre-reionization epoch. Our method relies on the splitting of the $F = 1$ hyperfine level of neutral hydrogen due to the quadrupole moment of the CMB. We show that unlike the Zeeman effect, where $M_{F}=\pm 1$ have opposite energy shifts, the CMB quadrupole shifts $M_{F}=\pm 1$ together relative to $M_{F}= 0$. This splitting leads to a small circular polarization of the emitted 21 cm radiation. In this paper (Paper I in a series on this effect), we present calculations on the microphysics behind this effect, accounting for all processes that affect the hyperfine transition. We conclude with an analytic formula for the circular polarization from the Dark Ages as a function of pre-reionization parameters and the value of the remote quadrupole of the CMB. We also calculate the splitting of the $F = 1$ hyperfine level due to other anisotropic radiation sources and show that they are not dominant. In a companion paper (Paper II) we make forecasts for measuring the tensor-to-scalar ratio $r$ using future radio arrays.

The intergalactic medium is expected to be at its coldest point before the formation of the first stars in the universe. Motivated by recent results from the EDGES experiment, we revisit the standard calculation of the kinetic temperature of the neutral gas through this period. When the first ultraviolet (UV) sources turn on, photons redshift into the Lyman lines of neutral hydrogen and repeatedly scatter within the Lyman-$\alpha$ line. They heat the gas via atomic recoils, and, through the Wouthuysen-Field effect, set the spin temperature of the 21-cm hyperfine (spin-flip) line of atomic hydrogen in competition with the resonant cosmic microwave background (CMB) photons. We show that the Lyman-$\alpha$ photons also mediate energy transfer between the CMB photons and the thermal motions of the hydrogen atoms. In the absence of X-ray heating, this new mechanism is the major correction to the temperature of the adiabatically cooling gas ($\sim 10 \%$ at $z=17$), and is several times the size of the heating rate found in previous calculations. We also find that the effect is more dramatic in non-standard scenarios that either enhance the radio background above the CMB or invoke new physics to cool the gas in order to explain the EDGES results. The coupling with the radio background can reduce the depth of the 21-cm absorption feature by almost a factor of two relative to the case with no sources of heating, and prevent the feature from developing a flattened bottom. As an inevitable consequence of the UV background that generates the absorption feature, this heating should be accounted for in any theoretical prediction.

We propose a method of measuring extremely weak magnetic fields in the intergalactic medium prior to and during the epoch of cosmic reionization. The method utilizes the Larmor precession of spin-polarized neutral hydrogen in the triplet state of the hyperfine transition. This precession leads to a systematic change in the brightness temperature fluctuations of the 21-cm line from the high-redshift universe, and thus the statistics of these fluctuations encode information about the magnetic field the atoms are immersed in. The method is most suited to probing fields that are coherent on large scales; in this paper, we consider a homogenous magnetic field over the scale of the 21-cm fluctuations. Due to the long lifetime of the triplet state of the 21-cm transition, this technique is naturally sensitive to extremely weak field strengths, of order $10^{-19}$ G at a reference redshift of $\sim 20$ (or $10^{-21}$ G if scaled to the present day). Therefore, this might open up the possibility of probing primordial magnetic fields just prior to reionization. If the magnetic fields are much stronger, it is still possible to use this method to infer their direction, and place a lower limit on their strength. In this paper (Paper I in a series on this effect), we perform detailed calculations of the microphysics behind this effect, and take into account all the processes that affect the hyperfine transition, including radiative decays, collisions, and optical pumping by Lyman-$\alpha$ photons. We conclude with an analytic formula for the brightness temperature of linear-regime fluctuations in the presence of a magnetic field, and discuss its limiting behavior for weak and strong fields.

Recent observations of lensed galaxies at cosmological distances have detected individual stars that are extremely magnified when crossing the caustics of lensing clusters. In idealized cluster lenses with smooth mass distributions, two images of a star of radius $R$ approaching a caustic brighten as $t^{-1/2}$ and reach a peak magnification $\sim 10^{6}\, (10\, R_{\odot}/R)^{1/2}$ before merging on the critical curve. We show that a mass fraction ($\kappa_\star \gtrsim \, 10^{-4.5}$) in microlenses inevitably disrupts the smooth caustic into a network of corrugated microcaustics, and produces light curves with numerous peaks. Using analytical calculations and numerical simulations, we derive the characteristic width of the network, caustic-crossing frequencies, and peak magnifications. For the lens parameters of a recent detection and a population of intracluster stars with $\kappa_\star \sim 0.01$, we find a source-plane width of $\sim 20 \, {\rm pc}$ for the caustic network, which spans $0.2 \, {\rm arcsec}$ on the image plane. A source star takes $\sim 2\times 10^4$ years to cross this width, with a total of $\sim 6 \times 10^4$ crossings, each one lasting for $\sim 5\,{\rm hr}\,(R/10\,R_\odot)$ with typical peak magnifications of $\sim 10^{4} \left( R/ 10\,R_\odot \right)^{-1/2}$. The exquisite sensitivity of caustic-crossing events to the granularity of the lens-mass distribution makes them ideal probes of dark matter components, such as compact halo objects and ultralight axion dark matter.

Strong lensing by intervening galaxies can produce multiple images of gravitational waves from sources at cosmological distances. These images acquire additional phase-shifts as the over-focused wavefront passes through itself along the line of sight. Time domain waveforms of Type-II images (associated with saddle points of the time delay) exhibit a non-trivial distortion from the unlensed waveforms. This phenomenon is in addition to the usual frequency-independent magnification, and happens even in the geometric limit where the wavelength is much shorter than the deflector's gravitational length scale. Similarly, Type-III images preserve the original waveform's shape but exhibit a sign flip. We show that for non-precessing binaries undergoing circular inspiral and merger, these distortions are equivalent to rotating the line of sight about the normal to the orbital plane by $45^\circ$ (Type II) and $90^\circ$ (Type III). This effect will enable us to distinguish between the different topological types among a set of multiple images, and give us valuable insight into the lens model. Furthermore, we show that for eccentric binaries, the waveform of a Type-II image is distorted in a manner that is inequivalent to a change of the source's orbital parameters.

The recent detection of gravitational waves indicates that stellar-mass black hole binaries are likely to be a key population of sources for forthcoming observations. With future upgrades, ground-based detectors could detect merging black hole binaries out to cosmological distances. Gravitational wave bursts from high redshifts ($z \gtrsim 1$) can be strongly magnified by gravitational lensing due to intervening galaxies along the line of sight. In the absence of electromagnetic counterparts, the mergers' intrinsic mass scale and redshift are degenerate with the unknown magnification factor $\mu$. Hence, strongly magnified low-mass mergers from high redshifts appear as higher-mass mergers from lower redshifts. We assess the impact of this degeneracy on the mass-redshift distribution of observable events for generic models of binary black hole formation from normal stellar evolution, Pop III star remnants, or a primordial black hole population. We find that strong magnification ($\mu \gtrsim 3$) generally creates a heavy tail of apparently massive mergers in the event distribution from a given detector. For LIGO and its future upgrades, this tail may dominate the population of intrinsically massive, but unlensed mergers in binary black hole formation models involving normal stellar evolution or primordial black holes. Modeling the statistics of lensing magnification can help account for this magnification bias when testing astrophysical scenarios of black hole binary formation and evolution.

The sky-averaged, or global, background of redshifted $21$ cm radiation is expected to be a rich source of information on cosmological reheating and reionizaton. However, measuring the signal is technically challenging: one must extract a small, frequency-dependent signal from under much brighter spectrally smooth foregrounds. Traditional approaches to study the global signal have used single antennas, which require one to calibrate out the frequency-dependent structure in the overall system gain (due to internal reflections, for example) as well as remove the noise bias from auto-correlating a single amplifier output. This has motivated proposals to measure the signal using cross-correlations in interferometric setups, where additional calibration techniques are available. In this paper we focus on the general principles driving the sensitivity of the interferometric setups to the global signal. We prove that this sensitivity is directly related to two characteristics of the setup: the cross-talk between readout channels (i.e. the signal picked up at one antenna when the other one is driven) and the correlated noise due to thermal fluctuations of lossy elements (e.g. absorbers or the ground) radiating into both channels. Thus in an interferometric setup, one cannot suppress cross-talk and correlated thermal noise without reducing sensitivity to the global signal by the same factor -- instead, the challenge is to characterize these effects and their frequency dependence. We illustrate our general theorem by explicit calculations within toy setups consisting of two short dipole antennas in free space and above a perfectly reflecting ground surface, as well as two well-separated identical lossless antennas arranged to achieve zero cross-talk.

In the first paper of this series, we proposed a novel method to probe large-scale intergalactic magnetic fields during the cosmic Dark Ages, using 21-cm tomography. This method relies on the effect of spin alignment of hydrogen atoms in a cosmological setting, and on the effect of magnetic precession of the atoms on the statistics of the 21-cm brightness-temperature fluctuations. In this paper, we forecast the sensitivity of future tomographic surveys to detecting magnetic fields using this method. For this purpose, we develop a minimum-variance estimator formalism to capture the characteristic anisotropy signal using the two-point statistics of the brightness-temperature fluctuations. We find that, depending on the reionization history, and subject to the control of systematics from foreground subtraction, an array of dipole antennas in a compact-grid configuration with a collecting area slightly exceeding one square kilometer can achieve a $1\sigma$ detection of $\sim$$10^{-21}$ Gauss comoving (scaled to present-day value) within three years of observation. Using this method, tomographic 21-cm surveys could thus probe ten orders of magnitude below current CMB constraints on primordial magnetic fields, and provide exquisite sensitivity to large-scale magnetic fields in situ at high redshift.

In this paper, we study small-scale fluctuations (baryon pressure sound waves) in the baryon fluid during recombination. In particular, we look at their evolution in the presence of relative velocities between baryons and photons on large scales ($k \sim 10^{-1} \ {\rm Mpc}^{-1}$), which are naturally present during the era of decoupling. Previous work concluded that the fluctuations grow due to an instability of sound waves in a recombining plasma, but that the growth factor is small for typical cosmological models. These analyses model recombination in an inhomogenous universe as a perturbation to the parameters of the homogenous solution. We show that for relevant wavenumbers $k\gtrsim 10^3\ {\rm Mpc}^{-1}$ the dynamics are significantly altered by the transport of both ionizing continuum ($h\nu>13.6$ eV) and Lyman-$\alpha$ photons between crests and troughs of the density perturbations. We solve the radiative transfer of photons in both these frequency ranges and incorporate the results in a perturbed three-level atom model. We conclude that the instability persists at intermediate scales. We use the results to estimate a distribution of growth rates in $10^{7}$ random realizations of large-scale relative velocities. Our results indicate that there is no appreciable growth; out of these $10^7$ realizations, the maximum growth factor we find is less than $\approx 1.2$ at wavenumbers of $k \approx 10^{3} \ {\rm Mpc}^{-1}$. The instability's low growth factors are due to the relatively short duration of the recombination epoch during which the electrons and photons are coupled.

We perform a detailed study of the weak interactions of standard model neutrinos with the primordial plasma and their effect on the resonant production of sterile neutrino dark matter. Motivated by issues in cosmological structure formation on small scales, and reported X-ray signals that could be due to sterile neutrino decay, we consider $7$ keV-scale sterile neutrinos. Oscillation-driven production of such sterile neutrinos occurs at temperatures $T \gtrsim 100$ MeV, where we study two significant effects of weakly charged species in the primordial plasma: (1) the redistribution of an input lepton asymmetry; (2) the opacity for active neutrinos. We calculate the redistribution analytically above and below the quark-hadron transition, and match with lattice QCD calculations through the transition. We estimate opacities due to tree level processes involving leptons and quarks above the quark-hadron transition, and the most important mesons below the transition. We report final sterile neutrino dark matter phase space densities that are significantly influenced by these effects, and yet relatively robust to remaining uncertainties in the nature of the quark-hadron transition. We also provide transfer functions for cosmological density fluctuations with cutoffs at $k \simeq 10 \ h \ {\rm Mpc}^{-1}$, that are relevant to galactic structure formation.

SPHEREx (Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer) ( http://spherex.caltech.edu ) is a proposed all-sky spectroscopic survey satellite designed to address all three science goals in NASA's Astrophysics Division: probe the origin and destiny of our Universe; explore whether planets around other stars could harbor life; and explore the origin and evolution of galaxies. SPHEREx will scan a series of Linear Variable Filters systematically across the entire sky. The SPHEREx data set will contain R=40 spectra fir 0.75$<\lambda<$4.1$\mu$m and R=150 spectra for 4.1$<\lambda<$4.8$\mu$m for every 6.2 arc second pixel over the entire-sky. In this paper, we detail the extra-galactic and cosmological studies SPHEREx will enable and present detailed systematic effect evaluations. We also outline the Ice and Galaxy Evolution Investigations.

It has recently been suggested that the tidal deformation of a neutron star excites daughter p- and g-modes to large amplitudes via a quasi-static instability. This would remove energy from the tidal bulge, resulting in dissipation and possibly affecting the phase evolution of inspiralling binary neutron stars and hence the extraction of binary parameters from gravitational wave observations. This instability appears to arise because of a large three-mode interaction among the tidal mode and high-order p- and g-modes of similar radial wavenumber. We show that additional four-mode interactions enter into the analysis at the same order as the three-mode terms previously considered. We compute these four-mode couplings by finding a volume-preserving coordinate transformation that relates the energy of a tidally deformed star to that of a radially perturbed spherical star. Using this method, we relate the four-mode coupling to three-mode couplings and show that there is a near-exact cancellation between the destabilizing effect of the three-mode interactions and the stabilizing effect of the four-mode interaction. We then show that the equilibrium tide is stable against the quasi-static decay into daughter p- and g-modes to leading order. The leading deviation from the quasi-static approximation due to orbital motion of the binary is considered; while it may slightly spoil the near-cancellation, any resulting instability timescale is at least of order the gravitational-wave inspiral time. We conclude that the p-/g-mode coupling does not lead to a quasi-static instability, and does not impact the phase evolution of gravitational waves from binary neutron stars.