Source author record

Rory Barnes

Rory Barnes 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

50works
5topics
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

50 published item(s)

preprint2022arXiv

The Contribution of M-Dwarf Flares to the Thermal Escape of Potentially Habitable Planet Atmospheres

The habitability of planets around M dwarfs ($\lesssim 0.5 M_\odot$) can be affected by the XUV (X rays + extreme UV) emission of these stars, with flares occasionally increasing the XUV flux by more than 2 orders of magnitude above quiescent levels. This wavelength range can warm and ionize terrestrial planets' upper atmospheres, which expands the planetary radius and promotes atmospheric loss. In this work, we study the contribution of the XUV flux due to flares on the atmospheric escape of Earth-like planets orbiting M dwarfs through numerical simulations. We considered the first Gyr of planets with initial surface water abundances between 1 and 10 terrestrial oceans (TO), a small primordial hydrogen envelope ($\le$ $10^{-3}$~$M_{\oplus}$), and with host star masses between 0.2 and 0.6 $M_{\odot}$. In this parameter range, we find that flares can remove up to two TO more than nonflaring stars, which, in some cases, translates to a doubling of the total water loss. We also find that flaring can increase atmospheric oxygen partial pressures by hundreds of bars in some cases. These results were obtained by adding a new module for flares to the \vplanet software package and upgrading its atmospheric escape module to account for Roche lobe overflow and radiation/recombination-limited escape.

preprint2020arXiv

A Coupled Analysis of Atmospheric Mass Loss and Tidal Evolution in XUV Irradiated Exoplanets: the TRAPPIST-1 Case Study

Exoplanets residing close to their stars can experience evolution of both their physical structures and their orbits due to the influence of their host stars. In this work, we present a coupled analysis of dynamical tidal dissipation and atmospheric mass loss for exoplanets in XUV irradiated environments. As our primary application, we use this model to study the TRAPPIST-1 system, and place constraints on the interior structure and orbital evolution of the planets. We start by reporting on a UV continuum flux measurement (centered around $\sim1900$ Angstroms) for the star TRAPPIST-1, based on 300 ks of Neil Gehrels Swift Observatory data, and which enables an estimate of the XUV-driven thermal escape arising from XUV photo-dissociation for each planet. We find that the X-ray flaring luminosity, measured from our X-ray detections, of TRAPPIST-1 is 5.6 $\times$10$^{-4} L_{*}$, while the full flux including non-flaring periods is 6.1 $\times$10$^{-5} L_{*}$, when $L_{*}$ is TRAPPIST-1's bolometric luminosity. We then construct a model that includes both atmospheric mass-loss and tidal evolution, and requires the planets to attain their present-day orbital elements during this coupled evolution. We use this model to constrain the ratio $Q'=3Q/2k_{2}$ for each planet. Finally, we use additional numerical models implemented with the Virtual Planet Simulator \texttt{VPLanet} to study ocean retention for these planets using our derived system parameters.

preprint2020arXiv

On The XUV Luminosity Evolution of TRAPPIST-1

We model the long-term XUV luminosity of TRAPPIST-1 to constrain the evolving high-energy radiation environment experienced by its planetary system. Using Markov Chain Monte Carlo (MCMC), we derive probabilistic constraints for TRAPPIST-1's stellar and XUV evolution that account for observational uncertainties, degeneracies between model parameters, and empirical data of low-mass stars. We constrain TRAPPIST-1's mass to $m_{\star} = 0.089 \pm{0.001}$ M$_{\odot}$ and find that its early XUV luminosity likely saturated at $\log_{10}(L_{XUV}/L_{bol}) = -3.03^{+0.23}_{-0.12}$. From the posterior distribution, we infer that there is a ${\sim}40\%$ chance that TRAPPIST-1 is still in the saturated phase today, suggesting that TRAPPIST-1 has maintained high activity and $L_{XUV}/L_{bol} \approx 10^{-3}$ for several Gyrs. TRAPPIST-1's planetary system therefore likely experienced a persistent and extreme XUV flux environment, potentially driving significant atmospheric erosion and volatile loss. The inner planets likely received XUV fluxes ${\sim}10^3 - 10^4\times$ that of the modern Earth during TRAPPIST-1's ${\sim}1$ Gyr-long pre-main sequence phase. Deriving these constraints via MCMC is computationally non-trivial, so scaling our methods to constrain the XUV evolution of a larger number of M dwarfs that harbor terrestrial exoplanets would incur significant computational expenses. We demonstrate that approxposterior, an open source Python machine learning package for approximate Bayesian inference using Gaussian processes, accurately and efficiently replicates our analysis using $980\times$ less computational time and $1330\times$ fewer simulations than MCMC sampling using emcee. We find that approxposterior derives constraints with mean errors on the best fit values and $1σ$ uncertainties of $0.61\%$ and $5.5\%$, respectively, relative to emcee.

preprint2019arXiv

VPLanet: The Virtual Planet Simulator

We describe a software package called VPLanet that simulates fundamental aspects of planetary system evolution over Gyr timescales, with a focus on investigating habitable worlds. In this initial release, eleven physics modules are included that model internal, atmospheric, rotational, orbital, stellar, and galactic processes. Many of these modules can be coupled simultaneously to simulate the evolution of terrestrial planets, gaseous planets, and stars. The code is validated by reproducing a selection of observations and past results. VPLanet is written in C and designed so that the user can choose the physics modules to apply to an individual object at runtime without recompiling, i.e., a single executable can simulate the diverse phenomena that are relevant to a wide range of planetary and stellar systems. This feature is enabled by matrices and vectors of function pointers that are dynamically allocated and populated based on user input. The speed and modularity of VPLanet enables large parameter sweeps and the versatility to add/remove physical phenomena to assess their importance. VPLanet is publicly available from a repository that contains extensive documentation, numerous examples, Python scripts for plotting and data management, and infrastructure for community input and future development.

preprint2016arXiv

EVEREST: Pixel Level Decorrelation of K2 Light curves

We present EVEREST, an open-source pipeline for removing instrumental noise from K2 light curves. EVEREST employs a variant of pixel level decorrelation (PLD) to remove systematics introduced by the spacecraft's pointing error and a Gaussian process (GP) to capture astrophysical variability. We apply EVEREST to all K2 targets in campaigns 0-7, yielding light curves with precision comparable to that of the original Kepler mission for stars brighter than $K_p \approx 13$, and within a factor of two of the Kepler precision for fainter targets. We perform cross-validation and transit injection and recovery tests to validate the pipeline, and compare our light curves to the other de-trended light curves available for download at the MAST High Level Science Products archive. We find that EVEREST achieves the highest average precision of any of these pipelines for unsaturated K2 stars. The improved precision of these light curves will aid in exoplanet detection and characterization, investigations of stellar variability, asteroseismology, and other photometric studies. The EVEREST pipeline can also easily be applied to future surveys, such as the TESS mission, to correct for instrumental systematics and enable the detection of low signal-to-noise transiting exoplanets. The EVEREST light curves and the source code used to generate them are freely available online.

preprint2016arXiv

Identifying Planetary Biosignature Impostors: Spectral Features of CO and O4 Resulting from Abiotic O2/O3 Production

O2 and O3 have been long considered the most robust individual biosignature gases in a planetary atmosphere, yet multiple mechanisms that may produce them in the absence of life have been described. However, these abiotic planetary mechanisms modify the environment in potentially identifiable ways. Here we briefly discuss two of the most detectable spectral discriminants for abiotic O2/O3: CO and O4. We produce the first explicit self-consistent simulations of these spectral discriminants as they may be seen by JWST. If JWST-NIRISS and/or NIRSpec observe CO (2.35, 4.6 um) in conjunction with CO2 (1.6, 2.0, 4.3 um) in the transmission spectrum of a terrestrial planet it could indicate robust CO2 photolysis and suggest that a future detection of O2 or O3 might not be biogenic. Strong O4 bands seen in transmission at 1.06 and 1.27 um could be diagnostic of a post-runaway O2-dominated atmosphere from massive H-escape. We find that for these false positive scenarios, CO at 2.35 um, CO2 at 2.0 and 4.3 um, and O4 at 1.27 um are all stronger features in transmission than O2/O3 and could be detected with SNRs $\gtrsim$ 3 for an Earth-size planet orbiting a nearby M dwarf star with as few as 10 transits, assuming photon-limited noise. O4 bands could also be sought in UV/VIS/NIR reflected light (at 0.345, 0.36, 0.38, 0.445, 0.475, 0.53, 0.57, 0.63, 1.06, and 1.27 um) by a next generation direct-imaging telescope such as LUVOIR/HDST or HabEx and would indicate an oxygen atmosphere too massive to be biologically produced.

