Quark nuggets are theoretical objects composed of approximately equal numbers of up, down, and strange quarks. They are also called strangelets, nuclearites, AQNs, slets, Macros, SQNs, and MQNs. Quark nuggets are a candidate for dark matter, which has been a mystery for decades despite constituting ${\sim}85\%$ of mass in universe. Most models of quark nuggets assume no intrinsic magnetic field; however, Tatsumi found that a magnetar core may be a quark-nugget ferromagnetic liquid with a magnetic field B\textsubscript{S} between $10^{ 11}$ T and $10^{ 13}$ T. Applying that result to quark-nugget dark matter indicates magnetized quark nuggets $(MQNs)$ magnetically aggregated in the early universe before they could decay by the weak interaction and formed a broad and magnetically-stabilized mass distribution. These MQNs satisfy the requirements for dark-matter even though they are Standard Model baryons. They interact with normal matter through a magnetopause and can form non-meteorite impact craters, which are reported approximately annually. We report results from excavating such a crater. Hydrodynamic simulations indicate fractured granite below the crater and under 4.7 m of soft sediments is the first observational evidence of MQN dark matter. The results also constrain values of B\textsubscript{S} between $4x10^{ 11}$ T and $3x10^{ 12}$ T.

Main sequence stars exhibit a clear rotation-activity relationship, in which rapidly rotating stars drive strong chromospheric/coronal ultraviolet and X-ray emission. While the vast majority of red giant stars are inactive, a few percent exhibit strong ultraviolet emission. Here we use a sample of 133 red giant stars observed by SDSS APOGEE and GALEX to demonstrate an empirical relationship between NUV excess and rotational velocity (vsini). Beyond this simple relationship, we find that NUV excess also correlates with rotation period and with Rossby number in a manner that shares broadly similar trends to those found in M dwarfs, including activity saturation among rapid rotators. Our data also suggest that the most extremely rapidly rotating giants may exhibit so-called "super-saturation", which could be caused by centrifugal stripping of these stars rotating at a high fraction of breakup speed. As an example application of our empirical rotation-activity relation, we demonstrate that the NUV emission observed from a recently reported system comprising a red giant with a black hole companion is fully consistent with arising from the rapidly rotating red giant in that system. Most fundamentally, our findings suggest a common origin of chromospheric activity in rotation and convection for cool stars from main sequence to red giant stages of evolution.

Citizen science has helped astronomers comb through large data sets to identify patterns and objects that are not easily found through automated processes. The Milky Way Project (MWP), a citizen science initiative on the Zooniverse platform, presents internet users with infrared (IR) images from Spitzer Space Telescope Galactic plane surveys. MWP volunteers make classification drawings on the images to identify targeted classes of astronomical objects. We present the MWP second data release (DR2) and an updated data reduction pipeline written in Python. We aggregate ${\sim}3$ million classifications made by MWP volunteers during the years 2012-2017 to produce the DR2 catalogue, which contains 2600 IR bubbles and 599 candidate bow-shock driving stars. The reliability of bubble identifications, as assessed by comparison to visual identifications by trained experts and scoring by a machine-learning algorithm, is found to be a significant improvement over DR1. We assess the reliability of IR bow shocks via comparison to expert identifications and the colours of candidate bow-shock driving stars in the 2MASS point-source catalogue. We hence identify highly-reliable subsets of 1394 DR2 bubbles and 453 bow-shock driving stars. Uncertainties on object coordinates and bubble size/shape parameters are included in the DR2 catalog. Compared with DR1, the DR2 bubbles catalogue provides more accurate shapes and sizes. The DR2 catalogue identifies 311 new bow shock driving star candidates, including three associated with the giant HII regions NGC 3603 and RCW 49.

Arcuate infrared nebulae are ubiquitous throughout the Galactic Plane and are candidates for partial shells, bubbles, or bowshocks produced by massive runaway stars. We tabulate infrared photometry for 709 such objects using images from the Spitzer Space Telescope (SST), Wide-Field Infrared Explorer (WISE), and Herschel Space Observatory (HSO). Of the 709 objects identified at 24 or 22 microns, 422 are detected at the HSO 70 micron bandpass. Of these, only 39 are detected at HSO 160 microns. The 70 micron peak surface brightnesses are 0.5 to 2.5 Jy/square arcminute. Color temperatures calculated from the 24 micron to 70 micron ratios range from 80 K to 400 K. Color temperatures from 70 micron to 160 micron ratios are systematically lower, 40 K to 200 K. Both of these temperature are, on average, 75% higher than the nominal temperatures derived by assuming that dust is in steady-state radiative equilibrium. This may be evidence of stellar wind bowshocks sweeping up and heating --- possibly fragmenting but not destroying --- interstellar dust. Infrared luminosity correlates with standoff distance, R_0, as predicted by published hydrodynamical models. Infrared spectral energy distributions are consistent with interstellar dust exposed to a either single radiant energy density, U=10^3 to 10^5 (in more than half of the objects) or a range of radiant energy densities U_min=25 to U_max=10^3 to 10^5 times the mean interstellar value for the remainder. Hence, the central OB stars dominate the energetics, making these enticing laboratories for testing dust models in constrained radiation environments. SEDs are consistent with PAH fractions q_PAH <1% in most objects.

