Source author record

Gregory Laughlin

Gregory Laughlin 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

58works
8topics
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

58 published item(s)

preprint2022arXiv

Identifying Interstellar Object Impact Craters

The discoveries of two Interstellar Objects (ISOs) in recent years has generated significant interest in constraining their physical properties and the mechanisms behind their formation. However, their ephemeral passages through our Solar System permitted only incomplete characterization. We investigate avenues for identifying craters that may have been produced by ISOs impacting terrestrial Solar System bodies, with particular attention towards the Moon. A distinctive feature of ISOs is their relatively high encounter velocity compared to asteroids and comets. Local stellar kinematics indicate that terrestrial Solar System bodies should have experienced of order unity ISO impacts exceeding 100 km/s. By running hydrodynamical simulations for projectiles of different masses and impact velocities, up to 100 km/s, we show how late-stage equivalence dictates that transient crater dimensions are alone insufficient for inferring the projectile's velocity. On the other hand, the melt volume within craters of a fixed diameter may be a potential route for identifying ISO craters, as faster impacts produce more melt. This method requires that the melt volume scales with the energy of the projectile, while crater diameter scales with the point-source limit (sub-energy). Given that there are probably only a few ISO craters in the Solar System at best, and that transient crater dimensions are not a distinguishing feature for impact velocities at least up to 100 km/s, identification of an ISO crater proves a challenging task. Melt volume and high-pressure petrology may be diagnostic features once large volumes of material can be analyzed in situ.

preprint2022arXiv

New Constraints on Gliese 876 -- Exemplar of Mean-Motion Resonance

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.

preprint2022arXiv

Origins of Hot Jupiters from the Stellar Obliquity Distribution

The obliquity of a star, or the angle between its spin axis and the average orbit normal of its companion planets, provides a unique constraint on that system's evolutionary history. Unlike the Solar System, where the Sun's equator is nearly aligned with its companion planets, many hot Jupiter systems have been discovered with large spin-orbit misalignments, hosting planets on polar or retrograde orbits. We demonstrate that, in contrast to stars harboring hot Jupiters on circular orbits, those with eccentric companions follow no population-wide obliquity trend with stellar temperature. This finding can be naturally explained through a combination of high-eccentricity migration and tidal damping. Furthermore, we show that the joint obliquity and eccentricity distributions observed today are consistent with the outcomes of high-eccentricity migration, with no strict requirement to invoke the other hot Jupiter formation mechanisms of disk migration or in-situ formation. At a population-wide level, high-eccentricity migration can consistently shape the dynamical evolution of hot Jupiter systems.

preprint2022arXiv

Revisiting the Full Sets of Orbital Parameters for the XO-3 System: No evidence for Temporal Variation of the Spin-Orbit Angle

We present 12 new transit light curves and 16 new out-of-transit radial velocity measurements for the XO-3 system. By modelling our newly collected measurements together with archival photometric and Doppler velocimetric data, we confirmed the unusual configuration of the XO-3 system, which contains a massive planet ($M_P=11.92^{+0.59}_{-0.63} M_J$) on a relatively eccentric ($e=0.2853^{+0.0027}_{-0.0026}$) and short-period ($3.19152 \pm 0.00145\,$day) orbit around a massive star ($M_*=1.219^{+0.090}_{-0.095} M_{\odot}$). Furthermore, we find no strong evidence for a temporal change of either $V\sin i_{*}$ (and by extension, the stellar spin vector of XO-3), or the transit profile (and thus orbital angular momentum vector of XO-3b). We conclude that the discrepancy in previous Rossiter-McLaughlin measurements ($70.0^{\circ} \pm 15.0^{\circ}$ (Hebrard et al. 2008); $37.3^{\circ} \pm 3.7^{\circ}$ (Winn et al. 2009); $37.3^{\circ} \pm 3.0^{\circ}$ (Hirano et al. 2011)) may have stemmed from systematic noise sources.

preprint2022arXiv

Stacked Periodograms as a Probe of Exoplanetary Populations

Ongoing, extreme-precision Doppler radial velocity surveys seek planets with masses less than several M$_\oplus$; population-level studies to determine the distribution of planetary masses, however, remain difficult due to the required observational time investment, as well as challenges associated with robustly detecting the lowest mass planets. We outline a novel approach that leverages extensive, existing RV datasets to constrain masses of exoplanet populations: stacking periodograms of RV timeseries across many targets. We show that an exoplanet population may be statistically identifiable in the stacked periodogram, even when individual planets do not pass the threshold of detection. We discuss analytical, statistical properties of the stacked periodogram, perform simulations to demonstrate the efficacy of the method, and investigate the influence of semi-structured window functions and stellar activity. Analysis of the Lick-Carnegie Exoplanet Survey data set reveals a marginally significant ($1.6σ$) signal consistent with a population of exoplanets occupying $3-7$ day periods with typical $K$ between $1.6-5.1$ m s$^{-1}$. More detailed investigation of signals associated with stellar activity and yearly systematics may be necessary to confirm this result or detect other underlying Keplerian contributions.

preprint2022arXiv

The Aligned Orbit of WASP-148b, the Only Known Hot Jupiter with a Nearby Warm Jupiter Companion, from NEID and HIRES

We present spectroscopic measurements of the Rossiter-McLaughlin effect for WASP-148b, the only known hot Jupiter with a nearby warm-Jupiter companion, from the WIYN/NEID and Keck/HIRES instruments. This is one of the first scientific results reported from the newly commissioned NEID spectrograph, as well as the second obliquity constraint for a hot Jupiter system with a close-in companion, after WASP-47. WASP-148b is consistent with being in alignment with the sky-projected spin axis of the host star, with $λ=-8^{\circ}.2^{{+8^{\circ}.7}}_{-9^{\circ}.7}$. The low obliquity observed in the WASP-148 system is consistent with the orderly-alignment configuration of most compact multi-planet systems around cool stars with obliquity constraints, including our solar system, and may point to an early history for these well-organized systems in which migration and accretion occurred in isolation, with relatively little disturbance. By contrast, previous results have indicated that high-mass and hot stars appear to more commonly host a wide range of misaligned planets: not only single hot Jupiters, but also compact systems with multiple super-Earths. We suggest that, to account for the high rate of spin-orbit misalignments in both compact multi-planet and isolated-hot-Jupiter systems orbiting high-mass and hot stars, spin-orbit misalignments may be caused by distant giant planet perturbers, which are most common around these stellar types.

preprint2022arXiv

The Prospects for Hurricane-like Vortices in Protoplanetary Disks

When ice on the surface of dust grains in protoplanetary disk sublimates, it adds its latent heat of water sublimation to the surrounding flow. Drawing on the analogy provided by tropical cyclones on Earth, we investigate if this energy source is sufficient to sustain or magnify anticyclonic disk vortices that would otherwise fall victim to viscous dissipation. An analytical treatment, supported by exploratory two-dimensional simulations, suggests that even modestly under-saturated flows can extend the lifetime of vortices, potentially to a degree sufficient to aid particle trapping and planetesimal formation. We expect the best conditions for this mechanism to occur will be found near the disk's water ice line if turbulent motions displace gas parcels out of thermodynamic equilibrium with the dust mid-plane.

preprint2020arXiv

Energy Optimization in Extrasolar Planetary Systems: The Transition from Peas-in-a-Pod to Runaway Growth

Motivated by the trends found in the observed sample of extrasolar planets, this paper determines tidal equilibrium states for forming planetary systems --- subject to conservation of angular momentum, constant total mass, and fixed orbital spacing. In the low-mass limit, valid for superearth-class planets with masses of order $m_{\rm p}\sim10M_\oplus$, previous work showed that energy optimization leads to nearly equal mass planets, with circular orbits confined to a plane. The present treatment generalizes previous results by including the self-gravity of the planetary bodies. For systems with sufficiently large total mass $m_{\scriptstyle T}$ in planets, the optimized energy state switches over from the case of nearly equal mass planets to a configuration where one planet contains most of the material. This transition occurs for a critical mass threshold of approximately $m_{\scriptstyle T} \ge m_{\scriptstyle C}\sim40M_\oplus$ (where the value depends on the semimajor axes of the planetary orbits, the stellar mass, and other system properties). These considerations of energy optimization apply over a wide range of mass scales, from binary stars to planetary systems to the collection of moons orbiting the giant planets in our solar system.

preprint2020arXiv