preprint2016arXiv

The Effect of Orbital Configuration on the Possible Climates and Habitability of Kepler-62f

As lower-mass stars often host multiple rocky planets, gravitational interactions among planets can have significant effects on climate and habitability over long timescales. Here we explore a specific case, Kepler-62f, a potentially habitable planet in a five-planet system with a K2V host star. N-body integrations reveal the stable range of initial eccentricities for Kepler-62f is $0.00 \leqslant e \leqslant 0.32$, absent the effect of additional, undetected planets. We simulate the tidal evolution of Kepler-62f in this range and find that, for certain assumptions, the planet can be locked in a synchronous rotation state. Simulations using LMD Generic GCM indicate that with 3 bars of CO$_2$ in its atmosphere, Kepler-62f would only be warm enough for surface liquid water at the upper limit of this eccentricity range, providing it has a high planetary obliquity (between 60$^\circ$ and 90$^\circ$). A climate similar to modern-day Earth is possible for the entire range of stable eccentricities if atmospheric CO$_2$ is increased to 5-bar levels. In a low-CO$_2$ case, simulations with CCSM4 and LMD Generic GCM indicate that increases in planetary obliquity and orbital eccentricity coupled with an orbital configuration that places the summer solstice at or near pericenter permit regions of the planet with above-freezing surface temperatures. This may melt ice sheets formed during colder seasons. If Kepler-62f is synchronously rotating and has an ocean, CO$_2$ levels above 3 bars would be required to distribute enough heat to the night side of the planet to avoid atmospheric freeze-out and permit a large enough region of open water at the planet's substellar point to remain stable. Overall, we find multiple plausible combinations of orbital and atmospheric properties that permit surface liquid water on Kepler-62f.

preprint2015arXiv

Comparative Habitability of Transiting Exoplanets

Exoplanet habitability is traditionally assessed by comparing a planet's semi-major axis to the location of its host star's "habitable zone," the shell around a star for which Earth-like planets can possess liquid surface water. The Kepler space telescope has discovered numerous planet candidates near the habitable zone, and many more are expected from missions such as K2, TESS and PLATO. These candidates often require significant follow-up observations for validation, so prioritizing planets for habitability from transit data has become an important aspect of the search for life in the universe. We propose a method to compare transiting planets for their potential to support life based on transit data, stellar properties and previously reported limits on planetary emitted flux. For a planet in radiative equilibrium, the emitted flux increases with eccentricity, but decreases with albedo. As these parameters are often unconstrained, there is an "eccentricity-albedo degeneracy" for the habitability of transiting exoplanets. Our method mitigates this degeneracy, includes a penalty for large-radius planets, uses terrestrial mass-radius relationships, and, when available, constraints on eccentricity to compute a number we call the "habitability index for transiting exoplanets" that represents the relative probability that an exoplanet could support liquid surface water. We calculate it for Kepler Objects of Interest and find that planets that receive between 60-90% of the Earth's incident radiation, assuming circular orbits, are most likely to be habitable. Finally, we make predictions for the upcoming TESS and JWST missions.

preprint2015arXiv

Habitable Evaporated Cores: Transforming Mini-Neptunes into Super-Earths in the Habitable Zones of M Dwarfs

We show that photoevaporation of small gaseous exoplanets ("mini-Neptunes") in the habitable zones of M dwarfs can remove several Earth masses of hydrogen and helium from these planets and transform them into potentially habitable worlds. We couple X-ray/extreme ultraviolet (XUV)-driven escape, thermal evolution, tidal evolution and orbital migration to explore the types of systems that may harbor such "habitable evaporated cores" (HECs). We find that HECs are most likely to form from planets with $\sim 1 M_\oplus$ solid cores with up to about 50% H/He by mass, though whether or not a given mini-Neptune forms a HEC is highly dependent on the early XUV evolution of the host star. As terrestrial planet formation around M dwarfs by accumulation of local material is likely to form planets that are small and dry, evaporation of small migrating mini-Neptunes could be one of the dominant formation mechanisms for volatile-rich Earths around these stars.

preprint2015arXiv

Long-lived Chaotic Orbital Evolution of Exoplanets in Mean Motion Resonances with Mutual Inclinations

We present N-body simulations of resonant planets with inclined orbits that show chaotically evolving eccentricities and inclinations that can persist for at least 10 Gyr. A wide range of behavior is possible, from fast, low amplitude variations to systems in which eccentricities reach 0.9999 and inclinations 179.9 degrees. While the orbital elements evolve chaotically, at least one resonant argument always librates. We show that the HD 73526, HD 45364 and HD 60532 systems may be in chaotically-evolving resonances. Chaotic evolution is apparent in the 2:1, 3:1 and 3:2 resonances, and for planetary masses from lunar- to Jupiter-mass. In some cases, orbital disruption occurs after several Gyr, implying the mechanism is not rigorously stable, just long-lived relative to the main sequence lifetimes of solar-type stars. Planet-planet scattering appears to yield planets in inclined resonances that evolve chaotically in about 0.5% of cases. These results suggest that 1) approximate methods for identifying unstable orbital architectures may have limited applicability, 2) the observed close-in exoplanets may be produced during the high eccentricity phases induced by inclined resonances, 3) those exoplanets' orbital planes may be misaligned with the host star's spin axis, 4) systems with resonances may be systematically younger than those without, 5) the distribution of period ratios of adjacent planets detected via transit may be skewed due to inclined resonances, and 6) potentially habitable planets in resonances may have dramatically different climatic evolution than the Earth. The GAIA spacecraft is capable of discovering giant planets in these types of orbits.

preprint2015arXiv

Spin-Orbit Misalignment of Two-Planet-System KOI-89 Via Gravity Darkening

We constrain the true spin-orbit alignment of the KOI-89 system by numerically fitting the two \emph{Kepler} photometric lightcurves produced by transiting planets KOI-89.01 and KOI-89.02. The two planets have periods of 84.69 days and 207.58 days, respectively. We find that the two bodies are low-density giant planets with radii $0.45 \pm 0.03~\mathrm{R_{jup}}$ and $0.43 \pm 0.05~\mathrm{R_{jup}}$ and spin-orbit misalignments $72^{\circ} \pm 3^\circ$ and $73^{\circ+11}_{-5}$, respectively. Via dynamic stability tests we demonstrate the general trend of higher system stability with the two planets close to mutual alignment and estimate their coalignment angle to $20^\circ \pm 20^\circ$ -- i.e. the planets are misaligned with the star but may be aligned with each other. From these results, we limit KOI-89's misalignment mechanisms to star-disk-binary interactions, disk warping via planet-disk interactions, planet-planet scattering, Kozai resonance, or internal gravity waves.

preprint2015arXiv

Tidal heating of Earth-like exoplanets around M stars: Thermal, magnetic, and orbital evolutions

