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 {\Lambda}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.

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.

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.

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.

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.

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.

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_{\nu} = \nu^{\alpha}$ where $S_{\nu}$}. A possible origin of these radio sources is turbulence induced by sloshing of the gas in the cluster core.

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.

This is the final report submitted to NASA for a Probe-class concept study of the "Farside Array for Radio Science Investigations of the Dark ages and Exoplanets" (FARSIDE), a low radio frequency interferometric array on the farside of the Moon. The design study focused on the instrument, a deployment rover, the lander and base station, and delivered an architecture broadly consistent with the requirements for a Probe mission. This notional architecture consists of 128 dipole antennas deployed across a 10 km area by a rover, and tethered to a base station for central processing, power and data transmission to the Lunar Gateway, or an alternative relay satellite. FARSIDE would provide the capability to image the entire sky each minute in 1400 channels spanning frequencies from 150 kHz to 40 MHz, extending down two orders of magnitude below bands accessible to ground-based radio astronomy. The lunar farside can simultaneously provide isolation from terrestrial radio frequency interference, auroral kilometric radiation, and plasma noise from the solar wind. This would enable near-continuous monitoring of the nearest stellar systems in the search for the radio signatures of coronal mass ejections and energetic particle events, and would also detect the magnetospheres for the nearest candidate habitable exoplanets. Simultaneously, FARSIDE would be used to characterize similar activity in our own solar system, from the Sun to the outer planets. Through precision calibration via an orbiting beacon, and exquisite foreground characterization, FARSIDE would also measure the Dark Ages global 21-cm signal at redshifts from 50-100. It will also be a pathfinder for a larger 21-cm power spectrum instrument by carefully measuring the foreground with high dynamic range.

Neutral hydrogen pervades the infant Universe, and its redshifted 21-cm signal allows one to chart the Universe. This signal allows one to probe astrophysical processes such as the formation of the first stars, galaxies, (super)massive black holes and enrichment of the pristine gas from z~6 to z~30, as well as fundamental physics related to gravity, dark matter, dark energy and particle physics at redshifts beyond that. As one enters the Dark Ages (z>30), the Universe becomes pristine. Ground-based low-frequency radio telescopes aim to detect the spatial fluctuations of the 21-cm signal. Complementary, global 21-cm experiments aim to measure the sky-averaged 21-cm signal. Escaping RFI and the ionosphere has motivated space-based missions, such as the Dutch-Chinese NCLE instrument (currently in lunar L2), the proposed US-driven lunar or space-based instruments DAPPER and FARSIDE, the lunar-orbit interferometer DSL (China), and PRATUSH (India). To push beyond the current z~25 frontier, though, and measure both the global and spatial fluctuations (power-spectra/tomography) of the 21-cm signal, low-frequency (1-100MHz; BW~50MHz; z>13) space-based interferometers with vast scalable collecting areas (1-10-100 km2), large filling factors (~1) and large fields-of-view (4pi sr.) are needed over a mission lifetime of >5 years. In this ESA White Paper, we argue for the development of new technologies enabling interferometers to be deployed, in space (e.g. Earth-Sun L2) or in the lunar vicinity (e.g. surface, orbit or Earth-Moon L2), to target this 21-cm signal. This places them in a stable environment beyond the reach of most RFI from Earth and its ionospheric corruptions, enabling them to probe the Dark Ages as well as the Cosmic Dawn, and allowing one to investigate new (astro)physics that is inaccessible in any other way in the coming decades. [Abridged]

Detecting the cosmological sky-averaged (global) 21 cm signal as a function of observed frequency will provide a powerful tool to study the ionization and thermal history of the intergalactic medium (IGM) in the early Universe ($\sim$ 400 million years after the Big Bang). The greatest challenge in conventional total-power global 21 cm experiments is the removal of the foreground synchrotron emission ($\sim 10^3$-$10^4$ K) to uncover the weak cosmological signal (tens to hundreds of mK), especially since the intrinsic smoothness of the foreground spectrum is corrupted by instrumental effects. Although the EDGES team has recently reported an absorption profile at 78 MHz in the sky-averaged spectrum, it is necessary to confirm this detection with an independent approach. The projection effect from observing anisotropic foreground source emission with a wide-view antenna pointing at the North Celestial Pole (NCP) can induce a net polarization, referred as the Projection-Induced Polarization Effect (PIPE). Due to Earth's rotation, observation centered at the circumpolar region will impose a dynamic sky modulation on the net polarization's waveforms which is unique to the foreground component. In this study, we review the implementation practicality and underlying instrumental effects of this new polarimetry-based technique with detailed numerical simulation and a testbed instrument, the Cosmic Twilight Polarimeter (CTP). In addition, we explore an SVD-based analysis approach for separating the foreground and instrumental effects from the background global 21 cm signal using the sky-modulated PIPE.

Due to ionosphere absorption and the interference by natural and artificial radio emissions, ground observation of the sky at the decameter or longer is very difficult. This unexplored part of electromagnetic spectrum has the potential of great discoveries, notably in the study of cosmic dark ages and dawn, but also in heliophysics and space weather, planets, cosmic ray and neutrinos, pulsar and interstellar medium, extragalactic radio sources, and even SETI. At a forum organized by the International Space Science Institute-Beijing (ISSI-BJ), we discussed the prospect of opening up this window for astronomical observations by using a constellation of small or micro-satellites. We discussed the past experiments and the current ones such as the low frequency payload on Chang'e-4 mission lander, relay satellite and the Longjiang satellite, and also the future DSL mission, which is a linear array on lunar orbit which can make synthesized map of the whole sky as well as measure the global spectrum. We also discuss the synergy with other experiments, including ground global experiments such as EDGES, SARAS, SCI-HI and High-z, PRIZM/Albatros, ground imaging facillities such as LOFAR and MWA, and space experiments such as SUNRISE, DARE/DAPPER and PRATUSH. We also discussed some technical aspects of the DSL concept.

FARSIDE (Farside Array for Radio Science Investigations of the Dark ages and Exoplanets) is a Probe-class concept to place a low radio frequency interferometric array on the farside of the Moon. A NASA-funded design study, focused on the instrument, a deployment rover, the lander and base station, delivered an architecture broadly consistent with the requirements for a Probe mission. This notional architecture consists of 128 dual polarization antennas deployed across a 10 km area by a rover, and tethered to a base station for central processing, power and data transmission to the Lunar Gateway. FARSIDE would provide the capability to image the entire sky each minute in 1400 channels spanning frequencies from 100 kHz to 40 MHz, extending down two orders of magnitude below bands accessible to ground-based radio astronomy. The lunar farside can simultaneously provide isolation from terrestrial radio frequency interference, auroral kilometric radiation, and plasma noise from the solar wind. This would enable near-continuous monitoring of the nearest stellar systems in the search for the radio signatures of coronal mass ejections and energetic particle events, and would also detect the magnetospheres for the nearest candidate habitable exoplanets. Simultaneously, FARSIDE would be used to characterize similar activity in our own solar system, from the Sun to the outer planets, including the hypothetical Planet Nine. Through precision calibration via an orbiting beacon, and exquisite foreground characterization, FARSIDE would also measure the Dark Ages global 21-cm signal at redshifts z=50-100. The unique observational window offered by FARSIDE would enable an abundance of additional science ranging from sounding of the lunar subsurface to characterization of the interstellar medium in the solar system neighborhood.

ClusterPyXT is a new software pipeline to generate spectral temperature, X-ray surface brightness, pressure, and density maps from X-ray observations of galaxy clusters. These data products help elucidate the physics of processes occurring within clusters of galaxies, including turbulence, shock fronts, nonthermal phenomena, and the overall dynamics of cluster mergers. ClusterPyXT automates the creation of these data products with minimal user interaction, and allows for rapid analyses of archival data with user defined parameters and the ability to straightforwardly incorporate additional observations. In this paper, we describe in detail the use of this code and release it as an open source Python project on GitHub.

The epoch of reionization, when photons from early galaxies ionized the intergalactic medium about a billion years after the Big Bang, is the last major phase transition in the Universe's history. Measuring the characteristics of the transition is important for understanding early galaxies and the cosmic web and for modeling dwarf galaxies in the later Universe. But such measurements require probes of the intergalactic medium itself. Here we describe how the 21-cm line of neutral hydrogen provides a powerful probe of the reionization process and therefore important constraints on both the galaxies and intergalactic absorbers at that time. While existing experiments will make precise statistical measurements over the next decade, we argue that improved 21-cm analysis techniques - allowing imaging of the neutral gas itself - as well as improved theoretical models, are crucial for testing our understanding of this important era.

The early phases of galaxy formation constitute one of the most exciting frontiers in astrophysics. It is during this era that the first luminous sources reionize the intergalactic medium - the moment when structure formation affects every baryon in the Universe. Here we argue that we will obtain a complete picture of this era by combining observations of galaxies with direct measurements of the reionization process: the former will provide a detailed understanding of bright sources, while the latter will constrain the (substantial) faint source population. We further describe how optimizing the comparison of these two measurements requires near-infrared galaxy surveys covering large volumes and retaining redshift information and also improvements in 21-cm analysis, moving those experiments into the imaging regime.

In addition to being a probe of Cosmic Dawn and Epoch of Reionization astrophysics, the 21cm line at $z>6$ is also a powerful way to constrain cosmology. Its power derives from several unique capabilities. First, the 21cm line is sensitive to energy injections into the intergalactic medium at high redshifts. It also increases the number of measurable modes compared to existing cosmological probes by orders of magnitude. Many of these modes are on smaller scales than are accessible via the CMB, and moreover have the advantage of being firmly in the linear regime (making them easy to model theoretically). Finally, the 21cm line provides access to redshifts prior to the formation of luminous objects. Together, these features of 21cm cosmology at $z>6$ provide multiple pathways toward precise cosmological constraints. These include the "marginalizing out" of astrophysical effects, the utilization of redshift space distortions, the breaking of CMB degeneracies, the identification of signatures of relative velocities between baryons and dark matter, and the discovery of unexpected signs of physics beyond the $\Lambda$CDM paradigm at high redshifts.

The Dark Ages are the period between the last scattering of the cosmic microwave background and the appearance of the first luminous sources, spanning approximately 1100 < z < 30. The only known way to measure fluctuations in this era is through the 21-cm line of neutral hydrogen. Such observations have enormous potential for cosmology, because they span a large volume while the fluctuations remain linear even on small scales. Observations of 21-cm fluctuations during this era can therefore constrain fundamental aspects of our Universe, including inflation and any exotic physics of dark matter. While the observational challenges to these low-frequency 21-cm observations are enormous, especially from the terrestrial environment, they represent an important goal for cosmology.

The "cosmic dawn" refers to the period of the Universe's history when stars and black holes first formed and began heating and ionizing hydrogen in the intergalactic medium (IGM). Though exceedingly difficult to detect directly, the first stars and black holes can be constrained indirectly through measurements of the cosmic 21-cm background, which traces the ionization state and temperature of intergalactic hydrogen gas. In this white paper, we focus on the science case for such observations, in particular those targeting redshifts z $\gtrsim$ 10 when the IGM is expected to be mostly neutral. 21-cm observations provide a unique window into this epoch and are thus critical to advancing first star and black hole science in the next decade.

Cosmic dawn and the Epoch of Reionization (EoR) are among the least explored observational eras in cosmology: a time at which the first galaxies and supermassive black holes formed and reionized the cold, neutral Universe of the post-recombination era. With current instruments, only a handful of the brightest galaxies and quasars from that time are detectable as individual objects, due to their extreme distances. Fortunately, a multitude of multi-wavelength intensity mapping measurements, ranging from the redshifted 21 cm background in the radio to the unresolved X-ray background, contain a plethora of synergistic information about this elusive era. The coming decade will likely see direct detections of inhomogenous reionization with CMB and 21 cm observations, and a slew of other probes covering overlapping areas and complementary physical processes will provide crucial additional information and cross-validation. To maximize scientific discovery and return on investment, coordinated survey planning and joint data analysis should be a high priority, closely coupled to computational models and theoretical predictions.

The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.