Evidence that 1I/2017 U1 (`Oumuamua) was composed of molecular hydrogen ice

`Oumuamua (I1 2017) was the first macroscopic ($l\sim100\,{\rm m}$) body observed to traverse the inner solar system on an unbound hyperbolic orbit. Its light curve displayed strong periodic variation, and it showed no hint of a coma or emission from molecular outgassing. Astrometric measurements indicate that 'Oumuamua experienced non-gravitational acceleration on its outbound trajectory, but energy balance arguments indicate this acceleration is inconsistent with a water ice sublimation-driven jet of the type exhibited by solar system comets. We show that all of `Oumaumua's observed properties can be explained if it contained a significant fraction of molecular hydrogen (H$_{2}$) ice. H$_{2}$ sublimation at a rate proportional to the incident solar flux generates a surface-covering jet that reproduces the observed acceleration. Mass wasting from sublimation leads to monotonic increase in the body axis ratio, explaining `Oumuamua's shape. Back-tracing `Oumuamua's trajectory through the Solar System permits calculation of its mass and aspect ratio prior to encountering the Sun. We show that H$_{2}$-rich bodies plausibly form in the coldest dense cores of Giant Molecular Clouds, where number densities are of order $n\sim10^5$, and temperatures approach the $T=3\,{\rm K}$ background. Post-formation exposure to galactic cosmic rays implies a $τ\sim 100$ Myr age, explaining the kinematics of `Oumuamua's inbound trajectory.

preprint2020arXiv

EXPRES I. HD~3651 an Ideal RV Benchmark

The next generation of exoplanet-hunting spectrographs should deliver up to an order of magnitude improvement in radial velocity precision over the standard 1 m/s state of the art. This advance is critical for enabling the detection of Earth-mass planets around Sun-like stars. New calibration techniques such as laser frequency combs and stabilized etalons ensure that the instrumental stability is well characterized. However, additional sources of error include stellar noise, undetected short-period planets, and telluric contamination. To understand and ultimately mitigate error sources, the contributing terms in the error budget must be isolated to the greatest extent possible. Here, we introduce a new high cadence radial velocity program, the EXPRES 100 Earths program, which aims to identify rocky planets around bright, nearby G and K dwarfs. We also present a benchmark case: the 62-d orbit of a Saturn-mass planet orbiting the chromospherically quiet star, HD 3651. The combination of high eccentricity (0.6) and a moderately long orbital period, ensures significant dynamical clearing of any inner planets. Our Keplerian model for this planetary orbit has a residual RMS of 58 cm/s over a $\sim 6$ month time baseline. By eliminating significant contributors to the radial velocity error budget, HD 3651 serves as a standard for evaluating the long term precision of extreme precision radial velocity (EPRV) programs.

preprint2016arXiv

3.6 and 4.5 $μ$m ${\it Spitzer}$ Phase Curves of the Highly-Irradiated Hot Jupiters WASP-19b and HAT-P-7b

We analyze full-orbit phase curve observations of the transiting hot Jupiters WASP-19b and HAT-P-7b at 3.6 and 4.5 $μ$m obtained using the Spitzer Space Telescope. For WASP-19b, we measure secondary eclipse depths of $0.485\%\pm 0.024\%$ and $0.584\%\pm 0.029\%$ at 3.6 and 4.5 $μ$m, which are consistent with a single blackbody with effective temperature $2372 \pm 60$ K. The measured 3.6 and 4.5 $μ$m secondary eclipse depths for HAT-P-7b are $0.156\%\pm 0.009\%$ and $0.190\%\pm 0.006\%$, which are well-described by a single blackbody with effective temperature $2667\pm 57$ K. Comparing the phase curves to the predictions of one-dimensional and three-dimensional atmospheric models, we find that WASP-19b's dayside emission is consistent with a model atmosphere with no dayside thermal inversion and moderately efficient day-night circulation. We also detect an eastward-shifted hotspot, suggesting the presence of a superrotating equatorial jet. In contrast, HAT-P-7b's dayside emission suggests a dayside thermal inversion and relatively inefficient day-night circulation; no hotspot shift is detected. For both planets, these same models do not agree with the measured nightside emission. The discrepancies in the model-data comparisons for WASP-19b might be explained by high-altitude silicate clouds on the nightside and/or high atmospheric metallicity, while the very low 3.6 $μ$m nightside planetary brightness for HAT-P-7b may be indicative of an enhanced global C/O ratio. We compute Bond albedos of 0 ($<0.08$ at $1σ$) and $0.38\pm 0.06$ for WASP-19b and HAT-P-7b, respectively. In the context of other planets with thermal phase curve measurements, we show that WASP-19b and HAT-P-7b fit the general trend of decreasing day-night heat recirculation with increasing irradiation.

preprint2016arXiv

Direct Measure of Radiative and Dynamical Properties of an Exoplanet Atmosphere

Two decades after the discovery of 51 Peg b, the formation processes and atmospheres of short-period gas giants remain poorly understood. Observations of eccentric systems provide key insights on those topics as they can illuminate how a planet's atmosphere responds to changes in incident flux. We report here the analysis of multi-day multi-channel photometry of the eccentric (e~ 0.93) hot Jupiter HD 80606 b obtained with the Spitzer Space Telescope. The planet's extreme eccentricity combined with the long coverage and exquisite precision of new periastron-passage observations allow us to break the degeneracy between the radiative and dynamical timescales of HD 80606 b's atmosphere and constrain its global thermal response. Our analysis reveals that the atmospheric layers probed heat rapidly (~4 hr radiative timescale) from $\lt$500 to 1400 K as they absorb ~ 20% of the incoming stellar flux during the periastron passage, while the planet's rotation period is 93$\pm_{35}^{85}$ hr, which exceeds the predicted pseudo-synchronous period (40 hr).

preprint2016arXiv

On the Detection of Non-Transiting Hot Jupiters in Multiple-Planet Systems

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$.

preprint2016arXiv

The Hunt for Planet Nine: Atmosphere, Spectra, Evolution, and Detectability

We investigate the physical characteristics of the Solar System's proposed Planet Nine using modeling tools with a heritage in studying Uranus and Neptune. For a range of plausible masses and interior structures, we find upper limits on the intrinsic Teff, from ~35-50 K for masses of 5-20 M_Earth, and we also explore lower Teff values. Possible planetary radii could readily span from 3 to 6 R_Earth depending on the mass fraction of any H/He envelope. Given its cold temperature, the planet encounters significant methane condensation, which dramatically alters the atmosphere away from simple Neptune-like expectations. We find the atmosphere is strongly depleted in molecular absorption at visible wavelengths, suggesting a Rayleigh scattering atmosphere with a high geometric albedo approaching 0.75. We highlight two diagnostics for the atmosphere's temperature structure, the first being the value of the methane mixing ratio above the methane cloud. The second is the wavelength at which cloud scattering can be seen, which yields the cloud-top pressure. Surface reflection may be seen if the atmosphere is thin. Due to collision-induced opacity of H2 in the infrared, the planet would be extremely blue (instead of red) in the shortest wavelength WISE colors if methane is depleted, and would, in some cases, exist on the verge of detectability by WISE. For a range of models, thermal fluxes from ~3-5 microns are ~20 orders of magnitude larger than blackbody expectations. We report a search of the AllWISE Source Catalog for Planet Nine, but find no detection.

preprint2015arXiv

3.6 and 4.5 $μ$m Phase Curves of the Highly-Irradiated Eccentric Hot Jupiter WASP-14b

We present full-orbit phase curve observations of the eccentric ($e\sim 0.08$) transiting hot Jupiter WASP-14b obtained in the 3.6 and 4.5 $μ$m bands using the \textit{Spitzer Space Telescope}. We use two different methods for removing the intrapixel sensitivity effect and compare their efficacy in decoupling the instrumental noise. Our measured secondary eclipse depths of $0.1882\%\pm 0.0048\%$ and $0.2247\%\pm 0.0086\%$ at 3.6 and 4.5 $μ$m, respectively, are both consistent with a blackbody temperature of $2402\pm 35$ K. We place a $2σ$ upper limit on the nightside flux at 3.6 $μ$m and find it to be $9\%\pm 1\%$ of the dayside flux, corresponding to a brightness temperature of 1079 K. At 4.5 $μ$m, the minimum planet flux is $30\%\pm 5\%$ of the maximum flux, corresponding to a brightness temperature of $1380\pm 65$ K. We compare our measured phase curves to the predictions of one-dimensional radiative transfer and three-dimensional general circulation models. We find that WASP-14b's measured dayside emission is consistent with a model atmosphere with equilibrium chemistry and a moderate temperature inversion. These same models tend to over-predict the nightside emission at 3.6 $μ$m, while under-predicting the nightside emission at 4.5 $μ$m. We propose that this discrepancy might be explained by an enhanced global C/O ratio. In addition, we find that the phase curves of WASP-14b ($7.8 M_{\mathrm{Jup}}$) are consistent with a much lower albedo than those of other Jovian mass planets with thermal phase curve measurements, suggesting that it may be emitting detectable heat from the deep atmosphere or interior processes.

preprint2015arXiv

A Six-Planet System Orbiting HD 219134

We present new, high-precision Doppler radial velocity (RV) data sets for the nearby K3V star HD 219134. The data include 175 velocities obtained with the HIRES Spectrograph at the Keck I Telescope, and 101 velocities obtained with the Levy Spectrograph at the Automated Planet Finder Telescope (APF) at Lick Observatory. Our observations reveal six new planetary candidates, with orbital periods of P=3.1, 6.8, 22.8, 46.7, 94.2 and 2247 days, spanning masses of msini=3.8, 3.5, 8.9, 21.3, 10.8 and 108 M_earth respectively. Our analysis indicates that the outermost signal is unlikely to be an artifact induced by stellar activity. In addition, several years of precision photometry with the T10 0.8~m automatic photometric telescope (APT) at Fairborn Observatory demonstrated a lack of brightness variability to a limit of ~0.0002 mag, providing strong support for planetary-reflex motion as the source of the radial velocity variations. The HD 219134 system, with its bright (V=5.6) primary provides an excellent opportunity to obtain detailed orbital characterization (and potentially follow-up observations) of a planetary system that resembles many of the multiple-planet systems detected by Kepler, and which are expected to be detected by NASA's forthcoming TESS Mission and by ESA's forthcoming PLATO Mission.

preprint2015arXiv

Jupiter's Decisive Role in the Inner Solar System's Early Evolution

The statistics of extrasolar planetary systems indicate that the default mode of planet formation generates planets with orbital periods shorter than 100 days, and masses substantially exceeding that of the Earth. When viewed in this context, the Solar System is unusual. Here, we present simulations which show that a popular formation scenario for Jupiter and Saturn, in which Jupiter migrates inward from a > 5 AU to a ~ 1.5 AU before reversing direction, can explain the low overall mass of the Solar System's terrestrial planets, as well as the absence of planets with a < 0.4 AU. Jupiter's inward migration entrained s ~ 10-100 km planetesimals into low-order mean-motion resonances, shepherding and exciting their orbits. The resulting collisional cascade generated a planetesimal disk that, evolving under gas drag, would have driven any pre-existing short-period planets into the Sun. In this scenario, the Solar System's terrestrial planets formed from gas-starved mass-depleted debris that remained after the primary period of dynamical evolution.

preprint2015arXiv

The Lick-Carnegie Exoplanet Survey: HD32963 -- A New Jupiter Analog Orbiting a Sun-like Star

We present a set of 109 new, high-precision Keck/HIRES radial velocity (RV) observations for the solar-type star HD 32963. Our dataset reveals a candidate planetary signal with a period of 6.49 $\pm$ 0.07 years and a corresponding minimum mass of 0.7 $\pm$ 0.03 Jupiter masses. Given Jupiter's crucial role in shaping the evolution of the early Solar System, we emphasize the importance of long-term radial velocity surveys. Finally, using our complete set of Keck radial velocities and correcting for the relative detectability of synthetic planetary candidates orbiting each of the 1,122 stars in our sample, we estimate the frequency of Jupiter analogs across our survey at approximately 3%.

preprint2014arXiv

A 4-Planet System Orbiting the K0V Star HD 141399

We present precision radial velocity (RV) data sets from Keck-HIRES and from Lick Observatory's new Automated Planet Finder Telescope and Levy Spectrometer on Mt. Hamilton that reveal a multiple-planet system orbiting the nearby, slightly evolved, K-type star HD 141399. Our 91 observations over 10.5 years suggest the presence of four planets with orbital periods of 94.35, 202.08, 1070.35, and 3717.35 days and minimum masses of 0.46, 1.36, 1.22, and 0.69 Jupiter masses respectively. The orbital eccentricities of the three inner planets are small, and the phase curves are well sampled. The inner two planets lie just outside the 2:1 resonance, suggesting that the system may have experienced dissipative evolution during the protoplanetary disk phase. The fourth companion is a Jupiter-like planet with a Jupiter-like orbital period. Its orbital eccentricity is consistent with zero, but more data will be required for an accurate eccentricity determination.

preprint2014arXiv

Constraints on the Atmospheric Circulation and Variability of the Eccentric Hot Jupiter XO-3b

We report secondary eclipse photometry of the hot Jupiter XO-3b in the 4.5~$μ$m band taken with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We measure individual eclipse depths and center of eclipse times for a total of twelve secondary eclipses. We fit these data simultaneously with two transits observed in the same band in order to obtain a global best-fit secondary eclipse depth of $0.1580\pm 0.0036\%$ and a center of eclipse phase of $0.67004\pm 0.00013 $. We assess the relative magnitude of variations in the dayside brightness of the planet by measuring the size of the residuals during ingress and egress from fitting the combined eclipse light curve with a uniform disk model and place an upper limit of 0.05$\%$. The new secondary eclipse observations extend the total baseline from one and a half years to nearly three years, allowing us to place an upper limit on the periastron precession rate of $2.9\times 10^{-3}$ degrees/day the tightest constraint to date on the periastron precession rate of a hot Jupiter. We use the new transit observations to calculate improved estimates for the system properties, including an updated orbital ephemeris. We also use the large number of secondary eclipses to obtain the most stringent limits to date on the orbit-to-orbit variability of an eccentric hot Jupiter and demonstrate the consistency of multiple-epoch Spitzer observations.

preprint2014arXiv

The 4.5 $μ$m full-orbit phase curve of the hot Jupiter HD 209458b

The hot Jupiter HD 209458b is particularly amenable to detailed study as it is among the brightest transiting exoplanet systems currently known (V-mag = 7.65; K-mag = 6.308) and has a large planet-to-star contrast ratio. HD 209458b is predicted to be in synchronous rotation about its host star with a hot spot that is shifted eastward of the substellar point by superrotating equatorial winds. Here we present the first full-orbit observations of HD 209458b, in which its 4.5 $μ$m emission was recorded with $Spitzer$/IRAC. Our study revises the previous 4.5 $μ$m measurement of HD 209458b's secondary eclipse emission downward by $\sim$35% to $0.1391%^{+0.0072%}_{-0.0069%}$, changing our interpretation of the properties of its dayside atmosphere. We find that the hot spot on the planet's dayside is shifted eastward of the substellar point by $40.9^{\circ}\pm{6.0^{\circ}}$, in agreement with circulation models predicting equatorial superrotation. HD 209458b's dayside (T$_{bright}$ = 1499 $\pm$ 15 K) and nightside (T$_{bright}$ = 972 $\pm$ 44 K) emission indicates a day-to-night brightness temperature contrast smaller than that observed for more highly irradiated exoplanets, suggesting that the day-to-night temperature contrast may be partially a function of the incident stellar radiation. The observed phase curve shape deviates modestly from global circulation model predictions potentially due to disequilibrium chemistry or deficiencies in the current hot CH$_{4}$ line lists used in these models. Observations of the phase curve at additional wavelengths are needed in order to determine the possible presence and spatial extent of a dayside temperature inversion, as well as to improve our overall understanding of this planet's atmospheric circulation.

preprint2014arXiv

The Lick-Carnegie Exoplanet Survey: Gliese 687b: A Neptune-Mass Planet Orbiting a Nearby Red Dwarf

Precision radial velocities from the Automated Planet Finder and Keck/HIRES reveal an M*sin(i) =18 +/- 2 Earth mass planet orbiting the nearby M3V star GJ 687. This planet has an orbital period, P = 38.14 days, and a low orbital eccentricity. Our Stromgren b and y photometry of the host star suggests a stellar rotation signature with a period of P = 60 days. The star is somewhat chromospherically active, with a spot filling factor estimated to be several percent. The rotationally{induced 60-day signal, however, is well-separated from the period of the radial velocity variations, instilling confidence in the interpretation of a Keplerian origin for the observed velocity variations. Although GJ 687b produces relatively little specific interest in connection with its individual properties, a compelling case can be argued that it is worthy of remark as an eminently typical, yet at a distance of 4.52 pc, a very nearby representative of the galactic planetary census. The detection of GJ 687b indicates that the APF telescope is well suited to the discovery of low-mass planets orbiting low-mass stars in the as-yet relatively un-surveyed region of the sky near the north celestial pole.

preprint2014arXiv

The Random Walk of High Frequency Trading

This paper builds a model of high-frequency equity returns by separately modeling the dynamics of trade-time returns and trade arrivals. Our main contributions are threefold. First, we characterize the distributional behavior of high-frequency asset returns both in ordinary clock time and in trade time. We show that when controlling for pre-scheduled market news events, trade-time returns of the highly liquid near-month E-mini S&P 500 futures contract are well characterized by a Gaussian distribution at very fine time scales. Second, we develop a structured and parsimonious model of clock-time returns by subordinating a trade-time Gaussian distribution with a trade arrival process that is associated with a modified Markov-Switching Multifractal Duration (MSMD) model. This model provides an excellent characterization of high-frequency inter-trade durations. Over-dispersion in this distribution of inter-trade durations leads to leptokurtosis and volatility clustering in clock-time returns, even when trade-time returns are Gaussian. Finally, we use our model to extrapolate the empirical relationship between trade rate and volatility in an effort to understand conditions of market failure. Our model suggests that the 1,200 km physical separation of financial markets in Chicago and New York/New Jersey provides a natural ceiling on systemic volatility and may contribute to market stability during periods of extremely heavy trading.

preprint2014arXiv

The Transiting Exoplanet Survey Satellite

The Transiting Exoplanet Survey Satellite (TESS) will search for planets transiting bright and nearby stars. TESS has been selected by NASA for launch in 2017 as an Astrophysics Explorer mission. The spacecraft will be placed into a highly elliptical 13.7-day orbit around the Earth. During its two-year mission, TESS will employ four wide-field optical CCD cameras to monitor at least 200,000 main-sequence dwarf stars with I = 4-13 for temporary drops in brightness caused by planetary transits. Each star will be observed for an interval ranging from one month to one year, depending mainly on the star's ecliptic latitude. The longest observing intervals will be for stars near the ecliptic poles, which are the optimal locations for follow-up observations with the James Webb Space Telescope. Brightness measurements of preselected target stars will be recorded every 2 min, and full frame images will be recorded every 30 min. TESS stars will be 10-100 times brighter than those surveyed by the pioneering Kepler mission. This will make TESS planets easier to characterize with follow-up observations. TESS is expected to find more than a thousand planets smaller than Neptune, including dozens that are comparable in size to the Earth. Public data releases will occur every four months, inviting immediate community-wide efforts to study the new planets. The TESS legacy will be a catalog of the nearest and brightest stars hosting transiting planets, which will endure as highly favorable targets for detailed investigations.

preprint2013arXiv

Host Star Properties and Transit Exclusion for the HD 38529 Planetary System

The transit signature of exoplanets provides an avenue through which characterization of exoplanetary properties may be undertaken, such as studies of mean density, structure, and atmospheric composition. The Transit Ephemeris Refinement and Monitoring Survey (TERMS) is a program to expand the catalog of transiting planets around bright host stars by refining the orbits of known planets discovered with the radial velocity technique. Here we present results for the HD 38529 system. We determine fundamental properties of the host star through direct interferometric measurements of the radius and through spectroscopic analysis. We provide new radial velocity measurements that are used to improve the Keplerian solution for the two known planets, and we find no evidence for a previously postulated third planet. We also present 12 years of precision robotic photometry of HD 38529 that demonstrate the inner planet does not transit and the host star exhibits cyclic variations in seasonal mean brightness with a timescale of approximately six years.

preprint2013arXiv

Information Transmission Between Financial Markets in Chicago and New York

High frequency trading has led to widespread efforts to reduce information propagation delays between physically distant exchanges. Using relativistically correct millisecond-resolution tick data, we document a 3-millisecond decrease in one-way communication time between the Chicago and New York areas that has occurred from April 27th, 2010 to August 17th, 2012. We attribute the first segment of this decline to the introduction of a latency-optimized fiber optic connection in late 2010. A second phase of latency decrease can be attributed to line-of-sight microwave networks, operating primarily in the 6-11 GHz region of the spectrum, licensed during 2011 and 2012. Using publicly available information, we estimate these networks' latencies and bandwidths. We estimate the total infrastructure and 5-year operations costs associated with these latency improvements to exceed $500 million.

preprint2013arXiv

Orbital Phase Variations of the Eccentric Giant Planet HAT-P-2b

We present the first secondary eclipse and phase curve observations for the highly eccentric hot Jupiter HAT-P-2b in the 3.6, 4.5, 5.8, and 8.0 μm bands of the Spitzer Space Telescope. The 3.6 and 4.5 μm data sets span an entire orbital period of HAT-P-2b, making them the longest continuous phase curve observations obtained to date and the first full-orbit observations of a planet with an eccentricity exceeding 0.2. We present an improved non-parametric method for removing the intrapixel sensitivity variations in Spitzer data at 3.6 and 4.5 μm that robustly maps position-dependent flux variations. We find that the peak in planetary flux occurs at 4.39+/-0.28, 5.84+/-0.39, and 4.68+/-0.37 hours after periapse passage with corresponding maxima in the planet/star flux ratio of 0.1138%+/-0.0089%, 0.1162%+/-0.0080%, and 0.1888%+/-0.0072% in the 3.6, 4.5, and 8.0 μm bands respectively. We compare our measured secondary eclipse depths to the predictions from a one-dimensional radiative transfer model, which suggests the possible presence of a transient day side inversion in HAT-P-2b's atmosphere near periapse. We also derive improved estimates for the system parameters, including its mass, radius, and orbital ephemeris. Our simultaneous fit to the transit, secondary eclipse, and radial velocity data allows us to determine the eccentricity and argument of periapse of HAT-P-2b's orbit with a greater precision than has been achieved for any other eccentric extrasolar planet. We also find evidence for a long-term linear trend in the radial velocity data. This trend suggests the presence of another substellar companion in the HAT-P-2 system, which could have caused HAT-P-2b to migrate inward to its present-day orbit via the Kozai mechanism.

preprint2013arXiv

Searching for Terrestrial Planets Orbiting in the Habitable Zone of Ultra-Cool Stars and Brown Dwarfs

We propose to use Kepler in 2-wheel mode to conduct a detailed search for Earth-sized planets orbiting ultra-cool stars and brown dwarfs (spectral types from M7 to L3). This population of objects presents several advantages for exoplanet surveys. First, ultra-cool stars and brown dwarfs are small and thus result in favorable planet-to-star area ratios. Second, because of their low effective temperature, the inner edge of their habitable zone is extremely close (2 to 3 days only). Third, our targets are bright at infrared wavelengths, which will enable detailed follow-up studies. Our program therefore represents a unique opportunity to find a transiting Earth-size exoplanet for which atmospheric features (including biosignatures) could be detected with near-to-come facilities such as JWST. Such exoplanet has not been discovered yet. Kepler in 2-wheel mode provides the required stability and photometric precision to make this survey successful. Our initial target sample includes 60 ultra-cool stars and brown dwarfs from which we expect to detect at least one transiting planet. We propose to monitor each source for 4 days, resulting in a total program duration of ~240 days.

preprint2013arXiv

Secondary Eclipse Photometry of the Exoplanet WASP-5b with Warm Spitzer

We present photometry of the extrasolar planet WASP-5b in the 3.6 and 4.5 micron bands taken with the Spitzer Space Telescope's Infrared Array Camera as part of the extended warm mission. By examining the depth of the planet's secondary eclipse at these two wavelengths, we can place joint constraints on the planet's atmospheric pressure-temperature profile and chemistry. We measure secondary eclipse depths of 0.197% +/- 0.028% and 0.227% +/- 0.025% in the 3.6 micron and 4.5 micron bands, respectively. Our observations are best matched by models showing a hot dayside and, depending on our choice of model, a weak thermal inversion or no inversion at all. We measure a mean offset from the predicted center of eclipse of 0.078 +/- 0.032 hours, translating to ecos(omega) = 0.0031 +/- 0.0013 and consistent with a circular orbit. We see no evidence for any eclipse timing variations comparable to those reported in a previous transit study.

preprint2013arXiv

Warm Spitzer Photometry of Three Hot Jupiters: HAT-P-3b, HAT-P-4b and HAT-P-12b

We present Warm Spitzer/IRAC secondary eclipse time series photometry of three short-period transiting exoplanets, HAT-P-3b, HAT-P-4b and HAT-P-12b, in both the available 3.6 and 4.5 micron bands. HAT-P-3b and HAT-P-4b are Jupiter-mass, objects orbiting an early K and an early G dwarf stars, respectively. For HAT-P-3b we find eclipse depths of 0.112%+0.015%-0.030% (3.6 micron) and 0.094%+0.016%-0.009% (4.5 micron). The HAT-P-4b values are 0.142%+0.014%-0.016% (3.6 micron) and 0.122%+0.012%-0.014% (4.5micron). The two planets' photometry is consistent with inefficient heat redistribution from their day to night sides (and low albedos), but it is inconclusive about possible temperature inversions in their atmospheres. HAT-P-12b is a Saturn-mass planet and is one of the coolest planets ever observed during secondary eclipse, along with hot Neptune GJ 436b and hot Saturn WASP-29b. We are able to place 3$σ$ upper limits on the secondary eclipse depth of HAT-P-12b in both wavelengths: < 0.042% (3.6 micron) and <0.085% (4.5 micron). We discuss these results in the context of the {\it Spitzer} secondary eclipse measurements of GJ 436b and WASP-29b. It is possible that we do not detect the eclipses of HAT-P-12b due to high eccentricity, but find that weak planetary emission in these wavelengths is a more likely explanation. We place 3 sigma upper limits on the |e cos(omega)| quantity (where e is eccentricity and omega is the argument of periapsis) for HAT-P-3b (<0.0081) and HAT-P-4b (<0.0042), based on the secondary eclipse timings.

preprint2012arXiv

3.6 and 4.5 Micron Phase Curves and Evidence for Non-Equilibrium Chemistry in the Atmosphere of Extrasolar Planet HD 189733b

We present new, full-orbit observations of the infrared phase variations of the canonical hot Jupiter HD 189733b obtained in the 3.6 and 4.5 micron bands using the Spitzer Space Telescope. When combined with previous phase curve observations at 8.0 and 24 micron, these data allow us to characterize the exoplanet's emission spectrum as a function of planetary longitude. We utilize improved methods for removing the effects of intrapixel sensitivity variations and accounting for the presence of time-correlated noise in our data. We measure a phase curve amplitude of 0.1242% +/- 0.0061% in the 3.6 micron band and 0.0982% +/- 0.0089% in the 4.5 micron band. We find that the times of minimum and maximum flux occur several hours earlier than predicted for an atmosphere in radiative equilibrium, consistent with the eastward advection of gas by an equatorial super-rotating jet. The locations of the flux minima in our new data differ from our previous observations at 8 micron, and we present new evidence indicating that the flux minimum observed in the 8 micron is likely caused by an over-shooting effect in the 8 micron array. We obtain improved estimates for HD 189733b's dayside planet-star flux ratio of 0.1466% +/- 0.0040% at 3.6 micron and 0.1787% +/- 0.0038% at 4.5 micron; these are the most accurate secondary eclipse depths obtained to date for an extrasolar planet. We compare our new dayside and nightside spectra for HD 189733b to the predictions of models from Burrows et al. (2008) and Showman et al. (2009). We find that HD 189733b's 4.5 micron nightside flux is 3.3 sigma smaller than predicted by the Showman et al. models, which assume that the chemistry is in local thermal equilibrium. We conclude that this discrepancy is best-explained by vertical mixing, which should lead to an excess of CO and correspondingly enhanced 4.5 micron absorption in this region. [abridged]

preprint2012arXiv

The HD 192263 system: planetary orbital period and stellar variability disentangled

As part of the Transit Ephemeris Refinement and Monitoring Survey (TERMS), we present new radial velocities and photometry of the HD 192263 system. Our analysis of the already available Keck-HIRES and CORALIE radial velocity measurements together with the five new Keck measurements we report in this paper results in improved orbital parameters for the system. We derive constraints on the size and phase location of the transit window for HD 192263b, a Jupiter-mass planet with a period of 24.3587 \pm 0.0022 days. We use 10 years of Automated Photoelectric Telescope (APT) photometry to analyze the stellar variability and search for planetary transits. We find continuing evidence of spot activity with periods near 23.4 days. The shape of the corresponding photometric variations changes over time, giving rise to not one but several Fourier peaks near this value. However, none of these frequencies coincides with the planet's orbital period and thus we find no evidence of star-planet interactions in the system. We attribute the ~23-day variability to stellar rotation. There are also indications of spot variations on longer (8 years) timescales. Finally, we use the photometric data to exclude transits for a planet with the predicted radius of 1.09 RJ, and as small as 0.79 RJ.

preprint2011arXiv

A Search for the Transit of HD 168443b: Improved Orbital Parameters and Photometry

The discovery of transiting planets around bright stars holds the potential to greatly enhance our understanding of planetary atmospheres. In this work we present the search for transits of HD 168443b, a massive planet orbiting the bright star HD 168443 V=6.92 with a period of 58.11 days. The high eccentricity of the planetary orbit e=0.53 significantly enhances the a-priori transit probability beyond that expected for a circular orbit, making HD 168443 a candidate for our ongoing Transit Ephemeris Refinement and Monitoring Survey (TERMS). Using additional radial velocities from Keck-HIRES, we refined the orbital parameters of this multi-planet system and derived a new transit ephemeris for HD 168443b. The reduced uncertainties in the transit window make a photometric transit search practicable. Photometric observations acquired during predicted transit windows were obtained on three nights. CTIO 1.0 m photometry acquired on 2010 September 7 had the required precision to detect a transit but fell just outside of our final transit window. Nightly photometry from the T8 0.8 m Automated Photometric Telescope (APT) at Fairborn Observatory, acquired over a span of 109 nights, demonstrates that HD 168443 is constant on a time scale of weeks. Higher-cadence photometry on 2011 April 28 and June 25 shows no evidence of a transit. We are able to rule out a non-grazing transit of HD 168443b.

preprint2011arXiv

A Spitzer Transmission Spectrum for the Exoplanet GJ 436b, Evidence for Stellar Variability, and Constraints on Dayside Flux Variations

In this paper we describe a uniform analysis of eight transits and eleven secondary eclipses of the extrasolar planet GJ 436b obtained in the 3.6, 4.5, and 8.0 micron bands using the IRAC instrument on the Spitzer Space Telescope between UT 2007 June 29 and UT 2009 Feb 4. We find that the best-fit transit depths for visits in the same bandpass can vary by as much as 8% of the total (4.7 sigma significance) from one epoch to the next. Although we cannot entirely rule out residual detector effects or a time-varying, high-altitude cloud layer in the planet's atmosphere as the cause of these variations, we consider the occultation of active regions on the star in a subset of the transit observations to be the most likely explanation. We reconcile the presence of magnetically active regions with the lack of significant visible or infrared flux variations from the star by proposing that the star's spin axis is tilted with respect to our line of sight, and that the planet's orbit is therefore likely to be misaligned. These observations serve to illustrate the challenges associated with transmission spectroscopy of planets orbiting late-type stars; we expect that other systems, such as GJ 1214, may display comparably variable transit depths. Our measured 8 micron secondary eclipse depths are consistent with a constant value, and we place a 1 sigma upper limit of 17% on changes in the planet's dayside flux in this band. Averaging over the eleven visits gives us an improved estimate of 0.0452% +/- 0.0027% for the secondary eclipse depth. We combine timing information from our observations with previously published data to produce a refined orbital ephemeris, and determine that the best-fit transit and eclipse times are consistent with a constant orbital period. [ABRIDGED]

preprint2011arXiv

Improved Orbital Parameters and Transit Monitoring for HD 156846b

HD 156846b is a Jovian planet in a highly eccentric orbit (e = 0.85) with a period of 359.55 days. The pericenter passage at a distance of 0.16 AU is nearly aligned to our line of sight, offering an enhanced transit probability of 5.4% and a potentially rich probe of the dynamics of a cool planetary atmosphere impulsively heated during close approach to a bright star (V = 6.5). We present new radial velocity (RV) and photometric measurements of this star as part of the Transit Ephemeris Refinement and Monitoring Survey (TERMS). The RV measurements from Keck-HIRES reduce the predicted transit time uncertainty to 20 minutes, an order of magnitude improvement over the ephemeris from the discovery paper. We photometrically monitored a predicted transit window under relatively poor photometric conditions, from which our non-detection does not rule out a transiting geometry. We also present photometry that demonstrates stability at the millimag level over its rotational timescale.

preprint2011arXiv

Kepler Exoplanet Candidate Host Stars are Preferentially Metal Rich

We find that Kepler exoplanet candidate (EC) host stars are preferentially metal-rich, including the low-mass stellar hosts of small-radius ECs. The last observation confirms a tentative hint that there is a correlation between the metallicity of low-mass stars and the presence of low-mass and small-radius exoplanets. In particular, we compare the J-H--g-r color-color distribution of Kepler EC host stars with a control sample of dwarf stars selected from the ~150,000 stars observed during Q1 and Q2 of the Kepler mission but with no detected planets. We find that at J-H = 0.30 characteristic of solar-type stars, the average g-r color of stars that host giant ECs is 4-sigma redder than the average color of the stars in the control sample. At the same time, the average g-r color of solar-type stars that host small-radius ECs is indistinguishable from the average color of the stars in the control sample. In addition, we find that at J-H = 0.62 indicative of late K dwarfs, the average g-r color of stars that host small-radius ECs is 4-sigma redder than the average color of the stars in the control sample. These offsets are unlikely to be caused by differential reddening, age differences between the two populations, or the presence of giant stars in the control sample. Stellar models suggest that the first color offset is due to a 0.2 dex enhancement in [Fe/H] of the giant EC host population at M_star = 1 M_Sun, while Sloan photometry of M 67 and NGC 6791 suggests that the second color offset is due to a similar [Fe/H] enhancement of the small-radius EC host population at M_star = 0.7 M_Sun. These correlations are a natural consequence of the core-accretion model of planet formation.

preprint2011arXiv

On the Anomalous Radii of the Transiting Extrasolar Planets

We present a systematic evaluation of the agreement between the observed radii of 90 well-characterized transiting extrasolar giant planets and their corresponding model radii. Our model radii are drawn from previously published calculations of core-less giant planets that have attained their asymptotic radii, and which have been tabulated for a range of planet masses and equilibrium temperatures. (We report a two-dimensional polynomial fitting function that accurately represents the models). As expected, the model radii provide a statistically significant improvement over a null hypothesis that the sizes of giant planets are completely independent of mass and effective temperature. As is well known, however, fiducial models provide an insufficient explanation; the planetary radius anomalies are strongly correlated with planetary equilibrium temperature. We find that the radius anomalies have a best-fit dependence, ${\cal R}\propto T_{\rm eff}^α$, with $α=1.4\pm0.6$. Incorporating this relation into the model radii leads to substantially less scatter in the radius correlation. The extra temperature dependence represents an important constraint on theoretical models for Hot Jupiters. Using simple scaling arguments, we find support for the hypothesis of Batygin and Stevenson (2010) that this correlation can be attributed to a planetary heating mechanism that is mediated by magnetohydrodynamic coupling between the planetary magnetic field and near-surface flow that is accompanied by ohmic dissipation at adiabatic depth. Additionally, we find that the temperature dependence is likely too strong to admit kinetic heating as the primary source of anomalous energy generation within the majority of the observed transiting planets.

preprint2011arXiv

Resolving the Sin(I) degeneracy in Low-Mass Multi-Planet Systems

Long-term orbital evolution of multi-planet systems under tidal dissipation often converges to a stationary state, known as the tidal fixed point. The fixed point is characterized by a lack of oscillations in the eccentricities and apsidal alignment among the orbits. Quantitatively, the nature of the fixed point is dictated by mutual interactions among the planets as well as non-Keplerian effects. We show that if a roughly coplanar system hosts a hot, sub-Saturn mass planet, and is tidally relaxed, separation of planet-planet interactions and non-Keplerian effects in the equations of motion leads to a direct determination of the true masses of the planets. Consequently, a "snap-shot" observational determination of the orbital state resolves the sin(I) degeneracy, and opens up a direct avenue towards identification of the true lowest-mass exo-planets detected. We present an approximate, as well as a general, mathematical framework for computation of the line of sight inclination of secular systems, and apply our models illustratively to the 61 Vir system. We conclude by discussing the observability of planetary systems to which our method is applicable and we set our analysis into a broader context by presenting a current summary of the various possibilities for determining the physical properties of planets from observations of their orbital states.

preprint2011arXiv

Warm Spitzer Photometry of XO-4b, HAT-P-6b and HAT-P-8b

We have analyzed Warm Spitzer/IRAC observations of the secondary eclipses of three planets, XO-4b, HAT-P-6b and HAT-P-8b. We measure secondary eclipse amplitudes at 3.6μm and 4.5μm for each target. XO-4b exhibits a stronger eclipse depth at 4.5μm than at 3.6μm, which is consistent with the presence of a temperature inversion. HAT-P-8b shows a stronger eclipse amplitude at 3.6μm, and is best-described by models without a temperature inversion. The eclipse depths of HAT-P-6b can be fitted with models with a small or no temperature inversion. We consider our results in the context of a postulated relationship between stellar activity and temperature inversions and a relationship between irradiation level and planet dayside temperature, as discussed by Knutson et al. (2010) and Cowan & Agol (2011), respectively. Our results are consistent with these hypotheses, but do not significantly strengthen them. To measure accurate secondary eclipse central phases, we require accurate ephemerides. We obtain primary transit observations and supplement them with publicly available observations to update the orbital ephemerides of the three planets. Based on the secondary eclipse timing, we set upper boundaries for e cos(ω) for HAT-P-6b, HAT-P-8b and XO-4b and find that the values are consistent with circular orbits.

preprint2010arXiv

A Physically-Motivated Photometric Calibration of M Dwarf Metallicity

The location of M dwarfs in the V-K_s--M_Ks color-magnitude diagram (CMD) has been shown to correlate with metallicity. We demonstrate that previous empirical photometric calibrations of M dwarf metallicity exploiting this correlation systematically underestimate or overestimate metallicity at the extremes of their range. We improve upon previous calibrations in three ways. We use both a volume-limited and kinematically-matched sample of F and G dwarfs from the Geneva-Copehnagen Survey (GCS) to infer the mean metallicity of M dwarfs in the Solar Neighborhood, we use theoretical models of M dwarf interiors and atmospheres to determine the effect of metallicity on M dwarfs in the V-K_s--M_Ks CMD, and we base our final calibration purely on high-resolution spectroscopy of FGK primaries with M dwarf companions. As a result, we explain an order of magnitude more of the variance in the calibration sample than previous photometric calibrations. We non-parametrically quantify the significance of the observation that M dwarfs that host exoplanets are preferentially in a region of the V-K_s--M_Ks plane populated by metal-rich M dwarfs. We find that the probability p that planet-hosting M dwarfs are distributed across the V-K_s--M_Ks CMD in the same way as field M dwarfs is p = 0.06 +/- 0.008. Interestingly, the subsample of M dwarfs that host Neptune and sub-Neptune mass planets may also be preferentially located in the region of the V-K_s--M_Ks plane populated by high-metallicity M dwarfs. The probability of this occurrence by chance is p = 0.40 +/- 0.02, and this observation hints that low-mass planets may be more likely to be found around metal-rich M dwarfs. An increased rate of low-mass planet occurrence around metal-rich M dwarfs would be a natural consequence of the core-accretion model of planet formation. (abridged)

preprint2010arXiv

A Scientometric Prediction of the Discovery of the First Potentially Habitable Planet with a Mass Similar to Earth

The search for a habitable extrasolar planet has long interested scientists, but only recently have the tools become available to search for such planets. In the past decades, the number of known extrasolar planets has ballooned into the hundreds, and with it the expectation that the discovery of the first Earth-like extrasolar planet is not far off. Here we develop a novel metric of habitability for discovered planets, and use this to arrive at a prediction for when the first habitable planet will be discovered. Using a bootstrap analysis of currently discovered exoplanets, we predict the discovery of the first Earth-like planet to be announced in the first half of 2011, with the likeliest date being early May 2011. Our predictions, using only the properties of previously discovered exoplanets, accord well with external estimates for the discovery of the first potentially habitable extrasolar planet, and highlights the the usefulness of predictive scientometric techniques to understand the pace of scientific discovery in many fields.

preprint2010arXiv

Secondary Eclipse Photometry of WASP-4b with Warm Spitzer

We present photometry of the giant extrasolar planet WASP-4b at 3.6 and 4.5 micron taken with the Infrared Array Camera on board the Spitzer Space Telescope as part of Spitzer's extended warm mission. We find secondary eclipse depths of 0.319+/-0.031% and 0.343+/-0.027% for the 3.6 and 4.5 micron bands, respectively and show model emission spectra and pressure-temperature profiles for the planetary atmosphere. These eclipse depths are well fit by model emission spectra with water and other molecules in absorption, similar to those used for TrES-3 and HD 189733b. Depending on our choice of model, these results indicate that this planet has either a weak dayside temperature inversion or no inversion at all. The absence of a strong thermal inversion on this highly irradiated planet is contrary to the idea that highly irradiated planets are expected to have inversions, perhaps due the presence of an unknown absorber in the upper atmosphere. This result might be explained by the modestly enhanced activity level of WASP-4b's G7V host star, which could increase the amount of UV flux received by the planet, therefore reducing the abundance of the unknown stratospheric absorber in the planetary atmosphere as suggested in Knutson et al. (2010). We also find no evidence for an offset in the timing of the secondary eclipse and place a 2 sigma upper limit on |ecos(omega)| of 0.0024, which constrains the range of tidal heating models that could explain this planet's inflated radius.

preprint2010arXiv

Systemic: A Testbed for Characterizing the Detection of Extrasolar Planets. II. Numerical approaches to the Transit Timing Inverse Problem

Transit timing variations - deviations from strict periodicity between successive passages of a transiting planet - can be used to probe the structure and dynamics of multiple-planet systems. In this paper, we examine prospects for numerically solving the so-called inverse problem, the determination of the orbital elements of a perturbing body from the transit timing variations it induces. We assume that the planetary systems under examination have a limited number of Doppler velocity measurements, and show that a more extensive radial velocity characterization with precision comparable to the semiamplitude of the perturber may remove degeneracies in the solution. We examine several configurations of interest, including (1) a prototypical non-resonant system, modeled after HD40307 b and c, which contains multiple super-Earth mass planets, (2) a hypothetical system containing a transiting giant planet with a terrestrial-mass companion trapped in low-order mean motion resonance, and (3) the HAT-P-13 system, in which forced precession by an outer perturbing body that is well characterized by Doppler radial velocity measurements can give insight into the interior structure of a perturbing planet, and for which the determination of mutual inclination between the transiting planet and its perturber is a key issue.

preprint2010arXiv

The Lick-Carnegie Exoplanet Survey: A Uranus-mass Fourth Planet for GJ 876 in an Extrasolar Laplace Configuration

(Abreviated) Continued radial velocity monitoring of the nearby M4V red dwarf star GJ~876 with Keck/HIRES has revealed the presence of a Uranus-mass fourth planetary companion in the system. The new planet has a mean period of $P_e=126.6$ days (over the 12.6-year baseline of the radial velocity observations), and a minimum mass of $m_e\sin{i_e}=12.9\pm 1.7\,M_{\oplus}$. Self-consistent, N-body fits to the radial velocity data set show that the four-planet system has an invariable plane with an inclination relative to the plane of the sky of $i=59.5^{\circ}$. The fit is not significantly improved by the introduction of a mutual inclination between the planets ``b'' and ``c,'' but the new data do confirm a non-zero eccentricity, $e_d=0.207\pm0.055$ for the innermost planet, ``d.'' In our best-fit coplanar model, the mass of the new component is $m_e=14.6\pm1.7\,M_{\oplus}$. Our best-fitting model places the new planet in a 3-body resonance with the previously known giant planets (which have mean periods of $P_c=30.4$ and $P_b=61.1$ days). The critical argument, $φ_{\rm Laplace}=λ_c-3λ_b+2λ_e$, for the Laplace resonance librates with an amplitude of $Δφ_{\rm Laplace}=40\pm13^{\circ}$ about $φ_{\rm Laplace}=0^{\circ}$. Numerical integration indicates that the four-planet system is stable for at least a billion years (at least for the coplanar cases). This resonant configuration of three giant planets orbiting an M-dwarf primary differs from the well-known Laplace configuration of the three inner Galilean satellites of Jupiter, which are executing very small librations about $φ_{\rm Laplace}=180^{\circ}$, and which never experience triple conjunctions. The GJ~876 system, by contrast, comes close to a triple conjunction between the outer three planets once per every orbit of the outer planet, ``e.''

preprint2010arXiv

The Lick-Carnegie Survey: Four New Exoplanet Candidates

We present new precise HIRES radial velocity (RV) data sets of five nearby stars obtained at Keck Observatory. HD 31253, HD 218566, HD 177830, HD 99492 and HD 74156 are host stars of spectral classes F through K and show radial velocity variations consistent with new or additional planetary companions in Keplerian motion. The orbital parameters of the candidate planets in the five planetary systems span minimum masses of M sin i = 27.43 M_{earth} to M sin i = 8.28 M_{jup}, periods of 17.05 to 4696.95 days and eccentricities ranging from circular to extremely eccentric (e ~ 0.63). The 5th star, HD 74156, was known to have both a 52-day and a 2500-day planet, and was claimed to also harbor a 3rd planet at 336d, in apparent support of the "Packed Planetary System" hypothesis. Our greatly expanded data set for HD 74156 provides strong confirmation of both the 52-day and 2500-d planets, but strongly contradicts the existence of a 336-day planet, and offers no significant evidence for any other planets in the system.

preprint2010arXiv

Warm Spitzer Photometry of the Transiting Exoplanets CoRoT-1 and CoRoT-2 at Secondary Eclipse

We measure secondary eclipses of the hot giant exoplanets CoRoT-1 at 3.6 and 4.5 microns, and CoRoT-2 at 3.6 microns, both using Warm Spitzer. We find that the Warm Spitzer mission is working very well for exoplanet science. For consistency of our analysis we also re-analyze archival cryogenic Spitzer data for secondary eclipses of CoRoT-2 at 4.5 and 8 microns. We compare the total data for both planets, including optical eclipse measurements by the CoRoT mission, and ground-based eclipse measurements at 2 microns, to existing models. Both planets exhibit stronger eclipses at 4.5 than at 3.6 microns, which is often indicative of an atmospheric temperature inversion. The spectrum of CoRoT-1 is best reproduced by a 2460K blackbody, due either to a high altitude layer that strongly absorbs stellar irradiance, or an isothermal region in the planetary atmosphere. The spectrum of CoRoT-2 is unusual because the 8 micron contrast is anomalously low. Non-inverted atmospheres could potentially produce the CoRoT-2 spectrum if the planet exhibits line emission from CO at 4.5 microns, caused by tidal-induced mass loss. However, the viability of that hypothesis is questionable because the emitting region cannot be more than about 30-percent larger than the planet's transit radius, based on the ingress and egress times at eclipse. An alternative possibility to account for the spectrum of CoRoT-2 is an additional opacity source that acts strongly at wavelengths less than 5 microns, heating the upper atmosphere while allowing the deeper atmosphere seen at 8 microns to remain cooler. We obtain a similar result as Gillon et al. for the phase of the secondary eclipse of CoRoT-2, implying an eccentric orbit with e*cos(omega)=-0.0030 +/- 0.0004.

preprint2009arXiv

A Neptune-mass Planet Orbiting the Nearby G Dwarf HD16417

Precision Doppler measurements from an intensive 48 night "Rocky Planet Search" observing campaign on the Anglo-Australian Telescope (AAT) have revealed the presence of a low-mass exoplanet orbiting the G1 dwarf HD16417. Subsequent Doppler observations with the AAT, as well as independent observations obtained by the Keck Planet Search, have confirmed this initial detection and refine the orbital parameters to period 17.24+/-0.01 d, eccentricity 0.20+/-0.09, orbital semi-major axis 0.14+/-0.01 AU and minimum planet mass 22.1+/-2.0 Mearth. HD 16417 raises the number of published exoplanets with minimum masses of less than 25 Mearth to eighteen. Interestingly, the distribution of detected sub-25 Mearth planets over the spectral types G, K and M is almost uniform. The detection of HD 16417b by an intensive observing campaign clearly demonstrates the need for extended and contiguous observing campaigns when aiming to detect low-amplitude Doppler planets in short period orbits. Perhaps most critically it demonstrates that the search for low-mass Doppler planets will eventually require these traditional "bright-time" projects to extend throughout dark lunations.

preprint2009arXiv

A Quasi-Stationary Solution to Gliese 436b's Eccentricity

We investigate the possibility that the large orbital eccentricity of the transiting Neptune-mass planet Gliese 436b is maintained in the face of tidal dissipation by a second planet in the system. We find that the currently observed configuration can be understood if Gliese 436b and a putative companion have evolved to a quasi-stationary fixed point in which the planets' orbital apses are co-linear and in which secular variations in the orbital eccentricities of the two planets have been almost entirely damped out. We adopt an octopole-order secular theory based on a Legendre expansion in the semi-major axis ratio to delineate well-defined regions of (P_c, M_c, e_c) space that can be occupied by a perturbing companion. We incorporate the evolutionary effect of tidal dissipation into our secular model of the system, and solve the resulting initial value problems for a large sample of the allowed configurations. We then polish the stationary configurations derived from secular theory with full numerical integrations. We present our results in the form of candidate companion planets to Gliese 436b. For these candidates, radial velocity half-amplitudes, K_c, are of order 3 m/s, and the maximum amplitude of orbit-to-orbit transit timing variations are of order Delta t=1 s to Delta t=5s. For the particular example case of a perturber with orbital period, P_c=40 d, mass, M_c=8.5 M_Earth, and eccentricity, e_c=0.58, we confirm our semi-analytic calculations with a full numerical 3-body integration of the orbital decay that includes tidal damping and spin evolution. Additionally, we discuss the possibility of many-perturber stationary configurations, utilizing modified Laplace-Lagrange secular theory.

preprint2009arXiv

A Super-Earth and two Neptunes Orbiting the Nearby Sun-like star 61 Virginis

We present precision radial velocity (RV) data that reveal a multiple exoplanet system orbiting the bright nearby G5V star 61 Virginis. Our 4.6 years of combined Keck/HIRES and Anglo-Australian Telescope precision RVs indicate the hitherto unknown presence of at least three planets orbiting this well-studied star. These planets are all on low-eccentricity orbits with periods of 4.2, 38.0, and 124.0 days, and projected masses (M sini) of 5.1, 18.2, and 24.0 M_Earth, respectively. Test integrations of systems consistent with the RV data suggest that the configuration is dynamically stable. Depending on the effectiveness of tidal dissipation within the inner planet, the inner two planets may have evolved into an eccentricity fixed-point configuration in which the apsidal lines of all three planets corotate. This conjecture can be tested with additional observations. We present a 16-year time series of photometric observations of 61 Virginis, which comprise 1194 individual measurements, and indicate that it has excellent photometric stability. No significant photometric variations at the periods of the proposed planets have been detected. This new system is the first known example of a G-type Sun-like star hosting a Super-Earth mass planet. It joins HD 75732 (55 Cnc), HD 69830, GJ 581, HD 40307, and GJ 876 as a growing group of exoplanet systems that have multiple planets orbiting with periods less than an Earth-year. The ubiquity of such systems portends that space-based transit-search missions such as KEPLER and COROT will find many multi-transiting systems.

preprint2009arXiv

A Super-Earth Orbiting the Nearby Sun-like Star HD 1461

We present precision radial velocity data that reveal a Super-Earth mass planet and two probable additional planets orbiting the bright nearby G0V star HD 1461. Our 12.8 years of Keck HIRES precision radial velocities indicate the presence of a 7.4M_Earth planet on a 5.77-day orbit. The data also suggest, but cannot yet confirm, the presence of outer planets on low-eccentricity orbits with periods of 446.1 and 5017 days, and projected masses (M sin i) of 27.9 and 87.1M_Earth, respectively. Test integrations of systems consistent with the radial velocity data suggest that the configuration is dynamically stable. We present a 12.2-year time series of photometric observations of HD 1461, which comprise 799 individual measurements, and indicate that it has excellent long-term photometric stability. However, there are small amplitude variations with periods comparable to those of the suspected 2nd and 3rd signals in the radial velocities near 5000 and 446 days, thus casting some suspicion on those periodicities as Keplerian signals. If the 5.77-day companion has a Neptune-like composition, then its expected transit depth is of order ~0.5 millimags. The geometric a priori probability of transits is ~8%. Phase-folding of the ground-based photometry shows no indication that transits of the 5.77-day companion are occurring, but high-precision follow-up of HD 1461 during upcoming transit phase windows will be required to definitively rule out or confirm transits. This new system joins a growing list of solar-type stars in the immediate galactic neighborhood that are accompanied by at least one Neptune- (or lower) mass planets having orbital periods of 50 days or less.

preprint2009arXiv

Determination of the Interior Structure of Transiting Planets in Multiple-Planet Systems

Tidal dissipation within a short-period transiting extrasolar planet perturbed by a companion object can drive orbital evolution of the system to a so-called tidal fixed point, in which the apsidal lines of the transiting planet and its perturber are aligned, and for which variations in the orbital eccentricities of both planet and perturber are damped out. Significant contributions to the apsidal precession rate are made by the secular planet-planet interaction, by general relativity, and by the gravitational quadropole fields created by the transiting planet's tidal and rotational distortions. The fixed-point orbital eccentricity of the inner planet is therefore a strong function of the planet's interior structure. We illustrate these ideas in the specific context of the recently discovered HAT-P-13 exo-planetary system, and show that one can already glean important insights into the physical properties of the inner transiting planet. We present structural models of the planet, which indicate that its observed radius can be maintained for a one-parameter sequence of models that properly vary core mass and tidal energy dissipation in the interior. We use an octopole-order secular theory of the orbital dynamics to derive the dependence of the inner planet's eccentricity, on its tidal Love number. We find that the currently measured eccentricity, implies 0.116 < k2_{b} < 0.425, 0 M_{Earth}<M_{core}<120 M_{Earth}$, and Q_{b} < 300,000. Improved measurement of the eccentricity will soon allow for far tighter limits to be placed on all three of these quantities, and will provide an unprecedented probe into the interior structure of an extrasolar planet.

preprint2009arXiv

Refining Exoplanet Ephemerides and Transit Observing Strategies

Transiting planet discoveries have yielded a plethora of information regarding the internal structure and atmospheres of extra-solar planets. These discoveries have been restricted to the low-periastron distance regime due to the bias inherent in the geometric transit probability. Monitoring known radial velocity planets at predicted transit times is a proven method of detecting transits, and presents an avenue through which to explore the mass-radius relationship of exoplanets in new regions of period/periastron space. Here we describe transit window calculations for known radial velocity planets, techniques for refining their transit ephemerides, target selection criteria, and observational methods for obtaining maximum coverage of transit windows. These methods are currently being implemented by the Transit Ephemeris Refinement and Monitoring Survey (TERMS).

preprint2009arXiv

Systemic: A Testbed For Characterizing the Detection of Extrasolar Planets. I. The Systemic Console Package

We present the systemic Console, a new all-in-one, general-purpose software package for the analysis and combined multiparameter fitting of Doppler radial velocity (RV) and transit timing observations. We give an overview of the computational algorithms implemented in the Console, and describe the tools offered for streamlining the characterization of planetary systems. We illustrate the capabilities of the package by analyzing an updated radial velocity data set for the HD128311 planetary system. HD128311 harbors a pair of planets that appear to be participating in a 2:1 mean motion resonance. We show that the dynamical configuration cannot be fully determined from the current data. We find that if a planetary system like HD128311 is found to undergo transits, then self-consistent Newtonian fits to combined radial velocity data and a small number of timing measurements of transit midpoints can provide an immediate and vastly improved characterization of the planet's dynamical state.

preprint2007arXiv

Observational Consequences of Hydrodynamic Flows on Hot Jupiters

We use a grid-based shallow water model to simulate the atmospheric dynamics of the transiting hot Jupiter HD 209458b. Under the usual assumption that the planet is in synchronous rotation with zero obliquity, a steady state is reached with a well-localized cold spot centered 76 degrees east of the antistellar point. This represents a departure from predictions made by previous simulations in the literature that used the shallow water formalism; we find that the disagreement is explained by the factor of 30 shorter radiative timescale used in our model. We also examine the case that the planet is in Cassini state 2, in which the expected obliquity is ~90 degrees. Under these circumstances, a periodic equilibrium is reached, with the temperature slightly leading the solar forcing. Using these temperature distributions, we calculate disk-integrated bolometric infrared light curves from the planet. The light curves for the two models are surprisingly similar, despite large differences in temperature patterns in the two cases. In the zero-obliquity case, the intensity at the minimum is 66% of the maximum intensity, with the minimum occuring 72 degrees ahead of transit. In the high-obliquity case, the minimum occurs 54 degrees ahead of transit, with an intensity of 58% of the maximum.

preprint2005arXiv

Identifying Very Metal-Rich Stars with Low-Resolution Spectra: Finding Planet-Search Targets

We present empirical calibrations that estimate stellar metallicity, effective temperature and surface gravity as a function of Lick/IDS indices. These calibrations have been derived from a training set of 261 stars for which (1) high-precision measurements of [Fe/H], T_eff and log g have been made using spectral-synthesis analysis of HIRES spectra, and (2) Lick indices have also been measured. Our [Fe/H] calibration, which has precision 0.07 dex, has identified a number of bright (V < 9) metal-rich stars which are now being screened for hot Jupiter-type planets. Using the Yonsei-Yale stellar models, we show that the calibrations provide distance estimates accurate to 20% for nearby stars. This paper outlines the second tier of the screening of planet-search targets by the N2K Consortium, a project designed to identify the stars most likely to harbor extrasolar planets. Discoveries by the N2K Consortium include the transiting hot Saturn HD 149026 b (Sato et al. 2005, astro-ph/0507009) and HD 88133 b (Fischer et al. 2005). See Ammons et al. (2005, In Press) for a description of the first tier of N2K metallicity screening, calibrations using broadband photometry.

preprint2004arXiv

Core-Accretion Model Predicts Few Jovian-Mass Planets Orbiting Red Dwarfs

The favored theoretical explanation for giant planet formation -- in both our solar system and others -- is the core accretion model (although it still has some serious difficulties). In this scenario, planetesimals accumulate to build up planetary cores, which then accrete nebular gas. With current opacity estimates for protoplanetary envelopes, this model predicts the formation of Jupiter-mass planets in 2--3 Myr at 5 AU around solar-mass stars, provided that the surface density of solids is enhanced over that of the minimum-mass solar nebula (by a factor of a few). Working within the core-accretion paradigm, this paper presents theoretical calculations which show that the formation of Jupiter-mass planets orbiting M dwarf stars is seriously inhibited at all radial locations (in sharp contrast to solar-type stars). Planet detection programs sensitive to companions of M dwarfs will test this prediction in the near future.

preprint2002arXiv

An Improved uvby-Metallicity Calibration for Metal-Rich Stars

We present an improved uvby-metallicity relation calibrated for F, G, and early K dwarfs, and an analogous uvby-T(eff) relation, both derived using a Levenberg-Marquardt minimization scheme. Our calibrations are based on 1533 stars which appear in both the Cayrel de Strobel (2001) metallicity compilation, and in the Hauck-Mermilliod (1998) catalog of uvby photometry. We also examine the speculative possibility of using uvby photometry to produce a uvby-planeticity calibration. We conclude that while there is likely no strong photometric indicator of the presence or absence of short-period planets, the possibility of a spectroscopic indicator of planeticity is well worth examining.

preprint1997arXiv

Constraints on the Intergalactic Transport of Cosmic Rays

Motivated by recent experimental proposals to search for extragalactic cosmic rays (including anti-matter from distant galaxies), we study particle propagation through the intergalactic medium (IGM). We first use estimates of the magnetic field strength between galaxies to constrain the mean free path for diffusion of particles through the IGM. We then develop a simple analytic model to describe the diffusion of cosmic rays. Given the current age of galaxies, our results indicate that, in reasonable models, a completely negligible number of particles can enter our Galaxy from distances greater than $\sim 100$ Mpc for relatively low energies ($E$ $< 10^6$ GeV/n). We also find that particle destruction in galaxies along the diffusion path produces an exponential suppression of the possible flux of extragalactic cosmic rays. Finally, we use gamma ray constraints to argue that the distance to any hypothetical domains of anti-matter must be roughly comparable to the horizon scale.