Source author record

Peter Gao

Peter Gao appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

32works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

32 published item(s)

preprint2023arXiv

A First Look at the JWST MIRI/LRS Phase Curve of WASP-43b

We observed a full-orbit phase curve of the hot Jupiter WASP-43b with MIRI/LRS as part of the Transiting Exoplanet Community Early Release Science Program. Here we report preliminary findings for the instrument performance from the team's MIRI Working Group. Overall we find that MIRI's performance for phase curve observations is excellent, with a few minor caveats. The key takeaways for Cycle 2 planning with MIRI/LRS are: (1) long-duration observations (> 24 hours) have now been successfully executed; (2) for phase curves, we recommend including a one-hour burn-in period prior to taking science data to mitigate the effects of the ramp systematic; and (3) we do not yet recommend partial phase curve observations. In addition, we also find that: the position of the spectrum on the detector is stable to within 0.03 pixels over the full 26.5-hour observation; the light curves typically show a systematic downward ramp that is strongest for the first 30 minutes, but continues to decay for hours; from 10.6-11.8 microns, the ramp effect has remarkably different behavior, possibly due to a different illumination history for the affected region of the detector; after trimming the integrations most affected by the initial ramps and correcting the remaining systematics with analytic models, we obtain residuals to the light-curve fits that are typically within 25% of the photon noise limit for 0.5-micron spectroscopic bins; non-linearity correction is not a significant source of additional noise for WASP-43, though it may be an issue for brighter targets; the gain value of 5.5 electrons/DN currently on CRDS and JDox is known to be incorrect, and the current best estimate for the gain is approximately 3.1 electrons/DN; new reference files for the JWST calibration pipeline reflecting these findings are under development at STScI.

preprint2022arXiv

A bimodal distribution of haze in Pluto's atmosphere

Pluto, Titan, and Triton make up a unique class of solar system bodies, with icy surfaces and chemically reducing atmospheres rich in organic photochemistry and haze formation. Hazes play important roles in these atmospheres, with physical and chemical processes highly dependent on particle sizes, but the haze size distribution in reducing atmospheres is currently poorly understood. Here we report observational evidence that Pluto's haze particles are bimodally distributed, which successfully reproduces the full phase scattering observations from New Horizons. Combined with previous simulations of Titan's haze, this result suggests that haze particles in reducing atmospheres undergo rapid shape change near pressure levels ~0.5Pa and favors a photochemical rather than a dynamical origin for the formation of Titan's detached haze. It also demonstrates that both oxidizing and reducing atmospheres can produce multi-modal hazes, and encourages reanalysis of observations of hazes on Titan and Triton.

preprint2022arXiv

A close-in puffy Neptune with hidden friends: The enigma of TOI 620

We present the validation of a transiting low-density exoplanet orbiting the M2.5 dwarf TOI 620 discovered by the NASA TESS mission. We utilize photometric data from both TESS and ground-based follow-up observations to validate the ephemerides of the 5.09-day transiting signal and vet false positive scenarios. High-contrast imaging data are used to resolve the stellar host and exclude stellar companions at separations $\gtrsim 0.2''$. We obtain follow-up spectroscopy and corresponding precise radial velocities (RVs) with multiple PRV spectrographs to confirm the planetary nature of the transiting exoplanet. We calculate a 5$σ$ upper limit of $M_P < 7.1$ M$_\oplus$ and $ρ_P < 0.74$ g cm$^{-3}$, and we identify a non-transiting 17.7-day candidate. We also find evidence for a substellar (1-20 M$_{\rm J}$) companion with a projected separation $\lesssim 20$ au from a combined analysis of Gaia, AO imaging, and RVs. With the discovery of this outer companion, we carry out a detailed exploration of the possibilities that TOI 620 b might instead be a circum-secondary planet or a pair of eclipsing binary stars orbiting the host in a hierarchical triple system. We find, under scrutiny, that we can exclude both of these scenarios from the multi-wavelength transit photometry, thus validating TOI 620 b as a low-density exoplanet transiting the central star in this system. The low density of TOI 620 b makes it one of the most amenable exoplanets for atmospheric characterization, such as with JWST and Ariel, validated or confirmed by the TESS mission to date.

preprint2022arXiv

AU Microscopii in the FUV: Observations in Quiescence, During Flares, and Implications for AU Mic b and c

High energy X-ray and ultraviolet (UV) radiation from young stars impacts planetary atmospheric chemistry and mass loss. The active $\sim 22$ Myr M dwarf AU Mic hosts two exoplanets orbiting interior to its debris disk. Therefore, this system provides a unique opportunity to quantify the effects of stellar XUV irradiation on planetary atmospheres as a function of both age and orbital separation. In this paper we present over 5 hours of Far-UV (FUV) observations of AU Mic taken with the Cosmic Origins Spectrograph (COS; 1070-1360 Angstrom) on the Hubble Space Telescope (HST). We provide an itemization of $120$ emission features in the HST/COS FUV spectrum and quantify the flux contributions from formation temperatures ranging from $10^4-10^7$ K. We detect 13flares in the FUV white-light curve with energies ranging from $10^{29} - 10^{31}$ ergs. The majority of the energy in each of these flares is released from the transition region between the chromosphere and the corona. There is a 100$\times$ increase in flux at continuum wavelengths $λ< 1100$ Angstrom in each flare which may be caused by thermal Bremsstrahlung emission. We calculate that the baseline atmospheric mass-loss rate for AU Mic b is $\sim 10^8$ g s$^{-1}$, although this rate can be as high as $\sim 10^{14}$ g s$^{-1}$ during flares with $L_\textrm{flare} \simeq 10^{33}$ erg s$^{-1}$. Finally, we model the transmission spectra for AU Mic b and c with a new panchromatic spectrum of AU Mic c and motivate future JWST observations of these planets.