The internal thermal and magnetic evolution of rocky exoplanets is critical to their habitability. We focus on the thermal-orbital evolution of Earth-mass planets around low mass M stars whose radiative habitable zone overlaps with the "tidal zone". We develop a thermal-orbital evolution model calibrated to Earth that couples tidal dissipation, with a temperature-dependent Maxwell rheology, to orbital circularization and migration. We illustrate thermal-orbital steady states where surface heat flow is balanced by tidal dissipation and cooling can be stalled for billions of years until circularization occurs. Orbital energy dissipated as tidal heat in the interior drives both inward migration and circularization, with a circularization time that is inversely proportional to the dissipation rate. We identify a peak in the internal dissipation rate as the mantle passes through a visco-elastic state at mantle temperatures near 1800 K. Planets orbiting a 0.1 solar-mass star within $0.07$ AU circularize before 10 Gyr, independent of initial eccentricity. Once circular, these planets cool monotonically and maintain dynamos similar to Earth. Planets forced into eccentric orbits can experience a super-cooling of the core and rapid core solidification, inhibiting dynamo action for planets in the habitable zone. We find that tidal heating is insignificant in the habitable zone around $0.45$ (or larger) solar mass stars because tidal dissipation is a stronger function of orbital distance than stellar mass, and the habitable zone is further from larger stars. Suppression of the planetary magnetic field exposes the atmosphere to stellar wind erosion and the surface to harmful radiation. In addition to weak magnetic fields, massive melt eruption rates and prolonged magma oceans may render eccentric planets in the habitable zone of low mass stars inhospitable for life.

preprint2014arXiv

Extreme Water Loss and Abiotic O$_2$ Buildup On Planets Throughout the Habitable Zones of M Dwarfs

We show that terrestrial planets in the habitable zones of M dwarfs older than $\sim$ 1 Gyr could have been in runaway greenhouses for several hundred Myr following their formation due to the star's extended pre-main sequence phase, provided they form with abundant surface water. Such prolonged runaway greenhouses can lead to planetary evolution divergent from that of Earth. During this early runaway phase, photolysis of water vapor and hydrogen/oxygen escape to space can lead to the loss of several Earth oceans of water from planets throughout the habitable zone, regardless of whether the escape is energy-limited or diffusion-limited. We find that the amount of water lost scales with the planet mass, since the diffusion-limited hydrogen escape flux is proportional to the planet surface gravity. In addition to undergoing potential desiccation, planets with inefficient oxygen sinks at the surface may build up hundreds to thousands of bars of abiotically produced O$_2$, resulting in potential false positives for life. The amount of O$_2$ that builds up also scales with the planet mass; we find that O$_2$ builds up at a constant rate that is controlled by diffusion: $\sim$ 5 bars/Myr on Earth-mass planets and up to $\sim$ 25 bars/Myr on super-Earths. As a result, some recently discovered super-Earths in the habitable zone such as GJ 667Cc could have built up as many as 2000 bars of O$_2$ due to the loss of up to 10 Earth oceans of water. The fate of a given planet strongly depends on the extreme ultraviolet flux, the duration of the runaway regime, the initial water content, and the rate at which oxygen is absorbed by the surface. In general, we find that the initial phase of high luminosity may compromise the habitability of many terrestrial planets orbiting low mass stars.

preprint2014arXiv

The 3-dimensional architecture of the Upsilon Andromedae planetary system

The Upsilon Andromedae system is the first exoplanetary system to have the relative inclination of two planets' orbital planes directly measured, and therefore offers our first window into the 3-dimensional configurations of planetary systems. We present, for the first time, full 3-dimensional, dynamically stable configurations for the 3 planets of the system consistent with all observational constraints. While the outer 2 planets, c and d, are inclined by about 30 degrees, the inner planet's orbital plane has not been detected. We use N-body simulations to search for stable 3-planet configurations that are consistent with the combined radial velocity and astrometric solution. We find that only 10 trials out of 1000 are robustly stable on 100 Myr timescales, or about 8 billion orbits of planet b. Planet b's orbit must lie near the invariable plane of planets c and d, but can be either prograde or retrograde. These solutions predict b's mass is in the range 2 - 9 $M_{Jup}$ and has an inclination angle from the sky plane of less than 25 degrees. Combined with brightness variations in the combined star/planet light curve ("phase curve"), our results imply that planet b's radius is about 1.8 $R_{Jup}$, relatively large for a planet of its age. However, the eccentricity of b in several of our stable solutions reaches values greater than 0.1, generating upwards of $10^{19}$ watts in the interior of the planet via tidal dissipation, possibly inflating the radius to an amount consistent with phase curve observations.

preprint2014arXiv

Tides, planetary companions, and habitability: Habitability in the habitable zone of low-mass stars

Earth-scale planets in the classical habitable zone (HZ) are more likely to be habitable if they possess active geophysics. Without a constant internal energy source, planets cool as they age, eventually terminating tectonic activity and rendering the planet sterile to life. However, for planets orbiting low-mass stars, the presence of an outer companion could generate enough tidal heat in the HZ planet to prevent such cooling. The range of mass and orbital parameters for the companion that give adequate long-term heating of the inner HZ planet, while avoiding very early total desiccation, is probably substantial. We locate the ideal location for the outer of a pair of planets, under the assumption that the inner planet has the same incident flux as Earth, orbiting example stars: a generic late M dwarf ($T_{eff}=2670 K$) and the M9V/L0 dwarf DEN1048. Thus discoveries of Earth-scale planets in the HZ zone of old small stars should be followed by searches for outer companion planets that might be essential for current habitability.

preprint2013arXiv

A dynamically-packed planetary system around GJ667C with three super-Earths in its habitable zone

Since low-mass stars have low luminosities, orbits at which liquid water can exist on Earth-sized planets are relatively close-in, which produces Doppler signals that are detectable using state-of-the-art Doppler spectroscopy. GJ 667C is already known to be orbited by two super-Earth candidates. We investigate whether the data supports the presence of additional companions. We obtain new Doppler measurements from HARPS extracted spectra and combined them with those obtained from the PFS and HIRES spectrographs. We used Bayesian and periodogram-based methods to re-assess the number of candidates and evaluated the confidence of each detection. Among other tests, we validated the planet candidates by analyzing correlations of each Doppler signal activity indices and investigate quasi-periodicity. Doppler measurements of GJ 667C are described better by six Keplerian-like signals: the two known candidates (b and c); three additional few-Earth mass candidates with periods of 92, 62 and 39 days (d, e and f); a cold super-Earth in a 260-day orbit (g) and tantalizing evidence of a $\sim$ 1 M$_\oplus$ object in a close-in orbit of 17 days (h). We explore whether long-term stable orbits are compatible with the data by integrating 8$\times 10^4$ solutions derived from the Bayesian samplings. The system consisting of six planets is compatible with dynamically stable configurations. As for the solar system, the most stable solutions do not contain mean-motion resonances and are described well by analytic Laplace-Lagrange solutions. The presence of a seventh planet (h) is supported by the fact that it appears squarely centered on the only island of stability left in the six-planet solution. Habitability assessments accounting for the stellar flux, as well as tidal dissipation effects, indicate that three (maybe four) planets are potentially habitable...

preprint2013arXiv

A Method to Identify the Boundary Between Rocky and Gaseous Exoplanets from Tidal Theory and Transit Durations

The determination of an exoplanet as rocky is critical for the assessment of planetary habitability. Observationally, the number of small-radius, transiting planets with accompanying mass measurements is insufficient for a robust determination of the transitional mass or radius. Theoretically, models predict that rocky planets can grow large enough to become gas giants when they reach ~10 Earth-masses, but the transitional mass remains unknown. Here I show how transit data, interpreted in the context of tidal theory, can reveal the critical radius that separates rocky and gaseous exoplanets. Standard tidal models predict that rocky exoplanets' orbits are tidally circularized much more rapidly than gaseous bodies', suggesting the former will tend to be found on circular orbits at larger semi-major axes than the latter. Well-sampled transits can provide a minimum eccentricity of the orbit, allowing a measurement of this differential circularization. I show that this effect should be present in the data from the Kepler spacecraft, but is not apparent. Instead, it appears that there is no evidence of tidal circularization at any planetary radius, probably because the publicly-available data, particularly the impact parameters, are not accurate enough. I also review the bias in the transit duration toward values that are smaller than that of planets on circular orbits, stressing that the azimuthal velocity of the planet determines the transit duration. The ensemble of Kepler planet candidates may be able to determine the critical radius between rocky and gaseous exoplanets, tidal dissipation as a function of planetary radius, and discriminate between tidal models.

preprint2013arXiv

APOSTLE: Longterm Transit Monitoring and Stability Analysis of XO-2b

The Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) observed 10 transits of XO-2b over a period of three years. We present measurements which confirm previous estimates of system parameters like the normalized semi-major axis (a/R_{*}), stellar density (ρ_{*}), impact parameter (b) and orbital inclination (i_{orb}). Our errors on system parameters like a/R_{*} and ρ_{*} have improved by ~40% compared to previous best ground-based measurements. Our study of the transit times show no evidence for transit timing variations and we are able to rule out co-planar companions with masses \ge 0.20 \mearth\ in low order mean motion resonance with XO-2b. We also explored the stability of the XO-2 system given various orbital configurations of a hypothetical planet near the 2:1 mean motion resonance. We find that a wide range of orbits (including Earth-mass perturbers) are both dynamically stable and produce observable TTVs. We find that up to 51% of our stable simulations show TTVs that are smaller than the typical transit timing errors (~20 sec) measured for XO-2b, and hence remain undetectable.

preprint2013arXiv

Exomoon habitability constrained by illumination and tidal heating

The detection of moons orbiting extrasolar planets ("exomoons") has now become feasible. Once they are discovered in the circumstellar habitable zone, questions about their habitability will emerge. Exomoons are likely to be tidally locked to their planet and hence experience days much shorter than their orbital period around the star and have seasons, all of which works in favor of habitability. These satellites can receive more illumination per area than their host planets, as the planet reflects stellar light and emits thermal photons. On the contrary, eclipses can significantly alter local climates on exomoons by reducing stellar illumination. In addition to radiative heating, tidal heating can be very large on exomoons, possibly even large enough for sterilization. We identify combinations of physical and orbital parameters for which radiative and tidal heating are strong enough to trigger a runaway greenhouse. By analogy with the circumstellar habitable zone, these constraints define a circumplanetary "habitable edge". We apply our model to hypothetical moons around the recently discovered exoplanet Kepler-22b and the giant planet candidate KOI211.01 and describe, for the first time, the orbits of habitable exomoons. If either planet hosted a satellite at a distance greater than 10 planetary radii, then this could indicate the presence of a habitable moon.

preprint2013arXiv

Hot Moons and Cool Stars

The exquisite photometric precision of the Kepler space telescope now puts the detection of extrasolar moons at the horizon. Here, we firstly review observational and analytical techniques that have recently been proposed to find exomoons. Secondly, we discuss the prospects of characterizing potentially habitable extrasolar satellites. With moons being much more numerous than planets in the solar system and with most exoplanets found in the stellar habitable zone being gas giants, habitable moons could be as abundant as habitable planets. However, satellites orbiting planets in the habitable zones of cool stars will encounter strong tidal heating and likely appear as hot moons.

preprint2013arXiv

Observations of the WASP-2 System by the APOSTLE Program

We present transit observations of the WASP-2 exoplanet system by the Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) program. Model fitting to these data allows us to improve measurements of the hot-Jupiter exoplanet WASP-2b and its orbital parameters by a factor of ~2 over prior studies; we do not find evidence for transit depth variations. We do find reduced chi^2 values greater than 1.0 in the observed minus computed transit times. A sinusoidal fit to the residuals yields a timing semi-amplitude of 32 seconds and a period of 389 days. However, random rearrangements of the data provide similar quality fits, and we cannot with certainty ascribe the timing variations to mutual exoplanet interactions. This inconclusive result is consistent with the lack of incontrovertible transit timing variations (TTVs) observed in other hot-Jupiter systems. This outcome emphasizes that unique recognition of TTVs requires dense sampling of the libration cycle (e.g. continuous observations from space-based platforms). However, even in systems observed with the Kepler spacecraft, there is a noted lack of transiting companions and TTVs in hot-Jupiter systems. This result is more meaningful, and indicates that hot-Jupiter systems, while they are easily observable from the ground, do not appear to be currently configured in a manner favorable to the detection of TTVs. The future of ground-based TTV studies may reside in resolving secular trends, and/or implementation at extreme quality observing sites to minimize atmospheric red noise.

preprint2013arXiv

Runaway greenhouse effect on exomoons due to irradiation from hot, young giant planets

The Kepler space telescope has detected transits of objects as small as the Earth's Moon, and moons as small as 0.2 Earth masses can be detected in the Kepler data by transit timing and transit duration variations of their host planets. Such massive moons around giant planets in the stellar habitable zone (HZ) could serve as habitats for extraterrestrial life. We here assess the danger of exomoons to be in a runaway greenhouse (RG) state due to extensive heating from the planet. We apply pre-computed evolution tracks for giant planets to calculate the incident planetary radiation on the moon as a function of time. The total energy budget of stellar flux, illumination from the planet, and tidal heating in the satellite is compared to the critical flux for the moon to experience an RG effect. Irradiation from a 13-Jupiter-mass planet onto an Earth-sized moon at a distance of ten Jupiter radii can drive an RG state on the moon for about 200 Myr. If stellar illumination equivalent to that received by Earth from the Sun is added, the RG holds for about 500 Myr. After 1000 Myr, the planet's habitable edge has moved inward to about 6 planetary radii. Exomoons in orbits with eccentricities of 0.1 experience strong tidal heating; they must orbit a 13-Jupiter-mass host beyond 29 or 18 Jupiter radii after 100 Myr (at the inner and outer boundaries of the stellar HZ, respectively), and beyond 13 Jupiter radii (in both cases) after 1000 Myr to be habitable. If a roughly Earth-sized moon would be detected in orbit around a giant planet, and if this binary would orbit in the stellar HZ, then it will be crucial to recover the moon's orbital history. If such a moon around a 13-Jupiter-mass planet has been closer than 20 Jupiter radii to its host during the first few hundred million years at least, then it might have lost substantial amounts of its initial water reservoir and be uninhabitable today.

preprint2013arXiv

Secular Behavior of Exoplanets: Self-Consistency and Comparisons with the Planet-Planet Scattering Hypothesis

If mutual gravitational scattering among exoplanets occurs, then it may produce unique orbital properties. For example, two-planet systems that lie near the boundary between circulation and libration of their periapses could result if planet-planet scattering ejected a former third planet quickly, leaving one planet on an eccentric orbit and the other on a circular orbit. We first improve upon previous work that examined the apsidal behavior of known multiplanet systems by doubling the sample size and including observational uncertainties. This analysis recovers previous results that demonstrated that many systems lay on the apsidal boundary between libration and circulation. We then performed over 12,000 three-dimensional N-body simulations of hypothetical three-body systems that are unstable, but stabilize to two-body systems after an ejection. Using these synthetic two-planet systems, we test the planet-planet scattering hypothesis by comparing their apsidal behavior, over a range of viewing angles, to that of the observed systems and find that they are statistically consistent regardless of the multiplicity of the observed systems. Finally, we combine our results with previous studies to show that, from the sampled cases, the most likely planetary mass function prior to planet-planet scattering follows a power law with index -1.1. We find that this pre-scattering mass function predicts a mutual inclination frequency distribution that follows an exponential function with an index between -0.06 and -0.1.

preprint2013arXiv

Spin-driven tidal pumping: Tidally driven changes in planetary spin coupled with secular interactions between planets

In a multiplanet system, tides acting on the inner planet can significantly affect the orbital evolution of the entire system. While tides usually damp eccentricities, a novel mechanism identified by Correia et al. (2012) tends to raise eccentricities as a result of the tides' effect on the inner planet's rotation. Our analytical description of this spin-driven tidal (SDT) effect shows that, while the inner planet's eccentricity undergoes pumping, the process is more completely described by an exchange of strength between the two eigenmodes of the dynamical system. Our analysis allows derivation of criteria for two-planet coplanar systems where the SDT effect can reverse tidal damping, and may preclude the effect's being significant for realistic systems. For the specific case quantified by Correia et al., the effect is strong because of the large adopted tidal time lag, which may not be appropriate for the assumed Saturn-like inner planet. On the other hand, the effective Q for any given planet in exotic circumstances is very uncertain, so the SDT effect could play a role in planetary evolution.

preprint2013arXiv

Very Low Mass Stellar and Substellar Companions to Solar-Like Stars From MARVELS V: A Low Eccentricity Brown Dwarf from the Driest Part of the Desert, MARVELS-6b

