We develop an improved sky background estimator which employs optimal filters for both spatial and pixel intensity distributions. It incorporates growth of masks around detected objects and a statistical estimate of the flux from undetected faint galaxies in the remaining sky pixels. We test this algorithm for underlying sky estimation and compare its performance with commonly used sky estimation codes on realistic simulations which include detected galaxies, faint undetected galaxies, and sky noise. We then test galaxy surface brightness recovery using GALFIT 3, a galaxy surface brightness profile fitting optimizer, yielding fits to S\'{e}rsic profiles. This enables robust sky background estimates accurate at the 4 parts-per-million level. This background sky estimator is more accurate and is less affected by surface brightness profiles of galaxies and the local image environment compared with other methods.

We evaluate the optical/near-infrared (OIR) color variability of 3C 279 in both gamma-ray flaring and non-flaring states over 7-year timescales using the Small and Medium Aperture Research Telescope System (SMARTS) in Cerro Tololo, Chile and gamma-ray fluxes obtained from the Fermi Gamma-ray Space Telescope. This observing strategy differs from previous blazar color variability studies in two key ways: 1) the reported color variability is assessed across optical through near-infrared wavelengths, and 2) the color variability is assessed over timescales significantly longer than an individual flare or ground-based observing season. We highlight 3C 279 because of its complex color variability, which is difficult to reconcile with the simple "redder when brighter" behavior often associated with Flat Spectrum Radio Quasar (FSRQ) color variability. We suggest that the observed OIR color changes depend on a combination of the jet and disk emission. We parameterize this behavior in terms of a single variable, $\zeta^m_n$, representing a smooth transition from disk-dominated, to a mixed contribution, to a jet-dominated system, which provides an explanation of the long-term OIR color variability in the same blazar over time. This suggests a general scheme that could apply to OIR color variability in other blazars.

The optical counterpart of the black-hole soft X-ray transient Nova Muscae 1991 has brightened by $\Delta{V}\approx0.8$ mag since its return to quiescence 23 years ago. We present the first clear evidence that the brightening of soft X-ray transients in quiescence occurs at a nearly linear rate. This discovery, and our precise determination of the disk component of emission obtained using our $simultaneous$ photometric and spectroscopic data, have allowed us to identify and accurately model archival ellipsoidal light curves of the highest quality. The simultaneity, and the strong constraint it provides on the component of disk emission, is a key element of our work. Based on our analysis of the light curves, and our earlier measurements of the mass function and mass ratio, we have obtained for Nova Muscae 1991 the first accurate estimates of its systemic inclination $i=43.2^{+2.1}_{-2.7}$ deg, and black hole mass $M=11.0^{+2.1}_{-1.4}\ M_\odot$. Based on our determination of the radius of the secondary, we estimate the distance to be $D=4.95^{+0.69}_{-0.65}$ kpc. We discuss the implications of our work for future dynamical studies of black-hole soft X-ray transients.

We report the results of a multi-band observing campaign on the famous blazar 3C 279 conducted during a phase of increased activity from 2013 December to 2014 April, including first observations of it with NuSTAR. The $\gamma$-ray emission of the source measured by Fermi-LAT showed multiple distinct flares reaching the highest flux level measured in this object since the beginning of the Fermi mission, with $F(E > 100\,{\rm MeV})$ of $10^{-5}$ photons cm$^{-2}$ s$^{-1}$, and with a flux doubling time scale as short as 2 hours. The $\gamma$-ray spectrum during one of the flares was very hard, with an index of $\Gamma_\gamma = 1.7 \pm 0.1$, which is rarely seen in flat spectrum radio quasars. The lack of concurrent optical variability implies a very high Compton dominance parameter $L_\gamma/L_{\rm syn} > 300$. Two 1-day NuSTAR observations with accompanying Swift pointings were separated by 2 weeks, probing different levels of source activity. While the 0.5$-$70 keV X-ray spectrum obtained during the first pointing, and fitted jointly with Swift-XRT is well-described by a simple power law, the second joint observation showed an unusual spectral structure: the spectrum softens by $\Delta\Gamma_{\rm X} \simeq 0.4$ at $\sim$4 keV. Modeling the broad-band SED during this flare with the standard synchrotron plus inverse Compton model requires: (1) the location of the $\gamma$-ray emitting region is comparable with the broad line region radius, (2) a very hard electron energy distribution index $p \simeq 1$, (3) total jet power significantly exceeding the accretion disk luminosity $L_{\rm j}/L_{\rm d} \gtrsim 10$, and (4) extremely low jet magnetization with $L_{\rm B}/L_{\rm j} \lesssim 10^{-4}$. We also find that single-zone models that match the observed $\gamma$-ray and optical spectra cannot satisfactorily explain the production of X-ray emission.

