We present a catalog of 316 trans-Neptunian bodies detected by the Dark Energy Survey (DES). These objects include 245 discoveries by DES (139 not previously published) detected in $\approx 60,000$ exposures from the first four seasons of the survey ("Y4" data). The survey covers a contiguous 5000 deg$^2$ of the southern sky in the $grizY$ optical/NIR filter set, with a typical TNO in this part of the sky being targeted by $25-30$ Y4 exposures. We describe the processes for detection of transient sources and the linkage into TNO orbits, which are made challenging by the absence of the few-hour repeat observations employed by TNO-optimized surveys. We also describe the procedures for determining detection efficiencies vs. magnitude and estimating rates of false-positive linkages. This work presents all TNOs which were detected on $\ge 6$ unique nights in the Y4 data and pass a "sub-threshold confirmation" test wherein we demand the the object be detectable in a stack of the individual images in which the orbit indicates an object should be present, but was not detected. This eliminates false positives and yields TNO detections complete to $r\lesssim 23.3$ mag with virtually no dependence on orbital properties for bound TNOs at distance $30\,{\rm AU}<d<2500\,{\rm AU}.$ The final DES TNO catalog is expected to yield $>0.3$ mag more depth, and arcs of $>4$ years for nearly all detections.

The outer Solar System contains a large number of small bodies (known as trans-Neptunian objects or TNOs) that exhibit diverse types of dynamical behavior. The classification of bodies in this distant region into dynamical classes -- sub-populations that experience similar orbital evolution -- aids in our understanding of the structure and formation of the Solar System. In this work, we propose an updated dynamical classification scheme for the outer Solar System. This approach includes the construction of a new (automated) method for identifying mean-motion resonances. We apply this algorithm to the current dataset of TNOs observed by the Dark Energy Survey (DES) and present a working classification for all of the DES TNOs detected to date. Our classification scheme yields 1 inner centaur, 19 outer centaurs, 21 scattering disk objects, 47 detached TNOs, 48 securely resonant objects, 7 resonant candidates, and 97 classical belt objects. Among the scattering and detached objects, we detect 8 TNOs with semi-major axes greater than 150 AU.

In the current model of early Solar System evolution, the stable members of the Jovian and Neptunian Trojan populations were captured into resonance from the leftover reservoir of planetesimals during the outward migration of the giant planets. As a result, both Jovian and Neptunian Trojans share a common origin with the primordial disk population, whose other surviving members constitute today's trans-Neptunian object (TNO) populations. The cold classical TNOs are ultra-red, while the dynamically excited "hot" population of TNOs contains a mixture of ultra-red and blue objects. In contrast, Jovian and Neptunian Trojans are observed to be blue. While the absence of ultra-red Jovian Trojans can be readily explained by the sublimation of volatile material from their surfaces due to the high flux of solar radiation at 5AU, the lack of ultra-red Neptunian Trojans presents both a puzzle and a challenge to formation models. In this work we report the discovery by the Dark Energy Survey (DES) of two new dynamically stable L4 Neptunian Trojans,2013 VX30 and 2014 UU240, both with inclinations i >30 degrees, making them the highest-inclination known stable Neptunian Trojans. We have measured the colors of these and three other dynamically stable Neptunian Trojans previously observed by DES, and find that 2013 VX30 is ultra-red, the first such Neptunian Trojan in its class. As such, 2013 VX30 may be a "missing link" between the Trojan and TNO populations. Using a simulation of the DES TNO detection efficiency, we find that there are 162 +/- 73 Trojans with Hr < 10 at the L4 Lagrange point of Neptune. Moreover, the blue-to-red Neptunian Trojan population ratio should be higher than 17:1. Based on this result, we discuss the possible origin of the ultra-red Neptunian Trojan population and its implications for the formation history of Neptunian Trojans.

