We report the detection of a transiting, dense Neptune planet candidate orbiting the bright ($V=8.6$) K0.5V star HD 95338. Detection of the 55-day periodic signal comes from the analysis of precision radial velocities from the Planet Finder Spectrograph on the Magellan II Telescope. Follow-up observations with HARPS also confirm the presence of the periodic signal in the combined data. HD 95338 was also observed by the Transiting Exoplanet Survey Satellite ({\it TESS}) where we identify a clear single transit in the photometry. A Markov Chain Monte Carlo period search on the velocities allows strong constraints on the expected transit time, matching well the epoch calculated from \tess{} data, confirming both signals describe the same companion. A joint fit model yields an absolute mass of 42.44$^{+2.22}_{-2.08} M_{\oplus}$ and a radius of 3.89$^{+0.19}_{-0.20}$ $R_{\oplus}$ which translates to a density of 3.98$^{+0.62}_{-0.64}$ \gcm\, for the planet. Given the planet mass and radius, structure models suggest it is composed of a mixture of ammonia, water, and methane. HD 95338\,b is one of the most dense Neptune planets yet detected, indicating a heavy element enrichment of $\sim$90\% ($\sim38\, M_{\oplus}$). This system presents a unique opportunity for future follow-up observations that can further constrain structure models of cool gas giant planets.

We present a transmission spectrum for the Neptune-size exoplanet HD 106315 c from optical to infrared wavelengths based on transit observations from the Hubble Space Telescope/Wide Field Camera 3, K2, and Spitzer. The spectrum shows tentative evidence for a water absorption feature in the $1.1 - 1.7\mu$m wavelength range with a small amplitude of 30 ppm (corresponding to just $0.8 \pm 0.04$ atmospheric scale heights). Based on an atmospheric retrieval analysis, the presence of water vapor is tentatively favored with a Bayes factor of 1.7 - 2.6 (depending on prior assumptions). The spectrum is most consistent with either enhanced metallicity, high altitude condensates, or both. Cloud-free solar composition atmospheres are ruled out at $>5\sigma$ confidence. We compare the spectrum to grids of cloudy and hazy forward models and find that the spectrum is fit well by models with moderate cloud lofting or haze formation efficiency, over a wide range of metallicities ($1 - 100\times$ solar). We combine the constraints on the envelope composition with an interior structure model and estimate that the core mass fraction is $\gtrsim0.3$. With a bulk composition reminiscent of that of Neptune and an orbital distance of 0.15 AU, HD 106315 c hints that planets may form out of broadly similar material and arrive at vastly different orbits later in their evolution.

Aerosols are common in the atmospheres of exoplanets across a wide swath of temperatures, masses, and ages. These aerosols strongly impact observations of transmitted, reflected, and emitted light from exoplanets, obfuscating our understanding of exoplanet thermal structure and composition. Knowing the dominant aerosol composition would facilitate interpretations of exoplanet observations and theoretical understanding of their atmospheres. A variety of compositions have been proposed, including metal oxides and sulphides, iron, chromium, sulphur, and hydrocarbons. However, the relative contributions of these species to exoplanet aerosol opacity is unknown. Here we show that the aerosol composition of giant exoplanets observed in transmission is dominated by silicates and hydrocarbons. By constraining an aerosol microphysics model with trends in giant exoplanet transmission spectra, we find that silicates dominate aerosol opacity above planetary equilibrium temperatures of 950 K due to low nucleation energy barriers and high elemental abundances, while hydrocarbon aerosols dominate below 950 K due to an increase in methane abundance. Our results are robust to variations in planet gravity and atmospheric metallicity within the range of most giant transiting exoplanets. We predict that spectral signatures of condensed silicates in the mid-infrared are most prominent for hot (>1600 K), low-gravity (<10 m s$^{-2}$) objects.

