Source author record

Jack O. Burns

Jack O. Burns appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

37works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

37 published item(s)

preprint2022arXiv

Constraining Warm Dark Matter and Pop III stars with the Global 21-cm Signal

Upcoming ground and space-based experiments may have sufficient accuracy to place significant constraints upon high-redshift star formation, Reionization, and dark matter (DM) using the global 21-cm signal of the intergalactic medium. In the early universe, when the relative abundance of low-mass DM halos is important, measuring the global signal would place constraints on the damping of structure formation caused by DM having a higher relic velocity (warm dark matter, or WDM) than in cold dark matter (CDM). Such damping, however, can be mimicked by altering the star formation efficiency (SFE) and difficult to detect because of the presence of Pop III stars with unknown properties. We study these various cases and their degeneracies with the WDM mass parameter $m_X$ using a Fisher matrix analysis. We study the $m_X = 7$ keV case and a star-formation model that parametrizes the SFE as a strong function of halo mass and include several variations of this model along with three different input noise levels for the likelihood; we also use a minimum halo virial temperature for collapse near the molecular cooling threshold. We find that when the likelihood includes only Pop II stars, $m_X$ is constrained to an uncertainty of $\sim 0.4$ keV for all models and noise levels at 68$\%$ CI. When the likelihood includes weak Pop III stars, $m_X \sim 0.3$ keV, and if Pop III star formation is relatively efficient, $m_X \sim 0.1$ keV uncertainty, with tight Pop III star-formation parameter constraints. Our results show that the global 21-cm signal is a promising test-bed for WDM models, even in the presence of strong degeneracies with astrophysical parameters.

preprint2022arXiv

Virtual Reality Digital Twin and Environment for Troubleshooting Lunar-based Infrastructure Assembly Failures

Humans and robots will need to collaborate in order to create a sustainable human lunar presence by the end of the 2020s. This includes cases in which a human will be required to teleoperate an autonomous rover that has encountered an instrument assembly failure. To aid teleoperators in the troubleshooting process, we propose a virtual reality digital twin placed in a simulated environment. Here, the operator can virtually interact with a digital version of the rover and mechanical arm that uses the same controls and kinematic model. The user can also adopt the egocentric (a first person view through using stereoscopic passthrough) and exocentric (a third person view where the operator can virtually walk around the environment and rover as if they were on site) view. We also discuss our metrics for evaluating the differences between our digital and physical robot, as well as the experimental concept based on real and applicable missions, and future work that would compare our platform to traditional troubleshooting methods.

preprint2021arXiv

Low Radio Frequency Observations from the Moon Enabled by NASA Landed Payload Missions

A new era of exploration of the low radio frequency Universe from the Moon will soon be underway with landed payload missions facilitated by NASA's Commercial Lunar Payload Services (CLPS) program. CLPS landers are scheduled to deliver two radio science experiments, ROLSES to the nearside and LuSEE to the farside, beginning in 2021. These instruments would be pathfinders for a 10-km diameter interferometric array, FARSIDE, composed of 128 pairs of dipole antennas proposed to be delivered to the lunar surface later in the decade. ROLSES and LuSEE, operating at frequencies from 100 kHz to a few tens of MHz, will investigate the plasma environment above the lunar surface and measure the fidelity of radio spectra on the surface. Both use electrically-short, spiral-tube deployable antennas and radio spectrometers based upon previous flight models. ROLSES will measure the photoelectron sheath density to better understand the charging of the lunar surface via photoionization and impacts from the solar wind, charged dust, and current anthropogenic radio frequency interference. LuSEE will measure the local magnetic field and exo-ionospheric density, interplanetary radio bursts, Jovian and terrestrial natural radio emission, and the galactic synchrotron spectrum. FARSIDE, and its precursor risk-reduction six antenna-node array PRIME, would be the first radio interferometers on the Moon. FARSIDE would break new ground by imaging radio emission from Coronal Mass Ejections (CME) beyond 2 solar radii, monitor auroral radiation from the B-fields of Uranus and Neptune (not observed since Voyager), and detect radio emission from stellar CMEs and the magnetic fields of nearby potentially habitable exoplanets.

preprint2020arXiv

Characterizing the Radio Quiet Region Behind the Lunar Farside for Low Radio Frequency Experiments

Low radio frequency experiments performed on Earth are contaminated by both ionospheric effects and radio frequency interference (RFI) from Earth-based sources. The lunar farside provides a unique environment above the ionosphere where RFI is heavily attenuated by the presence of the Moon. We present electrodynamics simulations of the propagation of radio waves around and through the Moon in order to characterize the level of attenuation on the farside. The simulations are performed for a range of frequencies up to 100 kHz, assuming a spherical lunar shape with an average, constant density. Additionally, we investigate the role of the topography and density profile of the Moon in the propagation of radio waves and find only small effects on the intensity of RFI. Due to the computational demands of performing simulations at higher frequencies, we propose a model for extrapolating the width of the quiet region above 100 kHz that also takes into account height above the lunar surface as well as the intensity threshold chosen to define the quiet region. This model, which we make publicly available through a Python package, allows the size of the radio quiet region to be easily calculated both in orbit or on the surface, making it directly applicable for lunar satellites as well as surface missions.

preprint2020arXiv

Diffuse radio emission in the galaxy cluster SPT-CL J2031-4037: a steep spectrum intermediate radio halo?