preprint2022arXiv

Depletion of gaseous CO in protoplanetary disks by surface-energy-regulated ice formation

Empirical constraints of fundamental properties of protoplanetary disks are essential for understanding planet formation and planetary properties (1,2). Carbon monoxide (CO) gas is often used to constrain disk properties (3). However, estimates show that the CO gas abundance in disks is depleted relative to expected values (4,5,6,7), and models of various disk processes impacting the CO abundance could not explain this depletion on observed 1Myr timescales (8,9,10,11,12,13,14). Here we demonstrate that surface energy effects on particles in disks, such as the Kelvin effect, that arise when ice heterogeneously nucleates onto an existing particle can efficiently trap CO in its ice phase. In previous ice formation models, CO gas was released when small ice-coated particles were lofted to warmed disk layers. Our model can reproduce the observed abundance, distribution and time evolution of gaseous CO in the four most studied protoplanetary disks (7). We constrain the solid and gaseous CO inventory at the midplane and disk diffusivities and resolve inconsistencies in estimates of the disk mass -- three crucial parameters that control planetary formation.

preprint2022arXiv

Microphysics of Water Clouds in the Atmospheres of Y Dwarfs and Temperate Giant Planets

Water clouds are expected to form on Y dwarfs and giant planets with equilibrium temperatures near or below that of Earth, drastically altering their atmospheric compositions and their albedos and thermal emission spectra. Here we use the 1D Community Aerosol and Radiation Model for Atmospheres (CARMA) to investigate the microphysics of water clouds on cool substellar worlds to constrain their typical particle sizes and vertical extent, taking into consideration nucleation and condensation, which have not been considered in detail for water clouds in H/He atmospheres. We compute a small grid of Y dwarf and temperate giant exoplanet atmosphere models with water clouds forming through homogeneous nucleation and heterogeneous nucleation on cloud condensation nuclei composed of meteoritic dust, organic photochemical hazes, and upwelled potassium chloride cloud particles. We present comparisons with the Ackerman & Marley parameterization of cloud physics to extract the optimal sedimentation efficiency parameter (f$_{sed}$) using Virga. We find that no Virga model replicates the CARMA water clouds exactly and that a transition in f$_{sed}$ occurs from the base of the cloud to the cloud top. Furthermore, we generate simulated thermal emission and geometric albedo spectra and find large, wavelength-dependent differences between the CARMA and Virga models, with different gas absorption bands reacting differently to the different cloud distributions and particularly large differences in the M band. Therefore, constraining the vertically-dependent properties of water clouds will be essential to estimating the gas abundances in these atmospheres.

preprint2022arXiv

Patchy nightside clouds on ultra-hot Jupiters: General Circulation Model simulations with radiatively active cloud tracers

The atmospheres of ultra-hot Jupiters have been characterized in detail through recent phase curve and low- and high-resolution emission and transmission spectroscopic observations. Previous numerical studies have analyzed the effect of the localized recombination of hydrogen on the atmospheric dynamics and heat transport of ultra-hot Jupiters, finding that hydrogen dissociation and recombination lead to a reduction in the day-to-night contrasts of ultra-hot Jupiters relative to previous expectations. In this work, we add to previous efforts by also considering the localized condensation of clouds in the atmospheres of ultra-hot Jupiters, their resulting transport by the atmospheric circulation, and the radiative feedback of clouds on the atmospheric dynamics. To do so, we include radiatively active cloud tracers into the existing MITgcm framework for simulating the atmospheric dynamics of ultra-hot Jupiters. We take cloud condensate properties appropriate for the high-temperature condensate corundum from CARMA cloud microphysics models. We conduct a suite of GCM simulations with varying cloud microphysical and radiative properties, and we find that partial cloud coverage is a ubiquitous outcome of our simulations. This patchy cloud distribution is inherently set by atmospheric dynamics in addition to equilibrium cloud condensation, and causes a cloud greenhouse effect that warms the atmosphere below the cloud deck. Nightside clouds are further sequestered at depth due to a dynamically induced high-altitude thermal inversion. We post-process our GCMs with the Monte Carlo radiative transfer code gCMCRT and find that the patchy clouds on ultra-hot Jupiters do not significantly impact transmission spectra but can affect their phase-dependent emission spectra.

preprint2022arXiv

The First Near-Infrared Transmission Spectrum of HIP 41378 f, a Low-Mass Temperate Jovian World in a Multi-Planet System

We present a near-infrared transmission spectrum of the long period (P=542 days), temperate ($T_{eq}$=294 K) giant planet HIP 41378 f obtained with the Wide-Field Camera 3 (WFC3) instrument aboard the Hubble Space Telescope (HST). With a measured mass of 12 $\pm$ 3 $M_{\oplus}$ and a radius of 9.2 $\pm$ 0.1 $R_{\oplus}$, HIP 41378 f has an extremely low bulk density (0.09 $\pm$ 0.02 g/cm$^{3}$). We measure the transit depth with a median precision of 84 ppm in 30 spectrophotometric channels with uniformly-sized widths of 0.018 microns. Within this level of precision, the spectrum shows no evidence of absorption from gaseous molecular features between 1.1-1.7 microns. Comparing the observed transmission spectrum to a suite of 1D radiative-convective-thermochemical-equilibrium forward models, we rule out clear, low-metallicity atmospheres and find that the data prefer high-metallicity atmospheres or models with an additional opacity source such as high-altitude hazes and/or circumplanetary rings. We explore the ringed scenario for this planet further by jointly fitting the K2 and HST light curves to constrain the properties of putative rings. We also assess the possibility of distinguishing between hazy, ringed, and high-metallicity scenarios at longer wavelengths with JWST. HIP 41378 f provides a rare opportunity to probe the atmospheric composition of a cool giant planet spanning the gap between the Solar System giants, directly imaged planets, and the highly-irradiated hot Jupiters traditionally studied via transit spectroscopy.