Understanding the anomalous radii of many transiting hot gas giant planets is a fundamental problem of planetary science. Recent detections of re-inflated warm Jupiters orbiting post-main-sequence stars and the re-inflation of hot Jupiters while their host stars evolve on the main-sequence may help constrain models for the anomalous radii of hot Jupiters. In this work, we present evolution models studying the re-inflation of gas giants to determine how varying the depth and intensity of deposited heating affects both main-sequence re-inflation of hot Jupiters and post-main-sequence re-inflation of warm Jupiters. We find that deeper heating is required to re-inflate hot Jupiters than is needed to suppress their cooling, and that the timescale of re-inflation decreases with increasing heating rate and depth. We find a strong degeneracy between heating rate and depth, with either strong shallow heating or weak deep heating providing an explanation for main-sequence re-inflation of hot Jupiters. This degeneracy between heating rate and depth can be broken in the case of post-main-sequence re-inflation of warm Jupiters, as the inflation must be rapid to occur within post-main-sequence evolution timescales. We also show that the dependence of heating rate on incident stellar flux inferred from the sample of hot Jupiters can explain re-inflation of both warm and hot Jupiters. TESS will obtain a large sample of warm Jupiters orbiting post-main-sequence stars, which will help to constrain the mechanism(s) causing the anomalous radii of gas giant planets.

In giant planet atmosphere modelling, the intrinsic temperature $T_\mathrm{int}$ and radiative-convective boundary (RCB) are important lower boundary conditions. Often in one-dimensional radiative-convective models and in three-dimensional general circulation models it is assumed that $T_\mathrm{int}$ is similar to that of Jupiter itself, around 100 K, which yields a RCB around 1 kbar for hot Jupiters. In this work, we show that the inflated radii, and hence high specific entropy interiors, of hot Jupiters suggest much higher $T_\mathrm{int}$ values. Assuming the effect is primarily due to current heating (rather than delayed cooling), we derive an equilibrium relation between $T_\mathrm{eq}$ and $T_\mathrm{int}$, showing that the latter can take values as high as 700 K. In response, the RCB moves upward in the atmosphere. Using one-dimensional radiative-convective atmosphere models, we find RCBs of only a few bars, rather than the kilobar typically supposed. This much shallower RCB has important implications for the atmospheric structure, vertical and horizontal circulation, interpretations of phase curves, and the effect of deep cold traps on cloud formation.

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.

The gravity field of a giant planet is typically our best window into its interior structure and composition. Through comparison of a model planet's calculated gravitational potential with the observed potential, inferences can be made about interior quantities, including possible composition and the existence of a core. Necessarily, a host of assumptions go into such calculations, making every inference about a giant planet's structure strongly model dependent. In this work we present a more general picture by setting Saturn's gravity field, as measured during the \emph{Cassini} Grand Finale, as a likelihood function driving a Markov-Chain Monte Carlo exploration of the possible interior density profiles. The result is a posterior distribution of the interior structure that is not tied to assumed composition, thermal state, or material equations of state. Constraints on interior structure derived in this Bayesian framework are necessarily less informative, but are also less biased and more general. These empirical and probabilistic constraints on the density structure are our main data product which we archive for continued analysis. We find that the outer half of Saturn's radius is relatively well constrained, and we interpret our findings as suggesting a significant metal enrichment, in line with atmospheric abundances from remote sensing. As expected, the inner half of Saturn's radius is less well-constrained by gravity, but we generally find solutions that include a significant density enhancement, which can be interpreted as a core, although this core is often lower in density and larger in radial extent than typically found by standard models. This is consistent with a dilute core and/or composition gradients.

The relationship between the compositions of giant planets and their host stars is of fundamental interest in understanding planet formation. The solar system giant planets are enhanced above solar composition in metals, both in their visible atmospheres and bulk compositions. A key question is whether the metal enrichment of giant exoplanets is correlated with that of their host stars. Thorngren et al. (2016) showed that in cool (Teq < 1000 K) giant exoplanets, the total heavy-element mass increases with total Mp and the heavy element enrichment relative to the parent star decreases with total Mp. In their work, the host star metallicity was derived from literature [Fe/H] measurements. Here we conduct a more detailed and uniform study to determine whether different host star metals (C, O, Mg, Si, Fe, and Ni) correlate with the bulk metallicity of their planets, using correlation tests and Bayesian linear fits. We present new host star abundances of 19 cool giant planet systems, and combine these with existing host star data for a total of 22 cool giant planet systems (24 planets). Surprisingly, we find no clear correlation between stellar metallicity and planetary residual metallicity (the relative amount of metal versus that expected from the planet mass alone), which is in conflict with common predictions from formation models. We also find a potential correlation between residual planet metals and stellar volatile-to-refractory element ratios. These results provide intriguing new relationships between giant planet and host star compositions for future modeling studies of planet formation.