Most of the current Global 21-cm experiments include ground screens that help moderate effects from the Earth. In this paper, we report on a possible systematic artifact within the ground plane that may produce broad absorption features in the spectra observed by these experiments. Using analytical approximations and numerical modeling, the origin of the artifact and its impact on the sky-averaged spectrum are described. The publicly released EDGES dataset, from which a 78 MHz absorption feature was recently suggested, is used to probe for the potential presence of ground plane resonances. While the lack of a noise level for the EDGES spectrum makes traditional goodness-of-fit statistics unattainable, the rms residual can be used to assess the relative goodness of fits performed under similar circumstances. The fit to the EDGES spectrum using a model with a simple 2-term foreground and three cavity-mode resonances is compared to a fit to the same spectrum with a model used by the EDGES team consisting of a 5-term foreground and a flattened Gaussian signal. The fits with the physically motivated resonance and empirical flattened Gaussian models have rms residuals of 20.8 mK (11 parameters) and 24.5 mK (9 parameters), respectively, allowing us to conclude that ground plane resonances constitute another plausible explanation for the EDGES data.

The Dark Ages, probed by the redshifted 21-cm signal, is the ideal epoch for a new rigorous test of the standard LCDM cosmological model. Divergences from that model would indicate new physics, such as dark matter decay (heating) or baryonic cooling beyond that expected from adiabatic expansion of the Universe. In the early Universe, most of the baryonic matter was in the form of neutral hydrogen (HI), detectable via its ground state's spin-flip transition. A measurement of the redshifted 21-cm spectrum maps the history of the HI gas through the Dark Ages and Cosmic Dawn and up to the Epoch of Reionization (EoR). The Experiment to Detect the Global EoR Signature (EDGES) recently reported an absorption trough at 78 MHz (redshift z of 17), similar in frequency to expectations for Cosmic Dawn, but about 3 times deeper than was thought possible from standard cosmology and adiabatic cooling of HI. Interactions between baryons and slightly-charged dark matter particles with electron-like mass provide a potential explanation of this difference but other cooling mechanisms are also being investigated to explain these results. The Cosmic Dawn trough is affected by cosmology and the complex astrophysical history of the first luminous objects. Another trough is expected during the Dark Ages, prior to the formation of the first stars and thus determined entirely by cosmological phenomena (including dark matter). Observations on or in orbit above the Moon's farside can investigate this pristine epoch (15-40 MHz; z=100-35), which is inaccessible from Earth. A single cross-dipole antenna or compact array can measure the amplitude of the 21-cm spectrum to the level required to distinguish (at >5$\sigma$}) the standard cosmological model from that of additional cooling derived from current EDGES results. This observation constitutes a powerful, clean probe of exotic physics in the Dark Ages.

The reduced chi-squared statistic is a commonly used goodness-of-fit measure, but it cannot easily detect features near the noise level, even when a large amount of data is available. In this paper, we introduce a new goodness-of-fit measure that we name the reduced psi-squared statistic. It probes the two-point correlations in the residuals of a fit, whereas chi-squared accounts for only the absolute values of each residual point, not considering the relationship between these points. The new statistic maintains sensitivity to individual outliers, but is superior to chi-squared in detecting wide, low level features in the presence of a large number of noisy data points. After presenting this new statistic, we show an instance of its use in the context of analyzing radio spectroscopic data for 21-cm cosmology experiments. We perform fits to simulated data with four components: foreground emission, the global 21-cm signal, an idealized instrument systematic, and noise. This example is particularly timely given the ongoing efforts to confirm the first observational result for this signal, where this work found its original motivation. In addition, we release a Python script dubbed $\texttt{psipy}$ which allows for quick, efficient calculation of the reduced psi-squared statistic on arbitrary data arrays, to be applied in any field of study.

NASA and ESA are preparing a series of human exploration missions using the four-person Orion crew vehicle, launched by NASA's Space Launch System, and a Lunar Orbiting Platform-Gateway (LOP-G) that enable long duration (>30 days) operations in cis-lunar space. This will provide an opportunity for science and exploration from the lunar surface facilitated by low latency surface telerobotics. We describe two precursor experiments, using the International Space Station (ISS) and a student-built teleoperated rover, which are laying the groundwork for remote operation of rovers on the Moon by astronauts aboard the LOP-G. Such missions will open the lunar far side, among other sites, for exploration and science. We describe examples of two high-priority, lunar science missions that can be conducted using low latency surface telerobotics including an astronaut-assisted far side sample return and the deployment/construction of a low frequency radio telescope array to observe the first stars and galaxies (Cosmic Dawn). The lessons learned from these lunar operations will feed-forward to future low latency telepresence missions on Mars.

We present results from an X-ray and radio study of the merging galaxy cluster Abell 115. We use the full set of 5 Chandra observations taken of A115 to date (360 ks total integration) to construct high-fidelity temperature and surface brightness maps. We also examine radio data from the Very Large Array at 1.5 GHz and the Giant Metrewave Radio Telescope at 0.6 GHz. We propose that the high X-ray spectral temperature between the subclusters results from the interaction of the bow shocks driven into the intracluster medium by the motion of the subclusters relative to one another. We have identified morphologically similar scenarios in Enzo numerical N-body/hydrodynamic simulations of galaxy clusters in a cosmological context. In addition, the giant radio relic feature in A115, with an arc-like structure and a relatively flat spectral index, is likely consistent with other shock-associated giant radio relics seen in other massive galaxy clusters. We suggest a dynamical scenario that is consistent with the structure of the X-ray gas, the hot region between the clusters, and the radio relic feature.

We document how the narrative and the policies of space exploration in the U.S. have changed over the past 50 years. We first examine the history of the U.S. space exploration program and also assess three current conditions of space exploration including: (1) the increasing role of the private sector, (2) the influence of global politics, and (3) the focus on a human mission to Mars. We identify five rhetorical themes: competition, prestige, collaboration, leadership, and a new paradigm. These themes are then used to analyze the content of forty documents from eight presidential administrations. The historical narrative and content analysis together suggest that space exploration has developed from a discourse about a bipolar world composed of the U.S. and U.S.S.R. into a complicated field that encompass many new players. We make three observations: (1) there is a disconnect between stated U.S. policy goals and the implementation of those goals, (2) the U.S. communicates mixed messages regarding its intent to be both the dominant leader in space exploration and also a committed participant in international collaborations, and (3) the U.S. cannot remain a true pioneer in space exploration if it does not embrace the realities of globalization and the changing dynamics within space exploration. We conclude with three suggestions: (1) the U.S. government and NASA should critically examine space exploration priorities and commit to implementing a program that will further realistic policy and goals, (2) the U.S. should re-examine its intention to play a dominant leadership role in space exploration and consider emphasizing a commitment toward active participation in international collaboration in space, and (3) the U.S. should fully embrace the new paradigm of space exploration by lowering barriers that hinder competitiveness.

NASA's SLS and Orion crew vehicle will launch humans to cislunar space to begin the new era of space exploration. NASA plans to use the Orion crew vehicle to transport humans between Earth and cislunar space where there will be a stationed habitat known as the Deep Space Gateway (DSG). The proximity to the lunar surface allows for direct communication between the DSG and surface assets, which enables low-latency telerobotic exploration. The operational constraints for telerobotics must be fully explored on Earth before being utilized on space exploration missions. We identified two constraints on space exploration using low-latency surface telerobotics and attempts to quantify these constraints. A constraint associated with low-latency surface telerobotics is the bandwidth available between the orbiting command station and the ground assets. The bandwidth available will vary during operation. As a result, it is critical to quantify the operational video conditions required for effective exploration. We designed an experiment to quantify the threshold frame rate required for effective exploration. The experiment simulated geological exploration via low-latency surface telerobotics using a COTS rover in a lunar analog environment. The results from this experiment indicate that humans should operate above a threshold frame rate of 5 frames per second. In a separate, but similar experiment, we introduced a 2.6 second delay in the video system. This delay recreated the latency conditions present when operating rovers on the lunar farside from an Earth-based command station. This time delay was compared to low-latency conditions for teleoperation at the DSG ($\leq$0.4 seconds). The results from this experiment show a 150% increase in exploration time when the latency is increased to 2.6 seconds. This indicates that such a delay significantly complicates real-time exploration strategies.

The sky-averaged (global) highly redshifted 21-cm spectrum from neutral hydrogen is expected to appear in the VHF range of $\sim20-200$ MHz and its spectral shape and strength are determined by the heating properties of the first stars and black holes, by the nature and duration of reionization, and by the presence or absence of exotic physics. Measurements of the global signal would therefore provide us with a wealth of astrophysical and cosmological knowledge. However, the signal has not yet been detected because it must be seen through strong foregrounds weighted by a large beam, instrumental calibration errors, and ionospheric, ground and radio-frequency-interference effects, which we collectively refer to as "systematics". Here, we present a signal extraction method for global signal experiments which uses Singular Value Decomposition (SVD) of "training sets" to produce systematics basis functions specifically suited to each observation. Instead of requiring precise absolute knowledge of the systematics, our method effectively requires precise knowledge of how the systematics can vary. After calculating eigenmodes for the signal and systematics, we perform a weighted least square fit of the corresponding coefficients and select the number of modes to include by minimizing an information criterion. We compare the performance of the signal extraction when minimizing various information criteria and find that minimizing the Deviance Information Criterion (DIC) most consistently yields unbiased fits. The methods used here are built into our widely applicable, publicly available Python package, $\texttt{pylinex}$, which analytically calculates constraints on signals and systematics from given data, errors, and training sets.

The viscous spreading of planetary rings is believed to be counteracted by satellite torques, either through an individual resonance or through overlapping resonances. For the A ring of Saturn, it has been commonly believed that the satellite Janus alone can prevent the ring from spreading via its 7:6 Lindblad resonance. We discuss this common misconception and show that, in reality, the A ring is confined by the contributions from the group of satellites Pan, Atlas, Prometheus, Pandora, Janus, Epimetheus, and Mimas, whose cumulative torques from various resonances gradually decrease the angular momentum flux transported outward through the ring via density and bending waves. We further argue that this decrease in angular momentum flux occurs through 'flux reversal'. Furthermore, we use the magnitude of the satellites' resonance torques to estimate the effective viscosity profile across the A ring, showing that it decreases with radius from ~50 cm2 s-1 to less than ~10 cm2 s-1. The gradual estimated decrease of the angular momentum flux and effective viscosity are roughly consistent with results obtained by balancing the shepherding torques from Pan and Daphnis with the viscous torque at the edges of the Encke and Keeler gaps, as well as the edge of the A ring. On the other hand, the Mimas 2:1 Lindblad resonance alone seems to be capable of confining the edge of the B ring, and contrary to the situation in the A ring, we show that the effective viscosity across the B ring is relatively constant at ~24-30 cm2 s-1.

Many obsects in the solar system are suspected to have experience reorientation of their spin axes. As their rotation rates are slow and their shapes are nearly spherical, the formation of mass anomalies, by either endogenic of exogenic processes, can change objects' moments of inertia. Therefore, the objects reorient to align their largest moment of inertia with their spin axis. Such phenomenon is called True Polar Wander (TPW). Here we report the discovery of a global series of topographic lows on Saturn's satellite Enceladus that we interpret to show that this synchronously locked moon has undergone TPW by ~55{\deg} about the tidal axis. We use improved topographic data from the spherical harmonic expansion of Cassini limb and stereogrammetric measurements to characterize regional topography over the surface of Enceladus. We identify a group of nearly antipodal basins orthogonal to a topographic basin chain tracing a non-equatorial circumglobal belt across Enceladus' surface. We argue that the belt and the antipodal regions are fossil remnants of an earlier equator and poles, respectively. We argue that these lows arise from isostasic compensation and that their pattern reflects spatial variations in internal dynamics of the ice shell. Our hypothesis is consistent with a variety of geological features visible in Cassini images.

We analyze several thousand Cassini ISS images in order to study the inner edge of the Keeler gap in Saturn's outer A ring. We find strong evidence for an m=32 perturbation with a mean amplitude of radial variation of 4.5 km. Phase analysis yields a pattern speed consistent with the mean motion of Prometheus, indicating that this pattern is generated by the 32:31 Inner Lindblad resonance with Prometheus. In addition, we find evidence of 18-lobed and 20-lobed patterns with amplitudes of ~1.5 km. These patterns, whose rotation rates correspond to resonance locations ~4 km interior to the gap edge, are believed to be normal modes. The former is probably related to the nearby 18:17 (m=18) resonance with Pandora. In addition to these resonant and normal mode patterns, we also observe multiple localized features that appear to move at the local keplerian rate and that persist for only a few months. One hypothesis is that different groups of ring particles at the inner edge of the gap may be reacting differently to the resonance with Prometheus, with local variations in the forced eccentricity and/or pericenter; an alternative hypothesis is the existence of several unseen objects embedded at or near the inner edge of the Keeler gap, similar to those suspected to exist at the outer edges of the A and B rings (Spitale and Porco 2009, 2010). In either case, observations of the ring edge at opposite ansae demonstrate that the localized features must be on eccentric orbits.