preprint2022arXiv

The Hubble PanCET Program: A Featureless Transmission Spectrum for WASP-29b and Evidence of Enhanced Atmospheric Metallicity on WASP-80b

We present a uniform analysis of transit observations from the Hubble Space Telescope and Spitzer Space Telescope of two warm gas giants orbiting K-type stars - WASP-29b and WASP-80b. The transmission spectra, which span 0.4-5.0 $μ$m, are interpreted using a suite of chemical equilibrium PLATON atmospheric retrievals. Both planets show evidence of significant aerosol opacity along the day-night terminator. The spectrum of WASP-29b is flat throughout the visible and near-infrared, suggesting the presence of condensate clouds extending to low pressures. The lack of spectral features hinders our ability to constrain the atmospheric metallicity and C/O ratio. In contrast, WASP-80b shows a discernible, albeit muted H$_2$O absorption feature at 1.4 $μ$m, as well as a steep optical spectral slope that is caused by fine-particle aerosols and/or contamination from unocculted spots on the variable host star. WASP-80b joins the small number of gas-giant exoplanets that show evidence for enhanced atmospheric metallicity: the transmission spectrum is consistent with metallicities ranging from $\sim$30-100 times solar in the case of cloudy limbs to a few hundred times solar in the cloud-free scenario. In addition to the detection of water, we infer the presence of CO$_2$ in the atmosphere of WASP-80b based on the enhanced transit depth in the Spitzer 4.5 $μ$m bandpass. From a complementary analysis of Spitzer secondary eclipses, we find that the dayside emission from WASP-29b and WASP-80b is consistent with brightness temperatures of $937 \pm 48$ and $851 \pm 14$ K, respectively, indicating relatively weak day-night heat transport and low Bond albedo.

preprint2022arXiv

Transit Timing Variations for AU Microscopii b & c

We explore the transit timing variations (TTVs) of the young (22 Myr) nearby AU Mic planetary system. For AU Mic b, we introduce three Spitzer (4.5 $μ$m) transits, five TESS transits, 11 LCO transits, one PEST transit, one Brierfield transit, and two transit timing measurements from Rossiter-McLaughlin observations; for AU Mic c, we introduce three TESS transits. We present two independent TTV analyses. First, we use EXOFASTv2 to jointly model the Spitzer and ground-based transits and to obtain the midpoint transit times. We then construct an O--C diagram and model the TTVs with Exo-Striker. Second, we reproduce our results with an independent photodynamical analysis. We recover a TTV mass for AU Mic c of 10.8$^{+2.3}_{-2.2}$ M$_{E}$. We compare the TTV-derived constraints to a recent radial-velocity (RV) mass determination. We also observe excess TTVs that do not appear to be consistent with the dynamical interactions of b and c alone, and do not appear to be due to spots or flares. Thus, we present a hypothetical non-transiting "middle-d" candidate exoplanet that is consistent with the observed TTVs, the candidate RV signal, and would establish the AU Mic system as a compact resonant multi-planet chain in a 4:6:9 period commensurability. These results demonstrate that the AU Mic planetary system is dynamically interacting producing detectable TTVs, and the implied orbital dynamics may inform the formation mechanisms for this young system. We recommend future RV and TTV observations of AU Mic b and c to further constrain the masses and to confirm the existence of possible additional planet(s).

preprint2021arXiv

A New Sedimentation Model for Greater Cloud Diversity in Giant Exoplanets and Brown Dwarfs

The observed atmospheric spectrum of exoplanets and brown dwarfs depends critically on the presence and distribution of atmospheric condensates. The Ackerman & Marley (2001) methodology for predicting the vertical distribution of condensate particles is widely used to study cloudy atmospheres and has recently been implemented in an open-source python package virga. The model relies upon input parameter $f_{\text{sed}}$, the sedimentation efficiency, which until now has been held constant. The relative simplicity of this model renders it useful for retrieval studies due to its rapidly attainable solutions. However, comparisons with more complex microphysical models such as CARMA have highlighted inconsistencies between the two approaches, namely that the cloud parameters needed for radiative transfer produced by virga are dissimilar to those produced by CARMA. To address these discrepancies, we have extended the original Ackerman and Marley methodology in virga to allow for non-constant $f_{\text{sed}}$ values, in particular those that vary with altitude. We discuss one such parameterization and compare the cloud mass mixing ratio produced by virga with constant and variable $f_{\text{sed}}$ profiles to that produced by CARMA. We find that the variable $f_{\text{sed}}$ formulation better captures the profile produced by CARMA with heterogeneous nucleation, yet performs comparatively to constant $f_{\text{sed}}$ for homogeneous nucleation. In general, virga has the capacity to handle any $f_{\text{sed}}$ with an explicit anti-derivative, permitting a plethora of alternative cloud profiles that are otherwise unattainable by constant $f_{\text{sed}}$ values. The ensuing flexibility has the potential to better agree with increasingly complex models and observed data.

preprint2021arXiv