We report the discovery of a compact multi-planet system orbiting the relatively nearby (78pc) and bright ($K=8.9$) K-star, K2-266 (EPIC248435473). We identify up to six possible planets orbiting K2-266 with estimated periods of P$_b$ = 0.66, P$_{.02}$ = 6.1, P$_c$ = 7.8, P$_d$ = 14.7, P$_e$ = 19.5, and P$_{.06}$ = 56.7 days and radii of R$_P$ = 3.3 R$_{\oplus}$, 0.646 R$_{\oplus}$, 0.705 R$_{\oplus}$, 2.93 R$_{\oplus}$, 2.73 R$_{\oplus}$, and 0.90 R$_{\oplus}$, respectively. We are able to confirm the planetary nature of two of these planets (d & e) from analyzing their transit timing variations ($m_d= 8.9_{-3.8}^{+5.7} M_\oplus$ and $m_e=14.3_{-5.0}^{+6.4} M_\oplus$), confidently validate the planetary nature of two other planets (b & c), and classify the last two as planetary candidates (K2-266.02 & .06). From a simultaneous fit of all 6 possible planets, we find that K2-266 b's orbit has an inclination of 75.32$^{\circ}$ while the other five planets have inclinations of 87-90$^{\circ}$. This observed mutual misalignment may indicate that K2-266 b formed differently from the other planets in the system. The brightness of the host star and the relatively large size of the sub-Neptune sized planets d and e make them well-suited for atmospheric characterization efforts with facilities like the Hubble Space Telescope and upcoming James Webb Space Telescope. We also identify an 8.5-day transiting planet candidate orbiting EPIC248435395, a co-moving companion to K2-266.

This paper reports the discovery and orbital characterization of two extreme trans-Neptunian objects (ETNOs), 2016 QV$_{89}$ and 2016 QU$_{89}$, which have orbits that appear similar to that of a previously known object, 2013 UH$_{15}$. All three ETNOs have semi-major axes $a\approx 172$ AU and eccentricities $e\approx0.77$. The angular elements $(i,\omega,\Omega)$ vary by 6, 15, and 49 deg, respectively between the three objects. The two new objects add to the small number of TNOs currently known to have semi-major axes between 150 and 250 AU, and serve as an interesting dynamical laboratory to study the outer realm of our Solar System. Using a large ensemble of numerical integrations, we find that the orbits are expected to reside in close proximity in the $(a,e)$ phase plane for roughly 100 Myr before diffusing to more separated values. We then explore other scenarios that could influence their orbits. With aphelion distances over 300 AU, the orbits of these ETNOs extend far beyond the classical Kuiper Belt, and an order of magnitude beyond Neptune. As a result, their orbital dynamics can be affected by the proposed new Solar System member, referred to as Planet Nine in this work. With perihelion distances of 35-40 AU, these orbits are also influenced by resonant interactions with Neptune. A full assessment of any possible, new Solar System planets must thus take into account this emerging class of TNOs.

We report the discovery and dynamical analysis of 2015 BP$_{519}$, an extreme Trans-Neptunian Object detected detected by the Dark Energy Survey at a heliocentric distance of 55 AU and absolute magnitude Hr= 4.3. The current orbit, determined from a 1110-day observational arc, has semi-major axis $a\approx$ 450 AU, eccentricity $e\approx$ 0.92 and inclination $i\approx$ 54 degrees. With these orbital elements, 2015 BP$_{519}$ is the most extreme TNO discovered to date, as quantified by the reduced Kozai action, which is is a conserved quantity at fixed semi-major axis $a$ for axisymmetric perturbations. We discuss the orbital stability and evolution of this object in the context of the known Solar System, and find that 2015 BP$_{519}$ displays rich dynamical behavior, including rapid diffusion in semi-major axis and more constrained variations in eccentricity and inclination. We also consider the long term orbital stability and evolutionary behavior within the context of the Planet Nine Hypothesis, and find that BP$_{519}$ adds to the circumstantial evidence for the existence of this proposed new member of the Solar System, as it would represent the first member of the population of high-i, $\varpi$-shepherded TNOs.

The observation that the orbits of long-period Kuiper Belt objects are anomalously clustered in physical space has recently prompted the Planet Nine hypothesis - the proposed existence of a distant and eccentric planetary member of our solar system. Within the framework of this model, a Neptune-like perturber sculpts the orbital distribution of distant Kuiper Belt objects through a complex interplay of resonant and secular effects, such that in addition to perihelion-circulating objects, the surviving orbits get organized into apsidally aligned and anti-aligned configurations with respect to Planet Nine's orbit. In this work, we investigate the role of Kuiper Belt initial conditions on the evolution of the outer solar system using numerical simulations. Intriguingly, we find that the final perihelion distance distribution depends strongly on the primordial state of the system, and demonstrate that a bimodal structure corresponding to the existence of both aligned and anti-aligned clusters is only reproduced if the initial perihelion distribution is assumed to extend well beyond $\sim 36$ AU. The bimodality in the final perihelion distance distribution is due to the existence of permanently stable objects, with the lower perihelion peak corresponding to the anti-aligned orbits and the higher perihelion peak corresponding to the aligned orbits. We identify the mechanisms which enable the persistent stability of these objects and locate the regions of phase space in which they reside. The obtained results contextualize the Planet Nine hypothesis within the broader narrative of solar system formation, and offer further insight into the observational search for Planet Nine.