We describe the discovery of a likely brown dwarf (BD) companion with a minimum mass of 31.7 +/- 2.0 M_Jup to GSC 03546-01452 from the MARVELS radial velocity survey, which we designate as MARVELS-6b. For reasonable priors, our analysis gives a probability of 72% that MARVELS-6b has a mass below the hydrogen-burning limit of 0.072 M_Sun, and thus it is a high-confidence BD companion. It has a moderately long orbital period of 47.8929 +0.0063/-0.0062 days with a low eccentricty of 0.1442 +0.0078/-0.0073, and a semi-amplitude of 1644 +12/-13 m/s. Moderate resolution spectroscopy of the host star has determined the following parameters: T_eff = 5598 +/- 63, log g = 4.44 +/- 0.17, and [Fe/H] = +0.40 +/- 0.09. Based upon these measurements, GSC 03546-01452 has a probable mass and radius of M_star = 1.11 +/- 0.11 M_Sun and R_star = 1.06 +/- 0.23 R_Sun with an age consistent with less than ~6 Gyr at a distance of 219 +/- 21 pc from the Sun. Although MARVELS-6b is not observed to transit, we cannot definitively rule out a transiting configuration based on our observations. There is a visual companion detected with Lucky Imaging at 7.7 arcsec from the host star, but our analysis shows that it is not bound to this system. The minimum mass of MARVELS-6b exists at the minimum of the mass functions for both stars and planets, making this a rare object even compared to other BDs.

preprint2012arXiv

APOSTLE: Eleven Transit Observations of TrES-3b

The Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) observed eleven transits of TrES-3b over two years in order to constrain system parameters and look for transit timing and depth variations. We describe an updated analysis protocol for APOSTLE data, including the reduction pipeline, transit model and Markov Chain Monte Carlo analyzer. Our estimates of the system parameters for TrESb are consistent with previous estimates to within the 2σ confidence level. We improved the errors (by 10--30%) on system parameters like the orbital inclination ($i_{\text{orb}}$), impact parameter (b) and stellar density (ρ$_{\star}$) compared to previous measurements. The near-grazing nature of the system, and incomplete sampling of some transits, limited our ability to place reliable uncertainties on individual transit depths and hence we do not report strong evidence for variability. Our analysis of the transit timing data show no evidence for transit timing variations and our timing measurements are able to rule out Super-Earth and Gas Giant companions in low order mean motion resonance with TrES-3b.

preprint2012arXiv

Constraints on the habitability of extrasolar moons

Detections of massive extrasolar moons are shown feasible with the Kepler space telescope. Kepler's findings of about 50 exoplanets in the stellar habitable zone naturally make us wonder about the habitability of their hypothetical moons. Illumination from the planet, eclipses, tidal heating, and tidal locking distinguish remote characterization of exomoons from that of exoplanets. We show how evaluation of an exomoon's habitability is possible based on the parameters accessible by current and near-future technology.

preprint2012arXiv

Habitable Planets Around White and Brown Dwarfs: The Perils of a Cooling Primary

White and brown dwarfs are astrophysical objects that are bright enough to support an insolation habitable zone (IHZ). Unlike hydrogen-burning stars, they cool and become less luminous with time, and hence their IHZ moves in with time. The inner edge of the IHZ is defined as the orbital radius at which a planet may enter a moist or runaway greenhouse, phenomena that can remove a planet's surface water forever. Thus, as the IHZ moves in, planets that enter it may no longer have any water, and are still uninhabitable. Additionally, the close proximity of the IHZ to the primary leads to concern that tidal heating may also be strong enough to trigger a runaway greenhouse, even for orbital eccentricities as small as 10^-6. Water loss occurs due to photolyzation by UV photons in the planetary stratosphere, followed by hydrogen escape. Young white dwarfs emit a large amount of these photons as their surface temperatures are over 10^4 K. The situation is less clear for brown dwarfs, as observational data do not constrain their early activity and UV emission very well. Nonetheless, both types of planets are at risk of never achieving habitable conditions, but planets orbiting white dwarfs may be less likely to sustain life than those orbiting brown dwarfs. We consider the future habitability of the planet candidates KOI 55.01 and 55.02 in these terms and find they are unlikely to become habitable.

preprint2012arXiv

The Ninth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Baryon Oscillation Spectroscopic Survey

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

preprint2012arXiv

Tidal Venuses: Triggering a Climate Catastrophe via Tidal Heating

Traditionally stellar radiation has been the only heat source considered capable of determining global climate on long timescales. Here we show that terrestrial exoplanets orbiting low-mass stars may be tidally heated at high enough levels to induce a runaway greenhouse for a long enough duration for all the hydrogen to escape. Without hydrogen, the planet no longer has water and cannot support life. We call these planets "Tidal Venuses," and the phenomenon a "tidal greenhouse." Tidal effects also circularize the orbit, which decreases tidal heating. Hence, some planets may form with large eccentricity, with its accompanying large tidal heating, and lose their water, but eventually settle into nearly circular orbits (i.e. with negligible tidal heating) in the habitable zone (HZ). However, these planets are not habitable as past tidal heating desiccated them, and hence should not be ranked highly for detailed follow-up observations aimed at detecting biosignatures. Planets orbiting stars with masses <0.3 solar masses may be in danger of desiccation via tidal heating. We apply these concepts to Gl 667C c, a ~4.5 Earth-mass planet orbiting a 0.3 solar mass star at 0.12 AU. We find that it probably did not lose its water via tidal heating as orbital stability is unlikely for the high eccentricities required for the tidal greenhouse. As the inner edge of the HZ is defined by the onset of a runaway or moist greenhouse powered by radiation, our results represent a fundamental revision to the HZ for non-circular orbits. In the appendices we review a) the moist and runaway greenhouses, b) hydrogen escape, c) stellar mass-radius and mass-luminosity relations, d) terrestrial planet mass-radius relations, and e) linear tidal theories. [abridged]

preprint2012arXiv

Very Low-Mass Stellar and Substellar Companions to Solar-Like Stars from MARVELS I: A Low Mass Ratio Stellar Companion to TYC 4110-01037-1 in a 79-day Orbit

TYC 4110-01037-1 has a low-mass stellar companion, whose small mass ratio and short orbital period are atypical amongst solar-like (Teff ~< 6000 K) binary systems. Our analysis of TYC 4110-01037-1 reveals it to be a moderately aged (~<5 Gyr) solar-like star having a mass of 1.07 +/- 0.08 MSun and radius of 0.99 +/- 0.18 RSun. We analyze 32 radial velocity measurements from the SDSS-III MARVELS survey as well as 6 supporting radial velocity measurements from the SARG spectrograph on the 3.6m TNG telescope obtained over a period of ~2 years. The best Keplerian orbital fit parameters were found to have a period of 78.994 +/- 0.012 days, an eccentricity of 0.1095 +/- 0.0023, and a semi-amplitude of 4199 +/- 11 m/s. We determine the minimum companion mass (if sin i = 1) to be 97.7 +/- 5.8 MJup. The system's companion to host star mass ratio, >0.087 +/- 0.003, places it at the lowest end of observed values for short period stellar companions to solar-like (Teff ~< 6000 K) stars. One possible way to create such a system would be if a triple-component stellar multiple broke up into a short period, low q binary during the cluster dispersal phase of its lifetime. A candidate tertiary body has been identified in the system via single-epoch, high contrast imagery. If this object is confirmed to be co-moving, we estimate it would be a dM4 star. We present these results in the context of our larger-scale effort to constrain the statistics of low mass stellar and brown dwarf companions to FGK-type stars via the MARVELS survey.

preprint2012arXiv

Very Low-mass Stellar and Substellar Companions to Solar-like Stars from MARVELS II: A Short-period Companion Orbiting an F Star with Evidence of a Stellar Tertiary And Significant Mutual Inclination