Spatially Resolved Modeling of Optical Albedos for a Sample of Six Hot Jupiters

Optical secondary eclipse measurements made by \emph{Kepler} reveal a diverse set of geometric albedos for hot Jupiters with equilibrium temperatures between $1550-1700$ K. The presence or absence of high altitude condensates, such as Mg$_2$SiO$_4$, Fe, Al$_2$O$_3$, and TiO$_2$, can significantly alter optical albedos, but these clouds are expected to be confined to localized regions in the atmospheres of these tidally locked planets. Here, we present 3D general circulation models and corresponding cloud and albedo maps for six hot Jupiters with measured optical albedos in this temperature range. We find that the observed optical albedos of K2-31b and K2-107b are best matched by either cloud free models or models with relatively compact cloud layers, while Kepler-8b and Kepler-17b's optical albedos can be matched by moderately extended ($f_{sed}$ = 0.1) parametric cloud models. HATS-11b has a high optical albedo, corresponding to models with bright Mg$_2$SiO$_4$ clouds extending to very low pressures ($f_{sed}$ = 0.03). We are unable to reproduce Kepler-7b's high albedo, as our models predict that the dayside will be dominated by dark Al$_2$O$_3$ clouds at most longitudes. We compare our parametric cloud model with a two-zone microphysical cloud model (\texttt{CARMA}). We find that even after accounting for the 3D thermal structure, no single cloud model can explain the full range of observed albedos within the sample. We conclude that a better knowledge of the vertical mixing profiles, cloud radiative feedback, cloud condensate properties, and atmospheric metallicities is needed in order to explain the unexpected diversity of albedos in this temperature range.

preprint2020arXiv

A planet within the debris disk around the pre-main-sequence star AU Microscopii

AU Microscopii (AU Mic) is the second closest pre main sequence star, at a distance of 9.79 parsecs and with an age of 22 million years. AU Mic possesses a relatively rare and spatially resolved3 edge-on debris disk extending from about 35 to 210 astronomical units from the star, and with clumps exhibiting non-Keplerian motion. Detection of newly formed planets around such a star is challenged by the presence of spots, plage, flares and other manifestations of magnetic activity on the star. Here we report observations of a planet transiting AU Mic. The transiting planet, AU Mic b, has an orbital period of 8.46 days, an orbital distance of 0.07 astronomical units, a radius of 0.4 Jupiter radii, and a mass of less than 0.18 Jupiter masses at 3 sigma confidence. Our observations of a planet co-existing with a debris disk offer the opportunity to test the predictions of current models of planet formation and evolution.

preprint2020arXiv

Deflating Super-Puffs: Impact of Photochemical Hazes on the Observed Mass-Radius Relationship of Low Mass Planets

The observed mass-radius relationship of low-mass planets informs our understanding of their composition and evolution. Recent discoveries of low mass, large radii objects ("super-puffs") have challenged theories of planet formation and atmospheric loss, as their high inferred gas masses make them vulnerable to runaway accretion and hydrodynamic escape. Here we propose that high altitude photochemical hazes could enhance the observed radii of low-mass planets and explain the nature of super-puffs. We construct model atmospheres in radiative-convective equilibrium and compute rates of atmospheric escape and haze distributions, taking into account haze coagulation, sedimentation, diffusion, and advection by an outflow wind. We develop mass-radius diagrams that include atmospheric lifetimes and haze opacity, which is enhanced by the outflow, such that young (~0.1-1 Gyr), warm (T$_{eq}$ $\geq$ 500 K), low mass objects ($M_c$ < 4M$_{\rm Earth}$) should experience the most apparent radius enhancement due to hazes, reaching factors of three. This reconciles the densities and ages of the most extreme super-puffs. For Kepler-51b, the inclusion of hazes reduces its inferred gas mass fraction to <10%, similar to that of planets on the large radius side of the sub-Neptune radius gap. This suggests that Kepler-51b may be evolving towards that population, and that some warm sub-Neptunes may have evolved from super-puffs. Hazes also render transmission spectra of super-puffs and sub-Neptunes featureless, consistent with recent measurements. Our hypothesis can be tested by future observations of super-puffs' transmission spectra at mid-infrared wavelengths, where we predict that the planet radius will be half of that observed in the near-infrared.

preprint2020arXiv

EarthFinder Probe Mission Concept Study: Characterizing nearby stellar exoplanet systems with Earth-mass analogs for future direct imaging

EarthFinder is a NASA Astrophysics Probe mission concept selected for study as input to the 2020 Astrophysics National Academies Decadal Survey. The EarthFinder concept is based on a dramatic shift in our understanding of how PRV measurements should be made. We propose a new paradigm which brings the high precision, high cadence domain of transit photometry as demonstrated by Kepler and TESS to the challenges of PRV measurements at the cm/s level. This new paradigm takes advantage of: 1) broad wavelength coverage from the UV to NIR which is only possible from space to minimize the effects of stellar activity; 2) extremely compact, highly stable, highly efficient spectrometers (R>150,000) which require the diffraction-limited imaging possible only from space over a broad wavelength range; 3) the revolution in laser-based wavelength standards to ensure cm/s precision over many years; 4) a high cadence observing program which minimizes sampling-induced period aliases; 5) exploiting the absolute flux stability from space for continuum normalization for unprecedented line-by-line analysis not possible from the ground; and 6) focusing on the bright stars which will be the targets of future imaging missions so that EarthFinder can use a ~1.5 m telescope.

preprint2020arXiv

Global Chemistry and Thermal Structure Models for the Hot Jupiter WASP-43b and Predictions for JWST

