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.

Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and analyze such data. On the positive side, it has become clear that the outstanding issues are purely technical in nature and can be solved with sufficient development activity. Such activity will enable science across redshifts, from early galaxy evolution in the pre-reionization era to dark energy evolution at low redshift.

Two groundbreaking new facilities will commence operations early in the 2020s and thereafter define much of the broad landscape of US optical-infrared astronomy in the remaining decade. The Large Synoptic Survey Telescope (LSST), perched atop Cerro Pachon in the Chilean Andes, will revolutionize the young field of Time Domain Astronomy through its wide-field, multi-band optical imaging survey. At the same time, the James Webb Space Telescope (JWST), orbiting at the Sun-Earth L2 Lagrange point, will provide stunningly high-resolution views of selected targets from the red end of the optical spectrum to the mid-infrared. However, the spatial resolution of the LSST observations will be limited by atmospheric seeing, while JWST will be limited in its time-domain capabilities. This paper highlights the scientific opportunities lying between these two landmark missions, i.e., science enabled by systems capable of astronomical observations with both high cadence in the time domain and high resolution in the spatial domain. The opportunities range from constraining the late phases of stellar evolution in nearby resolved populations to constraining dark matter distributions and cosmology using lensed transient sources. We describe a system that can deliver the required capabilities.

We present a non-perturbative calculation of the 1-point probability distribution function (PDF) for the spherically-averaged matter density field. The PDF is represented as a path integral and is evaluated using the saddle-point method. It factorizes into an exponent given by a spherically symmetric saddle-point solution and a prefactor produced by fluctuations. The exponent encodes the leading sensitivity of the PDF to the dynamics of gravitational clustering and statistics of the initial conditions. In contrast, the prefactor has only a weak dependence on cosmology. It splits into a monopole contribution which is evaluated exactly, and a factor corresponding to aspherical fluctuations. The latter is crucial for the consistency of the calculation: neglecting it would make the PDF incompatible with translational invariance. We compute the aspherical prefactor using a combination of analytic and numerical techniques. We demonstrate the factorization of spurious enhanced contributions of large bulk flows and their cancellation due the equivalence principle. We also identify the sensitivity to the short-scale physics and argue that it must be properly renormalized. The uncertainty associated with the renormalization procedure gives an estimate of the theoretical error. For zero redshift, the precision varies from sub percent for moderate density contrasts to tens of percent at the tails of the distribution. It improves at higher redshifts. We compare our results with N-body simulation data and find an excellent agreement.

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.

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 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.

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.

The galaxies with photometric redshifts observed in a close angular proximity might be either projection coincidences, strongly lensed images of the same galaxy, or separate galaxies that are in a stage of merging. We search for the groups of galaxies in the Hubble Ultra Deep Field (HUDF09) in $z\sim7$ and $z\sim8$ drop-out samples. We find no close pairs among 50 galaxies in the $z\sim7$ sample, while in the $z\sim8$ sample we find that 6 out of 22 galaxies have a companion within $\sim1'$ (3 pairs). Adopting a numerical simulation and performing forward modeling we show that even though mergers are unlikely to have such a high fraction, the projection coincidences and the strong lensing are even less likely mechanisms to account for all of three pairs. Alternatively, there is a possibility of the contamination in the drop-out catalog from lower redshifts, which potentially can account for all of the groups. Finally, we make projection on the sensitivity to mergers of the James Webb Space Telescope, and discuss the possible applications of the high-redshift merging galaxies for decreasing cosmic variance effect on the luminosity function and for improving the accuracy of photometric redshifts in general.

The 21cm neutral hydrogen line is likely to be a key probe for studying the epoch of reionization and comic dawn in the forthcoming decades. This prospect stimulates the development of the theoretical basis for simulating the power spectrum of this line. Because of the beam size of the upcoming radio telescopes at high redshifts, most of the theoretical models are focused on the inhomogeneities on scales above few comoving megaparsecs. Therefore, smaller scales are often neglected and modeled with approximated sub-grid models. In this study we explore whether the perturbations on small scales ($\lesssim 1h^{-1}\mathrm{Mpc}$) can affect the 21cm signal on larger scales. Two distinct mechanism are discussed. First, we show that during the cosmic dawn small scale perturbations regulate the formation time of the first Lyman-alpha Emitters (LAE), and consequently the coupling timing of spin and kinetic temperatures. Due to the low number density of LAE, the inhomogeneity of coupling includes the shot noise and manifests itself in the observed 21cm power spectrum. Second mechanism works during the reionization when the ionized bubbles actively grow and overlap. Small scales perturbations affect the galactic properties and merger histories, and consequently the number of ionizing photons produced by each galaxy. The ionizing photons bring the perturbations from the galactic scales to the scales of ionizing fronts, affecting the 21cm power spectrum. We conclude that these two effects introduce stochasticity in the potentially observed 21cm power spectrum and, moreover, might give another perspective into the physics of the first galaxies.

