-
Influence of turbulence on Lyman-alpha scattering
Authors:
Vadim R. Munirov,
Alexander A. Kaurov
Abstract:
We develop a Monte Carlo radiative transfer code to study the effect of turbulence with a finite correlation length on scattering of Lyman-alpha (Ly$α$) photons propagating through neutral atomic hydrogen gas. We investigate how the effective mean free path, the emergent spectrum, and the average number of scatterings that Ly$α$ photons experience change in the presence of turbulence. We find that…
▽ More
We develop a Monte Carlo radiative transfer code to study the effect of turbulence with a finite correlation length on scattering of Lyman-alpha (Ly$α$) photons propagating through neutral atomic hydrogen gas. We investigate how the effective mean free path, the emergent spectrum, and the average number of scatterings that Ly$α$ photons experience change in the presence of turbulence. We find that the correlation length is an important and sensitive parameter that has an influence on physically relevant properties of Ly$α$ radiative transfer. In particular, it can significantly, by orders of magnitude, reduce the number of scattering events that the average Ly$α$ photon undergoes before it escapes the turbulent cloud.
△ Less
Submitted 3 May, 2023; v1 submitted 27 August, 2022;
originally announced August 2022.
-
Detecting the non-Gaussianity of the 21-cm signal during reionisation with the Wavelet Scattering Transform
Authors:
Bradley Greig,
Yuan-Sen Ting,
Alexander A. Kaurov
Abstract:
Detecting the 21-cm hyperfine transition from neutral hydrogen in the intergalactic medium is our best probe for understanding the astrophysical processes driving the Epoch of Reionisation (EoR). The primary means for a detection of this 21-cm signal is through a statistical measurement of the spatial fluctuations using the 21-cm power spectrum (PS). However, the 21-cm signal is non-Gaussian meani…
▽ More
Detecting the 21-cm hyperfine transition from neutral hydrogen in the intergalactic medium is our best probe for understanding the astrophysical processes driving the Epoch of Reionisation (EoR). The primary means for a detection of this 21-cm signal is through a statistical measurement of the spatial fluctuations using the 21-cm power spectrum (PS). However, the 21-cm signal is non-Gaussian meaning the PS, which only measures the Gaussian fluctuations, is sub-optimal for characterising all of the available information. The upcoming Square Kilometre Array (SKA) will perform a deep, 1000 hr observation over 100 deg$.^{2}$ specifically designed to recover direct images of the 21-cm signal. In this work, we use the Wavelet Scattering Transform (WST) to extract the non-Gaussian information directly from these two-dimensional images of the 21-cm signal. The key advantage of the WST is its stability with respect to statistical noise for measuring non-Gaussian information, unlike the bispectrum whose statistical noise diverges. We introduce a novel method to isolate this non-Gaussian information from mock 21-cm images and demonstrate its detection at 150 (177)~MHz ($z\sim8.5$ and $\sim7$) for a fiducial model with signal-to-noise of $\sim$5~(8) assuming perfect foreground removal and $\sim2$~(3) assuming foreground wedge avoidance.
△ Less
Submitted 19 July, 2022;
originally announced July 2022.
-
Exploring the cosmic 21-cm signal from the Epoch of Reionisation using the Wavelet Scattering Transform
Authors:
Bradley Greig,
Yuan-Sen Ting,
Alexander A. Kaurov
Abstract:
Detecting the cosmic 21-cm signal during the Epoch of Reionisation and Cosmic Dawn will reveal insights into the properties of the first galaxies and advance cosmological parameter estimation. Until recently, the primary focus for astrophysical parameter inference from the 21-cm signal centred on the power spectrum (PS). However, the cosmic 21-cm signal is highly non-Gaussian rendering the PS sub-…
▽ More
Detecting the cosmic 21-cm signal during the Epoch of Reionisation and Cosmic Dawn will reveal insights into the properties of the first galaxies and advance cosmological parameter estimation. Until recently, the primary focus for astrophysical parameter inference from the 21-cm signal centred on the power spectrum (PS). However, the cosmic 21-cm signal is highly non-Gaussian rendering the PS sub-optimal for characterising the cosmic signal. In this work, we introduce a new technique to analyse the non-Gaussian information in images of the 21-cm signal called the Wavelet Scattering Transform (WST). This approach closely mirrors that of convolutional neural networks with the added advantage of not requiring tuning or training of a neural network. Instead, it compresses the 2D spatial information into a set of coefficients making it easier to interpret while also providing a robust statistical description of the non-Gaussian information contained in the cosmic 21-cm signal. First, we explore the application of the WST to mock 21-cm images to gain valuable physical insights by comparing to the known behaviour from the 21-cm PS. Then we quantitatively explore the WST applied to the 21-cm signal by extracting astrophysical parameter constraints using Fisher Matrices from a realistic 1000 hr mock observation with the Square Kilometre Array. We find that: (i) the WST applied only to 2D images can outperform the 3D spherically averaged 21-cm PS, (ii) the excision of foreground contaminated modes can degrade the constraining power by a factor of ~1.5-2 with the WST and (iii) higher cadences between the 21-cm images can further improve the constraining power.
△ Less
Submitted 5 April, 2022;
originally announced April 2022.
-
On the optical properties of Resonant Drag Instabilities: Variability of Asymptotic Giant Branch and R Coronae Borealis stars
Authors:
Ulrich P. Steinwandel,
Alexander A. Kaurov,
Philip F. Hopkins,
Jonathan Squire
Abstract:
In dusty cool-star outflow or ejection events around AGB or RCB-like stars, dust is accelerated by radiation from the star and coupled to the gas via collisional drag forces. But it has recently been shown that such dust-gas mixtures are unstable to a super-class of instabilities called the resonant drag instabilities (RDIs), which promote dust clustering. We therefore consider idealized simulatio…
▽ More
In dusty cool-star outflow or ejection events around AGB or RCB-like stars, dust is accelerated by radiation from the star and coupled to the gas via collisional drag forces. But it has recently been shown that such dust-gas mixtures are unstable to a super-class of instabilities called the resonant drag instabilities (RDIs), which promote dust clustering. We therefore consider idealized simulations of the RDIs operating on a spectrum of dust grain sizes subject to radiative acceleration (allowing for different grain optical properties), coupled to the gas with a realistic drag law, including or excluding the effects of magnetic fields and charged grains, and calculate for the first time how the RDIs could contribute to observed variability. We show that the RDIs naturally produce significant variations ($\sim 10-20\%$ $1σ$-level) in the extinction, corresponding to $\sim 0.1-1\,$mag level in the stellar types above, on timescales of order months to a year. The fluctuations are surprisingly robust to finite source-size effects as they are dominated by large-scale modes, which also means their spatial structure could be resolved in some nearby systems. We also quantify how this produces variations in the line-of-sight grain size-distribution. All of these variations are similar to those observed, suggesting that the RDIs may play a key role driving observed variability in dust extinction within dusty outflow/ejection events around cool stars. We further propose that the measured variations in grain sizes could directly be used to identify the presence of the RDIs in close by systems with observations.
△ Less
Submitted 17 November, 2021;
originally announced November 2021.
-
Simulations of $^{60}$Fe entrained in ejecta from a near-Earth supernova: Effects of observer motion
Authors:
Evgenii Chaikin,
Alexander A. Kaurov,
Brian D. Fields,
Camila A. Correa
Abstract:
Recent studies have shown that live (not decayed) radioactive $^{60}$Fe is present in deep-ocean samples, Antarctic snow, lunar regolith, and cosmic rays. $^{60}$Fe represents supernova (SN) ejecta deposited in the Solar System around $3 \, \rm Myr$ ago, and recently an earlier pulse $\approx 7 \ \rm Myr $ ago has been found. These data point to one or multiple near-Earth SN explosions that presum…
▽ More
Recent studies have shown that live (not decayed) radioactive $^{60}$Fe is present in deep-ocean samples, Antarctic snow, lunar regolith, and cosmic rays. $^{60}$Fe represents supernova (SN) ejecta deposited in the Solar System around $3 \, \rm Myr$ ago, and recently an earlier pulse $\approx 7 \ \rm Myr $ ago has been found. These data point to one or multiple near-Earth SN explosions that presumably participated in the formation of the Local Bubble. We explore this theory using 3D high-resolution smooth-particle hydrodynamical simulations of isolated supernovae with ejecta tracers in a uniform interstellar medium (ISM). The simulation allows us to trace the supernova ejecta in gas form and those eject in dust grains that are entrained with the gas. We consider two cases of diffused ejecta: when the ejecta are well-mixed in the shock and when they are not. In the latter case, we find that these ejecta remain far behind the forward shock, limiting the distance to which entrained ejecta can be delivered to $\approx 100$ pc in an ISM with $n_\mathrm{H}=0.1\; \rm cm^{-3}$ mean hydrogen density. We show that the intensity and the duration of $^{60}$Fe accretion depend on the ISM density and the trajectory of the Solar System. Furthermore, we show the possibility of reproducing the two observed peaks in $^{60}$Fe concentration with this model by assuming two linear trajectories for the Solar System with $30$-km s$^{-1}$ velocity. The fact that we can reproduce the two observed peaks further supports the theory that the $^{60}$Fe signal was originated from near-Earth SNe.
△ Less
Submitted 6 February, 2022; v1 submitted 23 September, 2021;
originally announced September 2021.
-
Asymmetric Surface Brightness Structure of Caustic Crossing Arc in SDSS J1226+2152: A Case for Dark Matter Substructure
Authors:
Liang Dai,
Alexander A. Kaurov,
Keren Sharon,
Michael K. Florian,
Jordi Miralda-Escudé,
Tejaswi Venumadhav,
Brenda Frye,
Jane R. Rigby,
Matthew Bayliss
Abstract:
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 t…
▽ More
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.
△ Less
Submitted 12 May, 2020; v1 submitted 1 January, 2020;
originally announced January 2020.
-
Research and Development for HI Intensity Mapping
Authors:
Zeeshan Ahmed,
David Alonso,
Mustafa A. Amin,
Réza Ansari,
Evan J. Arena,
Kevin Bandura,
Adam Beardsley,
Philip Bull,
Emanuele Castorina,
Tzu-Ching Chang,
Romeel Davé,
Joshua S. Dillon,
Alexander van Engelen,
Aaron Ewall-Wice,
Simone Ferraro,
Simon Foreman,
Josef Frisch,
Daniel Green,
Gilbert Holder,
Daniel Jacobs,
Dionysios Karagiannis,
Alexander A. Kaurov,
Lloyd Knox,
Emily Kuhn,
Adrian Liu
, et al. (29 additional authors not shown)
Abstract:
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 an…
▽ More
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.
△ Less
Submitted 29 July, 2019;
originally announced July 2019.
-
Probing the Time Domain with High Spatial Resolution
Authors:
J. P. Blakeslee,
S. A. Rodney,
J. M. Lotz,
G. Sivo,
S. Sivanandam,
M. Andersen,
R. Carrasco,
L. Ferrarese,
R. J. Foley,
S. Goodsell,
P. Hirst,
J. B. Jensen,
P. L. Kelly,
A. A. Kaurov,
M. Lemoine-Busserolle,
B. W. Miller,
J. O'Meara,
H. Roe,
M. E. Schwamb,
J. Scharwächter
Abstract:
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…
▽ More
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.
△ Less
Submitted 19 March, 2019;
originally announced March 2019.
-
Highly Magnified Stars in Lensing Clusters: New Evidence in a Galaxy Lensed by MACS J0416.1-2403
Authors:
Alexander A. Kaurov,
Liang Dai,
Tejaswi Venumadhav,
Jordi Miralda-Escudé,
Brenda Frye
Abstract:
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,…
▽ More
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 $μ\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.
△ Less
Submitted 26 February, 2019;
originally announced February 2019.
-
Non-perturbative probability distribution function for cosmological counts in cells
Authors:
Mikhail M. Ivanov,
Alexander A. Kaurov,
Sergey Sibiryakov
Abstract:
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 s…
▽ More
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.
△ Less
Submitted 18 March, 2019; v1 submitted 19 November, 2018;
originally announced November 2018.
-
Implication of the Shape of the EDGES Signal for the 21 cm Power Spectrum
Authors:
Alexander A. Kaurov,
Tejaswi Venumadhav,
Liang Dai,
Matias Zaldarriaga
Abstract:
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 ph…
▽ More
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-$α$ 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.
△ Less
Submitted 8 May, 2018;
originally announced May 2018.
-
Probing Dark Matter Subhalos in Galaxy Clusters Using Highly Magnified Stars
Authors:
Liang Dai,
Tejaswi Venumadhav,
Alexander A. Kaurov,
Jordi Miralda-Escudé
Abstract:
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…
▽ More
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\,μ{\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.
△ Less
Submitted 13 October, 2018; v1 submitted 9 April, 2018;
originally announced April 2018.
-
Observing galaxy mergers at the epoch of reionization
Authors:
Evgenii A. Chaikin,
Nadezda V. Tyulneva,
Alexander A. Kaurov
Abstract:
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,…
▽ More
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.
△ Less
Submitted 11 December, 2017;
originally announced December 2017.
-
Stochasticity in the 21cm power spectrum at the epoch of reionization and cosmic dawn
Authors:
Alexander A. Kaurov
Abstract:
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…
▽ More
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.
△ Less
Submitted 13 September, 2017;
originally announced September 2017.
-
On improving analytical models of cosmic reionization for matching numerical simulation
Authors:
Alexander A. Kaurov
Abstract:
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 fram…
▽ More
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.
△ Less
Submitted 8 September, 2016; v1 submitted 3 December, 2015;
originally announced December 2015.
-
The Effects of Dark Matter Annihilation on Cosmic Reionization
Authors:
Alexander A. Kaurov,
Dan Hooper,
Nickolay Y. Gnedin
Abstract:
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…
▽ More
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, $τ$. 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.
△ Less
Submitted 1 December, 2015;
originally announced December 2015.
-
Cosmic Reionization On Computers. Mean and Fluctuating Redshifted 21 cm Signal
Authors:
Alexander A. Kaurov,
Nickolay Y. Gnedin
Abstract:
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ΔT_B\rangle\sim-20\rm{mK}$, in agreement with prior simulatio…
▽ More
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Δ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.
△ Less
Submitted 18 April, 2016; v1 submitted 29 October, 2015;
originally announced October 2015.
-
Energy dissipation of energetic electrons in the inhomogeneous intergalactic medium during the epoch of reionization
Authors:
Alexander A. Kaurov
Abstract:
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 a…
▽ More
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.
△ Less
Submitted 18 April, 2016; v1 submitted 28 October, 2015;
originally announced October 2015.
-
Central compact objects in Kes 79 and RCW 103 as "hidden" magnetars with crustal activity
Authors:
S. B. Popov,
A. A. Kaurov,
A. D. Kaminker
Abstract:
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…
▽ More
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.
△ Less
Submitted 13 April, 2015;
originally announced April 2015.
-
Cosmic Reionization On Computers III. The Clumping Factor
Authors:
Alexander A. Kaurov,
Nickolay Y. Gnedin
Abstract:
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 s…
▽ More
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.
△ Less
Submitted 13 July, 2015; v1 submitted 17 December, 2014;
originally announced December 2014.
-
Thermal emission of neutron stars with internal heaters
Authors:
A. D. Kaminker,
A. A. Kaurov,
A. Y. Potekhin,
D. G. Yakovlev
Abstract:
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…
▽ More
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.
△ Less
Submitted 3 June, 2014;
originally announced June 2014.
-
Cosmic Reionization On Computers II. Reionization History and Its Back Reaction on Early Galaxies
Authors:
Nickolay Y. Gnedin,
Alexander A. Kaurov
Abstract:
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 independen…
▽ More
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.
△ Less
Submitted 4 August, 2014; v1 submitted 17 March, 2014;
originally announced March 2014.
-
Recombination clumping factor during cosmic reionization
Authors:
Alexander A. Kaurov,
Nickolay Y. Gnedin
Abstract:
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 pre…
▽ More
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.
△ Less
Submitted 14 April, 2014; v1 submitted 11 November, 2013;
originally announced November 2013.
-
Quasi-periodical features in the distribution of Luminous Red Galaxies
Authors:
Andrei I. Ryabinkov,
Alexander A. Kaurov,
Alexander D. Kaminker
Abstract:
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 $η$, dimensionless line-of-sight comoving distance calculated for the $Λ$C…
▽ More
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 $η$, dimensionless line-of-sight comoving distance calculated for the $Λ$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.
△ Less
Submitted 13 December, 2012;
originally announced December 2012.
-
Effect of halo bias and Lyman Limit Systems on the history of cosmic reionization
Authors:
Alexander A. Kaurov,
Nickolay Y. Gnedin
Abstract:
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-α forest at z<6. We also find that, for a wide range of values for the…
▽ More
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-α 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.
△ Less
Submitted 6 December, 2012;
originally announced December 2012.