The advent of sensitive low frequency radio observations has revealed a number of diffuse radio objects with peculiar properties that are challenging our understanding about the physics of the intracluster medium. Here, we report the discovery of a steep spectrum radio halo surrounding the central Brightest Cluster Galaxy (BCG) in the galaxy cluster SPT-CL J2031-4037. This cluster is morphologically disturbed yet has a weak cool core, an example of cool core/non-cool core transition system, which harbours a radio halo of $\sim 0.7$ Mpc in size. The halo emission detected at 1.7 GHz is less extended compared to that in the 325 MHz observation, and the spectral index of the part of the halo visible at 325 MHz to 1.7 GHz frequencies was found to be $-1.35 \pm 0.07$. Also, $P_{1.4\ \mathrm{GHz}}$ was found to be $0.77 \times 10^{24}$ W Hz$^{-1}$ which falls in the region where radio mini-halos, halo upper limits and ultra-steep spectrum (USS) halos are found in the $P_{1.4\ \mathrm{GHz}} - L_\mathrm{X}$ plane. Additionally, simulations presented in the paper provide support to the scenario of the steep spectrum. The diffuse radio emission found in this cluster may be a steep spectrum "intermediate" or "hybrid" radio halo which is transitioning into a mini-halo.

preprint2020arXiv

Formulating and critically examining the assumptions of global 21-cm signal analyses: How to avoid the false troughs that can appear in single spectrum fits

The assumptions inherent to global 21-cm signal analyses are rarely delineated. In this paper, we formulate a general list of suppositions underlying a given claimed detection of the global 21-cm signal. Then, we specify the form of these assumptions for two different analyses: 1) the one performed by the EDGES team showing an absorption trough in brightness temperature that they modeled separately from the sky foreground and 2) a new, so-called Minimum Assumption Analysis (MAA), that makes the most conservative assumptions possible for the signal. We show fits using the EDGES analysis on various beam-weighted foreground simulations from the EDGES latitude with no signal added. Depending on the beam used, these simulations produce large false troughs due to the invalidity of the foreground model to describe the combination of beam chromaticity and the shape of the Galactic plane in the sky, the residuals of which are captured by the ad hoc flattened Gaussian signal model. On the other hand, the MAA provides robust fits by including many spectra at different time bins and allowing any possible 21-cm spectrum to be modeled exactly. We present uncertainty levels and example signal reconstructions found with the MAA for different numbers of time bins. With enough time bins, one can determine the true 21-cm signal with the MAA to $<10$ times the noise level.

preprint2020arXiv

Global 21-cm signal extraction from foreground and instrumental effects III: Utilizing drift-scan time dependence and full Stokes measurements

When using valid foreground and signal models, the uncertainties on extracted signals in global 21-cm signal experiments depend principally on the overlap between signal and foreground models. In this paper, we investigate two strategies for decreasing this overlap: (i) utilizing time dependence by fitting multiple drift-scan spectra simultaneously and (ii) measuring all four Stokes parameters instead of only the total power, Stokes I. Although measuring polarization requires different instruments than are used in most existing experiments, all existing experiments can utilize drift-scan measurements merely by averaging their data differently. In order to evaluate the increase in constraining power from using these two techniques, we define a method for connecting Root-Mean-Square (RMS) uncertainties to probabilistic confidence levels. Employing simulations, we find that fitting only one total power spectrum leads to RMS uncertainties at the few K level, while fitting multiple time-binned, drift-scan spectra yields uncertainties at the $\lesssim 10$ mK level. This significant improvement only appears if the spectra are modeled with one set of basis vectors, instead of using multiple sets of basis vectors that independently model each spectrum. Assuming that they are simulated accurately, measuring all four Stokes parameters also leads to lower uncertainties. These two strategies can be employed simultaneously and fitting multiple time bins of all four Stokes parameters yields the best precision measurements of the 21-cm signal, approaching the noise level in the data.

preprint2020arXiv

Modelling the Galactic Foreground and Beam Chromaticity for Global 21-cm Cosmology

In order to characterize and model the beam-weighted foreground for global 21-cm signal experiments, we present a methodology for generating basis eigenvectors that combines analytical and observational models of both the galactic spectral index and sky brightness temperature with simulations of beams having various angular and spectral dependencies and pointings. Each combination creates a unique beam-weighted foreground. By generating eigenvectors to fit each foreground model using Singular Value Decomposition (SVD), we examine the effects of varying the components of the beam-weighted foreground. We find that the eigenvectors for modelling an achromatic, isotropic beam -- the ideal case -- are nearly identical regardless of the unweighted foreground model used, and are practicably indistinguishable from polynomial-based models. When anisotropic, chromatic beams weight the foreground, however, a coupling is introduced between the spatial and spectral structure of the foreground which distorts the eigenvectors away from the polynomial models and induces a dependence of the basis upon the exact features of the beam (chromaticity, pattern, pointing) and foreground (spectral index, sky brightness temperature map). We find that the beam has a greater impact upon the eigenvectors than foreground models. Any model which does not account for its distortion may produce RMS uncertainties on the order of $\sim 10$ - $10^3$ Kelvin for six-parameter, single spectrum fits. If the beam is incorporated directly using SVD and training sets, however, the resultant eigenvectors yield milli-Kelvin level uncertainties. Given a sufficiently detailed description of the sky, our methodology can be applied to any particular experiment with a suitably characterized beam for the purpose of generating accurate beam-weighted foreground models.

preprint2020arXiv

The Case for Probe-class NASA Astrophysics Missions

Astrophysics spans an enormous range of questions on scales from individual planets to the entire cosmos. To address the richness of 21st century astrophysics requires a corresponding richness of telescopes spanning all bands and all messengers. Much scientific benefit comes from having the multi-wavelength capability available at the same time. Most of these bands,or measurement sensitivities, require space-based missions. Historically, NASA has addressed this need for breadth with a small number of flagship-class missions and a larger number of Explorer missions. While the Explorer program continues to flourish, there is a large gap between Explorers and strategic missions. A fortunate combination of new astrophysics technologies with new, high capacity, low dollar-per-kg to orbit launchers, and new satellite buses allow for cheaper missions with capabilities approaching strategic mission levels. NASA has recognized these developments by calling for Probe-class mission ideas for mission studies, spanning most of the electromagnetic spectrum from GeV gamma-rays to the far infrared, and the new messengers of neutrinos and ultra-high energy cosmic rays. The key insight from the Probes exercise is that order-of-magnitude advances in science performance metrics are possible across the board for initial total cost estimates in the range 500M-1B dollars.