A review of the jovian ring system

The Global Exploration Roadmap indicates the need for increased human exploration of under-sampled regions of our solar system in order to make new scientific discoveries. The high costs and dangers of sending humans deeper into our solar system necessitates the use of human-robotic partnerships, especially in transitioning from low-Earth orbit to deep-space operations. Low-latency planetary surface exploration is an example of a human-robotic partnership that provides an exciting option for effective, low-cost exploration of our solar system. However, low-latency telerobotic exploration is a new concept for space exploration and needs to be tested for its limits and effectiveness. This paper focuses on a human operator's ability to identify exploration targets in an unfamiliar environment using real-time low-latency telerobotics under various frame rate conditions. This relationship was investigated using a Telerobotic Simulation System (TSS). The frame rates were varied and the order of the exploration tasks were randomized for each operator. The rover operated at peak speeds of one meter per second with a video stream resolution of 640x480 and colorscale of 24 bits. The results from this experiment indicate that 5 frames per second is the minimum necessary frame rate for effective exploration.

The redshifted 21-cm monopole is expected to be a powerful probe of the epoch of the first stars and galaxies ($10<z<35$). The global 21-cm signal is sensitive to the thermal and ionization state of hydrogen gas and thus provides a tracer of sources of energetic photons -- primarily hot stars and accreting black holes -- which ionize and heat the high redshift intergalactic medium (IGM). This paper presents a strategy for observations of the global spectrum with a realizable instrument placed in a low altitude lunar orbit, performing night-time 40-120 MHz spectral observations, while on the farside to avoid terrestrial radio frequency interference, ionospheric corruption, and solar radio emissions. The frequency structure, uniformity over large scales, and unpolarized state of the redshifted 21-cm spectrum are distinct from the spectrally featureless, spatially-varying, and polarized emission from the bright foregrounds. This allows a clean separation between the primordial signal and foregrounds. For signal extraction, we model the foreground, instrument, and 21-cm spectrum with eigenmodes calculated via Singular Value Decomposition analyses. Using a Markov Chain Monte Carlo algorithm to explore the parameter space defined by the coefficients associated with these modes, we illustrate how the spectrum can be measured and how astrophysical parameters (e.g., IGM properties, first star characteristics) can be constrained in the presence of foregrounds using the Dark Ages Radio Explorer (DARE).

The Dark Ages Radio Experiment (DARE) seeks to study the cosmic Dark Ages approximately 80 to 420 million years after the Big Bang. Observations require truly quiet radio conditions, shielded from Sun and Earth electromagnetic (EM) emissions, on the far side of the Moon. DARE's science orbit is a frozen orbit with respect to lunar gravitational perturbations. The altitude and orientation of the orbit remain nearly fixed indefinitely, maximizing science time without the need for maintenance. DARE's observation targets avoid the galactic center and enable investigation of the universe's first stars and galaxies.

The cosmological global (sky-averaged) 21-cm signal is a powerful tool to probe the evolution of the intergalactic medium (IGM) in high-redshift Universe ($z \leq 6$). One of the biggest observational challenges is to remove the foreground spectrum which is at least four orders of magnitude brighter than the cosmological 21-cm emission. Conventional global 21-cm experiments rely on the spectral smoothness of the foreground synchrotron emission to separate it from the unique 21-cm spectral structures in a single total-power spectrum. However, frequency-dependent instrumental and observational effects are known to corrupt such smoothness and complicates the foreground subtraction. We introduce a polarimetric approach to measure the projection-induced polarization of the anisotropic foreground onto a stationary dual-polarized antenna. Due to Earth rotation, when pointing the antenna at a celestial pole, the revolving foreground will modulate this polarization with a unique frequency-dependent sinusoidal signature as a function of time. In our simulations, by harmonic decomposing this dynamic polarization, our technique produces two separate spectra in parallel from the same observation: (i) a total sky power consisting both the foreground and the 21-cm background, (ii) a model-independent measurement of the foreground spectrum at a harmonic consistent to twice the sky rotation rate. In the absence of any instrumental effects, by scaling and subtracting the latter from the former, we recover the injected global 21-cm model within assumed uncertainty. We further discuss several limiting factors and potential remedies for future implementation.

Additional physics beyond standard hydrodynamics is needed to fully model the intracluster medium (ICM); however, as we move to more sophisticated models, it is important to consider the role of magnetic fields and the way the fluid approximation breaks down. This paper represents a first step towards developing a self-consistent model of the ICM by characterizing the statistical properties of magnetic fields in cosmological simulations of galaxy clusters. We find that plasma conditions are largely homogeneous across a range of cluster masses and relaxation states. We also find that the magnetic field length scales are resolution dependent and not based on any particular physical process. Energy transfer mechanisms and scales are also identified, and imply the existence of small scale dynamo action. The scales of the small scale dynamo are resolution limited and driven by numerical resistivity and viscosity.

We present observations at optical wavelengths with the Cassini Spacecraft's Imaging Science System of the Phoebe ring, a vast debris disk around Saturn that seems to be collisionally generated by its irregular satellites. The analysis reveals a radial profile from 80-260 Saturn radii ($R_S$) that changes behavior interior to $\approx 110 R_S$, which we attribute to either the moon Iapetus sweeping up small particles, or to orbital instabilities that cause the ring to flare up vertically. Our study yields an integrated I/F at 0.635 $\mu$m along Saturn's shadow in the Phoebe ring's midplane from 80-250 $R_S$ of $2.7^{+0.9}_{-0.3} \times 10^{-9}$. We develop an analytical model for the size-dependent secular dynamics of retrograde Phoebe ring grains, and compare this model to the observations. This analysis implies that 1) the "Phoebe" ring is partially sourced by debris from irregular satellites beyond Phoebe's orbit and 2) the scattered light signal is dominated by small grains ($\lesssim 20\mu$m in size). If we assume that the Phoebe ring is generated through steady-state micrometeoroid bombardment, this implies a power-law size distribution with index $> 4$, which is unusually steep among solar system rings. This suggests either a steep size distribution of ejecta when material is initially released, or a subsequent process that preferentially breaks up large grains.

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 $\nu^{-2}$ (where $\nu$ is the frequency of observation) in the presence of a bright galactic foreground ($\propto \nu^{-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.

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.

The James Webb Space Telescope (JWST) will provide unprecedented opportunities to observe the rings and small satellites in our solar system, accomplishing three primary objectives: 1) discovering new rings and moons, 2) unprecedented spectroscopy, and 3) time-domain observations. We give details on these science objectives and describe requirements that JWST must fulfill in order to accomplish the science objectives

We present a study of the behavior of Saturn's A ring outer edge, using images and occultation data obtained by the Cassini spacecraft over a period of 8 years from 2006 to 2014. More than 5000 images and 170 occultations of the A ring outer edge are analyzed. Our fits confirm the expected response to the Janus 7:6 Inner Lindblad resonance (ILR) between 2006 and 2010, when Janus was on the inner leg of its regular orbit swap with Epimetheus. During this period, the edge exhibits a regular 7-lobed pattern with an amplitude of 12.8 km and one minimum aligned with the orbital longitude of Janus, as has been found by previous investigators. However, between 2010 and 2014, the Janus/Epimetheus orbit swap moves the Janus 7:6 LR away from the A ring outer edge, and the 7-lobed pattern disappears. In addition to several smaller-amplitudes modes, indeed, we found a variety of pattern speeds with different azimuthal wave numbers, and many of them may arise from resonant cavities between the ILR and the ring edge; also we found some other signatures consistent with tesseral resonances that could be associated with inhomogeneities in Saturn's gravity field. Moreover, these signatures do not have a fixed pattern speed. We present an analysis of these data and suggest a possible dynamical model for the behavior of the A ring's outer edge after 2010.

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-$\alpha$, 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 \nu / \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.

Several planetary satellites apparently have subsurface seas that are of great interest for, among other reasons, their possible habitability. The geologically diverse Saturnian satellite Enceladus vigorously vents liquid water and vapor from fractures within a south polar depression and thus must have a liquid reservoir or active melting. However, the extent and location of any subsurface liquid region is not directly observable. We use measurements of control points across the surface of Enceladus accumulated over seven years of spacecraft observations to determine the satellite's precise rotation state, finding a forced physical libration of 0.120 $\pm$ 0.014{\deg} (2{\sigma}). This value is too large to be consistent with Enceladus's core being rigidly connected to its surface, and thus implies the presence of a global ocean rather than a localized polar sea. The maintenance of a global ocean within Enceladus is problematic according to many thermal models and so may constrain satellite properties or require a surprisingly dissipative Saturn.

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.

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.

The harsh environment on the lunar surface presents unique technological challenges for space exploration. This paper presents research on the design and development of the Tele- scope-deployment High-vacuum teleOperated Rover (THOR), currently being built and tested in the Lunar and Airless Bodies Simulator (LABS) facility at the University of Colorado Boulder. This rover is fabricated entirely out of cost-effective commercial off-the-shelf (COTS) components and materials. THOR can potentially survive for more than one simulated year in conditions similar to that of the lunar environment, demonstrating the successful initial results of a first phase research study on material and electronic survivability in an extreme environment such as the Moon.

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.

Previous investigations of Saturn's outer D ring (73,200-74,000 km from Saturn's center) identified periodic brightness variations whose radial wavenumber increased linearly over time. This pattern was attributed to a vertical corrugation, and its temporal variability implied that some event --possibly an impact with interplanetary debris-- caused the ring to become tilted out the planet's equatorial plane in 1983. This work examines these patterns in greater detail using a more extensive set of Cassini images in order to obtain additional insights into the 1983 event. These additional data reveal that the D ring is not only corrugated, but also contains a time-variable periodic modulation in its optical depth that probably represents organized eccentric motions of the D-ring's particles. This second pattern suggests that whatever event tilted the rings also disturbed the radial or azimuthal velocities of the ring particles. Furthermore, the relative amplitudes of the two patterns indicate that the vertical motions induced by the 1983 event were 2.3+/-0.5 times larger than the corresponding in-plane motions. If these structures were indeed produced by an impact, material would need to strike the ring at a steep angle (more than 60 degrees from the ring plane) to produce such motions. Meanwhile, the corrugation wavelengths in the D ring are about 0.7% shorter than one would predict based on extrapolations from similar structures in the nearby C ring. This could indicate that the D-ring was tilted/disturbed about 60 days before the C ring. Such a timing difference could be explained if the material that struck the ring was derived from debris released when some object broke up near Saturn some months earlier. This debris would need to have a substantial initial velocity dispersion and have its orbital properties perturbed by solar tides prior to its collision with the rings.

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.

X-ray observations were used to examine the complex temperature structure of Abell 85, a cool-core galaxy cluster. Temperature features can provide evidence of merging events which shock heat the intracluster gas. Temperature maps were made from both \textit{Chandra} and \textit{XMM-Newton} obervations. The combination of a new, long-exposure \textit{XMM} observation and an improved temperature map binning technique produced the highest fidelity temperature maps of A85 to date. Hot regions were detected near the subclusters to the South and Southwest in both the \textit{Chandra} and \textit{XMM} temperature maps. The presence of these structures implies A85 is not relaxed. The hot regions may indicate the presence of shocks. The Mach numbers were estimated to be $\sim$1.9 at the locations of the hot spots. Observational effects will tend to systematically reduce temperature jumps, so the measured Mach numbers are likely underestimated. Neither temperature map showed evidence for a shock in the vicinity of the presumed radio relic near the Southwest subcluster. However, the presence of a weak shock cannot be ruled out. There was tension between the temperatures measured by the two instruments.