The James Webb Space Telescope (JWST) is expected to revolutionize the field of exoplanets. The broad wavelength coverage and the high sensitivity of its instruments will allow characterization of exoplanetary atmospheres with unprecedented precision. Following the Call for the Cycle 1 Early Release Science Program, the Transiting Exoplanet Community was awarded time to observe several targets, including WASP-43b. The atmosphere of this hot Jupiter has been intensively observed but still harbors some mysteries, especially concerning the day-night temperature gradient, the efficiency of the atmospheric circulation, and the presence of nightside clouds. We will constrain these properties by observing a full orbit of the planet and extracting its spectroscopic phase curve in the 5--12 $μ$m range with JWST/MIRI. To prepare for these observations, we performed an extensive modeling work with various codes: radiative transfer, chemical kinetics, cloud microphysics, global circulation models, JWST simulators, and spectral retrieval. Our JWST simulations show that we should achieve a precision of 210 ppm per 0.1 $μ$m spectral bin on average, which will allow us to measure the variations of the spectrum in longitude and measure the night-side emission spectrum for the first time. If the atmosphere of WASP-43b is clear, our observations will permit us to determine if its atmosphere has an equilibrium or disequilibrium chemical composition, providing eventually the first conclusive evidence of chemical quenching in a hot Jupiter atmosphere. If the atmosphere is cloudy, a careful retrieval analysis will allow us to identify the cloud composition.

preprint2020arXiv

Keck/NIRC2 $L$'-Band Imaging of Jovian-Mass Accreting Protoplanets around PDS 70

We present $L$'-band imaging of the PDS 70 planetary system with Keck/NIRC2 using the new infrared pyramid wavefront sensor. We detected both PDS 70 b and c in our images, as well as the front rim of the circumstellar disk. After subtracting off a model of the disk, we measured the astrometry and photometry of both planets. Placing priors based on the dynamics of the system, we estimated PDS 70 b to have a semi-major axis of $20^{+3}_{-4}$~au and PDS 70 c to have a semi-major axis of $34^{+12}_{-6}$~au (95\% credible interval). We fit the spectral energy distribution (SED) of both planets. For PDS 70 b, we were able to place better constraints on the red half of its SED than previous studies and inferred the radius of the photosphere to be 2-3~$R_{Jup}$. The SED of PDS 70 c is less well constrained, with a range of total luminosities spanning an order of magnitude. With our inferred radii and luminosities, we used evolutionary models of accreting protoplanets to derive a mass of PDS 70 b between 2 and 4 $M_{\textrm{Jup}}$ and a mean mass accretion rate between $3 \times 10^{-7}$ and $8 \times 10^{-7}~M_{\textrm{Jup}}/\textrm{yr}$. For PDS 70 c, we computed a mass between 1 and 3 $M_{\textrm{Jup}}$ and mean mass accretion rate between $1 \times 10^{-7}$ and $5 \times~10^{-7} M_{\textrm{Jup}}/\textrm{yr}$. The mass accretion rates imply dust accretion timescales short enough to hide strong molecular absorption features in both planets' SEDs.

preprint2020arXiv

The Intrinsic Temperature and Radiative-Convective Boundary Depth in the Atmospheres of Hot Jupiters

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.

preprint2019arXiv

Transit Signatures of Inhomogeneous Clouds on Hot Jupiters: Insights From Microphysical Cloud Modeling

We determine the observability in transmission of inhomogeneous cloud cover on the limbs of hot Jupiters through post processing a general circulation model to include cloud distributions computed using a cloud microphysics model. We find that both the east and west limb often form clouds, but that the different properties of these clouds enhances the limb to limb differences compared to the clear case. Using JWST it should be possible to detect the presence of cloud inhomogeneities by comparing the shape of the transit lightcurve at multiple wavelengths because inhomogeneous clouds impart a characteristic, wavelength dependent signature. This method is statistically robust even with limited wavelength coverage, uncertainty on limb darkening coefficients, and imprecise transit times. We predict that the short wavelength slope varies strongly with temperature. The hot limb of the hottest planets form higher altitude clouds composed of smaller particles leading to a strong rayleigh slope. The near infrared spectral features of clouds are almost always detectable, even when no spectral slope is visible in the optical. In some of our models a spectral window between 5 and 9 microns can be used to probe through the clouds and detect chemical spectral features. Our cloud particle size distributions are not log-normal and differ from species to species. Using the area or mass weighted particle size significantly alters the relative strength of the cloud spectral features compared to using the predicted size distribution. Finally, the cloud content of a given planet is sensitive to a species' desorption energy and contact angle, two parameters that could be constrained experimentally in the future.

preprint2016arXiv

A High-Precision NIR Survey for RV Variable Low-Mass Stars

We present the results of a precise near-infrared (NIR) radial velocity (RV) survey of 32 low-mass stars with spectral types K2-M4 using CSHELL at the NASA IRTF in the $K$-band with an isotopologue methane gas cell to achieve wavelength calibration and a novel iterative RV extraction method. We surveyed 14 members of young ($\approx$ 25-150 Myr) moving groups, the young field star $\varepsilon$ Eridani as well as 18 nearby ($<$ 25 pc) low-mass stars and achieved typical single-measurement precisions of 8-15 m s$^{-1}$ with a long-term stability of 15-50 m s$^{-1}$. We obtain the best NIR RV constraints to date on 27 targets in our sample, 19 of which were never followed by high-precision RV surveys. Our results indicate that very active stars can display long-term RV variations as low as $\sim$ 25-50 m s$^{-1}$ at $\approx$ 2.3125 $μ$m, thus constraining the effect of jitter at these wavelengths. We provide the first multi-wavelength confirmation of GJ 876 bc and independently retrieve orbital parameters consistent with previous studies. We recovered RV variability for HD 160934 AB and GJ 725 AB that are consistent with their known binary orbits, and nine other targets are candidate RV variables with a statistical significance of 3-5$σ$. Our method combined with the new iSHELL spectrograph will yield long-term RV precisions of $\lesssim$ 5 m s$^{-1}$ in the NIR, which will allow the detection of Super-Earths near the habitable zone of mid-M dwarfs.

