We place the first constraints on the obliquity of a planetary-mass companion (PMC) outside of the Solar System. Our target is the directly imaged system 2MASS J01225093-2439505 (2M0122), which consists of a 120 Myr 0.4 M_sun star hosting a 12-27 M_J companion at 50 AU. We constrain all three of the system's angular momentum vectors: how the companion spin axis, the stellar spin axis, and the orbit normal are inclined relative to our line of sight. To accomplish this, we measure projected rotation rates (vsini) for both the star and the companion using new near-infrared high-resolution spectra with NIRSPEC at Keck Observatory. We combine these with a new stellar photometric rotation period from TESS and a published companion rotation period from HST to obtain spin axis inclinations for both objects. We also fitted multiple epochs of astrometry, including a new observation with NIRC2/Keck, to measure 2M0122b's orbital inclination. The three line-of-sight inclinations place limits on the true de-projected companion obliquity and stellar obliquity. We find that while the stellar obliquity marginally prefers alignment, the companion obliquity tentatively favors misalignment. We evaluate possible origin scenarios. While collisions, secular spin-orbit resonances, and Kozai-Lidov oscillations are unlikely, formation by gravitational instability in a gravito-turbulent disk - the scenario favored for brown dwarf companions to stars - appears promising.

While the Solar System contains no planets between the sizes of Uranus and Saturn, our current exoplanet census includes several dozen such planets with well-measured masses and radii. These sub-Saturns exhibit a diversity of bulk densities, ranging from ~$0.1-3\ \rm{g\ cm}^{-3}$. When modeled simply as hydrogen/helium envelopes atop rocky cores, this diversity in densities translates to a diversity in planetary envelope fractions, $f_\rm{env}=M_\rm{env}/M_p$ ranging from ~$10\%$ to ~$50\%$. Planets with $f_\rm{env}\sim50\%$ pose a challenge to traditional models of giant planet formation by core-nucleated accretion, which predict the onset of runaway gas accretion when $M_\rm{env}\sim M_\rm{core}$. Here we show that many of these apparent $f_\rm{env}\sim50\%$ planets are less envelope rich than they seem, after accounting for tidal heating. We present a new framework for modeling sub-Saturn interiors that incorporates envelope inflation due to tides, which are driven by the observed non-zero eccentricities, as well as potential obliquities. Consequently, when we apply our models to known sub-Saturns, we infer lower $f_\rm{env}$ than tides-free estimates. We present a case study of K2-19 b, a moderately eccentric sub-Saturn. Neglecting tides, K2-19 b appears to have $f_\rm{env}\sim50\%$, poised precariously near the runaway threshold; by including tides, we find $f_\rm{env}\sim10\%$, resolving the tension. Through a systematic analysis of $4-8\ R_{\oplus}$ planets, we find that most (but not all) of the similarly envelope-rich planets have more modest envelopes of $f_\rm{env}\sim10\%-20\%$. Thus, many sub-Saturns may be understood as sub-Neptunes that have undergone significant radius inflation, rather than a separate class of objects. Tidal radius inflation likely plays an important role in other size classes of planets including ultra-low-density Jupiter-size planets like WASP-107 b.

We report the discovery and confirmation of two new hot Jupiters discovered by the Transiting Exoplanet Survey Satellite (TESS): TOI 564 b and TOI 905 b. The transits of these two planets were initially observed by TESS with orbital periods of 1.651 d and 3.739 d, respectively. We conducted follow-up observations of each system from the ground, including photometry in multiple filters, speckle interferometry, and radial velocity measurements. For TOI 564 b, our global fitting revealed a classical hot Jupiter with a mass of $1.463^{+0.10}_{-0.096}\ M_J$ and a radius of $1.02^{+0.71}_{-0.29}\ R_J$. TOI 905 b is a classical hot Jupiter as well, with a mass of $0.667^{+0.042}_{-0.041}\ M_J$ and radius of $1.171^{+0.053}_{-0.051}\ R_J$. Both planets orbit Sun-like, moderately bright, mid-G dwarf stars with V ~ 11. While TOI 905 b fully transits its star, we found that TOI 564 b has a very high transit impact parameter of $0.994^{+0.083}_{-0.049}$, making it one of only ~20 known systems to exhibit a grazing transit and one of the brightest host stars among them. TOI 564 b is therefore one of the most attractive systems to search for additional non-transiting, smaller planets by exploiting the sensitivity of grazing transits to small changes in inclination and transit duration over the time scale of several years.