The existence of hot Jupiters has challenged theories of planetary formation since the first extrasolar planets were detected. Giant planets are generally believed to form far from their host stars, where volatile materials like water exist in their solid phase, making it easier for giant planet cores to accumulate. Several mechanisms have been proposed to explain how giant planets can migrate inward from their birth sites to short-period orbits. One such mechanism, called Kozai-Lidov migration, requires the presence of distant companions in orbits inclined by more than $\sim40$ degrees with respect to the plane of the hot Jupiter's orbit. The high occurrence rate of wide companions in hot Jupiter systems lends support to this theory for migration. However, the exact orbital inclinations of these detected planetary and stellar companions is not known, so it is not clear whether the mutual inclination of these companions is large enough for the Kozai-Lidov process to operate. This paper shows that in systems orbiting cool stars with convective outer layers, the orbits of most wide planetary companions to hot Jupiters must be well aligned with the orbits of the hot Jupiters and the spins of the host stars. For a variety of possible distributions for the inclination of the companion, the width of the distribution must be less than $\sim20$ degrees to recreate the observations with good fidelity. As a result, the companion orbits are likely well-aligned with those of the hot Jupiters, and the Kozai-Lidov mechanism does not enforce migration in these systems.

We evaluate the dynamical stability of a selection of outer solar system objects in the presence of the proposed new Solar System member Planet Nine. We use a Monte Carlo suite of numerical N-body integrations to construct a variety of orbital elements of the new planet and evaluate the dynamical stability of eight Trans-Neptunian objects (TNOs) in the presence of Planet Nine. These simulations show that some combinations of orbital elements ($a,e$) result in Planet Nine acting as a stabilizing influence on the TNOs, which can otherwise be destabilized by interactions with Neptune. These simulations also suggest that some TNOs transition between several different mean-motion resonances during their lifetimes while still retaining approximate apsidal anti-alignment with Planet Nine. This behavior suggests that remaining in one particular orbit is not a requirement for orbital stability. As one product of our simulations, we present an {\it a posteriori} probability distribution for the semi-major axis and eccentricity of the proposed Planet Nine based on TNO stability. This result thus provides additional evidence that supports the existence of this proposed planet. We also predict that TNOs can be grouped into multiple populations of objects that interact with Planet Nine in different ways: one population may contain objects like Sedna and 2012 VP$_{113}$, which do not migrate significantly in semi-major axis in the presence of Planet Nine and tend to stay in the same resonance; another population may contain objects like 2007 TG$_{422}$ and 2013 RF$_{98}$, which may both migrate and transition between different resonances.

We report the observation and physical characterization of the possible dwarf planet \UZ\ ("DeeDee"), a dynamically detached trans-Neptunian object discovered at 92 AU. This object is currently the second-most distant known trans-Neptunian object with reported orbital elements, surpassed in distance only by the dwarf planet Eris. The object was discovered with an $r$-band magnitude of 23.0 in data collected by the Dark Energy Survey between 2014 and 2016. Its 1140-year orbit has $(a,e,i) = (109~\mathrm{AU}, 0.65, 26.8^{\circ})$. It will reach its perihelion distance of 38 AU in the year 2142. Integrations of its orbit show it to be dynamically stable on Gyr timescales, with only weak interactions with Neptune. We have performed followup observations with ALMA, using 3 hours of on-source integration time to measure the object's thermal emission in the Rayleigh-Jeans tail. The signal is detected at 7$\sigma$ significance, from which we determine a $V$-band albedo of $13.1^{+3.3}_{-2.4}\mathrm{(stat)}^{+2.0}_{-1.4}\mathrm{(sys)}$ percent and a diameter of $635^{+57}_{-61}\mathrm{(stat)}^{+32}_{-39}\mathrm{(sys)}$~km, assuming a spherical body with uniform surface properties.