D68 is a narrow ringlet located only 67,627 km (1.12 planetary radii) from Saturn's spin axis. Images of this ringlet obtained by the Cassini spacecraft reveal that this ringlet exhibits persistent longitudinal brightness variations and a substantial eccentricity (ae=25+/-1 km). By comparing observations made at different times, we confirm that the brightness variations revolve around the planet at approximately the local orbital rate (1751.6 degrees/day), and that the ringlet's pericenter precesses at 38.243+/-0.008 degrees/day, consistent with the expected apsidal precession rate at this location due to Saturn's higher-order gravitational harmonics. Surprisingly, we also find that the ringlet's semi-major axis appears to be decreasing with time at a rate of 2.4+/-0.4 km/year between 2005 and 2013. A closer look at these measurements, along with a consideration of earlier Voyager observations of this same ringlet, suggests that the mean radius of D68 moves back and forth, perhaps with a period of around 15 Earth years or about half a Saturn year. These observations could place important constraints on both the ringlet's local dynamical environment and the planet's gravitational field.

The Phoebe ring, Saturn's largest and faintest ring, lies far beyond the planet's well-known main rings. It is primarily sourced by collisions with Saturn's largest irregular satellite Phoebe, perhaps through stochastic macroscopic collisions, or through more steady micrometeoroid bombardment. The ring was discovered with the Spitzer Space Telescope at 24 $\mu$m and has a normal optical depth of $\sim 2 \times 10^{-8}$ (Verbiscer et al. 2009). We report the first observations of sunlight scattered off the Phoebe ring using the Cassini spacecraft's ISS camera at optical wavelengths. We find that material between $\approx 130-210$ Saturnian radii ($R_S$) from the planet produces an I/F of $1.7 \pm 0.1 \times 10^{-11}$ per $R_S$ of the line-of-sight distance through the disk. Combining our measurements with the Spitzer infrared data, we can place constraints on the ring-particles' light-scattering properties. Depending on the particles' assumed phase function, the derived single-scattering albedo can match either photometric models of Phoebe's dark regolith or brighter sub-surface material excavated by macroscopic impacts on Phoebe.

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.

Dust grains generated by the Uranian irregular satellites will undergo chaotic large-amplitude eccentricity oscillations under the simultaneous action of radiation forces and the highly misaligned quadrupole potentials of the oblate planet and distant Sun. From a suite of orbital histories, we estimate collision proba- bilities of dust particles with the regular satellites and argue that this process may explain the observed hemispherical color asymmetries of the outermost four regular satellites of Uranus.

We present the design and performance of the multi-object fiber spectrographs for the Sloan Digital Sky Survey (SDSS) and their upgrade for the Baryon Oscillation Spectroscopic Survey (BOSS). Originally commissioned in Fall 1999 on the 2.5-m aperture Sloan Telescope at Apache Point Observatory, the spectrographs produced more than 1.5 million spectra for the SDSS and SDSS-II surveys, enabling a wide variety of Galactic and extra-galactic science including the first observation of baryon acoustic oscillations in 2005. The spectrographs were upgraded in 2009 and are currently in use for BOSS, the flagship survey of the third-generation SDSS-III project. BOSS will measure redshifts of 1.35 million massive galaxies to redshift 0.7 and Lyman-alpha absorption of 160,000 high redshift quasars over 10,000 square degrees of sky, making percent level measurements of the absolute cosmic distance scale of the Universe and placing tight constraints on the equation of state of dark energy. The twin multi-object fiber spectrographs utilize a simple optical layout with reflective collimators, gratings, all-refractive cameras, and state-of-the-art CCD detectors to produce hundreds of spectra simultaneously in two channels over a bandpass covering the near ultraviolet to the near infrared, with a resolving power R = \lambda/FWHM ~ 2000. Building on proven heritage, the spectrographs were upgraded for BOSS with volume-phase holographic gratings and modern CCD detectors, improving the peak throughput by nearly a factor of two, extending the bandpass to cover 360 < \lambda < 1000 nm, and increasing the number of fibers from 640 to 1000 per exposure. In this paper we describe the original SDSS spectrograph design and the upgrades implemented for BOSS, and document the predicted and measured performances.

The Iapetus -1:0 nodal bending wave, the first spiral wave ever described in Saturn's rings, has been seen again for the first time in 29 years. We demonstrate that it is in fact the nodal bending wave, not the 1:0 apsidal density wave as previously reported. We use wavelet analysis to determine the wavelength profile, thus deriving the surface density at every point in the region covered by the bending wave. This profile is consistent with surface densities measured from more localized spiral density waves in the outer Cassini Division and the inner and mid-A Ring, varying smoothly from the low values of the former to the higher values of the latter. Most remarkably, our analysis indicates that there is no significant change in surface density across the boundary between the outer Cassini Division and the inner-A ring, despite the very abrupt increase in optical depth and reflected brightness at this location. We consider anew the nature of the classically identified "inner edge of the A ring," given that it does not appear to be correlated with any abrupt increase in surface density. There is an abrupt increase in surface density at the Pandora 5:4 density wave, ~300 km outward of the A ring's inner edge. Further study is needed to robustly interpret our findings in terms of particle properties and abundances, much less to explain the origins of the implied structure.

The Encke Gap is a 320-km-wide opening in Saturn's outer A ring that contains the orbit of the small moon Pan and an array of dusty features composed of particles less than 100 microns across. In particular, there are three narrow ringlets in this region that are not longitudinally homogeneous, but instead contain series of bright clumps. Using images obtained by the Cassini spacecraft, we track the motions of these clumps and demonstrate that they do not follow the predicted trajectories of isolated ring particles moving under the influence of Saturn's and Pan's gravitational fields. We also examine the orbital properties of these ringlets by comparing images taken at different longitudes and times. We find evidence that the orbits of these particles have forced eccentricities induced by solar radiation pressure. In addition, the mean radial positions of the particles in these ringlets appear to vary with local co-rotating longitude, perhaps due to the combined action of drag forces, gravitational perturbations from Pan, and collisions among the ring particles. The dynamics of the dust within this gap therefore appears to be much more complex than previously appreciated.

We consider the Roche critical density (rho_Roche), the minimum density of an orbiting object that, at a given distance from its planet, is able to hold itself together by self-gravity. It is directly related to the more familiar "Roche limit," the distance from a planet at which a strengthless orbiting object of given density is pulled apart by tides. The presence of a substantial ring requires that transient clumps have an internal density less than rho_Roche. Conversely, in the presence of abundant material for accretion, an orbiting object with density greater than rho_Roche will grow. Comparing the rho_Roche values at which the Saturn and Uranus systems transition rapidly from disruption-dominated (rings) to accretion-dominated (moons), we infer that the material composing Uranus' rings is likely more rocky, as well as less porous, than that composing Saturn's rings. From the high values of rho_Roche at the innermost ring-moons of Jupiter and Neptune, we infer that those moons may be composed of denser material than expected, or more likely that they are interlopers that formed farther from their planets and have since migrated inward, now being held together by internal material strength. Finally, the "Portia group" of eight closely-packed Uranian moons has an overall surface density similar to that of Saturn's A ring. Thus, it can be seen as an accretion-dominated ring system, of similar character to the standard ring systems except that its material has a characteristic density greater than the local rho_Roche.

High-inclination circumplanetary orbits that are gravitationally perturbed by the central star can undergo Kozai oscillations---large-amplitude, coupled variations in the orbital eccentricity and inclination. We first study how this effect is modified by incorporating perturbations from the planetary oblateness. Tremaine et al. (2009) found that, for planets with obliquities > 68.875 degrees, orbits in the equilibrium local Laplace plane are unstable to eccentricity perturbations over a finite radial range, and execute large-amplitude chaotic oscillations in eccentricity and inclination. In the hope of making that treatment more easily understandable, we analyze the problem using orbital elements, confirming this threshold obliquity. Furthermore, we find that orbits inclined to the Laplace plane will be unstable over a broader radial range, and that such orbits can go unstable for obliquities less than 68.875 degrees. Finally, we analyze the added effects of radiation pressure, which are important for dust grains and provide a natural mechanism for particle semimajor axes to sweep via Poynting-Robertson drag through any unstable range. We find that generally the effect persists; however, the unstable radial range is shifted and small retrograde particles can avoid the instability altogether. We argue that this is occurs because radiation pressure modifies the equilibrium Laplace plane.

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.

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, \rho~10^(-28)-10^(-27) g/cm^3, with magnetic field strengths of 0.1-1.0 \mu 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.

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.