We identify 709 arc-shaped mid-infrared nebula in 24 micron Spitzer Space Telescope or 22 micron Wide Field Infrared Explorer surveys of the Galactic Plane as probable dusty interstellar bowshocks powered by early-type stars. About 20% are visible at 8 microns or shorter mid-infrared wavelengths as well. The vast majority (660) have no previous identification in the literature. These extended infrared sources are strongly concentrated near Galactic mid-Plane with an angular scale height of ~0.6 degrees. All host a symmetrically placed star implicated as the source of a stellar wind sweeping up interstellar material. These are candidate "runaway" stars potentially having high velocities in the reference frame of the local medium. Among the 286 objects with measured proper motions, we find an unambiguous excess having velocity vectors aligned with the infrared morphology --- kinematic evidence that many of these are "runaway" stars with large peculiar motions responsible for the bowshock signature. We discuss a population of "in-situ" bowshocks (103 objects) that face giant HII regions where the relative motions between the star and ISM may be caused by bulk outflows from an overpressured bubble. We also identify 58 objects that face 8 micron bright-rimmed clouds and apparently constitute a sub-class of in-situ bowshocks where the stellar wind interacts with a photo-evaporative flow from an eroding molecular cloud interface (i.e., "PEF bowshocks"). Orientations of the arcuate nebulae exhibit a correlation over small angular scales, indicating that external influences such as HII regions are responsible for producing some bowshock nebulae. However, the vast majority of this sample appear to be isolated (499 objects) from obvious external influences.

This paper deals with the problem of astrometric determination of the orbital elements of the outer planets, in particular by assessing the ability of astrometric observations to detect perturbations of the sort expected from the Pioneer effect or other small perturbations to gravity. We also show that while using simplified models of the dynamics can lead to some insights, one must be careful to not over-simplify the issues involved lest one be misled by the analysis onto false paths. Specifically, we show that the current ephemeris of Pluto does not preclude the existence of the Pioneer effect. We show that the orbit of Pluto is simply not well enough characterized at present to make such an assertion. A number of misunderstandings related to these topics have now propagated through the literature and have been used as a basis for drawing conclusions about the dynamics of the solar system. Thus, the objective of this paper is to address these issues. Finally, we offer some comments dealing with the complex topic of model selection and comparison.

The inverse square law of gravity is poorly probed by experimental tests at distances of ~ 10 AUs. Recent analysis of the trajectory of the Pioneer 10 and 11 spacecraft have shown an unmodeled acceleration directed toward the Sun which was not explained by any obvious spacecraft systematics, and occurred when at distances greater than 20 AUs from the Sun. If this acceleration represents a departure from Newtonian gravity or is indicative of an additional mass distribution in the outer solar system, it should be detectable in the orbits of Trans-Neptunian Objects (TNOs). To place limits on deviations from Newtonian gravity, we have selected a well observed sample of TNOs found orbiting between 20 and 100 AU from the Sun. By examining their orbits with modified orbital fitting software, we place tight limits on the perturbations of gravity that could exist in this region of the solar system.

The twin Pioneer spacecraft have been tracked for over thirty years as they headed out of the solar system. After passing 20 AU from the Sun, both exhibited a systematic error in their trajectories that can be interpreted as a constant acceleration towards the Sun. This Pioneer Effect is most likely explained by spacecraft systematics, but there have been no convincing arguments that that is the case. The alternative is that the Pioneer Effect represents a real phenomenon and perhaps new physics. What is lacking is a means of measuring the effect, its variation, its potential anisotropies, and its region of influence. We show that minor planets provide an observational vehicle for investigating the gravitational field in the outer solar system, and that a sustained observation campaign against properly chosen minor planets could confirm or refute the existence of the Pioneer Effect. Additionally, even if the Pioneer Effect does not represent a new physical phenomenon, minor planets can be used to probe the gravitational field in the outer Solar System and since there are very few intermediate range tests of gravity at the multiple AU distance scale, this is a worthwhile endeavor in its own right.