preprint2020arXiv

Transformative Science from the Lunar Farside: Observations of the Dark Ages and Exoplanetary Systems at Low Radio Frequencies

The farside of the Moon is a pristine, quiet platform to conduct low radio frequency observations of the early Universe's Dark Ages, as well as space weather and magnetospheres associated with habitable exoplanets. In this paper, the astrophysics associated with NASA-funded concept studies will be described including a lunar-orbiting spacecraft, DAPPER, that will measure the 21 cm global spectrum at redshifts 40-80, and an array of low frequency dipoles on the lunar farside surface, FARSIDE, that would detect exoplanet magnetic fields. DAPPER observations (17-38 MHz), using a single cross-dipole antenna, will measure the amplitude of the 21 cm spectrum to the level required to distinguish the standard ΛCDM cosmological model from those of additional cooling models possibly produced by exotic physics such as dark matter interactions. FARSIDE has a notional architecture consisting of 128 dipole antennas deployed across a 10 km area by a rover. FARSIDE would image the entire sky each minute in 1400 channels over 0.1-40 MHz. This would enable monitoring of the nearest stellar systems for the radio signatures of coronal mass ejections and energetic particle events, and would also detect the magnetospheres of the nearest candidate habitable exoplanets. In addition, FARSIDE would determine the Dark Ages global 21 cm signal at yet lower frequencies and provide a pathfinder for power spectrum measurements.

preprint2019arXiv

Global 21-cm Signal Extraction from Foreground and Instrumental Effects II: Efficient and Self-Consistent Technique for Constraining Nonlinear Signal Models

We present the completion of a data analysis pipeline that self-consistently separates global 21-cm signals from large systematics using a pattern recognition technique. In the first paper of this series, we obtain optimal basis vectors from signal and foreground training sets to linearly fit both components with the minimal number of terms that best extracts the signal given its overlap with the foreground. In this second paper, we utilize the spectral constraints derived in the first paper to calculate the full posterior probability distribution of any signal parameter space of choice. The spectral fit provides the starting point for a Markov Chain Monte Carlo (MCMC) engine that samples the signal without traversing the foreground parameter space. At each MCMC step, we marginalize over the weights of all linear foreground modes and suppress those with unimportant variations by applying priors gleaned from the training set. This method drastically reduces the number of MCMC parameters, augmenting the efficiency of exploration, circumvents the need for selecting a minimal number of foreground modes, and allows the complexity of the foreground model to be greatly increased to simultaneously describe many observed spectra without requiring extra MCMC parameters. Using two nonlinear signal models, one based on EDGES observations and the other on phenomenological frequencies and temperatures of theoretically expected extrema, we demonstrate the success of this methodology by recovering the input parameters from multiple randomly simulated signals at low radio frequencies (10-200 MHz), while rigorously accounting for realistically modeled beam-weighted foregrounds.

preprint2019arXiv

Probing the origin of diffuse radio emission in the cool-core of the Phoenix galaxy cluster