The methods for studying the epoch of cosmic reionization vary from full radiative transfer simulations to purely analytical models. While numerical approaches are computationally expensive and are not suitable for generating many mock catalogs, analytical methods are based on assumptions and approximations. We explore the interconnection between both methods. First, we ask how the analytical framework of excursion set formalism can be used for statistical analysis of numerical simulations and visual representation of the morphology of ionization fronts. Second, we explore the methods of training the analytical model on a given numerical simulation. We present a new code which emerged from this study. Its main application is to match the analytical model with a numerical simulation. Then, it allows one to generate mock reionization catalogs with volumes exceeding the original simulation quickly and computationally inexpensively, meanwhile reproducing large scale statistical properties. These mock catalogs are particularly useful for CMB polarization and 21cm experiments, where large volumes are required to simulate the observed signal.

We explore a time-dependent energy dissipation of the energetic electrons in the inhomogeneous intergalactic medium (IGM) during the epoch of cosmic reionization. In addition to the atomic processes we take into account the Inverse Compton (IC) scattering of the electrons on the comic microwave background (CMB) photons, which is the dominant channel of energy loss for the electrons with energies above a few MeV. We show that: (1) the effect on the IGM has both local (atomic processes) and non-local (IC radiation) components; (2) the energy distribution between Hydrogen and Helium ionizations depends on the initial electron energy; (3) the local baryon overdensity significantly affects the fractions of energy distributed in each channel; and (4) the relativistic effect of atomic cross section become important during the epoch of cosmic reionization. We release our code as open source for further modification by the community.

We explore the mean and fluctuating redshifted 21 cm signal in numerical simulations of cosmic reionization from the Cosmic Reionization On Computers (CROC) project. We find that the mean signal varies between about $\pm20\rm{mK}$. Most significantly, we find that the negative pre-reionization dip at $z\sim10-15$ only extends to $\langle\Delta T_B\rangle\sim-20\rm{mK}$, in agreement with prior simulation results and in significant contrast to Pritchard & Loeb analytical model, requiring substantially higher sensitivity from global signal experiments that operate in this redshift range (EDGES-II, LEDA, SCI-HI, and DARE). We also explore the role of dense substructure (filaments and embedded galaxies) in the formation of 21 cm power spectrum. We find that by neglecting the semi-neutral substructure inside ionized bubbles, the power spectrum can be mis-estimated by 25-50\% at scales $k\sim 0.1-1h\rm{Mpc}^{-1}$. This scale range is of a particular interest, because the upcoming 21 cm experiments (MWA, PAPER, HERA) are expected to be most sensitive within it.

We revisit the possibility of constraining the properties of dark matter (DM) by studying the epoch of cosmic reionization. Previous studies have shown that DM annihilation was unlikely to have provided a large fraction of the photons that ionized the universe, but instead played a subdominant role relative to stars and quasars. The DM, however, begins to efficiently annihilate with the formation of primordial microhalos at $z\sim100-200$, much earlier than the formation of the first stars. Therefore, if DM annihilation ionized the universe at even the percent level over the interval $z \sim 20-100$, it can leave a significant imprint on the global optical depth, $\tau$. Moreover, we show that cosmic microwave background (CMB) polarization data and future 21 cm measurements will enable us to more directly probe the DM contribution to the optical depth. In order to compute the annihilation rate throughout the epoch of reionization, we adopt the latest results from structure formation studies and explore the impact of various free parameters on our results. We show that future measurements could make it possible to place constraints on the dark matter's annihilation cross section that are at a level comparable to those obtained from the observations of dwarf galaxies, cosmic ray measurements, and studies of recombination.

We use fully self-consistent numerical simulations of cosmic reionization, completed under the Cosmic Reionization On Computers (CROC) project, to explore how well the recombinations in the ionized IGM can be quantified by the effective "clumping factor." The density distribution in the simulations (and, presumably, in a real universe) is highly inhomogeneous and more-or-less smoothly varying in space. However, even in highly complex and dynamic environments the concept of the IGM remains reasonably well-defined; the largest ambiguity comes from the unvirialized regions around galaxies that are over-ionized by the local enhancement in the radiation field ("proximity zones"). That ambiguity precludes computing the IGM clumping factor to better than about 20\%. We also discuss a "local clumping factor," defined over a particular spatial scale, and quantify its scatter on a given scale and its variation as a function of scale.

We propose that observations of "hidden" magnetars in central compact objects can be used to probe crustal activity of neutron stars with large internal magnetic fields. Estimates based on calculations by Perna \& Pons (2011), Pons \& Rea (2012) and Kaminker et al. (2014) suggest that central compact objects, which are proposed to be "hidden" magnetars, must demonstrate flux variations on the time scale of months-years. However, the most prominent candidate for the "hidden" magnetars --- CXO J1852.6+0040 in Kes 79, -- shows constant (within error bars) flux. This can be interpreted by lower variable crustal activity than in typical magnetars. Alternatively, CXO J1852.6+0040 can be in a high state of variable activity during the whole period of observations. Then we consider the source 1E161348-5055 in RCW103 as another candidate. Employing a simple 2D-modeling we argue that properties of the source can be explained by the crustal activity of the magnetar type. Thus, this object may be supplemented for the three known candidates for the "hidden" magnetars among central compact objects discussed in literature.