We report the discovery via radial velocity of a short-period (P = 2.430420 \pm 0.000006 days) companion to the F-type main sequence star TYC 2930-00872-1. A long-term trend in the radial velocities indicates the presence of a tertiary stellar companion with $P > 2000$ days. High-resolution spectroscopy of the host star yields T_eff = 6427 +/- 33 K, log(g) = 4.52 +/- 0.14, and [Fe/H]=-0.04 +/- 0.05. These parameters, combined with the broad-band spectral energy distribution and parallax, allow us to infer a mass and radius of the host star of M_1=1.21 +/- 0.08 M_\odot and R_1=1.09_{-0.13}^{+0.15} R_\odot. We are able to exclude transits of the inner companion with high confidence. The host star's spectrum exhibits clear Ca H and K core emission indicating stellar activity, but a lack of photometric variability and small v*sin(I) suggest the primary's spin axis is oriented in a pole-on configuration. The rotational period of the primary from an activity-rotation relation matches the orbital period of the inner companion to within 1.5 σ, suggesting they are tidally locked. If the inner companion's orbital angular momentum vector is aligned with the stellar spin axis, as expected through tidal evolution, then it has a stellar mass of M_2 ~ 0.3-0.4 M_\odot. Direct imaging limits the existence of stellar companions to projected separations < 30 AU. No set of spectral lines and no significant flux contribution to the spectral energy distribution from either companion are detected, which places individual upper mass limits of M < 1.0 M_\odot, provided they are not stellar remnants. If the tertiary is not a stellar remnant, then it likely has a mass of ~0.5-0.6 M_\odot, and its orbit is likely significantly inclined from that of the secondary, suggesting that the Kozai-Lidov mechanism may have driven the dynamical evolution of this system.

preprint2012arXiv

Very Low-mass Stellar and Substellar Companions to Solar-like Stars from Marvels III: A Short-Period Brown Dwarf Candidate Around An Active G0Iv Subgiant

We present an eccentric, short-period brown dwarf candidate orbiting the active, slightly evolved subgiant star TYC 2087-00255-1, which has effective temperature T_eff = 5903+/-42 K, surface gravity log (g) = 4.07+/-0.16 (cgs), and metallicity [Fe/H] = -0.23+/-0.07. This candidate was discovered using data from the first two years of the Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS), which is part of the third phase of Sloan Digital Sky Survey. From our 38 radial velocity measurements spread over a two-year time baseline, we derive a Keplerian orbital fit with semi-amplitude K=3.571+/-0.041 km/s, period P=9.0090+/-0.0004 days, and eccentricity e=0.226+/-0.011. Adopting a mass of 1.16+/-0.11 Msun for the subgiant host star, we infer that the companion has a minimum mass of 40.0+/-2.5 M_Jup. Assuming an edge-on orbit, the semimajor axis is 0.090+/-0.003 AU. The host star is photometrically variable at the \sim1% level with a period of \sim13.16+/-0.01 days, indicating that the host star spin and companion orbit are not synchronized. Through adaptive optics imaging we also found a point source 643+/-10 mas away from TYC 2087-00255-1, which would have a mass of 0.13 Msun if it is physically associated with TYC 2087-00255-1 and has the same age. Future proper motion observation should be able to resolve if this tertiary object is physically associated with TYC 2087-00255-1 and make TYC 2087-00255-1 a triple body system. Core Ca II H and K line emission indicate that the host is chromospherically active, at a level that is consistent with the inferred spin period and measured v_{rot}*sin i, but unusual for a subgiant of this T_eff. This activity could be explained by ongoing tidal spin-up of the host star by the companion.

preprint2011arXiv

Habitability of extrasolar planets and tidal spin evolution

Stellar radiation has conservatively been used as the key constraint to planetary habitability. We review here the effects of tides, exerted by the host star on the planet, on the evolution of the planetary spin. Tides initially drive the rotation period and the orientation of the rotation axis into an equilibrium state but do not necessarily lead to synchronous rotation. As tides also circularize the orbit, eventually the rotation period does equal the orbital period and one hemisphere will be permanently irradiated by the star. Furthermore, the rotational axis will become perpendicular to the orbit, i.e. the planetary surface will not experience seasonal variations of the insolation. We illustrate here how tides alter the spins of planets in the traditional habitable zone. As an example, we show that, neglecting perturbations due to other companions, the Super-Earth Gl581d performs two rotations per orbit and that any primordial obliquity has been eroded.

preprint2011arXiv

Tidal obliquity evolution of potentially habitable planets

Stellar insolation has been used as the main constraint on a planet's habitability. However, as more Earth-like planets are discovered around low-mass stars (LMSs), a re-examination of the role of tides on the habitability of exoplanets has begun. Those studies have yet to consider the misalignment between a planet's rotational axis and the orbital plane normal, i.e. the planetary obliquity. We apply two equilibrium tide theories to compute the obliquity evolution of terrestrial planets orbiting in the habitable zones around LMSs. The time for the obliquity to decrease from an Earth-like obliquity of 23.5 deg to 5 deg, the 'tilt erosion time', is compared to the traditional insolation habitable zone (IHZ) as a function of semi-major axis, eccentricity, and stellar mass. We also compute tidal heating and equilibrium rotation caused by obliquity tides. The Super-Earth Gl581d and the planet candidate Gl581g are studied as examples for tidal processes. Earth-like obliquities of terrestrial planets in the IHZ around <0.25 solar mass (M_sun) stars are eroded in <0.1 Gyr. Only terrestrial planets orbiting stars with masses >0.9 M_sun experience tilt erosion times larger than 1 Gyr throughout the IHZ. Terrestrial planets in the IHZ of stars with masses <0.25 M_sun undergo significant tidal heating due to obliquity tides. The predictions of the two tidal models diverge significantly for e>0.3. In our two-body simulations, Gl581d's obliquity is eroded to 0 and its rotation period reached its equilibrium state of half its orbital period in <0.1 Gyr. Tidal surface heating on the putative Gl581g is <150 mW/m^2 as long as its eccentricity is <0.3. Obliquity tides modify the concept of the habitable zone. Tilt erosion of terrestrial planets orbiting LMSs should be included by atmospheric modelers. Tidal heating needs to be considered by geologists.

preprint2010arXiv

Discovery of a Low-Mass Companion to a Metal-Rich F Star with the MARVELS Pilot Project

We report the discovery of a low-mass companion orbiting the metal-rich, main sequence F star TYC 2949-00557-1 during the MARVELS (Multi-object APO Radial Velocity Exoplanet Large-area Survey) Pilot Project. The host star has an effective temperature T_eff = 6135 +/- 40 K, log(g) = 4.4 +/- 0.1 and [Fe/H] = 0.32 +/- 0.01, indicating a mass of M = 1.25 +/- 0.09 M_\odot and R = 1.15 +/- 0.15 R_\odot. The companion has an orbital period of 5.69449 +/- 0.00023 days and straddles the hydrogen burning limit with a minimum mass of 64 M_J, and may thus be an example of the rare class of brown dwarfs orbiting at distances comparable to those of "Hot Jupiters." We present relative photometry that demonstrates the host star is photometrically stable at the few millimagnitude level on time scales of hours to years, and rules out transits for a companion of radius greater than 0.8 R_J at the 95% confidence level. Tidal analysis of the system suggests that the star and companion are likely in a double synchronous state where both rotational and orbital synchronization have been achieved. This is the first low-mass companion detected with a multi-object, dispersed, fixed-delay interferometer.

preprint2010arXiv

Habitability of Planets Orbiting Cool Stars

Terrestrial planets are more likely to be detected if they orbit M dwarfs due to the favorable planet/star size and mass ratios. However, M dwarf habitable zones are significantly closer to the star than the one around our Sun, which leads to different requirements for planetary habitability and its detection. We review 1) the current limits to detection, 2) the role of M dwarf spectral energy distributions on atmospheric chemistry, 3) tidal effects, stressing that tidal locking is not synonymous with synchronous rotation, 4) the role of atmospheric mass loss and propose that some habitable worlds may be the volatile-rich, evaporated cores of giant planets, and 5) the role of planetary rotation and magnetic field generation, emphasizing that slow rotation does not preclude strong magnetic fields and their shielding of the surface from stellar activity. Finally we present preliminary findings of the NASA Astrobiology Institute's workshop "Revisiting the Habitable Zone." We assess the recently-announced planet Gl 581 g and find no obvious barriers to habitability. We conclude that no known phenomenon completely precludes the habitability of terrestrial planets orbiting cool stars.