[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

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.

The Lunar University Network for Astrophysics Research (LUNAR) is a team of researchers and students at leading universities, NASA centers, and federal research laboratories undertaking investigations aimed at using the Moon as a platform for space science. LUNAR research includes Lunar Interior Physics & Gravitation using Lunar Laser Ranging (LLR), Low Frequency Cosmology and Astrophysics (LFCA), Planetary Science and the Lunar Ionosphere, Radio Heliophysics, and Exploration Science. The LUNAR team is exploring technologies that are likely to have a dual purpose, serving both exploration and science. There is a certain degree of commonality in much of LUNAR's research. Specifically, the technology development for a lunar radio telescope involves elements from LFCA, Heliophysics, Exploration Science, and Planetary Science; similarly the drilling technology developed for LLR applies broadly to both Exploration and Lunar Science.

Saturn's diffuse E ring consists of many tiny (micron and sub-micron) grains of water ice distributed between the orbits of Mimas and Titan. Various gravitational and non-gravitational forces perturb these particles' orbits, causing the ring's local particle density to vary noticeably with distance from the planet, height above the ring-plane, hour angle and time. Using remote-sensing data obtained by the Cassini spacecraft in 2005 and 2006, we investigate the E-ring's three-dimensional structure during a time when the Sun illuminated the rings from the south at high elevation angles (> 15 degrees). These observations show that the ring's vertical thickness grows with distance from Enceladus' orbit and its peak brightness density shifts from south to north of Saturn's equator plane with increasing distance from the planet. These data also reveal a localized depletion in particle density near Saturn's equatorial plane around Enceladus' semi-major axis. Finally, variations are detected in the radial brightness profile and the vertical thickness of the ring as a function of longitude relative to the Sun. Possible physical mechanisms and processes that may be responsible for some of these structures include solar radiation pressure, variations in the ambient plasma, and electromagnetic perturbations associated with Saturn's shadow.

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.

The Lunar University Network for Astrophysics Research (LUNAR) undertakes investigations across the full spectrum of science within the mission of the NASA Lunar Science Institute (NLSI), namely science of, on, and from the Moon. The LUNAR team's work on science of and on the Moon, which is the subject of this white paper, is conducted in the broader context of ascertaining the content, origin, and evolution of the solar system.

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.

The leading face of Saturn's moon Iapetus, Cassini Regio, has an albedo only one tenth that on its trailing side. The origin of this enigmatic dichotomy has been debated for over forty years, but with new data, a clearer picture is emerging. Motivated by Cassini radar and imaging observations, we investigate Soter's model of dark exogenous dust striking an originally brighter Iapetus by modeling the dynamics of the dark dust from the ring of the exterior retrograde satellite Phoebe under the relevant perturbations. In particular, we study the particles' probabilities of striking Iapetus, as well as their expected spatial distribution on the Iapetian surface. We find that, of the long-lived particles (greater than about 5 microns), most particle sizes (greater than about 10 microns) are virtually certain to strike Iapetus, and their calculated distribution on the surface matches up well with Cassini Regio's extent in its longitudinal span. The satellite's polar regions are observed to be bright, presumably because ice is deposited there. Thus, in the latitudinal direction we estimate polar dust deposition rates to help constrain models of thermal migration invoked to explain the bright poles (Spencer & Denk 2010). We also analyze dust originating from other irregular outer moons, determining that a significant fraction of that material will eventually coat Iapetus--perhaps explaining why the spectrum of Iapetus' dark material differs somewhat from that of Phoebe. Finally we track the dust particles that do not strike Iapetus, and find that most land on Titan, with a smaller fraction hitting Hyperion. As has been previously conjectured, such exogenous dust, coupled with Hyperion's chaotic rotation, could produce Hyperion's roughly isotropic, moderate-albedo surface.

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.

By volume, more than 99% of the solar system has not been imaged at radio frequencies. Almost all of this space (the solar wind) can be traversed by fast electrons producing radio emissions at frequencies lower than the terrestrial ionospheric cutoff, which prevents observation from the ground. To date, radio astronomy-capable space missions consist of one or a few satellites, typically far from each other, which measure total power from the radio sources, but cannot produce images with useful angular resolution. To produce such images, we require arrays of antennas distributed over many wavelengths (hundreds of meters to kilometers) to permit aperture synthesis imaging. Such arrays could be free-flying arrays of microsatellites or antennas laid out on the lunar surface. In this white paper, we present the lunar option. If such an array were in place by 2020, it would provide context for observations during Solar Probe Plus perihelion passes. Studies of the lunar ionosphere's density and time variability are also important goals. This white paper applies to the Solar and Heliospheric Physics study panel.

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.

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.

We have conducted an intensive search for any material that may orbit Rhea, using images obtained by the Cassini ISS narrow-angle camera. We find no evidence for any such material, contradicting an earlier and surprising inference that Rhea, the second-largest moon of Saturn, possesses a system of narrow rings embedded in a broad circum-satellite disk or cloud (Jones et al. 2008, Science). If one accepts the calculations of the previous authors, our results would require any narrow rings around Rhea to be composed of objects no smaller than 10 meters in radius, in order to be massive enough to account for the observed charged-particle absorptions but diffuse enough to evade detection in our images. Such a particle-size distribution is unrealistic given that erosion processes will break down 10-meter objects to smaller sizes, which would then have been seen in our images. Furthermore, Jones et al. (2008) assumed that the absorption of electrons by hypothetical circum-Rhea particles is proportional to the volume of the particles, even when they are much larger than the electron penetration depth (van Allen 1983, 1987, JGR). If we assume instead that large particles can only absorb electrons efficiently near their surfaces, then our results rule out any solid material orbiting Rhea as the explanation for the absorptions detected by Jones et al. (2008), regardless of particle size, by a margin of four (for narrow rings) or two (for a broad diffuse cloud) orders of magnitude. We conclude that the anomalous electron absorptions reported by Jones et al. (2008) are the signature of a new and hitherto unknown magnetospheric process.

We report the discovery of several large "propeller" moons in the outer part of Saturn's A ring, objects large enough to be followed over the 5-year duration of the Cassini mission. These are the first objects ever discovered that can be tracked as individual moons, but do not orbit in empty space. We infer sizes up to 1--2 km for the unseen moonlets at the center of the propeller-shaped structures, though many structural and photometric properties of propeller structures remain unclear. Finally, we demonstrate that some propellers undergo sustained non-keplerian orbit motion. (Note: This arXiv version of the paper contains supplementary tables that were left out of the ApJL version due to lack of space).

The so-called "Charming Ringlet" (R/2006 S3) is a low-optical-depth, dusty ringlet located in the Laplace gap in the Cassini Division. This ringlet is particularly interesting because its radial position varies systematically with longitude relative to the Sun in such a way that the ringlet's geometric center appears to be displaced away from Saturn's center in a direction roughly toward the Sun. In other words, the ringlet is always found at greater distances from the planet's center at longitudes near the sub-solar longitude than it is at longitudes near Saturn's shadow. This "heliotropic" behavior indicates that the dynamics of the particles in this ring are being influenced by solar radiation pressure. In order to investigate this phenomenon, which has been predicted theoretically but has never been observed this clearly, we analyze multiple image sequences of this ringlet obtained by Cassini in order to constrain its shape and orientation. These data can be fit reasonably well with a model in which both the eccentricity and the inclination of the ringlet have "forced" components (that maintain a fixed orientation relative to the Sun) as well as "free" components (that drift around the planet at steady rates determined by Saturn's oblateness). While the magnitude of the forced eccentricity is roughly consistent with theoretical expectations for radiation pressure acting on 10-to-100-micron-wide icy grains, the existence of significant free eccentricities and inclinations poses a significant challenge for models of low-optical-depth dusty rings.

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.

Aegaeon (Saturn LIII, S/2008 S1) is a small satellite of Saturn that orbits within a bright arc of material near the inner edge of Saturn's G ring. This object was observed in 21 images with Cassini's Narrow-Angle Camera between June 15 (DOY 166), 2007 and February 20 (DOY 51), 2009. If Aegaeon has similar surface scattering properties as other nearby small Saturnian satellites (Pallene, Methone and Anthe), then its diameter is approximately 500 m. Orbit models based on numerical integrations of the full equations of motion show that Aegaeon's orbital motion is strongly influenced by multiple resonances with Mimas. In particular, like the G-ring arc it inhabits, Aegaeon is trapped in the 7:6 corotation eccentricity resonance with Mimas. Aegaeon, Anthe and Methone therefore form a distinctive class of objects in the Saturn system: small moons in co-rotation eccentricity resonances with Mimas associated with arcs of debris. Comparisons among these different ring-arc systems reveal that Aegaeon's orbit is closer to the exact resonance than Anthe's and Methone's orbits are. This could indicate that Aegaeon has undergone significant orbital evolution via its interactions with the other objects in its arc, which would be consistent with the evidence that Aegaeon's mass is much smaller relative to the total mass in its arc than Anthe's and Methone's masses are.

We have developed a semi-analytic method of parameterizing N-body simulations of self-gravity wakes in Saturn's rings, describing their properties by means of only 6 numbers: 3 optical depths and 3 weighting factors. These numbers are obtained using a density-estimation procedure that finds the frequencies of various values of local density within a simulated ring patch. Application of our parameterization to a suite of N-body simulations implies that the distribution of local optical depths is trimodal, rather than bimodal as previous authors have assumed. Rings dominated by self-gravity wakes appear to be mostly empty space. The implications of this result for the analysis of occultation data are more conceptual than practical. The only adjustment needed is that the model parameter $\tau_{gap}$ should be interpreted as representing the area-weighted average optical depth within the gaps (or inter-wake regions). The most significant consequence of our results applies to the question of why "propeller" structures observed in the mid-A ring are seen as relative-bright features, even though the most prominent features of simulated propellers are regions of relatively low density. We find preliminary quantitative support for the hypothesis that propellers would be bright if they involve a local and temporary disruption of self-gravity wakes, flooding the region with more "photometrically active" material (i.e., material that can contribute to the rings' local optical depth) even though the overall density is lower, and thus raising their apparent brightnesses in agreement with observations. We suggest that this mechanism be tested by future detailed numerical models.

The Moon is a unique platform for fundamental astrophysical measurements of gravitation, the Sun, and the Universe. Lacking a permanent ionosphere and, on the farside, shielded from terrestrial radio emissions, a radio telescope on the Moon will be an unparalleled heliospheric and astrophysical observatory. Crucial stages in particle acceleration near the Sun can be imaged and tracked. The evolution of the Universe before and during the formation of the first stars will be traced, yielding high precision cosmological constraints. Lunar Laser Ranging of the Earth-Moon distance provides extremely high precision constraints on General Relativity and alternative models of gravity, and also reveals details about the interior structure of the Moon. With the aim of providing additional perspective on the Moon as a scientific platform, this white paper describes key research projects in these areas of astrophysics from the Moon that are being undertaken by the NLSI-funded LUNAR consortium. The NASA Lunar Science Institute (NLSI) recently funded 7 mostly university-based teams to study science of, on, and from the Moon. The LUNAR consortium was selected by the NLSI for astrophysical research and education that focuses on the key, unique instruments that most effectively take scientific advantage of sites on the lunar surface - low frequency heliophysics and cosmology, and lunar laser ranging. We are submitting this white paper to the Planetary Sciences Decadal Survey to provide additional perspective on the value of Moon for conducting cutting-edge research in astrophysics and gravitational physics by describing our key projects for LUNAR. This program of astrophysics from the Moon complements as well as takes advantage of expected scientific infrastructure on the Moon during the next few decades.

The feasibility of making highly redshifted HI 21-cm (rest frame) measurements from an early epoch of the Universe between the Dark Ages and Reionization (i.e., z>6 and nu<200 MHz) to probe the effects of feedback from the first stars and quasars is assessed in this paper. It may be possible to determine the distribution of hydrogen through the Universe and to constrain the birth of the first stars and black holes via HI tomography. Such observations may also place limits on the properties of Inflation and any exotic heating mechanisms before the first star formation begins (e.g., dark matter decay). The global (all-sky) HI signal after Recombination has distinct features at different frequencies between 30 and 200 MHz that changes as the relative balance between the CMB and spin temperatures changes due to the expansion of the Universe and the ignition of stars and/or black holes. A technology roadmap to approach these observations beginning with ground-based arrays and ending with a low frequency radio array on the lunar farside is described.

Epimetheus, a small moon of Saturn, has a rotational libration (an oscillation about synchronous rotation) of 5.9 +- 1.2 degrees, placing Epimetheus in the company of Earth's Moon and Mars' Phobos as the only natural satellites for which forced rotational libration has been detected. The forced libration is caused by the satellite's slightly eccentric orbit and non-spherical shape. Detection of a moon's forced libration allows us to probe its interior by comparing the measured amplitude to that predicted by a shape model assuming constant density. A discrepancy between the two would indicate internal density asymmetries. For Epimetheus, the uncertainties in the shape model are large enough to account for the measured libration amplitude. For Janus, on the other hand, although we cannot rule out synchronous rotation, a permanent offset of several degrees between Janus' minimum moment of inertia (long axis) and the equilibrium sub-Saturn point may indicate that Janus does have modest internal density asymmetries. The rotation states of Janus and Epimetheus experience a perturbation every four years, as the two moons "swap" orbits. The sudden change in the orbital periods produces a free libration about synchronous rotation that is subsequently damped by internal friction. We calculate that this free libration is small in amplitude (<0.1 degree) and decays quickly (a few weeks, at most), and is thus below the current limits for detection using Cassini images.

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.

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.

A moderate investment of observing time with the International X-ray Observatory to study high-redshift galaxy clusters detected in future large-scale surveys, will provide cosmological measurements of fundamental importance. IXO observations, combined with lensing follow-up, will measure the perturbation growth factor from z=0-2 with an accuracy comparable to, or possibly better than, that expected from observations of cosmic shear with JDEM, and redshift-space distortions with EUCLID. The growth of structure data derived from clusters will significantly improve our knowledge of the dark energy equation of state and will aid in constraining non-GR models for cosmic acceleration. IXO observations of the largest, dynamically relaxed clusters will provide a powerful, independent measurement of the cosmological expansion history using the apparent f_gas(z) trend. Systematic and statistical errors from this technique are competitive with SNIa and BAO studies, making the test extremely useful for improving the accuracy and reliability of the geometric cosmological measurements planned for LSST and JDEM. Only by employing a range of powerful, independent approaches, including those discussed here, can robust answers to puzzles as profound as the origin of cosmic acceleration be expected.

Given the impressive investment by the nation in observational Astronomy and Astrophysics facilities coming on line now and in the near future, we advocate for an increased investment in applied and fundamental research on Astrophysical and Celestial Dynamics (ACD). Specifically we call for a) continued and expanded support for applied research in ACD, b) creation of support for fundamental research in ACD and its subfields, and c) the creation of a unified program to help scientists coordinate and collaborate in their research in these fields. The benefits of this proposal are threefold. First, it will enable researchers to interpret and understand the implications of newly observed phenomena that will invariably arise from new facilities and surveys. Second, research on fundamentals will foster connections between specialists, leveraging advances found in one sub-field and making them available to others. Third, a coordinated approach for applied and fundamental research in ACD will help academic institutions in the United States to produce future researchers trained and knowledgeable in essential subfields such as Mathematical Celestial Mechanics and able to continue its advancement in conjunction with the increase in observations.

For over four decades, synchrotron-radiating sources have played a series of pathfinding roles in the study of galaxy clusters and large scale structure. Such sources are uniquely sensitive to the turbulence and shock structures of large-scale environments, and their cosmic rays and magnetic fields often play important dynamic and thermodynamic roles. They provide essential complements to studies at other wavebands. Over the next decade, they will fill essential gaps in both cluster astrophysics and the cosmological growth of structure in the universe, especially where the signatures of shocks and turbulence, or even the underlying thermal plasma itself, are otherwise undetectable. Simultaneously, synchrotron studies offer a unique tool for exploring the fundamental question of the origins of cosmic magnetic fields. This work will be based on the new generation of m/cm-wave radio telescopes now in construction, as well as major advances in the sophistication of 3-D MHD simulations.

We are learning much about how structure forms, in particular how clusters as nodes in the cosmic web evolve and accrete matter, and about the physical processes within these objects. In the next decade, the study of clusters will enable us to tackle important questions regarding the nature of Dark Matter and Dark Energy, how clusters co-evolve with super-massive black holes at their centers, and to advance our knowledge about fundamental plasma astrophysics. This science white paper outlines the key questions and research opportunities in cluster astrophysics that are emerging in the coming decade and beyond, and serves as an overview to other cluster related white papers.