Probabilistic relationships between the mass and radius of planets are useful for a variety of purposes, including estimating the yields of planet discovery efforts and the radii of discovered planets given their masses. Previous work on giant planets often do not make the important distinction between cool and (inflated) hot giants. This can lead to needlessly large uncertainties and overestimated radii for the cooler objects. Using hierarchical bayesian modeling and the AIC model selection criterion, we identify an empirical mass-radius relationship for cool giant planets above 15 Earth masses. It takes the form of a quadratic function in log mass, with a power-law dispersion. Parametric uncertainties from MCMC are provided, and we compare our results with the outputs of interior structure models and the radii of the solar system giant planets, which fit comfortably within our relation.

We present new 3.6 and 4.5 micron secondary eclipse measurements for five cool (less than approximately 1000 K) transiting gas giant planets: HAT-P-15b, HAT-P-17b, HAT-P-18b, HAT-P-26b, and WASP-69b. We detect eclipses in at least one bandpass for all planets except HAT-P-15b. We confirm and refine the orbital eccentricity of HAT-P-17b, which is also the only planet in our sample with a known outer companion. We compare our measured eclipse depths in these two bands, which are sensitive to the relative abundances of methane versus carbon monoxide and carbon dioxide, respectively, to predictions from 1D atmosphere models for each planet. For planets with hydrogen-dominated atmospheres and equilibrium temperatures cooler than approximately 1000 K, this ratio should vary as a function of both atmospheric metallicity and the carbon-to-oxygen ratio. For HAT-P-26b, our observations are in good agreement with the low atmospheric metallicity inferred from transmission spectroscopy. We find that all four of the planets with detected eclipses are best matched by models with relatively efficient circulation of energy to the nightside. We see no evidence for a solar-system-like correlation between planet mass and atmospheric metallicity, but instead identify a potential (1.9 sigma) correlation between the inferred methane/(carbon monoxide + carbon dioxide) ratio and stellar metallicity. Our ability to characterize this potential trend is limited by the relatively large uncertainties in the stellar metallicity values. Our observations provide a first look at the brightness of these planets at wavelengths accessible to the James Webb Space Telescope, which will be able to resolve individual methane, carbon monoxide, and carbon dioxide bands and provide much stronger constraints on their atmospheric compositions.

Atmospheric characterization through spectroscopic analysis, an essential tool of modern exoplanet science, can benefit significantly from the context provided by interior structure models. In particular, the planet's bulk metallicity $Z_p$ places an upper limit on potential atmospheric metallicity. Here we construct interior structure models to derive $Z_p$ and atmospheric metallicity upper limits for 403 known transiting giant exoplanets. These limits are low enough that they can usefully inform atmosphere models. Additionally, we argue that comparing $Z_p$ to the observed atmospheric metallicity gives a useful measure of how well-mixed metals are within the planet. This represents a new avenue for learning about planetary interiors. To aid in the future characterization of new planet discoveries we derive analytic prior predictions of atmosphere metallicity as a function of planet mass, and evaluate the effectiveness of our approach on Jupiter and Saturn. We we include log-linear fits for approximating the metallicities of planets not in our catalog.

The recently launched Transiting Exoplanet Survey Satellite (TESS) is expected to produce many new exoplanet discoveries which will be especially amenable to follow-up study. Assessments of the planet discovery yield of TESS, such as Sullivan et al. (2015) and Barclay et al. (2018), will be important for planning follow-up work. Analyzing these predicted planet samples, however, we find that giant planet radii derived from the current bulk transiting planet sample have been used at all potential orbits without accounting for the temperature dependence of radius inflation. These phantom inflated planets (PIPs) make up just over 1.1% of the Sullivan et al. (2015) predicted population and about 8% of the Barclay et al. (2018) sample. Similar population predictions for direct imaging studies should likewise take care not to include inflated planets at large separations as such planets will appear larger and, depending on assumed albedo, brighter and more easily detectable than physically possible. Despite their high false positive rate, giant planets are the testbeds for atmospheric characterization techniques and care should be taken to understand and account for potential contaminating factors in this population.