Recent work suggests that many short-period super-Earth and sub-Neptune planets may have significant spin axis tilts ("obliquities"). When planets are locked in high-obliquity states, the tidal dissipation rate may increase by several orders of magnitude. This intensified heat deposition within the planets' interiors should generate significant structural consequences, including atmospheric inflation leading to larger transit radii. Using up-to-date radius estimates from Gaia Data Release 2, we show evidence for $\sim50\%$ larger average radii of planets wide of first-order mean-motion resonances, a population of planets with a theorized frequent occurrence of high obliquities. We investigate whether this radius trend could be a signature of obliquity tides. Using an adaptation of the Modules for Experiments in Stellar Astrophysics (MESA) stellar evolution toolkit, we model the atmospheric evolution of sub-Neptune-mass planets in response to additional internal heat from obliquity tides. The degree of radius inflation predicted by the models is $\sim10\%-100\%$ for tidal luminosities $\gtrsim 10^{-5}$ of the incident stellar power; this degree of inflation is broadly consistent with the observations and can approximately be described by power law relationships. We present a few case studies of very low density "super-puff" planets -- Kepler-79 d, Kepler-31 c, and Kepler-27 b -- and show that they are strong candidates for potentially having undergone tidally-induced radius inflation. We also discuss how the discrepancy between the two populations of planets with masses derived from radial velocities and transit timing variations is connected to the radius distribution features we have identified. Altogether, the calculations in this work confirm that tidal dissipation has non-negligible consequences for the structural properties of short-period sub-Neptunes.

Recent work suggests that many short-period extrasolar planets may have spin obliquities that are significantly tilted with respect to their orbital planes. These large obliquities are a natural outcome of ``secular spin-orbit resonance'', a configuration in which the planetary spin precession frequency matches the frequency of orbit nodal regression, or a Fourier component thereof. While exoplanet spin obliquities have not yet been measured directly, they may be detectable indirectly through their signatures in various observations, such as photometric measurements across the full phase of a planet's orbit. In this work, we employ a thermal radiative model to explore how large polar tilts affect full-phase light curves, and we discuss the range of unique signatures that are expected to result. We show that the well-studied short-period planets HD 149026 b, WASP-12 b, and CoRoT-2 b all exhibit phase curve features that may arise from being in high-obliquity states. We also constrain the parameters and assess the detectability of hypothetical perturbing planets that could maintain the planets in these states. Among the three planets considered, CoRoT-2 b has the tightest constraints on its proposed obliquity ($45.8^{\circ} \pm 1.4^{\circ}$) and axial orientation. For HD 149026 b, we find no significant evidence for a non-zero obliquity, and the phase curve of WASP-12 b is too complicated by strong tidal distortions for a conclusive assessment.

The tilt of a planet's spin axis off its orbital axis ("obliquity") is a basic physical characteristic that plays a central role in determining the planet's global circulation and energy redistribution. Moreover, recent studies have also highlighted the importance of obliquities in sculpting not only the physical features of exoplanets but also their orbital architectures. It is therefore of key importance to identify and characterize the dominant processes of excitation of non-zero axial tilts. Here we highlight a simple mechanism that operates early on and is likely fundamental for many extrasolar planets and perhaps even Solar System planets. While planets are still forming in the protoplanetary disk, the gravitational potential of the disk induces nodal recession of the orbits. The frequency of this recession decreases as the disk dissipates, and when it crosses the frequency of a planet's spin axis precession, large planetary obliquities may be excited through capture into a secular spin-orbit resonance. We study the conditions for encountering this resonance and calculate the resulting obliquity excitation over a wide range of parameter space. Planets with semi-major axes in the range $0.3 \ \mathrm{AU} \lesssim a \lesssim 2 \ \mathrm{AU}$ are the most readily affected, but large-$a$ planets can also be impacted. We present a case study of Uranus and Neptune and show that this mechanism likely cannot help explain their high obliquities. While it could have played a role if finely tuned and envisioned to operate in isolation, large-scale obliquity excitation was likely inhibited by gravitational planet-planet perturbations.

NASA's Kepler mission revealed that $\sim 30\%$ of Solar-type stars harbor planets with sizes between that of Earth and Neptune on nearly circular and co-planar orbits with periods less than 100 days. Such short-period compact systems are rarely found with planet pairs in mean-motion resonances (MMRs) -- configurations in which the planetary orbital periods exhibit a simple integer ratio -- but there is a significant overabundance of planet pairs lying just wide of the first-order resonances. Previous work suggests that tides raised on the planets by the host star may be responsible for forcing systems into these configurations by draining orbital energy to heat. Such tides, however, are insufficient unless there exists a substantial and as-yet unidentified source of extra dissipation. Here we show that this cryptic heat source may be linked to "obliquity tides" generated when a large axial tilt (obliquity) is maintained by secular resonance-driven spin-orbit coupling. We present evidence that typical compact, nearly-coplanar systems frequently experience this mechanism, and we highlight additional features in the planetary orbital period and radius distributions that may be its signatures. Extrasolar planets that maintain large obliquities will exhibit infrared light curve features that are detectable with forthcoming space missions. The observed period ratio distribution can be explained if typical tidal quality factors for super-Earths and sub-Neptunes are similar to those of Uranus and Neptune.