preprint2016arXiv

Aggregate Particles in the Plumes of Enceladus

Estimates of the total particulate mass of the plumes of Enceladus are important to constrain theories of particle formation and transport at the surface and interior of the satellite. We revisit the calculations of Ingersoll & Ewald (2011), who estimated the particulate mass of the Enceladus plumes from strongly forward scattered light in Cassini ISS images. We model the plume as a combination of spherical particles and irregular aggregates resulting from the coagulation of spherical monomers, the latter of which allows for plumes of lower particulate mass. Though a continuum of solutions are permitted by the model, the best fits to the ISS data consist either of low mass plumes composed entirely of small aggregates or high mass plumes composed of mostly spheres. The high particulate mass plumes have total particulate masses of (166 $\pm$ 42) $\times$ 10$^3$ kg, consistent with the results of Ingersoll & Ewald (2011). The low particulate mass plumes have masses of (25 $\pm$ 4) $\times$ 10$^3$ kg, leading to a solid to vapor mass ratio of 0.07 $\pm$ 0.01 for the plume. If indeed the plumes are made of such aggregates, then a vapor-based origin for the plume particles cannot be ruled out. Finally, we show that the residence time of the monomers inside the plume vents is sufficiently long for Brownian coagulation to form the aggregates before they are ejected to space.

preprint2016arXiv

Constraints on the Microphysics of Pluto's Photochemical Haze from New Horizons Observations

The New Horizons flyby of Pluto confirmed the existence of hazes in its atmosphere. Observations of a large high- to low- phase brightness ratio, combined with the blue color of the haze, suggest that the haze particles are fractal aggregates, analogous to the photochemical hazes on Titan. Therefore, studying the Pluto hazes can shed light on the similarities and differences between the Pluto and Titan atmospheres. We model the haze distribution using the Community Aerosol and Radiation Model for Atmospheres assuming that the distribution is shaped by sedimentation and coagulation of particles originating from photochemistry. Hazes composed of both purely spherical and purely fractal aggregate particles are considered. Agreement between model results and occultation observations is obtained with aggregate particles when the downward flux of photochemical products is equal to the column-integrated methane destruction rate ~1.2 $\times$ 10$^{-14}$ g cm$^{-2}$ s$^{-1}$, while for spherical particles the mass flux must be 2-3 times greater. This flux is nearly identical to the haze production flux of Titan previously obtained by comparing microphysical model results to Cassini observations. The aggregate particle radius is sensitive to particle charging, and a particle charge to radius ratio of 30 e-/μm is necessary to produce ~0.1-0.2 μm aggregates near Pluto's surface, in accordance with forward scattering measurements. Such a particle charge to radius ratio is 2-4 times higher than those previously obtained for Titan. Hazes composed of spheres with the same particle charge to radius ratio have particles that are 4 times smaller. These results further suggest that the haze particles are fractal aggregates. We also consider the effect of condensation of HCN, and C$_{2}$-hydrocarbons on the haze particles, which may play an important role in shaping their distributions.

preprint2016arXiv

Hypotheses for near-surface exchange of methane on Mars

The Curiosity rover recently detected a background of 0.7 ppb and spikes of 7 ppb of methane on Mars. This in situ measurement reorients our understanding of the Martian environment and its potential for life, as the current theories do not entail any geological source or sink of methane that varies sub-annually. In particular, the 10-fold elevation during the southern winter indicates episodic sources of methane that are yet to be discovered. Here we suggest a near-surface reservoir could explain this variability. Using the temperature and humidity measurements from the rover, we find that perchlorate salts in the regolith deliquesce to form liquid solutions, and deliquescence progresses to deeper subsurface in the season of the methane spikes. We therefore formulate the following three testable hypotheses. The first scenario is that the regolith in Gale Crater adsorbs methane when dry and releases this methane to the atmosphere upon deliquescence. The adsorption energy needs to be 36 kJ/mol to explain the magnitude of the methane spikes, higher than existing laboratory measurements. The second scenario is that microorganisms convert organic matter in the soil to methane when they are in liquid solutions. This scenario does not require regolith adsorption, but entails extant life on Mars. The third scenario is that deep subsurface aquifers produce the bursts of methane. Continued in situ measurements of methane and water, as well as laboratory studies of adsorption and deliquescence, will test these hypotheses and inform the existence of the near-surface reservoir and its exchange with the atmosphere.

preprint2016arXiv

Precise Near-Infrared Radial Velocities

We present the results of two 2.3 micron near-infrared radial velocity surveys to detect exoplanets around 36 nearby and young M dwarfs. We use the CSHELL spectrograph (R ~46,000) at the NASA InfraRed Telescope Facility, combined with an isotopic methane absorption gas cell for common optical path relative wavelength calibration. We have developed a sophisticated RV forward modeling code that accounts for fringing and other instrumental artifacts present in the spectra. With a spectral grasp of only 5 nm, we are able to reach long-term radial velocity dispersions of ~20-30 m/s on our survey targets.