preprint2010arXiv

MARVELS-1b: A Short-Period, Brown Dwarf Desert Candidate from the SDSS-III MARVELS Planet Search

We present a new short-period brown dwarf candidate around the star TYC 1240-00945-1. This candidate was discovered in the first year of the Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS), which is part of the third phase of the Sloan Digital Sky Survey (SDSS-III), and we designate the brown dwarf as MARVELS-1b. MARVELS uses the technique of dispersed fixed-delay interferometery to simultaneously obtain radial velocity measurements for 60 objects per field using a single, custom-built instrument that is fiber fed from the SDSS 2.5-m telescope. From our 20 radial velocity measurements spread over a ~370 d time baseline, we derive a Keplerian orbital fit with semi-amplitude K=2.533+/-0.025 km/s, period P=5.8953+/-0.0004 d, and eccentricity consistent with circular. Independent follow-up radial velocity data confirm the orbit. Adopting a mass of 1.37+/-0.11 M_Sun for the slightly evolved F9 host star, we infer that the companion has a minimum mass of 28.0+/-1.5 M_Jup, a semimajor axis 0.071+/-0.002 AU assuming an edge-on orbit, and is probably tidally synchronized. We find no evidence for coherent instrinsic variability of the host star at the period of the companion at levels greater than a few millimagnitudes. The companion has an a priori transit probability of ~14%. Although we find no evidence for transits, we cannot definitively rule them out for companion radii ~<1 R_Jup.

preprint2010arXiv

Origin and Dynamics of the Mutually Inclined Orbits of Upsilon Andromedae c and d

We evaluate the orbital evolution and several plausible origins scenarios for the mutually inclined orbits of Upsilon Andromedae c and d. These two planets have orbital elements that oscillate with large amplitudes and lie close to the stability boundary. This configuration, and in particular the observed mutual inclination, demands an explanation. The planetary system may be influenced by a nearby low-mass star, Upsilon And B, which could perturb the planetary orbits, but we find it cannot modify two coplanar orbits into the observed mutual inclination of ~30 deg. However, it could incite ejections or collisions between planetary companions that subsequently raise the mutual inclination to >30 deg. Our simulated systems with large mutual inclinations tend to be further from the stability boundary than Upsilon And, but we are able to produce similar systems. We conclude that scattering is a plausible mechanism to explain the observed orbits of Upsilon And c and d, but we cannot determine whether the scattering was caused by instabilities among the planets themselves or by perturbations from Upsilon And B. We also develop a procedure to quantitatively compare numerous properties of the observed system to our numerical models. Although we only implement this procedure to Upsilon And, it may be applied to any exoplanetary system.

preprint2010arXiv

Stability analysis of single planet systems and their habitable zones

We study the dynamical stability of planetary systems consisting of one hypothetical terrestrial mass planet ($1 $ or $10 \mearth$) and one massive planet ($10 \mearth - 10 \mjup$). We consider masses and orbits that cover the range of observed planetary system architectures (including non-zero initial eccentricities), determine the stability limit through N-body simulations, and compare it to the analytic Hill stability boundary. We show that for given masses and orbits of a two planet system, a single parameter, which can be calculated analytically, describes the Lagrange stability boundary (no ejections or exchanges) but which diverges significantly from the Hill stability boundary. However, we do find that the actual boundary is fractal, and therefore we also identify a second parameter which demarcates the transition from stable to unstable evolution. We show the portions of the habitable zones of $ρ$ CrB, HD 164922, GJ 674, and HD 7924 which can support a terrestrial planet. These analyses clarify the stability boundaries in exoplanetary systems and demonstrate that, for most exoplanetary systems, numerical simulations of the stability of potentially habitable planets are only necessary over a narrow region of parameter space. Finally we also identify and provide a catalog of known systems which can host terrestrial planets in their habitable zones.

preprint2010arXiv

The Mass of HD 38529 c from Hubble Space Telescope Astrometry and High-Precision Radial Velocities

(Abridged) Hubble Space Telescope (HST) Fine Guidance Sensor astrometric observations of the G4 IV star HD 38529 are combined with the results of the analysis of extensive ground-based radial velocity data to determine the mass of the outermost of two previously known companions. Our new radial velocities obtained with the Hobby-Eberly Telescope and velocities from the Carnegie-California group now span over eleven years. With these data we obtain improved RV orbital elements for both the inner companion, HD 38529 b and the outer companion, HD 38529 c. We identify a rotational period of HD 38529 (P_{rot}=31.65 +/- 0.17 d) with FGS photometry. We model the combined astrometric and RV measurements to obtain the parallax, proper motion, perturbation period, perturbation inclination, and perturbation size due to HD 38529 c. For HD 38529 c we find P = 2136.1 +/- 0.3 d, perturbation semi-major axis α=1.05 +/-0.06$ mas, and inclination $i$ = 48.3 deg +/- 4 deg. Assuming a primary mass M_* = 1.48 M_{sun}, we obtain a companion mass M_c = 17.6 ^{+1.5}_{-1.2} M_{Jup}, 3-sigma above a 13 M_{Jup} deuterium burning, brown dwarf lower limit. Dynamical simulations incorporating this accurate mass for HD 38529 c indicate that a near-Saturn mass planet could exist between the two known companions. We find weak evidence of an additional low amplitude signal that can be modeled as a planetary-mass (~0.17 M$_{Jup}) companion at P~194 days. Additional observations (radial velocities and/or Gaia astrometry) are required to validate an interpretation of HD 38529 d as a planetary-mass companion. If confirmed, the resulting HD 38529 planetary system may be an example of a "Packed Planetary System".

preprint2010arXiv

The Roles of Tidal Evolution and Evaporative Mass Loss in the Origin of CoRoT-7 b

CoRoT-7 b is the first confirmed rocky exoplanet, but, with an orbital semi-major axis of 0.0172 AU, its origins may be unlike any rocky planet in our solar system. In this study, we consider the roles of tidal evolution and evaporative mass loss in CoRoT-7 b's history, which together have modified the planet's mass and orbit. If CoRoT-7 b has always been a rocky body, evaporation may have driven off almost half its original mass, but the mass loss may depend sensitively on the extent of tidal decay of its orbit. As tides caused CoRoT-7 b's orbit to decay, they brought the planet closer to its host star, thereby enhancing the mass loss rate. Such a large mass loss also suggests the possibility that CoRoT-7 b began as a gas giant planet and had its original atmosphere completely evaporated. In this case, we find that CoRoT-7 b's original mass probably didn't exceed 200 Earth masses (about 2/3 of a Jupiter mass). Tides raised on the host star by the planet may have significantly reduced the orbital semi-major axis, perhaps causing the planet to migrate through mean-motion resonances with the other planet in the system, CoRoT-7 c. The coupling between tidal evolution and mass loss may be important not only for CoRoT-7 b but also for other close-in exoplanets, and future studies of mass loss and orbital evolution may provide insight into the origin and fate of close-in planets, both rocky and gaseous.

preprint2010arXiv

Tidal effects on brown dwarfs: Application to the eclipsing binary 2MASSJ05352184-0546085 - The anomalous temperature reversal in the context of tidal heating

2MASSJ05352184-0546085 (2M0535-05) is the only known eclipsing brown dwarf (BD) binary, and so may serve as an important benchmark for models of BD formation and evolution. However, theoretical predictions of the system's properties seem inconsistent with observations: i. The more massive (primary) component is observed to be cooler than the less massive (secondary) one. ii. The secondary is more luminous (by roughly 10^{24} W) than expected. We study the impact of tidal heating to the energy budget of both components. We also compare various plausible tidal models to determine a range of predicted properties. We apply two versions of two different, well-known models for tidal interaction, respectively, (i.) the 'constant-phase-lag' model and (ii.) the 'constant-time-lag' model, and incorporate the predicted tidal heating into a model of BD structure. We find that the contribution of heat from tides in 2M0535-05 alone may only be large enough to account for the discrepancies between observation and theory in an unlikely region of the parameter space. The tidal quality factor of BDs, Q_{BD}, would have to be 10^{3.5} and the secondary needs a spin-orbit misalignment greater than 50 degrees. However, tidal synchronization time scales for 2M0535-05 restrict the tidal dissipation function Q_{BD} to values greater than 10^{4.5} and rule out intense tidal heating in 2M0535-05. We provide the first constraint on Q_{BD}. Tidal heating alone is unlikely to be responsible for the surprising temperature reversal within 2M0535-05. But an evolutionary embedment of tidal effects and a coupled treatment with the structural evolution of the BDs is necessary to corroborate or refute this result.