In anticipation of the July 2015 flyby of the Pluto system by NASA's New Horizons mission, we propose naming conventions and example names for surface features on Pluto and its satellites (Charon, Nix, Hydra, Kerberos, Styx) and names for newly discovered satellites.

We present multi-epoch optical spectroscopy of seven southern Fermi-monitored blazars from 2008 - 2013 using the Small and Medium Aperture Research Telescope System (SMARTS), with supplemental spectroscopy and polarization data from the Steward Observatory. We find that the emission lines are much less variable than the continuum; 4 of 7 blazars had no detectable emission line variability over the 5 years. This is consistent with photoionization primarily by an accretion disk, allowing us to use the lines as a probe of disk activity. Comparing optical emission line flux with Fermi $\gamma$-ray flux and optical polarized flux, we investigate whether relativistic jet variability is related to the accretion flow. In general, we see no such dependence, suggesting the jet variability is likely caused by internal processes like turbulence or shock acceleration rather than a variable accretion rate. However, three sources showed statistically significant emission line flares in close temporal proximity to very large Fermi $\gamma$-ray flares. While we do not have sufficient emission line data to quantitatively assess their correlation with the $\gamma$-ray flux, it appears that in some cases, the jet might provide additional photoionizing flux to the broad line region, which implies some gamma-rays are produced within the broad line region, at least for these large flares.

A fraction of multiple planet candidate systems discovered from transits by the Kepler mission contain pairs of planet candidates that are in orbital resonance or are spaced slightly too far apart to be in resonance. We focus here on the four planet system, KOI 730, that has planet periods satisfying the ratios 8:6:4:3. By numerically integrating four planets initially in this resonant configuration in proximity to an initially exterior cold planetesimal disk, we find that of the order of a Mars mass of planet-orbit-crossing planetesimals is sufficient to pull this system out of resonance. Approximately one Earth mass of planet-orbit-crossing planetesimals increases the interplanetary spacings sufficiently to resemble the multiple planet candidate Kepler systems that lie just outside of resonance. This suggests that the closely spaced multiple planet Kepler systems, host only low mass debris disks or their debris disks have been extremely stable. We find that the planetary inclinations increase as a function of the mass in planetesimals that have crossed the orbits of the planets. If systems are left at zero inclination and in resonant chains after depletion of the gas disk then we would expect a correlation between distance to resonance and mutual planetary inclinations. This may make it possible to differentiate between dynamical mechanisms that account for the fraction of multiple planet systems just outside of resonance.

By logging encounters between planetesimals and planets we compute the distribution of encounters in a numerically integrated two planet system that is migrating due to interactions with an exterior planetesimal belt. Capture of an irregular satellite in orbit about a planet through an exchange reaction with a binary planetesimal is only likely when the binary planetesimal undergoes a slow and close encounter with the planet. In our simulations we find that close and slow encounters between planetesimals and a planet primarily occur with the outermost and not innermost planet. Taking care to consider where a planet orbit crossing binary planetesimal would first be tidally disrupted, we estimate the probability of both tidal disruption and irregular satellite capture. We estimate that the probability that the secondary of a binary planetesimal is captured and becomes an irregular satellite about a Neptune mass outer planet is about 1/100 for binaries with masses and separations similar to transneptunian planetesimal binaries. If young exoplanetary debris disks host a binary planetesimal population then outwards migrating outer planets should host captured irregular satellite populations. We discuss interpretation of emission associated with the exoplanet Fomalhaut b in terms of collisional evolution of a captured irregular satellite population that is replenished due to planetary migration.