Recently, we introduced PLanetary Atmospheric Tool for Observer Noobs (PLATON), a Python package that calculates model transmission spectra for exoplanets and retrieves atmospheric characteristics based on observed spectra. We now expand its capabilities to include the ability to compute secondary eclipse depths. We have also added the option to calculate models using the correlated-$k$ method for radiative transfer, which improves accuracy without sacrificing speed. Additionally, we update the opacities in PLATON--many of which were generated using old or proprietary line lists--using the most recent and complete public line lists. These opacities are made available at R=1000 and R=10,000 over the 0.3--30 um range, and at R=375,000 in select near IR bands, making it possible to utilize PLATON for ground-based high resolution cross correlation studies. To demonstrate PLATON's new capabilities, we perform a retrieval on published HST and Spitzer transmission and emission spectra of the archetypal hot Jupiter HD 189733b. This is the first joint transit and secondary eclipse retrieval for this planet in the literature, as well as the most comprehensive set of both transit and eclipse data assembled for a retrieval to date. We find that these high signal-to-noise data are well-matched by atmosphere models with a C/O ratio of $0.66_{-0.09}^{+0.05}$ and a metallicity of $12_{-5}^{+8}$ times solar where the terminator is dominated by extended nanometer-sized haze particles at optical wavelengths. These are among the smallest uncertainties reported to date for an exoplanet, demonstrating both the power and the limitations of HST and Spitzer exoplanet observations.

We present the transmission spectrum of HAT-P-12b through a joint analysis of data obtained from the Hubble Space Telescope Space Telescope Imaging Spectrograph (STIS) and Wide Field Camera 3 (WFC3) and Spitzer, covering the wavelength range 0.3-5.0 $\mu$m. We detect a muted water vapor absorption feature at 1.4 $\mu$m attenuated by clouds, as well as a Rayleigh scattering slope in the optical indicative of small particles. We interpret the transmission spectrum using both the state-of-the-art atmospheric retrieval code SCARLET and the aerosol microphysics model CARMA. These models indicate that the atmosphere of HAT-P-12b is consistent with a broad range of metallicities between several tens to a few hundred times solar, a roughly solar C/O ratio, and moderately efficient vertical mixing. Cloud models that include condensate clouds do not readily generate the sub-micron particles necessary to reproduce the observed Rayleigh scattering slope, while models that incorporate photochemical hazes composed of soot or tholins are able to match the full transmission spectrum. From a complementary analysis of secondary eclipses by Spitzer, we obtain measured depths of $0.042\%\pm0.013\%$ and $0.045\%\pm0.018\%$ at 3.6 and 4.5 $\mu$m, respectively, which are consistent with a blackbody temperature of $890^{+60}_{-70}$ K and indicate efficient day-night heat recirculation. HAT-P-12b joins the growing number of well-characterized warm planets that underscore the importance of clouds and hazes in our understanding of exoplanet atmospheres.

We present ground-based infrared transit observations for four dynamically interacting \textit{Kepler} planets, including Kepler-29b, Kepler-36c, KOI-1783.01, and Kepler-177c, obtained using the Wide-field Infrared Camera on the Hale 200" telescope at Palomar Observatory. By utilizing an engineered diffuser and custom guiding software, we mitigate time-correlated telluric and instrumental noise sources in these observations. We achieve an infrared photometric precision comparable to or better than that of space-based observatories such as the \textit{Spitzer Space Telescope}, and detect transits with greater than 3$\sigma$ significance for all planets. For Kepler-177c ($J=13.9$) our measurement uncertainties are only $1.2\times$ the photon noise limit and 1.9 times better than the predicted photometric precision for \textit{Spitzer} IRAC photometry of this same target. We find that a single transit observation obtained $4-5$ years after the end of the original \textit{Kepler} mission can reduce dynamical mass uncertainties by as much as a factor of three for these systems. Additionally, we combine our new observations of KOI-1783.01 with information from the literature to confirm the planetary nature of this system. We discuss the implications of our new mass and radius constraints in the context of known exoplanets with low incident fluxes, and we note that Kepler-177c may be a more massive analog to the currently known super-puffs given its core mass (3.8$\pm0.9M_\Earth$) and large gas-to-core ratio (2.8$\pm0.7$). Our demonstrated infrared photometric performance opens up new avenues for ground-based observations of transiting exoplanets previously thought to be restricted to space-based investigation.

We present the first comprehensive look at the $0.35-5$ $\mu$m transmission spectrum of the warm ($\sim 800$ K) Neptune HAT-P-11b derived from thirteen individual transits observed using the Hubble and Spitzer Space Telescopes. Along with the previously published molecular absorption feature in the $1.1-1.7$ $\mu$m bandpass, we detect a distinct absorption feature at 1.15 $\mu$m and a weak feature at 0.95 $\mu$m, indicating the presence of water and/or methane with a combined significance of 4.4 $\sigma$. We find that this planet's nearly flat optical transmission spectrum and attenuated near-infrared molecular absorption features are best-matched by models incorporating a high-altitude cloud layer. Atmospheric retrievals using the combined $0.35-1.7$ $\mu$m HST transmission spectrum yield strong constraints on atmospheric cloud-top pressure and metallicity, but we are unable to match the relatively shallow Spitzer transit depths without under-predicting the strength of the near-infrared molecular absorption bands. HAT-P-11b's HST transmission spectrum is well-matched by predictions from our microphysical cloud models. Both forward models and retrievals indicate that HAT-P-11b most likely has a relatively low atmospheric metallicity ($<4.6 \; Z_{\odot}$ and $<86 \; Z_{\odot}$ at the $2 \sigma$ and $3 \sigma$ levels respectively), in contrast to the expected trend based on the solar system planets. Our work also demonstrates that the wide wavelength coverage provided by the addition of the HST STIS data is critical for making these inferences.