Recent analyses have revealed a mystery. The orbital period of the highly inflated hot Jupiter, WASP-12b, is decreasing rapidly. The rate of inspiral, however, is too fast to be explained by either eccentricity tides or equilibrium stellar tides. While dynamical stellar tides are possible, they require a subgiant structure for the star, whereas stellar models point toward a main sequence host. Here, we show that these hitherto irreconcilable observations might be explained by planetary obliquity tides if planet b's spin vector is trapped in a high-obliquity state maintained by a secular spin-orbit resonance with an unseen exterior perturbing planet. We derive constraints on the obliquity ($\epsilon\gtrsim50^{\circ}$), reduced tidal quality factor ($Q^{\prime}\sim10^{6}-10^{7}$), and perturbing planet parameters ($M_{2}\sim10-20M_{\oplus}$, $a_2\lesssim0.04\,{\rm AU}$) required to generate the observed orbital decay. Direct N-body simulations that include tidal and spin dynamics reinforce the plausibility of the scenario. Furthermore, we show that the resonance could have been captured when planet b's obliquity was small, making the proposed sequence of events easy to explain. The hypothetical perturbing planet is within the limits of current radial velocity constraints on the system yet is also detectable. If it exists, it could provide evidence in favor of the in situ formation hypothesis for hot Jupiters.

We report the first confirmation of a hot Jupiter discovered by the Transiting Exoplanet Survey Satellite (TESS) mission: HD 202772A b. The transit signal was detected in the data from TESS Sector 1, and was confirmed to be of planetary origin through radial-velocity measurements. HD 202772A b is orbiting a mildly evolved star with a period of 3.3 days. With an apparent magnitude of V = 8.3, the star is among the brightest known to host a hot Jupiter. Based on the 27days of TESS photometry, and radial velocity data from the CHIRON and HARPS spectrographs, the planet has a mass of 1.008+/-0.074 M_J and radius of 1.562+/-0.053 R_J , making it an inflated gas giant. HD 202772A b is a rare example of a transiting hot Jupiter around a quickly evolving star. It is also one of the most strongly irradiated hot Jupiters currently known.

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.

Gliese 876 harbors one of the most dynamically rich and well-studied exoplanetary systems. The nearby M4V dwarf hosts four known planets, the outer three of which are trapped in a Laplace mean-motion resonance. A thorough characterization of the complex resonant perturbations exhibited by the orbiting planets, and the chaotic dynamics therein, is key to a complete picture of the system's formation and evolutionary history. Here we present a reanalysis of the system using six years of new radial velocity (RV) data from four instruments. This new data augments and more than doubles the size of the decades-long collection of existing velocity measurements. We provide updated estimates of the system parameters by employing a computationally efficient Wisdom-Holman N-body symplectic integrator, coupled with a Gaussian Process (GP) regression model to account for correlated stellar noise. Experiments with synthetic RV data show that the dynamical characterization of the system can differ depending on whether a white noise or correlated noise model is adopted. Despite there being a region of stability for an additional planet in the resonant chain, we find no evidence for one. Our new parameter estimates place the system even deeper into resonance than previously thought and suggest that the system might be in a low energy, quasi-regular double apsidal corotation resonance. This result and others will be used in a subsequent study on the primordial migration processes responsible for the formation of the resonant chain.

The widespread prevalence of close-in, nearly coplanar super-Earth- and sub-Neptune-sized planets in multiple-planet systems was one of the most surprising results from the Kepler mission. By studying a uniform sample of Kepler "multis" with mass measurements from transit timing variations (TTVs), we show that a given planetary system tends to harbor a characteristic type of planet. That is, planets in a system have both masses and radii that are far more similar than if the system were assembled randomly from planets in the population. This finding has two important ramifications. First, the large intrinsic compositional scatter in the planet mass-radius relation is dominated by system-to-system variance rather than intra-system variance. Second, if provided enough properties of the star and primordial protoplanetary disk, there may be a substantial degree of predictability in the outcome of the planet formation process. We show that stellar mass and metallicity account for of order $20\%$ of the variation in outcomes; the remainder is as-yet unknown.