preprint2010arXiv

Tidal Effects on the Habitability of Exoplanets: The Case of GJ 581 d

Tides may be crucial to the habitability of exoplanets. If such planets form around low-mass stars, then those in the circumstellar habitable zone will be close enough to their host stars to experience strong tidal forces. Tides may result in orbital decay and circularization, evolution toward zero obliquity, a fixed rotation rate (not necessarily synchronous), and substantial internal heating. Due to tidal effects, the range of habitable orbital locations may be quite different from that defined by the traditional concept of a habitable zone (HZ) based on stellar insolation, atmospheric effects, and liquid water on a planet's surface. Tidal heating may make locations within the traditional HZ too hot, while planets outside the traditional zone could be rendered quite habitable due to tides. Here we consider these effects on the exoplanet GJ 581 d.

preprint2009arXiv

CoRoT-7 b: Super-Earth or Super-Io?

CoRoT-7 b, a planet about 70% larger than the Earth orbiting a Sun-like star, is the first-discovered rocky exoplanet, and hence has been dubbed a "super-Earth". Some initial studies suggested that since the planet is so close to its host star, it receives enough insolation to partially melt its surface. However, these past studies failed to take into consideration the role that tides may play in this system. Even if the planet's eccentricity has always been zero, we show that tidal decay of semi-major axis could have been large enough that the planet formed on a wider orbit which received less insolation. Moreover, CoRoT-7 b could be tidally heated at a rate that dominates its geophysics and drives extreme volcanism. In this case, CoRoT-7 b is a "super-Io" that, like Jupiter's volcanic moon, is dominated by volcanism and rapid resurfacing. Such heating could occur with an eccentricity of just 10^-5. This small value could be driven by CoRoT-7 c if its own eccentricity is larger than ~10^-4. CoRoT-7 b may be the first of a class of planetary super-Ios likely to be revealed by the CoRoT and Kepler spacecraft.

preprint2009arXiv

N-Body Simulations of Growth from 1 km Planetesimals at 0.4 AU

We present N-body simulations of planetary accretion beginning with 1 km radius planetesimals in orbit about a 1 solar mass star at 0.4 AU. The initial disk of planetesimals contains too many bodies for any current N-body code to integrate; therefore, we model a sample patch of the disk. Although this greatly reduces the number of bodies, we still track in excess of 10^5 particles. We consider three initial velocity distributions and monitor the growth of the planetesimals. The masses of some particles increase by more than a factor of 100. Additionally, the escape speed of the largest particle grows considerably faster than the velocity dispersion of the particles, suggesting impending runaway growth, although no particle grows large enough to detach itself from the power law size-frequency distribution. These results are in general agreement with previous statistical and analytical results. We compute rotation rates by assuming conservation of angular momentum around the center of mass at impact and that merged planetesimals relax to spherical shapes. At the end of our simulations, the majority of bodies that have undergone at least one merger are rotating faster than the breakup frequency. This implies that the assumption of completely inelastic collisions (perfect accretion), which is made in most simulations of planetary growth at sizes 1 km and above, is inappropriate. Our simulations reveal that, subsequent to the number of particles in the patch having been decreased by mergers to half its initial value, the presence of larger bodies in neighboring regions of the disk may limit the validity of simulations employing the patch approximation.

preprint2009arXiv

Observational Evidence for Tidal Destruction of Exoplanets

The distribution of the orbits of close-in exoplanets shows evidence for on-going removal and destruction by tides. Tides raised on a planet's host star cause the planet's orbit to decay, even after the orbital eccentricity has dropped to zero. Comparison of the observed orbital distribution and predictions of tidal theory show good qualitative agreement, suggesting tidal destruction of close-in exoplanets is common. The process can explain the observed cut-off in small a-values, the clustering of orbital periods near three days, and the relative youth of transiting planets. Contrary to previous considerations, a mechanism to stop the inward migration of close-in planets at their current orbits is not necessarily required. Planets nearing tidal destruction may be found with extremely small a, possibly already stripped of any gaseous envelope. The recently discovered CoRoT-Exo-7 b may be an example of such a planet and will probably be destroyed by tides within the next few Gyrs. Also, where one or more planets have already been accreted, a star may exhibit an unusual composition and/or spin rate.

preprint2009arXiv

Planet-planet scattering leads to tightly packed planetary systems

The known extrasolar multiple-planet systems share a surprising dynamical attribute: they cluster just beyond the Hill stability boundary. Here we show that the planet-planet scattering model, which naturally explains the observed exoplanet eccentricity distribution, can reproduce the observed distribution of dynamical configurations. We calculated how each of our scattered systems would appear over an appropriate range of viewing geometries; as Hill stability is weakly dependent on the masses, the mass-inclination degeneracy does not significantly affect our results. We consider a wide range of initial planetary mass distributions and find that some are poor fits to the observed systems. In fact, many of our scattering experiments overproduce systems very close to the stability boundary. The distribution of dynamical configurations of two-planet systems actually may provide better discrimination between scattering models than the distribution of eccentricity. Our results imply that, at least in their inner regions which are weakly affected by gas or planetesimal disks, planetary systems should be "packed", with no large gaps between planets.

preprint2009arXiv

The HD 40307 Planetary System: Super-Earths or Mini-Neptunes?

Three planets with minimum masses less than 10 Earth masses orbit the star HD 40307, suggesting these planets may be rocky. However, with only radial velocity data, it is impossible to determine if these planets are rocky or gaseous. Here we exploit various dynamical features of the system in order to assess the physical properties of the planets. Observations allow for circular orbits, but a numerical integration shows that the eccentricities must be at least 0.0001. Also, planets b and c are so close to the star that tidal effects are significant. If planet b has tidal parameters similar to the terrestrial planets in the Solar System and a remnant eccentricity larger than 0.001, then, going back in time, the system would have been unstable within the lifetime of the star (which we estimate to be 6.1 +/- 1.6 Gyr). Moreover, if the eccentricities are that large and the inner planet is rocky, then its tidal heating may be an order of magnitude greater than extremely volcanic Io, on a per unit surface area basis. If planet b is not terrestrial, e.g. Neptune-like, these physical constraints would not apply. This analysis suggests the planets are not terrestrial-like, and are more like our giant planets. In either case, we find that the planets probably formed at larger radii and migrated early-on (via disk interactions) into their current orbits. This study demonstrates how the orbital and dynamical properties of exoplanet systems may be used to constrain the planets' physical properties.

preprint2008arXiv

The Successful Prediction of the Extrasolar Planet HD 74156 d

Most of the first-discovered extrasolar multi-planet systems were found to lie close to dynamically unstable configurations. However a few observed multi-planet systems (e.g. HD 74156) did not show this trait. Those systems could share this property if they contain an additional planet in between those that are known. Previous investigations identified the properties of hypothetical planets that would place these systems near instability. The hypothetical planet in HD 74156 was expected to have a mass about equal to that of Saturn, a semi-major axis between 0.9 and 1.4 AU, and an eccentricity less than 0.2. HD 74156 d, a planet with a mass of 1.3 Saturn masses at 1.04 AU with an eccentricity of 0.25, was recently announced. We have reanalyzed all published data on this system in order to place tighter constraints on the properties of the new planet. We find two possible orbits for this planet, one close to that already identified and another (with a slightly better fit to the data) at ~0.89 AU. We also review the current status of other planet predictions, discuss the observed single planet systems, and suggest other systems which may contain planets in between those that are already known. The confirmation of the existence of HD 74156 d suggests that planet formation is an efficient process, and planetary systems should typically contain many planets.