One of the next decade's most exciting prospects is to explore the cosmic "dark ages," during which the first stars in the Universe formed, with the 21 cm line of neutral hydrogen. At z>6, this light redshifts into the low-frequency radio regime and can be observed with new generations of wide-field arrays. These experiments have the potential to observe enormous volumes and to improve both cosmological and astrophysical constraints dramatically. In this white paper, we describe how the next decade will see the first steps toward unlocking this potential, allowing us to answer two key questions about the fundamental physics of our Universe and the intergalactic medium: Does the standard cosmological model describe the Universe during the "dark ages?" How did the intergalactic medium evolve during this important time, ending with the reionization of hydrogen?

The cosmic dark ages and the epoch of reionization, during which the first generations of stars and galaxies formed, are among the most compelling frontiers of extragalactic astrophysics and cosmology. Here we describe an exciting new probe of these eras: the 21 cm line of neutral hydrogen, which will allow us to map the neutral intergalactic medium throughout the era of first galaxy formation. In this white paper, we describe how these studies can answer two key questions about galaxy formation: What were the properties of high-redshift galaxies? How did they affect the Universe around them?

We present new results characterizing cosmological shocks within adaptive mesh refinement N-Body/hydrodynamic simulations that are used to predict non-thermal components of large-scale structure. This represents the first study of shocks using adaptive mesh refinement. We propose a modified algorithm for finding shocks from those used on unigrid simulations that reduces the shock frequency of low Mach number shocks by a factor of ~3. We then apply our new technique to a large, (512 Mpc/h)^3, cosmological volume and study the shock Mach number (M) distribution as a function of pre-shock temperature, density, and redshift. Because of the large volume of the simulation, we have superb statistics that results from having thousands of galaxy clusters. We find that the Mach number evolution can be interpreted as a method to visualize large-scale structure formation. Shocks with Mach<5 typically trace mergers and complex flows, while 5<Mach<20 and Mach>20 generally follow accretion onto filaments and galaxy clusters, respectively. By applying results from nonlinear diffusive shock acceleration models using the first-order Fermi process, we calculate the amount of kinetic energy that is converted into cosmic ray protons. The acceleration of cosmic ray protons is large enough that in order to use galaxy clusters as cosmological probes, the dynamic response of the gas to the cosmic rays must be included in future numerical simulations.

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.

We present an extensive data set of ~150 localized features from Cassini images of Saturn's Ring A, a third of which are demonstrated to be persistent by their appearance in multiple images, and half of which are resolved well enough to reveal a characteristic "propeller" shape. We interpret these features as the signatures of small moonlets embedded within the ring, with diameters between 40 and 500 meters. The lack of significant brightening at high phase angle indicates that they are likely composed primarily of macroscopic particles, rather than dust. With the exception of two features found exterior to the Encke Gap, these objects are concentrated entirely within three narrow (~1000 km) bands in the mid-A Ring that happen to be free from local disturbances from strong density waves. However, other nearby regions are similarly free of major disturbances but contain no propellers. It is unclear whether these bands are due to specific events in which a parent body or bodies broke up into the current moonlets, or whether a larger initial moonlet population has been sculpted into bands by other ring processes.

Flux-limited X-ray samples indicate that about half of rich galaxy clusters have cool cores. Why do only some clusters have cool cores while others do not? In this paper, cosmological N-body + Eulerian hydrodynamic simulations, including radiative cooling and heating, are used to address this question as we examine the formation and evolution of cool core (CC) and non-cool core (NCC) clusters. These adaptive mesh refinement simulations produce both CC and NCC clusters in the same volume. They have a peak resolution of 15.6 h^{-1} kpc within a (256 h^{-1} Mpc)^3 box. Our simulations suggest that there are important evolutionary differences between CC clusters and their NCC counterparts. Many of the numerical CC clusters accreted mass more slowly over time and grew enhanced cool cores via hierarchical mergers; when late major mergers occurred, the CC's survived the collisions. By contrast, NCC clusters experienced major mergers early in their evolution that destroyed embryonic cool cores and produced conditions that prevented CC re-formation. As a result, our simulations predict observationally testable distinctions in the properties of CC and NCC beyond the core regions in clusters. In particular, we find differences between CC versus NCC clusters in the shapes of X-ray surface brightness profiles, between the temperatures and hardness ratios beyond the cores, between the distribution of masses, and between their supercluster environs. It also appears that CC clusters are no closer to hydrostatic equilibrium than NCC clusters, an issue important for precision cosmology measurements.

We present the first results from a new generation of 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 simulated a very large (512^3h^{-3}Mpc^3) volume with unprecedented dynamic range. We have generated simulated light cones to match the resolution and sensitivity of current and future SZE instruments. Unlike many previous studies of this type, our simulation includes unbound gas, where an appreciable fraction of the baryons in the universe reside. We have found that cluster line-of-sight overlap may be a significant issue in upcoming single-dish SZE surveys. Smaller beam surveys (~1 arcmin) have more than one massive cluster within a beam diameter 5-10% of the time, and a larger beam experiment like Planck has multiple clusters per beam 60% of the time. We explore the contribution of unresolved halos and unbound gas to the SZE signature at the maximum decrement. We find that there is a contribution from gas outside clusters of ~16% per object on average for upcoming surveys. This adds both bias and scatter to the deduced value of the integrated SZE, adding difficulty in accurately calibrating a cluster Y-M relationship. Finally, we find that in images where objects with M > 5x10^{13} M_{\odot} have had their SZE signatures removed, roughly a third of the total SZE flux still remains. This gas exists at least partially in the Warm Hot Intergalactic Medium (WHIM), and will possibly be detectable with the upcoming generation of SZE surveys.

We have analyzed a large sample of numerically simulated clusters to demonstrate the adverse effects resulting from use of X-ray fitted beta-model parameters with Sunyaev-Zeldovich effect (SZE) data. There is a fundamental incompatibility between model fits to X-ray surface brightness profiles and those done with SZE profiles. Since observational SZE radial profiles are in short supply, the X-ray parameters are often used in SZE analysis. We show that this leads to biased estimates of the integrated Compton y-parameter inside r_{500} calculated from clusters. We suggest a simple correction of the method, using a non-isothermal beta-model modified by a universal temperature profile, which brings these calculated quantities into closer agreement with the true values.

Observations at radio wavelengths address key problems in astrophysics, astrobiology, and lunar structure including the first light in the Universe (the Epoch of Reionization), the presence of magnetic fields around extrasolar planets, particle acceleration mechanisms, and the structure of the lunar ionosphere. Moreover, achieving the performance needed to address these scientific questions demands observations at wavelengths longer than those that penetrate the Earth's ionosphere, observations in extremely "radio quiet" locations such as the Moon's far side, or both. We describe a series of lunar-based radio wavelength interferometers of increasing capability. The Radio Observatory for Lunar Sortie Science (ROLSS) is an array designed to be deployed during the first lunar sorties (or even before via robotic rovers) and addressing particle acceleration and the lunar ionosphere. Future arrays would be larger, more capable, and deployed as experience is gained in working on the lunar surface.

We describe a powerful signal processing method, the continuous wavelet transform, and use it to analyze radial structure in Cassini ISS images of Saturn's rings. Wavelet analysis locally separates signal components in frequency space, causing many structures to become evident that are difficult to observe with the naked eye. Density waves, generated at resonances with saturnian satellites orbiting outside (or within) the rings, are particularly amenable to such analysis. We identify a number of previously unobserved weak waves, and demonstrate the wavelet transform's ability to isolate multiple waves superimposed on top of one another. We also present two wave-like structures that we are unable to conclusively identify. In a multi-step semi-automated process, we recover four parameters from clearly observed weak spiral density waves: the local ring surface density, the local ring viscosity, the precise resonance location (useful for pointing images, and potentially for refining saturnian astrometry), and the wave amplitude (potentially providing new constraints upon the masses of the perturbing moons). Our derived surface densities have less scatter than previous measurements that were derived from stronger non-linear waves, and suggest a gentle linear increase in surface density from the inner to the mid-A Ring. We show that ring viscosity consistently increases from the Cassini Division outward to the Encke Gap. Meaningful upper limits on ring thickness can be placed on the Cassini Division (3.0 m at r~118,800 km, 4.5 m at r~120,700 km) and the inner A Ring (10 to 15 m for r<127,000 km).

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.

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.

We describe a model that accounts for the complex morphology of spiral density waves raised in Saturn's rings by the co-orbital satellites, Janus and Epimetheus. Our model may be corroborated by future Cassini observations of these time-variable wave patterns.

Withdrawn at the current time due to errors.

We critically analyze the measurement of galaxy cluster gas masses, which is central to cosmological studies that rely on the galaxy cluster gas mass fraction. Using synthetic observations of numerically simulated clusters viewed through their X-ray emission and thermal Sunyaev-Zeldovich effect (SZE), we reduce the observations to obtain measurements of the cluster gas mass. We are thus able to quantify the possible sources of uncertainty and systematic bias associated with the common simplifying assumptions used in reducing real cluster observations including isothermality and hydrostatic equilibrium. We find that intrinsic variations in clusters limit the precision of observational gas mass estimation to ~10% to 1 sigma confidence excluding instrumental effects. Gas mass estimates performed via all methods surprisingly show little or no trending in the scatter as a function of cluster redshift. For the full cluster sample, methods that use SZE profiles out to roughly the virial radius are the simplest, most accurate, and unbiased way to estimate cluster mass. X-ray methods are systematically more precise mass estimators than are SZE methods if merger and cool core systems are removed, but X-ray methods slightly overestimate (5-10%) the cluster gas mass on average. We find that cool core clusters in our samples are particularly poor candidates for observational mass estimation, even when excluding emission from the core region. The effects of cooling in the cluster gas alter the radial profile of the X-ray and SZE surface brightness outside the cool core region, leading to poor gas mass estimates in cool core clusters. Finally, we find that methods using a universal temperature profile estimate cluster masses to higher precision than those assuming isothermality.

We investigate empirical scaling relations between the thermal Sunyaev-Zeldovich effect (SZE) and cluster mass in simulated clusters of galaxies. The simulated clusters have been compiled from four different samples that differ only in their assumed baryonic physics. We show that the strength of the thermal SZE integrated over a significant fraction of the virialized region of the clusters is relatively insensitive to the detailed heating and cooling processes in the cores of clusters by demonstrating that the derived scaling relations are nearly identical between the four cluster samples considered. For our synthetic images, the central Comptonization parameter shows significant boosting during transient merging events, but the integrated SZE appears to be relatively insensitive to these events. Most importantly, the integrated SZE closely tracks the underlying cluster mass. Observations through the thermal SZE allow a strikingly accurate mass estimation from relatively simple measurements that do not require either parametric modeling or geometric deprojection and thus avoid assumptions regarding the physics of the ICM or the symmetry of the cluster. This result offers significant promise for precision cosmology using clusters of galaxies.

Although analytical studies on the secular motion of the irregular satellites have been published recently, these theories have not yet been satisfactorily reconciled with the results of direct numerical integrations. These discrepancies occur because in secular theories the disturbing function is averaged over orbital motions, whereas instead one should take into account some large periodic terms, most notably the so-called ``evection''. We demonstrate that such terms can be incorporated into the Kozai formalism, and that our synthetic approach produces much better agreement with results from symplectic integrations. Using this method, we plot the locations of secular resonances in the orbital-element space, and we note that the distribution of irregular satellite clusters appears to be non-random. We find that the large majority of irregular-satellite groups cluster close to the secular resonances, with several objects having practically stationary pericenters. None of the largest satellites belong to this class, so we argue that this dichotomy implies that the smaller near-resonant satellites might have been captured differently than the largest irregulars.