The optical, full-phase photometric variations of a short-period planet provide a unique view of the planet's atmospheric composition and dynamics. The number of planets with optical phase curve detections, however, is currently too small to study them as an aggregate population, motivating an extension of the search to non-transiting planets. Here we present an algorithm for the detection of non-transiting, short-period giant planets in the Kepler field. The procedure uses the phase curves themselves as evidence for the planets' existence. We employ a supervised learning algorithm to recognize the salient time-dependent properties of synthetic phase curves; we then search for detections of signals that match these properties. After demonstrating the algorithm's capabilities, we classify 142,630 FGK Kepler stars without confirmed planets or KOIs and, for each one, assign a probability of a phase curve of a non-transiting planet being present. We identify 60 high-probability non-transiting hot Jupiter candidates. We also derive constraints on the candidates' albedos and offsets of the phase curve maxima. These targets are strong candidates for follow-up radial velocity confirmation and characterization. Once confirmed, the atmospheric information content in the phase curves may be studied in yet greater detail.

A number of authors have proposed that the statistically significant orbital alignment of the most distant Kuiper Belt Objects (KBOs) is evidence of an as-yet undetected planet in the outer solar system, now referred to colloquially a "Planet Nine". Dynamical simulations by Batygin & Brown (2016) have provided constraints on the range of the planet's possible orbits and sky locations. We extend these investigations by exploring the suggestion of Malhotra et al. (2016) that Planet Nine is in small integer ratio mean-motion resonances (MMRs) with several of the most distant KBOs. We show that the observed KBO semi-major axes present a set of commensurabilities with an unseen planet at $\sim 654 \ \mathrm{AU}$ ($P \sim 16,725 \ \mathrm{yr}$) that has a greater than $98\%$ chance of stemming from a sequence of MMRs rather than from a random distribution. We describe and implement a Monte-Carlo optimization scheme that drives billion-year dynamical integrations of the outer solar system to pinpoint the orbital properties of perturbers that are capable of maintaining the KBOs' apsidal alignment. This optimization exercise suggests that the unseen planet is most consistently represented with mass, $m \sim 6-12 M_{\oplus}$, semi-major axis, $a \sim 654 \ \mathrm{AU}$, eccentricity, $e\sim0.45$, inclination, $i \sim 30^{\circ}$, argument of periastron, $\omega \sim 150^{\circ}$, longitude of ascending node, $\Omega \sim 50^{\circ}$, and mean anomaly, $M \sim 180^{\circ}$. A range of sky locations relative to this fiducial ephemeris are possible. We find that the region $30^{\circ}\, \lesssim \mathrm{RA} \lesssim 50^{\circ}$, $-20^{\circ}\, \lesssim \mathrm{Dec} \lesssim 20^{\circ}$ is promising.

We outline a photometric method for detecting the presence of a non-transiting short-period giant planet in a planetary system harboring one or more longer period transiting planets. Within a prospective system of the type that we consider, a hot Jupiter on an interior orbit inclined to the line-of-sight signals its presence through approximately sinusoidal full-phase photometric variations in the stellar light curve, correlated with astrometrically induced transit timing variations for exterior transiting planets. Systems containing a hot Jupiter along with a low-mass outer planet or planets on inclined orbits are a predicted hallmark of in situ accretion for hot Jupiters, and their presence can thus be used to test planetary formation theories. We outline the prospects for detecting non-transiting hot Jupiters using photometric data from typical \textit{Kepler} objects of interest (KOIs). As a demonstration of the technique, we perform a brief assessment of \textit{Kepler} candidates and identify a potential non-transiting hot Jupiter in the KOI-1822 system. Candidate non-transiting hot Jupiters can be readily confirmed with a small number of Doppler velocity observations, even for stars with $V\gtrsim14$.

The standard model for eruptive flares has in the past few years been extended to 3D. It predicts typical J-shaped photospheric footprints of the coronal current layer, forming at similar locations as the Quasi-Separatrix Layers (QSLs). Such a morphology is also found for flare ribbons observed in the EUV band, as well as in non-linear force-free field (NLFFF) magnetic field extrapolations and models. We study the evolution of the photospheric traces of the current density and flare ribbons, both obtained with the SDO instruments. We investigate the photospheric current evolution during the 6 September 2011 X-class flare (SOL2011-09-06T22:20) from observational data of the magnetic field obtained with HMI. This evolution is compared with that of the flare ribbons observed in the EUV filters of the AIA. We also compare the observed electric current density and the flare ribbon morphology with that of the QSLs computed from the flux rope insertion method/NLFFF model. The NLFFF model shows the presence of a fan-spine configuration of overlying field lines, due to the presence of a parasitic polarity, embedding an elongated flux rope that appears in the observations as two parts of a filament. The QSLs, evolved via a magnetofrictional method, also show similar morphology and evolution as both the current ribbons and the EUV flare ribbons obtained at several times during the flare. For the first time, we propose a combined analysis of the photospheric traces of an eruptive flare, in a complex topology, with direct measurements of electric currents and QSLs from observational data and a magnetic field model. The results, obtained by two different and independent approaches, 1) confirm previous results of current increase during the impulsive phase of the flare, 2) show how NLFFF models can capture the essential physical signatures of flares even in a complex magnetic field topology.