The super-Neptune exoplanet WASP-107b is an exciting target for atmosphere characterization. It has an unusually large atmospheric scale height and a small, bright host star, raising the possibility of precise constraints on its current nature and formation history. We report the first atmospheric study of WASP-107b, a Hubble Space Telescope measurement of its near-infrared transmission spectrum. We determined the planet's composition with two techniques: atmospheric retrieval based on the transmission spectrum and interior structure modeling based on the observed mass and radius. The interior structure models set a $3\,\sigma$ upper limit on the atmospheric metallicity of $30\times$ solar. The transmission spectrum shows strong evidence for water absorption ($6.5\,\sigma$ confidence), and the retrieved water abundance is consistent with expectations for a solar abundance pattern. The inferred carbon-to-oxygen ratio is subsolar at $2.7\,\sigma$ confidence, which we attribute to possible methane depletion in the atmosphere. The spectral features are smaller than predicted for a cloud-free composition, crossing less than one scale height. A thick condensate layer at high altitudes (0.1 - 3 mbar) is needed to match the observations. We find that physically motivated cloud models with moderate sedimentation efficiency ($f_\mathrm{sed} = 0.3$) or hazes with a particle size of 0.3 $\mu$m reproduce the observed spectral feature amplitude. Taken together, these findings serve as an illustration of the diversity and complexity of exoplanet atmospheres. The community can look forward to more such results with the high precision and wide spectral coverage afforded by future observing facilities.

The cause of hot Jupiter radius inflation, where giant planets with $T_{\rm eq}$ $>1000$ K are significantly larger than expected, is an open question and the subject of many proposed explanations. Rather than examine these models individually, this work seeks to characterize the anomalous heating as a function of incident flux, $\epsilon(F)$, needed to inflate the population of planets to their observed sizes. We then compare that result to theoretical predictions for various models. We examine the population of about 300 giant planets with well-determined masses and radii and apply thermal evolution and Bayesian statistical models to infer the anomalous power as a function of incident flux that best reproduces the observed radii. First, we observe that the inflation of planets below about $M=0.5 \;\rm{M}_\rm{J}$ appears very different than their higher mass counterparts, perhaps as the result of mass loss or an inefficient heating mechanism. As such, we exclude planets below this threshold. Next, we show with strong significance that $\epsilon(F)$ increases with $T_{\rm{eq}}$ towards a maximum of $\sim 2.5\%$ at $T_{\rm{eq}} \approx 1500$ K, and then decreases as temperatures increase further, falling to $\sim0.2\%$ at $T_\rm{eff}= 2500$ K. This high-flux decrease in inflation efficiency was predicted by the Ohmic dissipation model of giant planet inflation but not other models. We also explicitly check the thermal tides model and find that it predicts far more variance in radii than is observed. Thus, our results provide evidence for the Ohmic dissipation model and a functional form for $\epsilon(F)$ that any future theories of hot Jupiter radii can be tested against.

A large fraction of known Jupiter like exoplanets are inflated as compared to Jupiter. These "hot" Jupiters orbit close to their parent star and are bombarded with intense starlight. Many theories have been proposed to explain their radius inflation and several suggest that a small fraction of the incident starlight is injected in to the planetary interior which helps to puff up the planet. How will such energy injection affect the planetary dynamo? In this Letter, we estimate the surface magnetic field strength of hot Jupiters using scaling arguments that relate energy available in planetary interiors to the dynamo generated magnetic fields. We find that if we take into account the energy injected in the planetary interior that is sufficient to inflate hot Jupiters to observed radii, then the resulting dynamo should be able generate magnetic fields that are more than an order of magnitude stronger than the Jovian values. Our analysis highlights the potential fundamental role of the stellar light in setting the field strength in hot Jupiters.

We predict the carbon-to-oxygen (C/O) ratios in the hydrogen-helium envelope and atmospheres of a sample of nearly 50 relatively cool ($T_{\mathrm eq}<$ 1000 K) transiting gas giant planets. The method involves planetary envelope metallicity estimates that use the structure models of Thorngren et al. (2016) and the disk and planetary accretion model of \"Oberg et al. (2011). We find that nearly all of these planets are strongly metal-enriched which, coupled with the fact that solid material is the main deliverer of metals in the protoplanetary disk, implies that the substellar C/O ratios of their accreted solid material dominate compared to the enhanced C/O ratio of their accreted gaseous component. We predict that these planets will have atmospheres that are typically reduced in their C/O compared to parent star values independent of the assessed formation locations, with C/O $<1$ a nearly universal outcome within the framework of the model. We expect water vapor absorption features to be ubiquitous in the atmospheres of these planets, and by extension, other gas giants.