The presence of planets in binary systems poses interesting problems for planet formation theories, both in cases where planets must have formed in very compact discs around the individual stars and where they are located near the edge of the stable circumbinary region, where in situ formation is challenging. Dust dynamics is expected to play an important role in such systems, since dust trapping at the inner edge of circumbinary discs could aid in situ formation, but would simultaneously starve the circumstellar discs of the solid material needed to form planets. Here we investigate the dynamics of dust in binary systems using Smooth Particle Hydrodynamics. We find that all our simulations tend towards dust trapping in the circumbinary disc, but the timescale on which trapping begins depends on binary mass ratio ($q$) and eccentricity as well as the angular momentum of the infalling material. For $q \gtrsim 0.1$, we find that dust can initially accrete onto the circumstellar discs, but as the circumbinary cavity grows in radius, dust eventually becomes trapped in the circumbinary disc. For $q = 0.01$, we find that increasing the binary eccentricity increases the time required for dust trapping to begin. However, even this longer timescale is likely to be shorter than the planet formation timescale in the inner disc and is insufficient to explain the observed pre-transitional discs. This indicates that increase in companion eccentricity alone is not enough to allow significant transfer of solids from the outer to the inner disc.

We introduce PLanetary Atmospheric Transmission for Observer Noobs (PLATON), a Python package that calculates transmission spectra for exoplanets and retrieves atmospheric characteristics based on observed spectra. PLATON is easy to install and use, with common use cases taking no more than a few lines of code. It is also fast, with the forward model taking much less than one second and a typical retrieval finishing in minutes on an ordinary desktop. PLATON supports the most common atmospheric parameters, such as temperature, metallicity, C/O ratio, cloud-top pressure, and scattering slope. It also has less commonly included features, such as a Mie scattering cloud model and unocculted starspot corrections. The code is available online at https://github.com/ideasrule/platon under the open source GPL-3.0 license.

The higher order gravitational moments of a differentially rotating planet are greatly affected or even dominated by the near surface differential flow. Unlike the contribution from rigid body rotation, this part of the gravity field can be well estimated by linear theory when the differential rotation is on cylinders and the corresponding gravity field arises from the higher order moments directly introduced by that flow. In the context of an n = 1 polytrope, we derive approximate analytical formulas for the gravity moments. We find that $\Delta J_{2\mathrm{n}}$ is typically at most a few times $(-1)^{n+1} \; a \; q \; d^{5/2}$, where $a$ is the amplitude of the differential rotation (as a fraction of the background rigid body rotation), $q = \Omega^2 R^3 / G M$ is the usual dimensionless measure of rotation for the planet (mass $M$, radius $R$), and $d << 1$ is the characteristic depth of the flow as a fraction of the planetary radius. Applied to Saturn, with $a$ set by the observed surface wind amplitude, we find first that the observed signs of the $\Delta J_{2\mathrm{n}}$ are a trivial consequence of the definition of the corresponding Legendre polynomials, but the {\it Cassini} observations can not be explained by a simple exponentially decaying flow and instead require a substantial retrograde flow at depth and a larger $d$ than the simple scaling suggests. This is consistent with the results reported by \citet{Iess2019}. However, there is no fluid dynamical requirement that the flows observed in the atmosphere are a guide to the flows thousands of km deeper. We explore a wide range of flow depths and amplitudes which yield values for $\Delta J_{\mathrm{n}}$ that are acceptable within the error estimates and thus highlight the inherent non-uniqueness of inferences made from the higher order gravity moments.

Low-mass low-density planets discovered by Kepler in the super-Earth mass regime typically have large radii for their inferred masses, implying the presence of H$_2$-He atmospheres. These planets are vulnerable to atmospheric mass loss due to heating by the parent star's XUV flux. Models coupling atmospheric mass loss with thermal evolution predicted a bimodal distribution of planetary radii, which has gained observational support. However, a key component that has been ignored in previous studies is the dissolution of these gases into the molten core of rock and iron that constitute most of their mass. Such planets have high temperatures ($>$2000 K) and pressures ($\sim$kbars) at the core-envelope boundary, ensuring a molten surface and a subsurface reservoir of hydrogen that can be 5-10 times larger than the atmosphere. This study bridges this gap by coupling the thermal evolution of the planet and the mass loss of the atmosphere with the thermodynamic equilibrium between the dissolved H$_2$ and the atmospheric H$_2$ (Henry's law). Dissolution in the interior allows a planet to build a larger hydrogen repository during the planet formation stage. We show that the dissolved hydrogen outgasses to buffer atmospheric mass loss. The slow cooling of the planet also leads to outgassing because solubility decreases with decreasing temperature. Dissolution of hydrogen in the interior therefore increases the atmosphere retention ability of super-Earths. The study highlights the importance of including the temperature- and pressure-dependent solubility of gases in magma oceans and coupling outgassing to planetary evolution models.