We compare the results from several sets of cosmological simulations of cosmic reionization, produced under Cosmic Reionization On Computers (CROC) project, with existing observational data on the high-redshift Ly-alpha forest and the abundance of Ly-alpha emitters. We find good consistency with the observational measurements and the previous simulation work. By virtue of having several independent realizations for each set of numerical parameters, we are able to explore the effect of cosmic variance on observable quantities. One unexpected conclusion we are forced into is that cosmic variance is unusually large at z>6, with both our simulations and, most likely, observational measurements are still not fully converged for even such basic quantities as the average Gunn-Peterson optical depth or the volume-weighted neutral fraction. We also find that reionization has little effect on the early galaxies or on global cosmic star formation history, because galaxies whose gas content is affected by photoionization contain no molecular (i.e. star-forming) gas in the first place. In particular, measurements of the faint end of the galaxy luminosity function by JWST are unlikely to provide a useful constraint on reionization.

Using 1D and 2D cooling codes we study thermal emission from neutron stars with steady state internal heaters of various intensities and geometries (blobs or spherical layers) located at different depths in the crust. The generated heat tends to propagate radially, from the heater down to the stellar core and up to the surface; it is also emitted by neutrinos. In local regions near the heater the results are well described with the 1D code. The heater's region projects onto the stellar surface forming a hot spot. There are two heat propagation regimes. In the first, conduction outflow regime (realized at heat rates $H_0 \lesssim 10^{20}$ erg cm$^{-3}$ s$^{-1}$ or temperatures $T_\mathrm{h} \lesssim 10^9$ K in the heater) the thermal surface emission of the star depends on the heater's power and neutrino emission in the stellar core. In the second, neutrino outflow regime ($H_0 \gtrsim 10^{20}$ erg cm$^{-3}$ s$^{-1}$ or $T_\mathrm{h} \gtrsim 10^9$ K) the surface thermal emission becomes independent of heater's power and the physics of the core. The largest (a few per cent) fraction of heat power is carried to the surface if the heater is in the outer crust and the heat regime is intermediate. The results can be used for modeling young cooling neutron stars (prior to the end of internal thermal relaxation), neutron stars in X-ray transients, magnetars and high-$B$ pulsars, as well as merging neutron stars.

We discuss the role of recombinations in the IGM, and the related concept of the clumping factor, during cosmic reionization. The clumping factor is, in general, a local quantity that depends on the local over-density and the scale below which the baryon density field can be assumed smooth. That scale, called the filtering scale, is itself depended on over-density and local thermal history. We present a method for building a self-consistent analytical model of inhomogeneous reionization assuming the linear growth rate of the density fluctuation, which accounts for these effects simultaneously. We show that taking into account the local clumping factor introduces significant corrections to the total recombination rate, comparing to the model with a globally uniform clumping factor.

A statistical analysis of radial distributions of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey (SDSS DR7) catalogue within an interval $0.16 \leq z \leq 0.47$ is carried out. We found that the radial distribution of $\sim$ 106,000 LRGs incorporates a few quasi-periodical components relatively to a variable $\eta$, dimensionless line-of-sight comoving distance calculated for the $\Lambda$CDM cosmological model. The most significant peaks of the power spectra are obtained for two close periodicities corresponding to the spatial comoving scales $(135 \pm 12) h^{-1}$ Mpc and $(101 \pm 6)h^{-1}$ Mpc. The latter one is dominant and consistent with the characteristic scale of the baryon acoustic oscillations. We analyse also the radial distributions of two other selected LRG samples: $\sim$ 33,400 bright LRGs ($-23.2 < M \leq -21.8$) and $\sim$ 60,300 all LRGs within a rectangle region on the sky, and show differences of the quasi-periodical features characteristic for different samples. Being confirmed the results would allow to give preference of the spatial against temporal models which could explain the quasi-periodicities discussed here. As a caveat we show that estimations of the significance levels of the peaks strongly depend on a smoothed radial function (trend) as well as characteristics of random fluctuations.

We extend the existing analytical model of reionization by Furlanetto et al. (2004) to include the biasing of reionization sources and additional absorption by Lyman Limit systems. Our model is, by construction, consistent with the observed evolution of the galaxy luminosity function at z<8 and with the observed evolution of Ly-{\alpha} forest at z<6. We also find that, for a wide range of values for the relative escape fraction that we consider reasonable, and which are consistent with the observational constraints on the relative escape fraction from lower redshifts, our reionization model is consistent with the WMAP constraint on the Thompson optical depth and with the SPT and EDGES constraints on the duration of reionization. We, therefore, conclude that it is possible to develop physically realistic models of reionization that are consistent with all existing observational constraints.