By combining OSSE, SMM and TGRS observations of the galactic center region, Purcell et al. (1997) and Cheng et al. (1997) produced the first maps of galactic positron annihilation. That data-set has been augmented with additional data, both recent and archival, and re-analyzed to improve the spectral fitting. The improved spectral fitting has enabled the first maps of positronium continuum emission and the most extensive maps of 511 keV line emission. Bulge and disk combinations have been compared with the 511 keV line data, demonstrating that extended bulges are favored over a GC point source for every disk model tested. This result is independent of whether OSSE-only, OSSE/SMM, or OSSE/SMM/TGRS data-sets are used. The estimated bulge to disk ratio (and to a lesser extent the total flux) is shown to be dependent upon the assumption of bulge shape. A positive latitude enhancement is shown to have an effect upon the B/D ratio, but this effect is secondary to the choice of bulge shape.

We present a new algorithm, called Multiresolution Regularized Expectation Maximization (MREM), for the reconstruction of gamma-ray intensity maps from COMPTEL data. The algorithm is based on the iterative Richardson-Lucy scheme to which we added a wavelet thresholding step in order to eliminate image-noise in the reconstruction. The wavelet thresholding explicitly accounts for spatial correlations in the data, and adapts the angular resolution locally, depending on the significance of the signal in the data. We compare the performance of MREM to that of the maximum entropy and the Richardson-Lucy algorithms by means of Monte-Carlo simulations of COMPTEL 1.809 MeV gamma-ray line observations. The simulations demonstrate that the maximum entropy and Richardson-Lucy algorithms provide virtually identical reconstructions which are heavily disturbed by image noise. MREM largely suppresses this noise in the reconstructions, showing only the significant structures that are present in the data. Application of MREM to COMPTEL 1.8 MeV gamma-ray line data results in a 1.809 MeV sky map that is much smoother than the maximum entropy or Richardson-Lucy reconstructions presented previously. The essential features of this map are (1) an asymmetric galactic ridge emission reaching from l=45 deg to l=240 deg, (2) a bright localised emission feature in the Cygnus region around (l,b)=(80 deg,0 deg), (3) two emission spots at l=317 deg and l=332 deg situated in the galactic plane, and (4) an extended emission region around (l,b)=(160 deg, 0 deg). Comparison of the MREM map to the simulated reconstructions demonstrates that the 1.809 MeV emission is confined to the galactic plane.

Deconvolving COMPTEL gamma-ray data into images presents a major methodological challenge. We developed a new algorithm called Multiresolution Regularized Expectation Maximization (MREM), which explicitly accounts for spatial correlations in the image by using wavelets (Dixon et al. 1998). We demonstrate that MREM largely suppresses image noise in the reconstruction, showing only significant structures that are present in the data. Application to 1.8 MeV data results in a sky map that is much smoother than the maximum entropy reconstructions presented previously, but it shows the same characteristic emission features which have been established earlier.

Image processing is an increasingly important aspect for analysis of data from X and $\gamma$-ray astrophysics missions. In this paper, I review a method proposed by Kebede (L. W. Kebede 1994, ApJ, 423, 878), and point out an error in the derivation of this method. It turns out that this error is not debilitating -- the method still ``works'' in some sense -- but as published is rather sub-optimal, as demonstrated both on theoretical grounds and via a set of examples.

We present quantitative statistical evidence for a $\gamma$-ray emission halo surrounding the Galaxy. Maps of the emission are derived. EGRET data were analyzed in a wavelet-based non-parametric hypothesis testing framework, using a model of expected diffuse (Galactic + isotropic) emission as a null hypothesis. The results show a statistically significant large scale halo surrounding the center of the Milky Way as seen from Earth. The halo flux at high latitudes is somewhat smaller than the isotropic gamma-ray flux at the same energy, though of the same order (O(10^(-7)--10^(-6)) ph/cm^2/s/sr above 1 GeV).

The region near the Galactic center may have experienced recurrent episodes of injection of energy in excess of $\sim$ 10$^{55}$ ergs due to repeated starbursts involving more than $\sim$ 10$^4$ supernovae. This hypothesis can be tested by measurements of $\gamma$-ray lines produced by the decay of radioactive isotopes and positron annihilation, or by searches for pulsars produced during starbursts. Recent OSSE observations of 511 keV emission extending above the Galactic center led to the suggestion of a starburst driven fountain from the Galactic center. We present EGRET observations that might support this picture.

We present progress toward using scanned OSSE observations for mapping and sky survey work. To this end, we have developed a technique for detecting pointlike sources of unknown number and location, given that they appear in a background which is relatively featureless or which can be modeled. The technique, based on the newly developed concept and mean field annealing, is described, with sample reconstructions of data from the OSSE Virgo Survey. The results demonstrate the capability of reconstructing source information without any a priori information about the number and/or location of pointlike sources in the field-of-view.