The Neptune-mass GJ 436b is one of the most-studied transiting exoplanets with repeated measurements of both its thermal emission and transmission spectra. We build on previous studies to answer outstanding questions about this planet, including its potentially high metallicity and tidal heating of its interior. We present new observations of GJ 436b's thermal emission at 3.6 and 4.5 micron, which reduce uncertainties in estimates of GJ 436b's flux at those wavelengths and demonstrate consistency between Spitzer observations spanning more than 7 years. We analyze the Spitzer thermal emission photometry and Hubble WFC3 transmission spectrum in tandem. We use a powerful dual-pronged modeling approach, comparing these data to both self-consistent and retrieval models. We vary the metallicity, intrinsic luminosity from tidal heating, disequilibrium chemistry, and heat redistribution. We also study the effect of clouds and photochemical hazes on the spectra, but do not find strong evidence for either. The self-consistent and retrieval modeling combine to suggest that GJ 436b has a high atmospheric metallicity, with best fits at or above several hundred times solar metallicity, tidal heating warming its interior with best-fit intrinsic effective effective temperatures around 300--350 K, and disequilibrium chemistry. High metal-enrichments (>600x solar) can only occur from the accretion of rocky, rather than icy, material. Assuming Tint~300--350 K, we find that Q'~2x10^5--10^6, larger than Neptune's Q', and implying a long tidal circularization timescale for the planet's orbit. We suggest that Neptune-mass planets may be a more diverse class than previously imagined, with metal-enhancements potentially spanning several orders of magnitude, to perhaps over 1000x solar metallicity. High fidelity observations with instruments like JWST will be critical for characterizing this diversity.

Exoplanet discoveries of recent years have provided a great deal of new data for studying the bulk compositions of giant planets. Here we identify 47 transiting giant planets ($20 M_\oplus < M < 20 M_{\mathrm{J}}$) whose stellar insolation is low enough ($F_* < 2\times10^8\; \text{erg}\; \text{s}^{-1}\; \text{cm}^{-2}$, or roughly $T_\text{eff} < 1000$) that they are not affected by the hot Jupiter radius inflation mechanism(s). We compute a set of new thermal and structural evolution models and use these models in comparison with properties of the 47 transiting planets (mass, radius, age) to determine their heavy element masses. A clear correlation emerges between the planetary heavy element mass $M_z$ and the total planet mass, approximately of the form $M_z \propto \sqrt{M}$. This finding is consistent with the core accretion model of planet formation. We also study how stellar metallicity [Fe/H] affects planetary metal-enrichment and find a weaker correlation than has been previously reported from studies with smaller sample sizes. We confirm a strong relationship between the planetary metal-enrichment relative to the parent star $Z_{\rm planet}/Z_{\rm star}$ and the planetary mass, but see no relation in $Z_{\rm planet}/Z_{\rm star}$ with planet orbital properties or stellar mass. The large heavy element masses of many planets ($>50$ $M_{\oplus}$) suggest significant amounts of heavy elements in H/He envelopes, rather than cores, such that metal-enriched giant planet atmospheres should be the rule. We also discuss a model of core-accretion planet formation in a one-dimensional disk and show that it agrees well with our derived relation between mass and $Z_{\rm planet}/Z_{\rm star}$.

The Noble Element Simulation Technique (NEST) is an exhaustive collection of models explaining both the scintillation light and ionization yields of noble elements as a function of particle type (nuclear recoil, electron recoil, alphas), electric field, and incident energy or energy loss dE/dx. It is packaged as C++ code for Geant4 that implements said models, overriding the default model which does not account for certain complexities, such as the reduction in yields for nuclear recoils (NR) compared to electron recoils (ER). We present here improvements to the existing NEST models and updates to the code which make the package even more realistic and turn it into a more full-fledged Monte Carlo simulation. All available liquid xenon data on NR and ER to date have been taken into consideration in arriving at the current models. Furthermore, NEST addresses the question of the magnitude of the light and charge yields of nuclear recoils, including their electric field dependence, thereby shedding light on the possibility of detection or exclusion of a low-mass dark matter WIMP by liquid xenon detectors.