preprint2016arXiv

Retrieval of Precise Radial Velocities from Near-Infrared High Resolution Spectra of Low Mass Stars

Given that low-mass stars have intrinsically low luminosities at optical wavelengths and a propensity for stellar activity, it is advantageous for radial velocity (RV) surveys of these objects to use near-infrared (NIR) wavelengths. In this work we describe and test a novel RV extraction pipeline dedicated to retrieving RVs from low mass stars using NIR spectra taken by the CSHELL spectrograph at the NASA Infrared Telescope Facility, where a methane isotopologue gas cell is used for wavelength calibration. The pipeline minimizes the residuals between the observations and a spectral model composed of templates for the target star, the gas cell, and atmospheric telluric absorption; models of the line spread function, continuum curvature, and sinusoidal fringing; and a parameterization of the wavelength solution. The stellar template is derived iteratively from the science observations themselves without a need for separate observations dedicated to retrieving it. Despite limitations from CSHELL's narrow wavelength range and instrumental systematics, we are able to (1) obtain an RV precision of 35 m/s for the RV standard star GJ 15 A over a time baseline of 817 days, reaching the photon noise limit for our attained SNR, (2) achieve ~3 m/s RV precision for the M giant SV Peg over a baseline of several days and confirm its long-term RV trend due to stellar pulsations, as well as obtain nightly noise floors of ~2 - 6 m/s, and (3) show that our data are consistent with the known masses, periods, and orbital eccentricities of the two most massive planets orbiting GJ 876. Future applications of our pipeline to RV surveys using the next generation of NIR spectrographs, such as iSHELL, will enable the potential detection of Super-Earths and Mini-Neptunes in the habitable zones of M dwarfs.

preprint2015arXiv

DeepLogo: Hitting Logo Recognition with the Deep Neural Network Hammer

Recently, there has been a flurry of industrial activity around logo recognition, such as Ditto's service for marketers to track their brands in user-generated images, and LogoGrab's mobile app platform for logo recognition. However, relatively little academic or open-source logo recognition progress has been made in the last four years. Meanwhile, deep convolutional neural networks (DCNNs) have revolutionized a broad range of object recognition applications. In this work, we apply DCNNs to logo recognition. We propose several DCNN architectures, with which we surpass published state-of-art accuracy on a popular logo recognition dataset.

preprint2015arXiv

Stability of CO2 Atmospheres on Desiccated M Dwarf Exoplanets

We investigate the chemical stability of CO2-dominated atmospheres of desiccated M dwarf terrestrial exoplanets using a 1-dimensional photochemical model. Around Sun-like stars, CO2 photolysis by Far-UV (FUV) radiation is balanced by recombination reactions that depend on water abundance. Planets orbiting M dwarf stars experience more FUV radiation, and could be depleted in water due to M dwarfs' prolonged, high-luminosity pre-main sequences (Luger & Barnes 2015). We show that, for water-depleted M dwarf terrestrial planets, a catalytic cycle relying on H2O2 photolysis can maintain a CO2 atmosphere. However, this cycle breaks down for atmospheric hydrogen mixing ratios <1 ppm, resulting in ~40% of the atmospheric CO2 being converted to CO and O2 on a time scale of 1 Myr. The increased O2 abundance leads to high O3 concentrations, the photolysis of which forms another CO2-regenerating catalytic cycle. For atmospheres with <0.1 ppm hydrogen, CO2 is produced directly from the recombination of CO and O. These catalytic cycles place an upper limit of ~50% on the amount of CO2 that can be destroyed via photolysis, which is enough to generate Earth-like abundances of (abiotic) O2 and O3. The conditions that lead to such high oxygen levels could be widespread on planets in the habitable zones of M dwarfs. Discrimination between biological and abiotic O2 and O3 in this case can perhaps be accomplished by noting the lack of water features in the reflectance and emission spectra of these planets, which necessitates observations at wavelengths longer than 0.95 microns.

preprint2013arXiv

Bimodal Distribution of Sulfuric Acid Aerosols in the Upper Haze of Venus

The upper haze (UH) of Venus is variable on the order of days and it is populated by two particle modes. We use a 1D microphysics and vertical transport model based on the Community Aerosol and Radiation Model for Atmospheres to evaluate whether interaction of upwelled cloud particles and sulfuric acid particles nucleated in situ on meteoric dust are able to generate the two size modes and whether their observed variability are due to cloud top vertical transient winds. Nucleation of photochemically produced sulfuric acid onto polysulfur condensation nuclei generates mode 1 cloud droplets that then diffuse upwards into the UH. Droplets generated in the UH from nucleation of sulfuric acid onto meteoric dust coagulate with the upwelled cloud particles and cannot reproduce the observed bimodal size distribution. The mass transport enabled by cloud top transient winds are able to generate a bimodal size distribution in a time scale consistent with observations. Sedimentation and convection in the middle and lower clouds causes the formation of large mode 2 and mode 3 particles. Evaporation of these particles below the clouds creates a local sulfuric acid vapor maximum that causes upwelling of sulfuric acid back into the clouds. If the polysulfur condensation nuclei are small and their production rate is high, coagulation of small droplets onto larger droplets in the middle cloud may result in sulfuric acid "rain" below the clouds once every few Earth months. Reduction of the polysulfur condensation nuclei production rate destroys this oscillation and reduces the mode 1 particle abundance in the middle cloud by two orders of magnitude, though it better reproduces the sulfur-to-sulfuric-acid mass ratio in the cloud and haze droplets. In general we find satisfactory agreement between our results and observations, though improvements could be made by incorporating sulfur microphysics.