Among the many new irregular satellites that have been discovered in the last five years, at least six are in the so-called Kozai resonance. Due to solar perturbations, the argument of pericenter of a satellite usually precesses from 0 to 360 degrees. However, at inclinations higher than 39.3 degrees and lower than 140.7 degrees a new kind of behavior occurs for which the argument of pericenter oscillates around +/-90 degrees. In this work we will concentrate on the orbital history of the saturnian satellite S/2000 S5 Kiviuq, one of the satellites currently known to be in such resonance Kiviuq's orbit is very close to the separatrix of the Kozai resonance. Due to perturbations from the other jovian planets, it is expected that orbits near the Kozai separatrix may show significant chaotic behavior. This is important because chaotic diffusion may transfer orbits from libration to circulation, and vice versa. To identify chaotic orbits we used two well-known methods: the Frequency Analysis Method (Laskar 1990) and Maximum Lyapunov Exponents (Benettin et al. 1980). Our results show that the Kozai resonance is crossed by a web of secondary resonances, whose arguments involve combinations of the argument of pericenter, the argument of the Great Inequality, longitude of the node, and other terms related to the secular frequencies g5, g6, and s6. Many test orbits whose precession period is close to the period of the Great Inequality (883 yrs), or some of its harmonics, are trapped by these secondary resonances, and show significant chaotic behavior. Planetary migration, by moving the locations of these secondary resonances, may have depleted an original population of Kozai resonators.

We present an analysis of the X-ray emission for a complete sample of 288 Abell clusters spanning the redshift range 0.016<= z <= 0.09 from the ROSAT All-Sky Survey. This sample is based on our 20cm VLA survey of nearby Abell clusters. We find an X-ray detection rate of 83%. We report cluster X-ray fluxes and luminosities and two different flux ratios indicative of the concentration and extent of the emission. We examine correlations between the X-ray luminosity, Abell Richness, and Bautz-Morgan and Rood-Sastry cluster morphologies. We find a strong correlation between Lx and cluster richness coupled to a dependence on the optical morphological type. These results are consistent with the observed scatter between X-ray luminosity and temperature and a large fraction of cooling flows. For each cluster field we also report the positions, peak X-ray fluxes, and flux-ratios of all X-ray peaks above 3-sigma significance within a box of 2x2 Mpc centered on Abell's position.

We present results from recent simulations of the formation and evolution of clusters of galaxies in a LambdaCDM cosmology. These simulations contain our most physically complete input physics to date including radiative cooling, star formation that transforms rapidly cooling material into aggregate star particles and we also model the thermal feedback from resulting supernovae in the star particles. We use an adaptive mesh refinement (AMR) Eulerian hydrodynamics scheme to obtain very high spatial resolution (~ 2 kpc) in a computational volume 256 Mpc on a side with mass resolution for dark matter and star particles of ~ 10^8 M_solar. We examine in detail the appearance and evolution of the core region of our simulated clusters.

We present a new model for the creation of cool cores in rich galaxy clusters within a LambdaCDM cosmological framework using the results from high spatial dynamic range, adaptive mesh hydro/N-body simulations. It is proposed that cores of cool gas first form in subclusters and these subclusters merge to create rich clusters with cool, central X-Ray excesses. The rich cool clusters do not possess ``cooling flows'' due to the presence of bulk velocities in the intracluster medium in excess of 1000 km/sec produced by on-going accretion of gas from supercluster filaments. This new model has several attractive features including the presence of substantial core substructure within the cool cores, and it predicts the appearance of cool bullets, cool fronts, and cool filaments all of which have been recently observed with X-Ray satellites. This hierarchical formation model is also consistent with the observation that cool cores in Abell clusters occur preferentially in dense supercluster environments. On the other hand, our simulations overproduce cool cores in virtually all of our numerical clusters, the central densities are high, and physical core temperatures are often below 1 keV (in contrast to recent observations). We will discuss additional preliminary simulations to ``soften'' the cool cores involving star formation and supernova feedback.

We present a new scenario for the formation of cool cores in rich galaxy clusters based on results from recent high spatial dynamic range, adaptive mesh Eulerian hydrodynamic simulations of large-scale structure formation. We find that cores of cool gas, material that would be identified as a classical cooling flow based on its X-ray luminosity excess and temperature profile, are built from the accretion of discrete, stable subclusters. Any ``cooling flow'' present is overwhelmed by the velocity field within the cluster - the bulk flow of gas through the cluster typically has speeds up to about 2,000 km s^-1 and significant rotation is frequently present in the cluster core. The inclusion of consistent initial cosmological conditions for the cluster within its surrounding supercluster environment is crucial when simulating the evolution of cool cores in rich galaxy clusters. This new model for the hierarchical assembly of cool gas naturally explains the high frequency of cool cores in rich galaxy clusters despite the fact that a majority of these clusters show evidence of substructure which is believed to arise from recent merger activity. Furthermore, our simulations generate complex cluster cores in concordance with recent X-ray observations of cool fronts, cool ``bullets'', and filaments in a number of galaxy clusters. Our simulations were computed with a coupled N-body, Eulerian, adaptive mesh refinement, hydrodynamics cosmology code that properly treats the effects of shocks and radiative cooling by the gas. We employ up to seven levels of refinement to attain a peak resolution of 15.6 h^-1 kpc within a volume 256 h-1 Mpc on a side and assume a standard LambdaCDM cosmology.

We investigate the predicted present-day temperature profiles of the hot, X-ray emitting gas in galaxy clusters for two cosmological models - a current best-guess LCDM model and standard cold dark matter (SCDM). Our numerically-simulated "catalogs" of clusters are derived from high-resolution (15/h kpc) simulations which make use of a sophisticated, Eulerian-based, Adaptive Mesh-Refinement (AMR) code that faithfully captures the shocks which are essential for correctly modelling cluster temperatures. We show that the temperature structure on Mpc-scales is highly complex and non-isothermal. However, the temperature profiles of the simulated LCDM and SCDM clusters are remarkably similar and drop-off as $T +AFw-propto (1+-r/a_x)^{-+AFw-delta}$ where $a_x +AFw-sim r_{vir}/1.5$ and $+AFw-delta +AFw-sim 1.6$. This decrease is in good agreement with the observational results of Markevitch et al.(1998) but diverges, primarily in the innermost regions, from their fit which assumes a polytropic equation of state. Our result is also in good agreement with a recent sample of clusters observed by BeppoSAX though there is some indication of missing physics at small radii ($r<0.2 r_{vir}$). We discuss the interpretation of our results and make predictions for new x-ray observations that will extend to larger radii than previously possible. Finally, we show that, for $r>0.2 r_{vir}$, our universal temperature profile is consistent with our most recent simulations which include both radiative cooling and supernovae feedback.

We present new redshifts and positions for 635 galaxies in nine rich clusters containing Wide-Angle Tailed (WAT) radio galaxies. Combined with existing data, we now have a sample of 18 WAT-containing clusters with more than 10 redshifts. This sample contains a substantial portion of the WAT clusters in the VLA 20 cm survey of Abell clusters, including 75% of WAT clusters in the complete survey (z<=0.09), and 20% of WAT clusters with z>0.09. It is a representative sample which should not contain biases other than selection by radio morphology. We graphically present the new data using histograms and sky maps. A semi-automated procedure is used to search for emission lines in the spectra in order to add and verify galaxy redshifts. We find that the average apparent fraction of emission line galaxies is about 9% in both the clusters and the field. We investigate the magnitude completeness of our redshift surveys with CCD data for a test case, Abell 690. This case indicates that our galaxy target lists are deeper than the detection limit of a typical MX exposure, and they are 82% complete down to R=19.0. The importance of the uniformity of the placement of fibers on targets is posited, and we evaluate this in our datasets. We find some cases of non-uniformities which may influence dynamical analyses. A second paper will use this database to look for correlations between the WAT radio morphology and the cluster's dynamical state.