Cool core galaxy clusters are considered to be dynamically relaxed clusters with regular morphology and highly X-ray luminous central region. However, cool core clusters can also be sites for merging events that exhibit cold fronts in X-ray and mini-halos in radio. We present recent radio/X-ray observations of the Phoenix Cluster or SPT-CL J2344-4243 at the redshift of $z=0.596$. Using archival {\it Chandra} X-ray observations, we detect spiraling cool gas around the cluster core as well as discover two cold fronts near the core. It is perhaps the most distant galaxy cluster to date known to host cold fronts. Also, we present JVLA\footnote{Jansky Very Large Array\\ \url{https://science.nrao.edu/facilities/vla}} 1.52 GHz observations of the minihalo, previously discovered at 610 MHz with GMRT\footnote{Giant Metrewave Radio Telescope\\ \url{http://www.gmrt.ncra.tifr.res.in}} observations in the center of the Phoenix galaxy cluster. The minihalo flux density at 1.52 GHz is $9.65 \pm 0.97$ mJy with the spectral index between 610 MHz and 1.52 GHz being $-0.98 \pm 0.16$\footnote{$S_ν = ν^α$ where $S_ν$}. A possible origin of these radio sources is turbulence induced by sloshing of the gas in the cluster core.

preprint2016arXiv

Effects Of The Ionosphere On Ground-Based Detection Of The Global 21 CM Signal From The Cosmic Dawn And The Dark Ages

Detection of the global HI 21 cm signal from Cosmic Dawn and Epoch of Reionization is the key science driver for several ongoing ground-based and future ground/space-based experiments. The crucial spectral features in the global 21 cm signal (turning points) occur at low radio frequencies <100 MHz. In addition to the human-generated RFI, Earth's ionosphere drastically corrupts low-frequency radio observations from the ground. In this paper, we examine the effects of time-varying ionospheric refraction, absorption and thermal emission at these low radio frequencies and their combined effect on any ground-based global 21 cm experiment. It should be noted that this is the first study of the effect of a dynamic ionosphere on global 21 cm experiments. The fluctuations in the ionosphere are influenced by solar activity with flicker noise characteristics. The same characteristics are reflected in the ionospheric corruption to any radio signal passing through the ionosphere. As a result, any ground based observations of the faint global 21 cm signal are corrupted by flicker noise (or "$1/f$" noise, where "$f$" is the dynamical frequency) which scales as $ν^{-2}$ (where $ν$ is the frequency of observation) in the presence of a bright galactic foreground ($\propto ν^{-s}$, where $s$ is radio spectral index). Hence, the calibration of the ionosphere for any such experiment is critical. Any attempt to calibrate the ionospheric effects will be subject to the inaccuracies in the current ionospheric measurements using GPS ionospheric measurements, riometer measurements, ionospheric soundings, etc. Even considering an optimistic improvement in the accuracy of GPS-TEC (Total Electron Content) measurements, we conclude that the detection of the global 21 cm signal below 100 MHz is best done from above the Earth's atmosphere in orbit of the Moon.

preprint2015arXiv

A New Lecture-Tutorial for Teaching about Molecular Excitations and Synchrotron Radiation

Light and spectroscopy are among the most important and frequently taught topics in introductory, college-level, general education astronomy courses. This is due to the fact that the vast majority of observational data studied by astronomers arrives at Earth in the form of light. While there are many processes by which matter can emit and absorb light, Astro 101 courses typically limit their instruction to the Bohr model of the atom and electron energy level transitions. In this paper, we report on the development of a new Lecture-Tutorial to help students learn about other processes that are responsible for the emission and absorption of light, namely molecular rotations, molecular vibrations, and the acceleration of charged particles by magnetic fields.

preprint2015arXiv

Interpreting the Global 21-cm Signal from High Redshifts. II. Parameter Estimation for Models of Galaxy Formation

Following our previous work, which related generic features in the sky-averaged (global) 21-cm signal to properties of the intergalactic medium, we now investigate the prospects for constraining a simple galaxy formation model with current and near-future experiments. Markov-Chain Monte Carlo fits to our synthetic dataset, which includes a realistic galactic foreground, a plausible model for the signal, and noise consistent with 100 hours of integration by an ideal instrument, suggest that a simple four-parameter model that links the production rate of Lyman-$α$, Lyman-continuum, and X-ray photons to the growth rate of dark matter halos can be well-constrained (to $\sim 0.1$ dex in each dimension) so long as all three spectral features expected to occur between $40 \lesssim ν/ \mathrm{MHz} \lesssim 120$ are detected. Several important conclusions follow naturally from this basic numerical result, namely that measurements of the global 21-cm signal can in principle (i) identify the characteristic halo mass threshold for star formation at all redshifts $z \gtrsim 15$, (ii) extend $z \lesssim 4$ upper limits on the normalization of the X-ray luminosity star-formation rate ($L_X$-SFR) relation out to $z \sim 20$, and (iii) provide joint constraints on stellar spectra and the escape fraction of ionizing radiation at $z \sim 12$. Though our approach is general, the importance of a broad-band measurement renders our findings most relevant to the proposed Dark Ages Radio Explorer, which will have a clean view of the global 21-cm signal from $\sim 40-120$ MHz from its vantage point above the radio-quiet, ionosphere-free lunar far-side.

preprint2015arXiv

Parametrizations of the global 21-cm signal and parameter estimation from single-dipole experiments

One approach to extracting the global 21-cm signal from total-power measurements at low radio frequencies is to parametrize the different contributions to the data and then fit for these parameters. We examine parametrizations of the 21-cm signal itself, and propose one based on modelling the Lyman-alpha background, IGM temperature and hydrogen ionized fraction using tanh functions. This captures the shape of the signal from a physical modelling code better than an earlier parametrization based on interpolating between maxima and minima of the signal, and imposes a greater level of physical plausibility. This allows less biased constraints on the turning points of the signal, even though these are not explicitly fit for. Biases can also be alleviated by discarding information which is less robustly described by the parametrization, for example by ignoring detailed shape information coming from the covariances between turning points or from the high-frequency parts of the signal, or by marginalizing over the high-frequency parts of the signal by fitting a more complex foreground model. The fits are sufficiently accurate to be usable for experiments gathering 1000 h of data, though in this case it may be important to choose observing windows which do not include the brightest areas of the foregrounds. Our assumption of pointed, single-antenna observations and very broad-band fitting makes these results particularly applicable to experiments such as the Dark Ages Radio Explorer, which would study the global 21-cm signal from the clean environment of a low lunar orbit, taking data from the far side.

preprint2015arXiv

Trajectory Design From GTO To Lunar Equatorial Orbit For The Dark Ages Radio Explorer (DARE) Spacecraft

The trajectory design for the Dark Ages Radio Explorer (DARE) mission con-cept involves launching the DARE spacecraft into a geosynchronous transfer orbit (GTO) as a secondary payload. From GTO, the spacecraft then transfers to a lunar orbit that is stable (i.e., no station-keeping maneuvers are required with minimum perilune altitude always above 40 km) and allows for more than 1,000 cumulative hours for science measurements in the radio-quiet region located on the lunar farside.

preprint2014arXiv

Dark Ages Radio Explorer Mission: Probing the Cosmic Dawn

The period between the creation of the cosmic microwave background at a redshift of ~1000 and the formation of the first stars and black holes that re-ionize the intergalactic medium at redshifts of 10-20 is currently unobservable. The baryonic component of the universe during this period is almost entirely neutral hydrogen, which falls into local regions of higher dark matter density. This seeds the formation of large-scale structures including the cosmic web that we see today in the filamentary distribution of galaxies and clusters of galaxies. The only detectable signal from these dark ages is the 21-cm spectral line of hydrogen, redshifted down to frequencies of approximately 10-100 MHz. Space-based observations of this signal will allow us to determine the formation epoch and physics of the first sources of ionizing radiation, and potentially detect evidence for the decay of dark matter particles. JPL is developing deployable low frequency antenna and receiver prototypes to enable both all-sky spectral measurements of neutral hydrogen and ultimately to map the spatial distribution of the signal as a function of redshift. Such observations must be done from space because of Earth's ionosphere and ubiquitous radio interference. A specific application of these technologies is the Dark Ages Radio Explorer (DARE) mission. This small Explorer class mission is designed to measure the sky-averaged hydrogen signal from the shielded region above the far side of the Moon. These data will complement ground-based radio observations of the final stages of intergalactic re-ionization at higher frequencies. DARE will also serve as a scientific percursor for space-based interferometry missions to image the distribution of hydrogen during the cosmic dark ages.

preprint2014arXiv

How Much Can We Learn From A Merging Cold Front Cluster? : Insights From X-ray Temperature and Radio Maps of Abell 3667

The galaxy cluster Abell 3667 is an ideal laboratory to study the plasma processes in the intracluster medium (ICM). High resolution Chandra X-ray observations show a cold front in Abell 3667. At radio wavelengths, Abell 3667 reveals a double radio-relic feature in the outskirts of the cluster. These suggest multiple merger events in this cluster. In this paper, we analyze the substantial archival X-ray observations of Abell 3667 from ChandraX-ray Observatory and compare these with existing radio observations as well as state-of-the-art AMR (Adaptive Mesh Refinement) MHD cosmological simulations using Enzo. We have used two temperature map making techniques, Weighted Voronoi Tessellation and Adaptive Circular Binning, to produce the high resolution and largest field-of-view temperature maps of Abell 3667. These high fidelity temperature maps allow us to study the X-ray shocks in the cluster using a new 2-dimensional shock-finding algorithm. We have also estimated the Mach numbers from the shocks inferred from previous ATCA radio observations. The combined shock statistics from the X-ray and radio data are in agreement with the shock statistics in a simulated MHD cluster. We have also studied the profiles of the thermodynamic properties across the cold front using 447 ksec from the combined Chandra observations on Abell 3667. Our results show that the stability of the cold front in Abell 3667 can be attributed to the suppression of the thermal conduction across the cold front by a factor of 100 - 700 compared to the classical Spitzer value.

preprint2013arXiv

Interpreting the Global 21-cm Signal from High Redshifts. I. Model Independent Constraints

The sky-averaged (global) 21-cm signal is a powerful probe of the intergalactic medium (IGM) prior to the completion of reionization. However, it has so far been unclear that even in the best case scenario, in which the signal is accurately extracted from the foregrounds, that it will provide more than crude estimates of when the Universe's first stars and black holes form. In contrast to previous work, which has focused on predicting the 21-cm signatures of the first luminous objects, we investigate an arbitrary realization of the signal, and attempt to translate its features to the physical properties of the IGM. Within a simplified global framework, the 21-cm signal yields quantitative constraints on the Lyman-alpha background intensity, net heat deposition, ionized fraction, and their time derivatives, without invoking models for the astrophysical sources themselves. The 21-cm absorption signal is most easily interpreted, setting strong limits on the heating rate density of the Universe with a measurement of its redshift alone, independent of the ionization history or details of the Lyman-alpha background evolution. In a companion paper we extend these results, focusing on the confidence with which one can infer source emissivities from IGM properties.

preprint2012arXiv

A Lunar L2-Farside Exploration and Science Mission Concept with the Orion Multi-Purpose Crew Vehicle and a Teleoperated Lander/Rover

A novel concept is presented in this paper for a human mission to the lunar L2 (Lagrange) point that would be a proving ground for future exploration missions to deep space while also overseeing scientifically important investigations. In an L2 halo orbit above the lunar farside, the astronauts aboard the Orion Crew Vehicle would travel 15% farther from Earth than did the Apollo astronauts and spend almost three times longer in deep space. Such a mission would serve as a first step beyond low Earth orbit and prove out operational spaceflight capabilities such as life support, communication, high speed re-entry, and radiation protection prior to more difficult human exploration missions. On this proposed mission, the crew would teleoperate landers and rovers on the unexplored lunar farside, which would obtain samples from the geologically interesting farside and deploy a low radio frequency telescope. Sampling the South Pole-Aitken basin, one of the oldest impact basins in the solar system, is a key science objective of the 2011 Planetary Science Decadal Survey. Observations at low radio frequencies to track the effects of the Universe's first stars/galaxies on the intergalactic medium are a priority of the 2010 Astronomy and Astrophysics Decadal Survey. Such telerobotic oversight would also demonstrate capability for human and robotic cooperation on future, more complex deep space missions such as exploring Mars.

preprint2012arXiv

Astrophysics Conducted by the Lunar University Network for Astrophysics Research (LUNAR) and the Center for Lunar Origins (CLOE)

[Abridged] The Moon is a unique platform from and on which to conduct astrophysical measurements. The Lunar University Network for Astrophysics Research (LUNAR) and the Center for Lunar Origins and Evolution (CLOE) teams within the NASA Lunar Science Institute (NLSI) are illustrating how the Moon can be used as a platform to advance important goals in astrophysics. Of relevance to Astrophysics and aligned with NASA strategic goals, all three of the primary research themes articulated by New Worlds, New Horizons in Astronomy & Astrophysics are being addressed by LUNAR and CLOE, namely Probing Cosmic Dawn, Understanding New Worlds, and Physics of the Universe

preprint2012arXiv

Cosmological MHD Simulations of Galaxy Cluster Radio Relics: Insights and Warnings for Observations

Non-thermal radio emission from cosmic ray electrons in the vicinity of merging galaxy clusters is an important tracer of cluster merger activity, and is the result of complex physical processes that involve magnetic fields, particle acceleration, gas dynamics, and radiation. In particular, objects known as radio relics are thought to be the result of shock-accelerated electrons that, when embedded in a magnetic field, emit synchrotron radiation in the radio wavelengths. In order to properly model this emission, we utilize the adaptive mesh refinement simulation of the magnetohydrodynamic evolution of a galaxy cluster from cosmological initial conditions. We locate shock fronts and apply models of cosmic ray electron acceleration that are then input into radio emission models. We have determined the thermodynamic properties of this radio-emitting plasma and constructed synthetic radio observations to compare to observed galaxy clusters. We find a significant dependence of the observed morphology and radio relic properties on the viewing angle of the cluster, raising concerns regarding the interpretation of observed radio features in clusters. We also find that a given shock should not be characterized by a single Mach number. We find that the bulk of the radio emission comes from gas with T>5x10^7, ρ~10^(-28)-10^(-27) g/cm^3, with magnetic field strengths of 0.1-1.0 μG and shock Mach numbers of M~3-6. We present an analysis of the radio spectral index which suggests that the spatial variation of the spectral index can mimic synchrotron aging. Finally, we examine the polarization fraction and position angle of the simulated radio features, and compare to observations.

preprint2012arXiv

On The Road To More Realistic Galaxy Cluster Simulations: The Effects of Radiative Cooling and Thermal Feedback Prescriptions on the Observational Properties of Simulated Galaxy Clusters

Flux limited X-ray surveys of galaxy clusters show that clusters come in two roughly equally proportioned varieties: "cool core" clusters (CCs) and non-"cool core" clusters (NCCs). In previous work, we have demonstrated using cosmological $N$-body + Eulerian hydrodynamic simulations that NCCs are often consistent with early major mergers events that destroy embryonic CCs. In this paper we extend those results and conduct a series of simulationsusing different methods of gas cooling, and of energy and metal feedback from supernovae, where we attempt to produce a population of clusters with realistic central cooling times, entropies, and temperatures. We find that the use of metallicity-dependent gas cooling is essential to prevent early overcooling,and that adjusting the amount of energy and metal feedback can have a significant impact on observable X-ray quantities of the gas. We are able to produce clusters with more realistic central observable quantities than have previously been attained. However, there are still significant discrepancies between the simulated clusters and observations, which indicates that a different approach to simulating galaxies in clusters is needed. We conclude by looking towards a promising subgrid method of modeling galaxy feedback in clusters which may help to ameliorate the discrepancies between simulations and observations.

preprint2012arXiv

Optimized Multi-Frequency Spectra for Applications in Radiative Feedback and Cosmological Reionization

The recent implementation of radiative transfer algorithms in numerous hydrodynamics codes has led to a dramatic improvement in studies of feedback in various astrophysical environments. However, because of methodological limitations and computational expense, the spectra of radiation sources are generally sampled at only a few evenly-spaced discrete emission frequencies. Using one-dimensional radiative transfer calculations, we investigate the discrepancies in gas properties surrounding model stars and accreting black holes that arise solely due to spectral discretization. We find that even in the idealized case of a static and uniform density field, commonly used discretization schemes induce errors in the neutral fraction and temperature by factors of two to three on average, and by over an order of magnitude in certain column density regimes. The consequences are most severe for radiative feedback operating on large scales, dense clumps of gas, and media consisting of multiple chemical species. We have developed a method for optimally constructing discrete spectra, and show that for two test cases of interest, carefully chosen four-bin spectra can eliminate errors associated with frequency resolution to high precision. Applying these findings to a fully three-dimensional radiation-hydrodynamic simulation of the early universe, we find that the HII region around a primordial star is substantially altered in both size and morphology, corroborating the one-dimensional prediction that discrete spectral energy distributions can lead to sizable inaccuracies in the physical properties of a medium, and as a result, the subsequent evolution and observable signatures of objects embedded within it.

preprint2011arXiv

An MCMC approach to extracting the global 21-cm signal during the cosmic dawn from sky-averaged radio observations

Efforts are being made to observe the 21-cm signal from the 'cosmic dawn' using sky-averaged observations with individual radio dipoles. In this paper, we develop a model of the observations accounting for the 21-cm signal, foregrounds, and several major instrumental effects. Given this model, we apply Markov Chain Monte Carlo techniques to demonstrate the ability of these instruments to separate the 21-cm signal from foregrounds and quantify their ability to constrain properties of the first galaxies. For concreteness, we investigate observations between 40 and 120 MHz with the proposed DARE mission in lunar orbit, showing its potential for science return.

preprint2011arXiv

Cosmological Numerical Simulations of Radio Relics in Galaxy Clusters: Insights for Future Observations

The acceleration of electrons at shock fronts is thought to be responsible for radio relics, extended radio features in the vicinity of merging galaxy clusters. By combining high resolution Adaptive Mesh Refinement Hydro/N-body cosmological simulations with an accurate shock-finding algorithm and a model for electron acceleration, we calculate the expected synchrotron emission resulting from cosmological structure formation. From these simulations, we produce radio, SZE and X-ray images for a large sample of galaxy clusters along with radio luminosity functions and scaling relationships. We find that with upcoming radio arrays, we expect to see an abundance of radio emission associated with merger shocks in the intracluster medium. By producing observationally motivated statistics, we provide predictions that can be compared with observations to further our understanding of electron shock acceleration and kinematic structure of galaxy clusters.

preprint2011arXiv

Galaxy Cluster Radio Relics in Adaptive Mesh Refinement Cosmological Simulations: Relic Properties and Scaling Relationships

Cosmological shocks are a critical part of large-scale structure formation, and are responsible for heating the intracluster medium in galaxy clusters. In addition, they are also capable of accelerating non-thermal electrons and protons. In this work, we focus on the acceleration of electrons at shock fronts, which is thought to be responsible for radio relics - extended radio features in the vicinity of merging galaxy clusters. By combining high resolution AMR/N-body cosmological simulations with an accurate shock finding algorithm and a model for electron acceleration, we calculate the expected synchrotron emission resulting from cosmological structure formation. We produce synthetic radio maps of a large sample of galaxy clusters and present luminosity functions and scaling relationships. With upcoming long wavelength radio telescopes, we expect to see an abundance of radio emission associated with merger shocks in the intracluster medium. By producing observationally motivated statistics, we provide predictions that can be compared with observations to further improve our understanding of magnetic fields and electron shock acceleration.

preprint2011arXiv

Probing the First Stars and Black Holes in the Early Universe with the Dark Ages Radio Explorer (DARE)

A concept for a new space-based cosmology mission called the Dark Ages Radio Explore (DARE) is presented in this paper. DARE's science objectives include (1) When did the first stars form? (2) When did the first accreting black holes form? (3) When did Reionization begin? (4) What surprises does the end of the Dark Ages hold (e.g., Dark Matter decay)? DARE will use the highly-redshifted hyperfine 21-cm transition from neutral hydrogen to track the formation of the first luminous objects by their impact on the intergalactic medium during the end of the Dark Ages and during Cosmic Dawn (redshifts z=11-35). It will measure the sky-averaged spin temperature of neutral hydrogen at the unexplored epoch 80-420 million years after the Big Bang, providing the first evidence of the earliest stars and galaxies to illuminate the cosmos and testing our models of galaxy formation. DARE's approach is to measure the expected spectral features in the sky-averaged, redshifted 21-cm signal over a radio bandpass of 40-120 MHz. DARE orbits the Moon for a mission lifetime of 3 years and takes data above the lunar farside, the only location in the inner solar system proven to be free of human-generated radio frequency interference and any significant ionosphere. The science instrument is composed of a three-element radiometer, including electrically-short, tapered, bi-conical dipole antennas, a receiver, and a digital spectrometer. The smooth frequency response of the antennas and the differential spectral calibration approach using a Markov Chain Monte Carlo technique will be applied to detect the weak cosmic 21-cm signal in the presence of the intense solar system and Galactic foreground emissions.

preprint2010arXiv

Galaxy Clusters at the Edge: Temperature, Entropy, and Gas Dynamics at the Virial Radius

Recently, Suzaku has produced temperature and entropy profiles, along with profiles of gas density, gas fraction, and mass, for multiple galaxy clusters out to ~r_200 (~= virial radius). In this paper, we compare these novel X-ray observations with results from N-body + hydrodynamic adaptive mesh refinement cosmological simulations using the Enzo code. There is excellent agreement in the temperature, density, and entropy profiles between a sample of 27 mostly substructure-free massive clusters in the simulated volume and the observed clusters. This supports our previous contention that clusters have "universal" outer temperature profiles. Furthermore, it appears that the simplest adiabatic gas physics used in these Enzo simulations is adequate to model the outer regions of these clusters without other mechanisms (e.g., non-gravitational heating, cooling, magnetic fields, or cosmic rays). However, the outskirts of these clusters are not in hydrostatic equilibrium. There is significant bulk flow and turbulence in the outer intracluster medium created by accretion from filaments. Thus, the gas is not fully supported by thermal pressure. The implications for mass estimation from X-ray data are discussed.

preprint2010arXiv

The Properties of X-ray Cold Fronts in a Statistical Sample of Simulated Galaxy Clusters

We examine the incidence of cold fronts in a large sample of galaxy clusters extracted from a (512h^-1 Mpc) hydrodynamic/N-body cosmological simulation with adiabatic gas physics computed with the Enzo adaptive mesh refinement code. This simulation contains a sample of roughly 4000 galaxy clusters with M > 10^14 M_sun at z=0. For each simulated galaxy cluster, we have created mock 0.3-8.0 keV X-ray observations and spectroscopic-like temperature maps. We have searched these maps with a new automated algorithm to identify the presence of cold fronts in projection. Using a threshold of a minimum of 10 cold front pixels in our images, corresponding to a total comoving length L_cf > 156h^-1 kpc, we find that roughly 10-12% of all projections in a mass-limited sample would be classified as cold front clusters. Interestingly, the fraction of clusters with extended cold front features in our synthetic maps of a mass-limited sample trends only weakly with redshift out to z=1.0. However, when using different selection functions, including a simulated flux limit, the trending with redshift changes significantly. The likelihood of finding cold fronts in the simulated clusters in our sample is a strong function of cluster mass. In clusters with M>7.5x10^14 M_sun the cold front fraction is 40-50%. We also show that the presence of cold fronts is strongly correlated with disturbed morphology as measured by quantitative structure measures. Finally, we find that the incidence of cold fronts in the simulated cluster images is strongly dependent on baryonic physics.

preprint2009arXiv

On the Origin of Cool Core Galaxy Clusters: Comparing X-Ray Observations with Numerical Simulations

To better constrain models of cool core galaxy cluster formation, we have used X-ray observations taken from the Chandra and ROSAT archives to examine the properties of cool core and non-cool core clusters, especially beyond the cluster cores. We produced X-ray images, surface brightness profiles, and hardness ratio maps of 30 nearby rich Abell clusters (17 cool cores and 13 non-cool cores). We show that the use of double beta-models with cool core surface brightness profiles and single beta-models for non-cool core profiles yield statistically significant differences in the slopes (i.e., beta values) of the outer surface brightness profiles, but similar cluster core radii, for the two types of clusters. Hardness ratio profiles as well as spectroscopically-fit temperatures suggest that non-cool core clusters are warmer than cool core clusters of comparable mass beyond the cluster cores. We compared the properties of these clusters with the results from analogously reduced simulations of 88 numerical clusters created by the AMR Enzo code. The simulated surface brightness profiles have steeper beta-model fits in the outer cluster regions for both cool cores and non-cool cores, suggesting additional ICM heating is required compared to observed cluster ICMs. Temperature and surface brightness profiles reveal that the simulated clusters are over-cooled in their cores. As in the observations, however, simulated hardness ratio and temperature profiles indicate that non-cool core clusters are warmer than cool core clusters of comparable mass far beyond the cluster cores. The general similarities between observations and simulations support a model described in Burns et al. 2008 suggesting that non-cool core clusters suffered early major mergers destroying nascent cool cores.

preprint2009arXiv

The Santa Fe Light Cone Simulation Project: II. The Prospects for Direct Detection of the WHIM with SZE Surveys

Detection of the Warm-Hot Intergalactic Medium (WHIM) using Sunyaev-Zeldovich effect (SZE) surveys is an intriguing possibility, and one that may allow observers to quantify the amount of "missing baryons" in the WHIM phase. We estimate the necessary sensitivity for detecting low density WHIM gas with the South Pole Telescope (SPT) and Planck Surveyor for a synthetic 100 square degree sky survey. This survey is generated from a very large, high dynamic range adaptive mesh refinement cosmological simulation performed with the Enzo code. We find that for a modest increase in the SPT survey sensitivity (a factor of 2-4), the WHIM gas makes a detectable contribution to the integrated sky signal. For a Planck-like satellite, similar detections are possible with a more significant increase in sensitivity (a factor of 8-10). We point out that for the WHIM gas, the kinematic SZE signal can sometimes dominate the thermal SZE where the thermal SZE decrement is maximal (150 GHz), and that using the combination of the two increases the chance of WHIM detection using SZE surveys. However, we find no evidence of unique features in the thermal SZE angular power spectrum that may aid in its detection. Interestingly, there are differences in the power spectrum of the kinematic SZE, which may not allow us to detect the WHIM directly, but could be an important contaminant in cosmological analyses of the kSZE-derived velocity field. Corrections derived from numerical simulations may be necessary to account for this contamination.

preprint2008arXiv

Cluster Structure in Cosmological Simulations I: Correlation to Observables, Mass Estimates, and Evolution

We use Enzo, a hybrid Eulerian AMR/N-body code including non-gravitational heating and cooling, to explore the morphology of the X-ray gas in clusters of galaxies and its evolution in current generation cosmological simulations. We employ and compare two observationally motivated structure measures: power ratios and centroid shift. Overall, the structure of our simulated clusters compares remarkably well to low-redshift observations, although some differences remain that may point to incomplete gas physics. We find no dependence on cluster structure in the mass-observable scaling relations, T_X-M and Y_X-M, when using the true cluster masses. However, estimates of the total mass based on the assumption of hydrostatic equilibrium, as assumed in observational studies, are systematically low. We show that the hydrostatic mass bias strongly correlates with cluster structure and, more weakly, with cluster mass. When the hydrostatic masses are used, the mass-observable scaling relations and gas mass fractions depend significantly on cluster morphology, and the true relations are not recovered even if the most relaxed clusters are used. We show that cluster structure, via the power ratios, can be used to effectively correct the hydrostatic mass estimates and mass-scaling relations, suggesting that we can calibrate for this systematic effect in cosmological studies. Similar to observational studies, we find that cluster structure, particularly centroid shift, evolves with redshift. This evolution is mild but will lead to additional errors at high redshift. Projection along the line of sight leads to significant uncertainty in the structure of individual clusters: less than 50% of clusters which appear relaxed in projection based on our structure measures are truly relaxed.

preprint2006arXiv

AMR Simulations of the Cosmological Light Cone: SZE Surveys of the Synthetic Universe

We present preliminary results from simulated large sky coverage (~100 square degrees) Sunyaev-Zeldovich effect (SZE) cluster surveys using the cosmological adaptive mesh refinement N-body/hydro code Enzo. We have generated simulated light cones to match the resolution and sensitivity of current and future SZE instruments. These simulations are the most advanced calculations of their kind. The simulated sky surveys allow a direct comparison of large N-body/hydro cosmological simulations to current and pending sky surveys. Our synthetic surveys provide an indispensable guide for observers in the interpretation of large area sky surveys, and will develop the tools necessary to discriminate between models for cluster baryonic physics, and to accurately determine cosmological parameters.

preprint2006arXiv

Beyond the Cool Core: The Formation of Cool Core Galaxy Clusters

Why do some clusters have cool cores while others do not? In this paper, cosmological simulations, including radiative cooling and heating, are used to examine the formation and evolution of cool core (CC) and non-cool core (NCC) clusters. Numerical CC clusters at z=0 accreted mass more slowly over time and grew enhanced cool cores via hierarchical mergers; when late major mergers occurred, the CCs survived the collisions. By contrast, NCC clusters of similar mass experienced major mergers early in their evolution that destroyed embryonic cool cores and produced conditions that prevent CC re-formation. We discuss observational consequences.

preprint1995arXiv

TWO-POINT ANGULAR CORRELATION FUNCTION FOR THE GREEN BANK 4.85 GHZ SKY SURVEY

This paper presents an angular correlation analysis of the Green Bank 4.85 GHz radio catalog (Gregory \& Condon 1991) of 54,579 sources (S \gsim 25 mJy). The Green Bank catalog is found to be complete to S $\geq$ 35 mJy over 20$^{\circ} \leq δ< 74^{\circ}$, 0$^h \leq α< 24^h$, and Galactic latitude $|b| \geq 10^{\circ}$. The 2-point angular correlation function shows evidence for the clustering of radio sources, with a power-law distribution consistent with a slope $γ= -0.8$. This may well provide the $first$ detection of an angular correlation in a large area, complete deep radio survey.