preprint2013arXiv

Nonhydrostatic Effects and the Determination of Icy Satellites' Moment of Inertia

We compare the moment of inertia (MOI) of a simple hydrostatic, two layer body as determined by the Radau-Darwin Approximation (RDA) to its exact hydrostatic MOI calculated to first order in the parameter q = w^2R^3/GM, where w, R, and M are the spin angular velocity, radius, and mass of the body, and G is the gravitational constant. RDA is in error by less than 1% for many configurations of core sizes and layer densities congruent with those of solid bodies in the Solar System. We determine the error in the MOI of icy satellites calculated with the RDA due to nonhydrostatic effects by using a simple model in which the core and outer shell have slight degree 2 distortions away from their expected hydrostatic shapes. Since the hydrostatic shape has an associated stress of order pw^2R^2 (where p is density) it follows that the importance of nonhydrostatic effects scales with the dimensionless number s/pw^2R^2, where s is the nonhydrostatic stress. This highlights the likely importance of this error for slowly rotating bodies (e.g., Titan and Callisto) and small bodies (e.g., Saturn moons other than Titan). We apply this model to Titan, Callisto, and Enceladus and find that the RDA-derived MOI can be 10% greater than the actual MOI for nonhydrostatic stresses as small as ~0.1 bars at the surface or ~1 bar at the core-mantle boundary, while for Ganymede the stresses necessary to produce the same MOI errors are an order of magnitude greater due to its faster rotation. If satellites can reorient to the lowest energy state then RDA will always give an overestimate of the true MOI. Observations have shown that small nonhydrostatic gravity anomalies exist on Ganymede and Titan. We conclude that nonhydrostatic effects could be present to an extent that allows Callisto and Titan to be fully differentiated.

preprint2013arXiv

Precision near-infrared radial velocity instrumentation I: Absorption Gas Cells

We have built and commissioned gas absorption cells for precision spectroscopic radial velocity measurements in the near-infrared in the H and K bands. We describe the construction and installation of three such cells filled with 13CH4, 12CH3D, and 14NH3 for the CSHELL spectrograph at the NASA Infrared Telescope Facility (IRTF). We have obtained their high-resolution laboratory Fourier Transform spectra, which can have other practical uses. We summarize the practical details involved in the construction of the three cells, and the thermal and mechanical control. In all cases, the construction of the cells is very affordable. We are carrying out a pilot survey with the 13CH4 methane gas cell on the CSHELL spectrograph at the IRTF to detect exoplanets around low mass and young stars. We discuss the current status of our survey, with the aim of photon-noise limited radial velocity precision. For adequately bright targets, we are able to probe a noise floor of 7 m/s with the gas cell with CSHELL at cassegrain focus. Our results demonstrate the feasibility of using a gas cell on the next generation of near-infrared spectrographs such as iSHELL on IRTF, iGRINS, and an upgraded NIRSPEC at Keck.

preprint2013arXiv

Precision near-infrared radial velocity instrumentation II: Non-Circular Core Fiber Scrambler

We have built and commissioned a prototype agitated non-circular core fiber scrambler for precision spectroscopic radial velocity measurements in the near-infrared H band. We have collected the first on-sky performance and modal noise tests of these novel fibers in the near-infrared at H and K bands using the CSHELL spectrograph at the NASA InfraRed Telescope Facility (IRTF). We discuss the design behind our novel reverse injection of a red laser for co-alignment of star-light with the fiber tip via a corner cube and visible camera. We summarize the practical details involved in the construction of the fiber scrambler, and the mechanical agitation of the fiber at the telescope. We present radial velocity measurements of a bright standard star taken with and without the fiber scrambler to quantify the relative improvement in the obtainable blaze function stability, the line spread function stability, and the resulting radial velocity precision. We assess the feasibility of applying this illumination stabilization technique to the next generation of near-infrared spectrographs such as iSHELL on IRTF and an upgraded NIRSPEC at Keck. Our results may also be applied in the visible for smaller core diameter fibers where fiber modal noise is a significant factor, such as behind an adaptive optics system or on a small < 1 meter class telescope such as is being pursued by the MINERVA and LCOGT collaborations.

preprint2012arXiv

Design and Construction of Absorption Cells for Precision Radial Velocities in the K Band using Methane Isotopologues

We present a method to optimize absorption cells for precise wavelength calibration in the near-infrared. We apply it to design and optimize methane isotopologue cells for precision radial velocity measurements in the K band. We also describe the construction and installation of two such cells for the CSHELL spectrograph at NASA's IRTF. We have obtained their high-resolution laboratory spectra, which we can then use in precision radial velocity measurements and which can also have other applications. In terms of obtainable RV precision methane should out-perform other proposed cells, such as the ammonia cell ($^{14}$NH$_{3}$) recently demonstrated on CRIRES/VLT. The laboratory spectra of Ammonia and the Methane cells show strong absorption features in the H band that could also be exploited for precision Doppler measurements. We present spectra and preliminary radial velocity measurements obtained during our first-light run. These initial results show that a precision down to 20-30 m s$^{-1}$ can be obtained using a wavelength interval of only 5 nm in the K band and S/N$\sim$150. This supports the prediction that a precision down to a few m s$^{-1}$ can be achieved on late M dwarfs using the new generation of NIR spectrographs, thus enabling the detection of terrestrial planets in their habitable zones. Doppler measurements in the NIR can also be used to mitigate the radial velocity jitter due to stellar activity enabling more efficient surveys on young active stars.