We report the results of recent numerical simulations of the head-on merger of a cooling flow cluster with an infalling subcluster of galaxies. These simulations examined the effects of different types of cluster mergers (with 16:1 and 4:1 mass ratios) on the evolution of cluster cooling flows (mass accretion rates of 100 and 400 M/year. The 2-dimensional simulations were performed with a Hydro/N-body code on a uniform grid with a resolution of 20 kpc. We find that the ram-pressure of the infalling gas is crucial in determining the fate of the cooling flow as disruption occurs when a substantial amount of subcluster gas reaches the primary's core. In such cases, the subcluster gas can increase the central cooling time by displacing the high-density cooling gas and by heating it via shocks and turbulent gas motions. Moreover, the fate of a merging cooling flow is also dependent on its initial cooling time. In cases where the initial cooling time is very small then, even if the flow is disrupted, the central cooling time will remain less than a Hubble time and the flow will likely re-establish itself. This has an important observational consequence because such clusters will be classified as cooling flows even though they have experienced a significant merger. In addition, there is a time delay (1-2 Gyrs) between core-crossing and the point at which the central cooling time of a disrupted flow becomes of order a Hubble time (abridged).

We have simulated the evolution of a large sample of X-ray clusters in a Lambda-CDM universe at high spatial resolution using adaptive mesh refinement and placed the results in an online archive for public access. The Simulated Cluster Archive website http://sca.ncsa.uiuc.edu provides tools for interactive 2D and 3D analysis of gas and dark matter fields, X-ray and SZ imaging, and data export. We encourage community use and solicit their feedback.

We present high resolution ROSAT x-ray and radio observations of three cooling flow clusters containing steep spectrum radio sources at their cores. All three systems exhibit strong signs of interaction between the radio plasma and the hot intracluster medium. Two clusters, A133 and A2626, show enhanced x-ray emission spatially coincident with the radio source whereas the third cluster, A2052, exhibits a large region of x-ray excess surrounding much of the radio source. Using 3-D numerical simulations, we show that a perturbed jet propagating through a cooling flow atmosphere can give rise to amorphous radio morphologies, particularly in the case where the jet was ``turned off'' and allowed to age passively. In addition, the simulated x-ray surface brightness produced both excesses and deficits as seen observationally.

New adaptive mesh refinement N-body + hydrodynamics numerical simulations are used to illustrate the complex and changing cluster environments in which many radio galaxies live and evolve. Groups and clusters of galaxies form at the intersections of filaments where they continue to accrete gas and dark matter to the present day. The accretion process produces shocks, turbulence, and transonic bulk flows forming a kind of stormy weather within the intracluster medium (ICM). Radio sources embedded within the stormy ICM form distorted, complex morphologies as observed in recent VLA cluster surveys. We show that the bending of wide-angle tailed radio sources can be understood as the result of recent cluster-subcluster mergers. We use new MHD simulations to illustrate how cluster radio halos can be formed by the shocks and turbulence produced during cluster mergers. Finally, we discuss new observations of distant Abell clusters that reveal a class of weak radio sources, probably starbursts, likely produced during the formation of the clusters as they accrete material from the supercluster environment.

Investigations of the cluster environment of radio sources have not shown a correlation between radio power and degree of clustering. However, it has been demonstrated that extended X-ray luminosity and galaxy clustering do exhibit a positive correlation. This study investigates a complete sample of 25 nearby (z less than 0.06) radio galaxies which are not cataloged members of Abell clusters. The environment of these radio galaxies is studied in both the X-ray and the optical by means of the ROSAT All-Sky Survey (RASS), ROSAT pointed observations, and the Palomar optical Digitized Sky Survey (DSS). X-ray luminosities and extents are determined from the RASS, and the DSS is used to quantify the degree of clustering via the spatial two-point correlation coefficient, Bgg. Of the 25 sources, 20 are greater than sigma detections in the X-ray and 11 possessed Bgg's significantly in excess of that expected for an isolated galaxy. Adding the criterion that the X-ray emission be resolved, 10 of the radio galaxies do appear to reside in poor clusters with extended X-ray emission suggestive of the presence of an intracluster medium. Eight of these galaxies also possess high spatial correlation coefficients. Taken together, these data suggest that the radio galaxies reside in a low richness extension of the Abell clusters. The unresolved X-ray emission from the other galaxies is most likely associated with AGN phenomena. Furthermore, although the sample size is small, it appears that the environments of FR I and FR II sources differ. FR I's tend to be more frequently associated with extended X-ray emission (10 of 18), whereas FR II's are typically point sources or non-detections in the X-ray (none of the 7 sources exhibit extended X-ray emission).

A catalog of 732 optically selected, nearby poor clusters of galaxies covering the entire sky north of $\rm -3^{\circ}$ declination is presented. The poor clusters, called WBL clusters, were identified as concentrations of 3 or more galaxies with photographic magnitudes brighter than 15.7, possessing a galaxy surface overdensity of $10^{4/3}$. These criteria are consistent with those used in the identification of the original Yerkes poor clusters, and this new catalog substantially increases the sample size of such objects. These poor clusters cover the entire range of galaxy associations up to and including Abell clusters, systematically including poor and rich galaxy systems spanning over three orders of magnitude in the cluster mass function. As a result, this new catalog contains a greater diversity of richness and structures than other group catalogs, such as the Hickson or Yerkes catalogs. The information on individual galaxies includes redshifts and cross-references to other galaxy catalogs. The entries for the clusters include redshift (where available) and cross-references to other group and cluster catalogs.

In this letter, we present a new determination of the local (z<0.09) X-ray luminosity function (XLF) using a large, statistical sample of 294 Abell clusters and the ROSAT All-Sky-Survey. Given our large sample size, we have reduced errors by a factor of two for L(X)(0.5-2keV)>10^43 ergs/sec. We combine our data with previous work in order to explore possible constraints imposed by the shape of the XLF on cosmological models. A set of currently viable cosmologies is used to construct theoretical XLFs assuming Lx is proportional to M^p and a sigma_8-Omega_0 constraint (from Viana & Liddle 1996) based on the local X-ray temperature function. We fit these models to our observed XLF and verify that the simplest adiabatic, analytic scaling relation (e.g. Kaiser 1986) disagrees strongly with observations. If we assume that clusters can be described by the pre-heated, constant core-entropy models of Evrard & Henry (1991) then the observed XLF is consistent only with 0.1 < Omega_0 < 0.4 if the energy per unit mass in galaxies is roughly equal to the gas energy (ie if beta=1). (abridged)

We present a numerical model for the extended steep-spectrum radio sources and the elongated X-ray structure in A3667 based on new 3-dimensional MHD/N-body simulations. The X-ray and optical analyses of A3667 indicate that it has undergone a recent subcluster merger event. We believe that the Mpc-scale radio sources identified in A3667 are also a consequence of the merger. Our previous numerical simulations show that mergers often produce large-scale shocks and turbulence capable of both magnetic field amplification and in-situ reacceleration of relativistic particles. Our model suggests that these radio structures, separated by ~2.6 h^{-1}_100 Mpc, are in fact causally linked via a slightly off-axis merger that occurred nearly in the plane of the sky approximately 1 Gyr ago with a subcluster having a total mass equal to ~20% of the primary cluster.

We present initial results from the first 3-dimensional numerical magnetohydrodynamical (MHD) simulations of magnetic field evolution in merging clusters of galaxies. Within the framework of idealized initial conditions similar to our previous work, we look at the gasdynamics and the magnetic field evolution during a major merger event in order to examine the suggestion that shocks and turbulence generated during a cluster/subcluster merger can produce magnetic field amplification and relativistic particle acceleration and, as such, may play a role in the formation and evolution of cluster-wide radio halos. The ICM, as represented by the equations of ideal MHD, is evolved self-consistently within a changing gravitational potential defined largely by the collisionless dark matter component represented by an N-body particle distribution. The MHD equations are solved by the Eulerian, finite-difference code, ZEUS. The particles are evolved by a standard particle-mesh (PM) code. We find significant evolution of the magnetic field structure and strength during two distinct epochs of the merger evolution.

Within Abell galaxy clusters containing wide-angle tailed radio sources, there is evidence of a ``prevailing wind'' which directs the WAT jets. We study the alignment of WAT jets and nearby clusters to test the idea that this wind may be a fossil of drainage along large-scale supercluster axes. We also test this idea with a study of the alignment of WAT jets and supercluster axes. Statistical test neighbours indicate no alignment of WAT jets towards nearest clusters, but do indicate approximately 98% confidence in alignment with the long axis of the supercluster in which the cluster lies. We find a preferred scale for such superclusters of order 25 Mpc $h^{-1}$.

We present ASCA X-ray observations of a sample of five groups selected from a cross-correlation of the ROSAT All-Sky Survey with the White et al. optical catalog of groups. These X-ray bright groups significantly increase the number of known systems with temperatures between 2 and 3 keV. They have element abundances of roughly 0.3 times solar that are typical of clusters, but their favored ratio of Si/Fe abundance is lower than the cluster value. Combining the ASCA results with ROSAT imaging data, we calculate total masses of a few to several times 10^13 M_sun, gas mass fractions of ~10%, and baryonic mass fractions of at least 15-20% within a radius of 0.5 Mpc. Upper limits for the ratios of gas to galaxy mass and of the iron mass to galaxy luminosity overlap with the range observed in rich clusters and extend to lower values, but not to such low values as seen in much poorer groups. These results support the idea that groups, unlike clusters, are subject to the loss of their primordial and processed gas, and show that this transition occurs at the mass scale of the 2-3 keV groups. A discussion of ASCA calibration issues and a comparison of ROSAT and ASCA temperatures are included in an Appendix.

We have examined the ROSAT PSPC X-ray properties of a sample of 15 Abell clusters containing 23 narrow-angle tailed (NAT) radio galaxies. We find that clusters with NATs show a significantly higher level of substructure than a similar sample of radio-quiet clusters, indicating that NAT radio sources are preferentially located in dynamically complex systems. Also, the velocity distribution of the NAT galaxies is similar to that of other cluster members; these velocities are inadequate for producing the ram pressure necessary to bend the radio jets. We therefore propose a new model for NAT formation, in which NATs are associated with dynamically complex clusters undergoing merger events. The U-shaped NAT morphology is produced in part by the merger-induced bulk motion of the ICM bending the jets.

We have conducted pencil-beam searches for outer solar system objects to a limiting magnitude of R ~ 26. Five new trans-neptunian objects were detected in these searches. Our combined data set provides an estimate of ~90 trans-neptunian objects per square degree brighter than ~ 25.9. This estimate is a factor of 3 above the expected number of objects based on an extrapolation of previous surveys with brighter limits, and appears consistent with the hypothesis of a single power-law luminosity function for the entire trans-neptunian region. Maximum likelihood fits to all self-consistent published surveys with published efficiency functions predicts a cumulative sky density Sigma(<R) obeying log10(Sigma) = 0.76(R-23.4) objects per square degree brighter than a given magnitude R.

We present the results of our multiwavelength study of two nearby Abell clusters that contain extended tailed radio sources. From our analysis of PSPC X-ray, VLA radio, and optical data, we find strong evidence that these clusters are not relaxed. Significant X-ray substructure is revealed by examining their X-ray surface brightness maps. We report 130 new redshifts from the fields of A578 and A1569 and detect the presence of substructure in the velocity and spatial distribution of galaxies. Furthermore, we find that the radio sources in A578 are at rest with respect to their subclusters and that the Wide-Angled Tailed (WAT) source in A1569 has a small peculiar motion with respect to its subcluster. The standard model of ram pressure induced jet curvature requires relative velocities of the radio sources with respect to the ICM of > 1000 km/s. Therefore, the apparent peculiar velocities of the radio galaxies can not explain the bending of the jets/tails for the Narrow-Angled Tailed (NAT) radio source in A578 and for the WAT in A1569. We suggest that a cluster- subcluster merger model provides a reasonable explanation for the observed X-ray and optical substructure. Numerical simulations of cluster-subcluster mergers predict a bulk flow of ICM gas that has sufficient dynamic pressure for bending the extended radio sources. Thus, a cluster-subcluster merger can provide the environment responsible for the observed cluster X-ray and optical morphologies and for the bending seen in the radio sources.

The radio continuum properties of galaxy clusters with cooling flows are reviewed along with the relationship to the X-ray environment. We find that 60-70% of cD galaxies in cooling flows are radio-loud, a much higher fraction than the 14% found for typical cluster ellipticals. New ROSAT HRI observations reveal a variety of interesting correlations and anticorrelations between the X-ray structure within the inner cooling flows and the radio morphologies. It appears that the radio plasma can have a strong effect on the inner structure of cluster cooling cores as virtually all cooling flow clusters with central radio sources have non-symmetric X-ray structure. Numerical simulations of radio jets in cooling flow atmospheres are presented. We also discuss the prospects for destroying cooling flows via cluster-cluster mergers, using new hydro/N-body simulations which incorporate radiative cooling.

We present an observational analysis of numerical simulations of galaxy cluster mergers. We identify several observational signatures of recent merger activity, and quantitatively assess the uncertainty introduced into cluster mass estimates when invoking the commonly held assumptions of hydrostatic equilibrium, virial equilibrium, spherical symmetry and isothermality. We find that mergers result in multiple X-ray peaks, long-lived elongation of the X-ray emission as well as isophotal twisting and centroid shifting to a degree consistent with recent observations. We also find an enlargement of the X-ray core relative to the dark matter core. Mergers result in non-isothermal clusters exhibiting observable inhomogeneities in the emission-weighted X-ray temperature of several keV on linear scales of less than 0.5 Mpc. The resulting gas dynamics are extremely complex, and we present an example of what might be observed by a high resolution X-ray spectrograph. We further speculate that the gas dynamics, via shocks, bulk flows and turbulence, play an important role in the evolution of cluster galaxies and associated radio sources, particularly wide-angle tailed (WAT) sources and radio halos. We find that X-ray based by cluster mass estimates made under equilibrium assumptions can be uncertain 50\% or more in the first 2 Gyrs after a merger and by up to 25\% after 2 Gyrs depending on the details of the analysis and projection effects. Uncertainties can be considerably larger if the temperature is not well constrained.

From the poor cluster catalog of White et al. (1996), we define a sample of 71 optically-selected poor galaxy clusters. The surface-density enhance- ment we require for our clusters falls between that of the loose associations of Turner and Gott (1976) and the Hickson compact groups (Hickson, 1982). We review the selection biases and determine the statistical comleteness of the sample. For this sample, we report new velocity measurements made with the ARC 3.5-m Dual-Imaging spectrograph and the 2.3-m Steward Observatory MX fiber spectrograph. Combining our own measurements with those from the literature, we examine the velocity distributions, velocity dispersions, and 1-d velocity substructure for our poor cluster sample, and compare our results to other poor cluster samples. We find that approximately half of the sample may have significant 1-d velocity substructure. The optical morphology, large-scale environment, and velocity field of many of these clusters is indicative of young, dynamically evolving systems. In future papers, we will use this sample to derive the poor cluster X-ray luminosity function and gas mass function (see astro-ph/9606120), and will examine the optical/X-ray properties of the clusters in more detail.

We present the first X-ray Luminosity Function (XLF) for an optically-selected sample of 49 nearby poor clusters of galaxies and a sample of 67 Abell clusters with z < 0.15. We have extended the measured cluster XLF by more than a factor of 10 in X-ray luminosity. Our poor cluster sample was drawn from an optical catalog of groups with 0.01<z<0.03 composed of Zwicky galaxies. The X-ray emission was measured from the ROSAT all-sky survey. About 45% of the poor clusters were detected with 0.5-2.0 keV luminosities from (1.7-65)x 10^{41} h^{-2} ergs/sec. These are among the X-ray brightest, optically-selected poor clusters in the northern hemisphere. For this sample, the poor cluster XLF was found to be smooth extrapolation of the rich cluster XLF. A new Hydro/N-body simulation of a Hot +Cold dark matter model with \Omega_{total}=1, \Omega_{\nu}=0.2, and a baryon fraction of 7.5% was used to model and understand our observational selection effects. We found that the observed cluster Gas Mass Function was consistent with our model.

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 \delta < 74^{\circ}$, 0$^h \leq \alpha < 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 $\gamma = -0.8$. This may well provide the $first$ detection of an angular correlation in a large area, complete deep radio survey.

We present the first self-consistently computed results on the evolution of X-ray properties of galaxy clusters in a Cold + Hot Dark Matter (CHDM) model. We have performed a hydrodynamic plus N-body simulation for the COBE-compatible CHDM model with standard mass components: Omega(hot) = 0.3, Omega(cold) = 0.6 and Omega(baryon) = 0.1 (h = 0.5). In contrast with the CDM model, which fails to reproduce the observed temperature distribution function dN/dT (Bryan et al. 1994b), the CHDM model fits the observational dN/dT quite well. Our results on X-ray luminosity are less firm but even more intriguing. We find that the resulting X-ray luminosity functions at redshifts z = 0.0, 0.2, 0.4, 0.7 are well fit by observations, where they overlap. The fact that both temperatures and luminosities provide a reasonable fit to the available observational data indicates that, unless we are missing some essential physics, there is neither room nor need for a large fraction of gas in rich clusters: 10% (or less) in baryons is sufficient to explain their X-ray properties. We also see a tight correlation between X-ray luminosity and gas temperature.