Source author record

Matthew J. Payne

Matthew J. Payne 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

26works
4topics
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

26 published item(s)

preprint2024arXiv

Giant Outer Transiting Exoplanet Mass (GOT 'EM) Survey. IV. Long-term Doppler Spectroscopy for 11 Stars Thought to Host Cool Giant Exoplanets

Discovering and characterizing exoplanets at the outer edge of the transit method's sensitivity has proven challenging owing to geometric biases and the practical difficulties associated with acquiring long observational baselines. Nonetheless, a sample of giant exoplanets on orbits longer than 100 days has been identified by transit hunting missions. We present long-term Doppler spectroscopy for 11 such systems with observation baselines spanning a few years to a decade. We model these radial velocity observations jointly with transit photometry to provide initial characterizations of these objects and the systems in which they exist. Specifically, we make new precise mass measurements for four long-period giant exoplanets (Kepler-111 c, Kepler-553 c, Kepler-849 b, and PH-2 b), we place new upper limits on mass for four others (Kepler-421 b, KOI-1431.01, Kepler-1513 b, and Kepler-952 b), and we show that several "confirmed" planets are in fact not planetary at all. We present these findings to complement similar efforts focused on closer-in short-period giant planets, and with the hope of inspiring future dedicated studies of cool giant exoplanets.

preprint2022arXiv

The Population of Interstellar Objects Detectable with the LSST and Accessible for $\textit{In Situ}$ Rendezvous with Various Mission Designs

The recently discovered population of interstellar objects presents us with the opportunity to characterize material from extrasolar planetary and stellar systems up close. The forthcoming Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will provide an unprecedented increase in sensitivity to these objects compared to the capabilities of currently operational observational facilities. We generate a synthetic population of `Oumuamua-like objects drawn from their galactic kinematics, and identify the distribution of impact parameters, eccentricities, hyperbolic velocities and sky locations of objects detectable with the LSST, assuming no cometary activity. This population is characterized by a clustering of trajectories in the direction of the solar apex and anti-apex, centered at orbital inclinations of $\sim90^\circ$. We identify the ecliptic or solar apex as the optimal sky locations to search for future interstellar objects as a function of survey limiting magnitude. Moreover, we identify the trajectories of detectable objects that will be reachable for \textit{in-situ} rendezvous with a dedicated mission with the capabilities of the forthcoming \textit{Comet Interceptor} or proposed \textit{Bridge} concept. By scaling our fractional population statistics with the inferred spatial number density, we estimate that the LSST will detect of order $\sim 15$ interstellar objects over the course of its $\sim10$ year observational campaign. Furthermore, we find that there should be $\sim1-3$ and $\sim0.0007-0.001$ reachable targets for missions with propulsion capabilities comparable to \textit{Bridge} and \textit{Comet Interceptor}, respectively. These numbers are lower limits, and will be readily updateable when the number density and size frequency distribution of interstellar objects are better constrained.

preprint2020arXiv

Modeling Radial Velocity Data of Resonant Planets to Infer Migration Histories

A number of giant planet pairs discovered by the radial velocity method with period ratios $\lesssim 2$ may reside in mean motion resonances. Convergent orbital migration and resonant capture at the time of formation would naturally explain the present-day resonant orbital configurations of these systems. Planets that experience smooth migration and eccentricity damping forces due to a proto-planetary disk should not only capture into mean motion resonances but also end up in a specific dynamical configuration within the resonance, sometimes referred to as apsidal corotation resonance (ACR). Here we develop a method for testing the hypothesis that a planet pair resides in an ACR by directly fitting radial velocity data. The ACR hypothesis strongly restricts the number of free parameters describing the radial velocity signal and we compare fits using this highly restricted model to fits using a more conventional two-planet RV model by using nested sampling simulations. We apply our method to HD 45364 and HD 33844, two systems hosting giant planet pairs in 3:2 and 5:3 resonances, respectively. We demonstrate that the observations of both systems support an ACR configuration and we use the results of our ACR model fits to constrain possible migration histories of these systems.

preprint2016arXiv

Discovery of A New Retrograde Trans-Neptunian Object: Hint of A Common Orbital Plane for Low Semi-Major Axis, High Inclination TNOs and Centaurs

Although the majority of Centaurs are thought to have originated in the scattered disk, with the high-inclination members coming from the Oort cloud, the origin of the high inclination component of trans-Neptunian objects (TNOs) remains uncertain. We report the discovery of a retrograde TNO, which we nickname "Niku", detected by the Pan-STARRS 1 Outer Solar System Survey. Our numerical integrations show that the orbital dynamics of Niku are very similar to that of 2008 KV$_{42}$ (Drac), with a half-life of $\sim 500$ Myr. Comparing similar high inclination TNOs and Centaurs ($q > 10$ AU, $a < 100$ AU and $i > 60^\circ$), we find that these objects exhibit a surprising clustering of ascending node, and occupy a common orbital plane. This orbital configuration has high statistical significance: 3.8-$σ$. An unknown mechanism is required to explain the observed clustering. This discovery may provide a pathway to investigate a possible reservoir of high-inclination objects.

preprint2016arXiv

Liberating exomoons in white dwarf planetary systems

Previous studies indicate that more than a quarter of all white dwarf (WD) atmospheres are polluted by remnant planetary material, with some WDs being observed to accrete the mass of Pluto in 10^6 years. The short sinking timescale for the pollutants indicate that the material must be frequently replenished. Moons may contribute decisively to this pollution process if they are liberated from their parent planets during the post-main-sequence evolution of the planetary systems. Here, we demonstrate that gravitational scattering events among planets in WD systems easily triggers moon ejection. Repeated close encounters within tenths of a planetary Hill radii are highly destructive to even the most massive, close-in moons. Consequently, scattering increases both the frequency of perturbing agents in WD systems, as well as the available mass of polluting material in those systems, thereby enhancing opportunities for collision and fragmentation and providing more dynamical pathways for smaller bodies to reach the WD. Moreover, during intense scattering, planets themselves have pericenters with respect to the WD of only a fraction of an AU, causing extreme Hill-sphere contraction, and the liberation of moons into WD-grazing orbits. Many of our results are directly applicable to exomoons orbiting planets around main sequence stars.

preprint2016arXiv

Observational Constraints on Planet Nine: Astrometry of Pluto and Other Trans-Neptunian Objects

We use astrometry of Pluto and other TNOs to constrain the sky location, distance, and mass of the possible additional planet (Planet Nine) hypothesized by Batygin and Brown (2016). We find that over broad regions of the sky, the inclusion of a massive, distant planet degrades the fits to the observations. However, in other regions, the fits are significantly improved by the addition of such a planet. Our best fits suggest a planet that is either more massive or closer than argued for by Batygin and Brown (2016) based on the orbital distribution of distant trans-neptunian objects (or by Fienga et al. (2016) based on range measured to the Cassini spacecraft). The trend to favor larger and closer perturbing planets is driven by the residuals to the astrometry of Pluto, remeasured from photographic plates using modern stellar catalogs (Buie and Folkner 2015), which show a clear trend in declination, over the course of two decades, that drive a preference for large perturbations. Although this trend may be the result of systematic errors of unknown origin in the observations, a possible resolution is that the declination trend may be due to perturbations from a body, additional to Planet Nine , that is closer to Pluto, but less massive than, Planet Nine .

preprint2016arXiv

The Fate of Exomoons in White Dwarf Planetary Systems

Roughly 1000 white dwarfs are known to be polluted with planetary material, and the progenitors of this material are typically assumed to be asteroids. The dynamical architectures which perturb asteroids into white dwarfs are still unknown, but may be crucially dependent on moons liberated from parent planets during post-main-sequence gravitational scattering. Here, we trace the fate of these exomoons, and show that they more easily achieve deep radial incursions towards the white dwarf than do scattered planets. Consequently, moons are likely to play a significant role in white dwarf pollution, and in some cases may be the progenitors of the pollution itself.

preprint2015arXiv

An Empirically Derived Three-Dimensional Laplace Resonance in the Gliese 876 Planetary System

We report constraints on the three-dimensional orbital architecture for all four planets known to orbit the nearby M dwarf Gliese 876 based solely on Doppler measurements and demanding long-term orbital stability. Our dataset incorporates publicly available radial velocities taken with the ELODIE and CORALIE spectrographs, HARPS, and Keck HIRES as well as previously unpublished HIRES velocities. We first quantitatively assess the validity of the planets thought to orbit GJ 876 by computing the Bayes factors for a variety of different coplanar models using an importance sampling algorithm. We find that a four-planet model is preferred over a three-planet model. Next, we apply a Newtonian MCMC algorithm to perform a Bayesian analysis of the planet masses and orbits using an n-body model in three-dimensional space. Based on the radial velocities alone, we find that a 99% credible interval provides upper limits on the mutual inclinations for the three resonant planets ($Φ_{cb}<6.20^\circ$ for the "c" and "b" pair and $Φ_{be}<28.5^\circ$ for the "b" and "e" pair). Subsequent dynamical integrations of our posterior sample find that the GJ 876 planets must be roughly coplanar ($Φ_{cb}<2.60^\circ$ and $Φ_{be}<7.87^\circ$), suggesting the amount of planet-planet scattering in the system has been low. We investigate the distribution of the respective resonant arguments of each planet pair and find that at least one argument for each planet pair and the Laplace argument librate. The libration amplitudes in our three-dimensional orbital model supports the idea of the outer-three planets having undergone significant past disk migration.

preprint2015arXiv

Kepler-432: a red giant interacting with one of its two long period giant planets

We report the discovery of Kepler-432b, a giant planet ($M_b = 5.41^{+0.32}_{-0.18} M_{\rm Jup}, R_b = 1.145^{+0.036}_{-0.039} R_{\rm Jup}$) transiting an evolved star $(M_\star = 1.32^{+0.10}_{-0.07} M_\odot, R_\star = 4.06^{+0.12}_{-0.08} R_\odot)$ with an orbital period of $P_b = 52.501129^{+0.000067}_{-0.000053}$ days. Radial velocities (RVs) reveal that Kepler-432b orbits its parent star with an eccentricity of $e = 0.5134^{+0.0098}_{-0.0089}$, which we also measure independently with asterodensity profiling (AP; $e=0.507^{+0.039}_{-0.114}$), thereby confirming the validity of AP on this particular evolved star. The well-determined planetary properties and unusually large mass also make this planet an important benchmark for theoretical models of super-Jupiter formation. Long-term RV monitoring detected the presence of a non-transiting outer planet (Kepler-432c; $M_c \sin{i_c} = 2.43^{+0.22}_{-0.24} M_{\rm Jup}, P_c = 406.2^{+3.9}_{-2.5}$ days), and adaptive optics imaging revealed a nearby (0\farcs87), faint companion (Kepler-432B) that is a physically bound M dwarf. The host star exhibits high signal-to-noise asteroseismic oscillations, which enable precise measurements of the stellar mass, radius and age. Analysis of the rotational splitting of the oscillation modes additionally reveals the stellar spin axis to be nearly edge-on, which suggests that the stellar spin is likely well-aligned with the orbit of the transiting planet. Despite its long period, the obliquity of the 52.5-day orbit may have been shaped by star-planet interaction in a manner similar to hot Jupiter systems, and we present observational and theoretical evidence to support this scenario. Finally, as a short-period outlier among giant planets orbiting giant stars, study of Kepler-432b may help explain the distribution of massive planets orbiting giant stars interior to 1 AU.

preprint2014arXiv

Formation of irregular and runaway moons/exomoons through moon-moon scattering

Gas giant planets in the Solar system host large satellite systems with multiple regular and irregular moons. Regular moons revolve around their host planet in circular, low inclination short period orbits, and are thought to form in-situ through coagulation processes. In contrast, irregular moons have highly inclined (and even retrograde), typically eccentric and long period orbits around their host planet. Irregular moons are therefore often thought to be unbound objects in helio-centric orbits that were later captured to their current orbits around the planet. Here we study the a different possibility in which regular moons form in-situ, and dynamically evolve through mutual moon-moon scattering. We find that such evolution can excite the satellites into high eccentricities and inclinations. We find that moons are either ejected from the host planet to become runaway moons, or stay bound and become prograde orbiting irregular moons with inclined and eccentric orbits around their host planet. Ejected moons, unbound to the planet, can later be temporarily re-captured by the host planet even at retrograde orbits. Such moons are eventually re-ejected from the system or collide with the planet, at least in the absence of dissipative processes (e.g. collisions with existing bound moons, a debris disk or through tidal interactions with the host planet), not currently modeled. Uncaptured runaway moons may eventually be ejected from the Solar system, or be captured into stable helio-centric orbits and contribute to the populations of asteroidal or trans-Neptunian objects. Such scenarios are potentially relevant both for the gas-giant satellites in the Solar system and for the dynamical evolution of exomoons.

preprint2014arXiv

Newly-Discovered Planets Orbiting HD~5319, HD~11506, HD~75784 and HD~10442 from the N2K Consortium

Initially designed to discover short-period planets, the N2K campaign has since evolved to discover new worlds at large separations from their host stars. Detecting such worlds will help determine the giant planet occurrence at semi-major axes beyond the ice line, where gas giants are thought to mostly form. Here we report four newly-discovered gas giant planets (with minimum masses ranging from 0.4 to 2.1 MJup) orbiting stars monitored as part of the N2K program. Two of these planets orbit stars already known to host planets: HD 5319 and HD 11506. The remaining discoveries reside in previously-unknown planetary systems: HD 10442 and HD 75784. The refined orbital period of the inner planet orbiting HD 5319 is 641 days. The newly-discovered outer planet orbits in 886 days. The large masses combined with the proximity to a 4:3 mean motion resonance make this system a challenge to explain with current formation and migration theories. HD 11506 has one confirmed planet, and here we confirm a second. The outer planet has an orbital period of 1627.5 days, and the newly-discovered inner planet orbits in 223.6 days. A planet has also been discovered orbiting HD 75784 with an orbital period of 341.7 days. There is evidence for a longer period signal; however, several more years of observations are needed to put tight constraints on the Keplerian parameters for the outer planet. Lastly, an additional planet has been detected orbiting HD 10442 with a period of 1043 days.

preprint2014arXiv

The 55 Cancri Planetary System: Fully Self-Consistent N-body Constraints and a Dynamical Analysis

We present an updated study of the planets known to orbit 55 Cancri A using 1,418 high-precision radial velocity observations from four observatories (Lick, Keck, Hobby-Eberly Telescope, Harlan J. Smith Telescope) and transit time/durations for the inner-most planet, 55 Cancri "e" (Winn et al. 2011). We provide the first posterior sample for the masses and orbital parameters based on self-consistent n-body orbital solutions for the 55 Cancri planets, all of which are dynamically stable (for at least $10^8$ years). We apply a GPU version of Radial velocity Using N-body Differential evolution Markov Chain Monte Carlo (RUN DMC; B. Nelson et al. 2014) to perform a Bayesian analysis of the radial velocity and transit observations. Each of the planets in this remarkable system has unique characteristics. Our investigation of high-cadence radial velocities and priors based on space-based photometry yields an updated mass estimate for planet "e" ($8.09\pm0.26$ M$_\oplus$), which affects its density ($5.51\pm^{1.32}_{1.00}$ g cm$^{-3}$) and inferred bulk composition. Dynamical stability dictates that the orbital plane of planet "e" must be aligned to within $60^o$ of the orbital plane of the outer planets (which we assume to be coplanar). The mutual interactions between the planets "b" and "c" may develop an apsidal lock about $180^o$. We find 36-45% of all our model systems librate about the anti-aligned configuration with an amplitude of $51^o\pm^{6^o}_{10^o}$. Other cases showed short-term perturbations in the libration of $\varpi_b-\varpi_c$, circulation, and nodding, but we find the planets are not in a 3:1 mean-motion resonance. A revised orbital period and eccentricity for planet "d" pushes it further toward the closest known Jupiter analog in the exoplanet population.

preprint2013arXiv

A Habitable Zone Census via Transit Timing and the Imperative for Continuing to Observe the Kepler Field

We propose a scientific program to complete a census of planets, characterizing their masses, orbital properties, and dynamical histories using continued observations of the Kepler field of view with the Kepler spacecraft in a two reaction wheel mission. Even with a significantly reduced photometric precision, extending time-domain observations of this field is uniquely capable of pursuing several critical science goals: 1) measuring the architectures of planetary systems by identifying non-transiting planets interleaved among known transiting planets, 2) establishing the mass-radius relationship for planets in the important transition region between small, gas-rich sub-Neptune planets and large, rocky super-Earths, and 3) uncovering dynamical evidence of the formation and evolution of the inner regions of planetary systems. To meet these objectives, the unique multi-object observing capabilities of Kepler will be used in a set of concurrent campaigns with specific motivations. These campaigns focus largely on the ability to interpret Transit Timing Variations (TTVs) that result from dynamical interactions among planets in a system and include: 1) observations of systems that exhibit large TTVs and are particularly rich in dynamical information, 2) observations of systems where additional transit times will yield mass measurements of the constituent planets, 3) observations of systems where the TTV signal evolves over very long timescales, and 4) observations of systems with long-period planet candidates where additional transits will remove orbital period ambiguities caused by gaps in the original Kepler data.

preprint2013arXiv

Characterizing the Orbital and Dynamical State of the HD 82943 Planetary System With Keck Radial Velocity Data

We present an updated analysis of radial velocity data of the HD 82943 planetary system based on 10 years of measurements obtained with the Keck telescope. Previous studies have shown that the HD 82943 system has two planets that are likely in 2:1 mean-motion resonance (MMR), with the orbital periods about 220 and 440 days (Lee et al. 2006). However, alternative fits that are qualitatively different have also been suggested, with two planets in a 1:1 resonance (Gozdziewski & Konacki 2006) or three planets in a Laplace 4:2:1 resonance (Beauge et al. 2008). Here we use \c{hi}2 minimization combined with parameter grid search to investigate the orbital parameters and dynamical states of the qualitatively different types of fits, and we compare the results to those obtained with the differential evolution Markov chain Monte Carlo method. Our results support the coplanar 2:1 MMR configuration for the HD 82943 system, and show no evidence for either the 1:1 or 3-planet Laplace resonance fits. The inclination of the system with respect to the sky plane is well constrained at about 20(+4.9 -5.5) degree, and the system contains two planets with masses of about 4.78 MJ and 4.80 MJ (where MJ is the mass of Jupiter) and orbital periods of about 219 and 442 days for the inner and outer planet, respectively. The best fit is dynamically stable with both eccentricity-type resonant angles θ1 and θ2 librating around 0 degree.

preprint2013arXiv

RUN DMC: An efficient, parallel code for analyzing Radial Velocity Observations using N-body Integrations and Differential Evolution Markov chain Monte Carlo

In the 20+ years of Doppler observations of stars, scientists have uncovered a diverse population of extrasolar multi-planet systems. A common technique for characterizing the orbital elements of these planets is Markov chain Monte Carlo (MCMC), using a Keplerian model with random walk proposals and paired with the Metropolis-Hastings algorithm. For approximately a couple of dozen planetary systems with Doppler observations, there are strong planet-planet interactions due to the system being in or near a mean-motion resonance (MMR). An N-body model is often required to accurately describe these systems. Further computational difficulties arise from exploring a high-dimensional parameter space ($\sim$7 x number of planets) that can have complex parameter correlations. To surmount these challenges, we introduce a differential evolution MCMC (DEMCMC) applied to radial velocity data while incorporating self-consistent N-body integrations. Our Radial velocity Using N-body DEMCMC (RUN DMC) algorithm improves upon the random walk proposal distribution of the traditional MCMC by using an ensemble of Markov chains to adaptively improve the proposal distribution. We describe the methodology behind the algorithm, along with results of tests for accuracy and performance. We find that most algorithm parameters have a modest effect on the rate of convergence. However, the size of the ensemble can have a strong effect on performance. We show that the optimal choice depends on the number of planets in a system, as well as the computer architecture used and the resulting extent of parallelization. While the exact choices of optimal algorithm parameters will inevitably vary due to the details of individual planetary systems, we offer recommendations for choosing the DEMCMC algorithm's algorithmic parameters that result in excellent performance for a wide variety of planetary systems.

preprint2013arXiv

Stability of Satellites in Closely Packed Planetary Systems

We perform numerical integrations of four-body (star, planet, planet, satellite) systems to investigate the stability of satellites in planetary Systems with Tightly-packed Inner Planets (STIPs). We find that the majority of closely-spaced stable two-planet systems can stably support satellites across a range of parameter-space which is only slightly decreased compared to that seen for the single-planet case. In particular, circular prograde satellites remain stable out to $\sim 0.4 R_H$ (where $R_H$ is the Hill Radius) as opposed to $\sim 0.5 R_H$ in the single-planet case. A similarly small restriction in the stable parameter-space for retrograde satellites is observed, where planetary close approaches in the range 2.5 to 4.5 mutual Hill radii destabilize most satellites orbits only if $a\sim 0.65 R_H$. In very close planetary pairs (e.g. the 12:11 resonance) the addition of a satellite frequently destabilizes the entire system, causing extreme close-approaches and the loss of satellites over a range of circumplanetary semi-major axes. The majority of systems investigated stably harbored satellites over a wide parameter-space, suggesting that STIPs can generally offer a dynamically stable home for satellites, albeit with a slightly smaller stable parameter-space than the single-planet case. As we demonstrate that multi-planet systems are not a priori poor candidates for hosting satellites, future measurements of satellite occurrence rates in multi-planet systems versus single-planet systems could be used to constrain either satellite formation or past periods of strong dynamical interaction between planets.

preprint2012arXiv

Constraining the Planetary System of Fomalhaut Using High-Resolution ALMA Observations

The dynamical evolution of planetary systems leaves observable signatures in debris disks. Optical images trace micron-sized grains, which are strongly affected by stellar radiation and need not coincide with their parent body population. Observations of mm-size grains accurately trace parent bodies, but previous images lack the resolution and sensitivity needed to characterize the ring's morphology. Here we present ALMA 350 GHz observations of the Fomalhaut debris ring. These observations demonstrate that the parent body population is 13-19 AU wide with a sharp inner and outer boundary. We discuss three possible origins for the ring, and suggest that debris confined by shepherd planets is the most consistent with the ring's morphology.

preprint2012arXiv

Interactions Between Moderate- and Long-Period Giant Planets: Scattering Experiments for Systems in Isolation and with Stellar Flybys

The chance that a planetary system will interact with another member of its host star's nascent cluster would be greatly increased if gas giant planets form in situ on wide orbits. In this paper, we explore the outcomes of planet-planet scattering for a distribution of multiplanet systems that all have one of the planets on an initial orbit of 100 AU. The scattering experiments are run with and without stellar flybys. We convolve the outcomes with distributions for protoplanetary disk and stellar cluster sizes to generalize the results where possible. We find that the frequencies of large mutual inclinations and high eccentricities are sensitive to the number of planets in a system, but not strongly to stellar flybys. However, flybys do play a role in changing the low and moderate portions of the mutual inclination distributions, and erase dynamically cold initial conditions on average. Wide-orbit planets can be mixed throughout the planetary system, and in some cases, can potentially become hot Jupiters, which we demonstrate using scattering experiments that include a tidal damping model. If planets form on wide orbits in situ, then there will be discernible differences in the proper motion distributions of a sample of wide-orbit planets compared with a pure scattering formation mechanism. Stellar flybys can enhance the frequency of ejections in planetary systems, but auto-ionization is likely to remain the dominant source of free-floating planets.

preprint2012arXiv

The Discovery of HD 37605c and a Dispositive Null Detection of Transits of HD 37605b

We report the radial-velocity discovery of a second planetary mass companion to the K0 V star HD 37605, which was already known to host an eccentric, P~55 days Jovian planet, HD 37605b. This second planet, HD 37605c, has a period of ~7.5 years with a low eccentricity and an Msini of ~3.4 MJup. Our discovery was made with the nearly 8 years of radial velocity follow-up at the Hobby-Eberly Telescope and Keck Observatory, including observations made as part of the Transit Ephemeris Refinement and Monitoring Survey (TERMS) effort to provide precise ephemerides to long-period planets for transit follow-up. With a total of 137 radial velocity observations covering almost eight years, we provide a good orbital solution of the HD 37605 system, and a precise transit ephemeris for HD 37605b. Our dynamic analysis reveals very minimal planet-planet interaction and an insignificant transit time variation. Using the predicted ephemeris, we performed a transit search for HD 37605b with the photometric data taken by the T12 0.8-m Automatic Photoelectric Telescope (APT) and the Microvariability and Oscillations of Stars (MOST) satellite. Though the APT photometry did not capture the transit window, it characterized the stellar activity of HD 37605, which is consistent of it being an old, inactive star, with a tentative rotation period of 57.67 days. The MOST photometry enabled us to report a dispositive null detection of a non-grazing transit for this planet. Within the predicted transit window, we exclude an edge-on predicted depth of 1.9% at >>10sigma, and exclude any transit with an impact parameter b>0.951 at greater than 5sigma. We present the BOOTTRAN package for calculating Keplerian orbital parameter uncertainties via bootstrapping. We found consistency between our orbital parameters calculated by the RVLIN package and error bars by BOOTTRAN with those produced by a Bayesian analysis using MCMC.

preprint2012arXiv

Traditional formation scenarios fail to explain 4:3 mean motion resonances

At least two multi-planetary systems in a 4:3 mean motion resonance have been found by radial velocity surveys. These planets are gas giants and the systems are only stable when protected by a resonance. Additionally the Kepler mission has detected at least 4 strong candidate planetary systems with a period ratio close to 4:3. This paper investigates traditional dynamical scenarios for the formation of these systems. We systematically study migration scenarios with both N-body and hydro-dynamic simulations. We investigate scenarios involving the in-situ formation of two planets in resonance. We look at the results from finely tuned planet-planet scattering simulations with gas disk damping. Finally, we investigate a formation scenario involving isolation-mass embryos. Although the combined planet-planet scattering and damping scenario seems promising, none of the above scenarios is successful in forming enough systems in 4:3 resonance with planetary masses similar to the observed ones. This is a negative result but it has important implications for planet formation. Previous studies were successful in forming 2:1 and 3:2 resonances. This is generally believed to be evidence of planet migration. We highlight the main differences between those studies and our failure in forming a 4:3 resonance. We also speculate on more exotic and complicated ideas. These results will guide future investigators toward exploring the above scenarios and alternative mechanisms in a more general framework.

preprint2011arXiv

An Analysis of Jitter and Transit Timing Variations in the HAT-P-13 System

If the two planets in the HAT-P-13 system are coplanar, the orbital states provide a probe of the internal planetary structure. Previous analyses of radial velocity and transit timing data of the system suggested that the observational constraints on the orbital states were rather small. We reanalyze the available data, treating the jitter as an unknown MCMC parameter, and find that a wide range of jitter values are plausible, hence the system parameters are less well constrained than previously suggested. For slightly increased levels of jitter ($\sim 4.5\,m\,s^{-1}$) the eccentricity of the inner planet can be in the range $0<e_{inner}<0.07$, the period and eccentricity of the outer planet can be $440<P_{outer}<470$ days and $0.55<e_{outer}<0.85$ respectively, while the relative pericenter alignment, $η$, of the planets can take essentially any value $-180^{\circ}<η<+180^{\circ}$. It is therefore difficult to determine whether $e_{inner}$ and $η$ have evolved to a fixed-point state or a limit cycle, or to use $e_{inner}$ to probe the internal planetary structure. We perform various transit timing variation (TTV) analyses, demonstrating that current constraints merely restrict $e_{outer}<0.85$, and rule out relative planetary inclinations within $\sim 2^{\circ}$ of $i_{rel}=90^{\circ}$, but that future observations could significantly tighten the restriction on both these parameters. We demonstrate that TTV profiles can readily distinguish the theoretically favored inclinations of $i_{rel}=0^{\circ}\,&\,45^{\circ}$, provided that sufficiently precise and frequent transit timing observations of HAT-P-13b can be made close to the pericenter passage of HAT-P-13c. We note the relatively high probability that HAT-P-13c transits and suggest observational dates and strategies.

preprint2011arXiv

The Heavy Element Composition of Disk Instability Planets Can Range From Sub- to Super-Nebular

Transit surveys combined with Doppler data have revealed a class of gas giant planets that are massive and highly enriched in heavy elements (e.g., HD149026b, GJ436b, and HAT-P-20b). It is tempting to consider these planets as validation of core accretion plus gas capture because it is often assumed that disk instability planets should be of nebular composition. We show in this paper, to the contrary, that gas giants that form by disk instability can have a variety of heavy element compositions, ranging from sub- to super-nebular values. High levels of enrichment can be achieved through one or multiple mechanisms, including enrichment at birth, planetesimal capture, and differentiation plus tidal stripping. As a result, the metallicity of an individual gas giant cannot be used to discriminate between gas giant formation modes.

preprint2010arXiv

From Dust To Planetesimal: The Snowball Phase ?

The standard model of planet formation considers an initial phase in which planetesimals form from a dust disk, followed by a phase of mutual planetesimal-planetesimal collisions, leading eventually to the formation of planetary embryos. However, there is a potential transition phase (which we call the "snowball phase"), between the formation of the first planetesimals and the onset of mutual collisions amongst them, which has often been either ignored or underestimated in previous studies. In this snowball phase, isolated planetesimals move on Keplerian orbits and grow solely via the direct accretion of sub-cm sized dust entrained with the gas in the protoplanetary disk. Using a simplified model in which planetesimals are progressively produced from the dust, we consider the expected sizes to which the planetesimals can grow before mutual collisions commence and derive the dependence of this size on a number of critical parameters, including the degree of disk turbulence, the planetesimal size at birth and the rate of planetesimal creation. For systems in which turbulence is weak and the planetesimals are created at a low rate and with relatively small birth size, we show that the snowball growth phase can be very important, allowing planetesimals to grow by a factor of 10^6 in mass before mutual collisions take over. In such cases, the snowball growth phase can be the dominant mode to transfer mass from the dust to planetesimals. Moreover, such growth can take place within the typical lifetime of a protoplanetary gas disk. A noteworthy result is that ... ...(see the paper). For the specific case of close binaries such as Alpha Centauri ... ... (see the paper). From a more general perspective, these preliminary results suggest that an efficient snowball growth phase provides a large amount of "room at the bottom" for theories of planet formation.

preprint2010arXiv

Quantifying the challenges of detecting unseen planetary companions with transit timing variations

Both ground and space-based transit observatories are poised to significantly increase the number of known transiting planets and the number of precisely measured transit times. The variation in a planet's transit times may be used to infer the presence of additional planets. Deducing the masses and orbital parameters of such planets from transit time variations (TTVs) alone is a rich and increasingly relevant dynamical problem. In this work, we evaluate the extent of the degeneracies in this process, systematically explore the dependence of TTV signals on several parameters and provide phase space plots that could aid observers in planning future observations. Our explorations are focused on a likely-to-be prevalent situation: a known transiting short-period Neptune or Jupiter-sized planet and a suspected external low-mass perturber on a nearly-coplanar orbit. Through approximately 10^7 N-body simulations, we demonstrate how TTV signal amplitudes may vary by orders of magnitude due to slight variations in any one orbital parameter (0.001 AU in semimajor axis, 0.005 in eccentricity, or a few degrees in orbital angles), and quantify the number of consecutive transit observations necessary in order to obtain a reasonable opportunity to characterize the unseen planet (approximately greater or equal to 50 observations). Planets in or near period commensurabilities of the form p:q, where p < 21 and q < 4, produce distinct TTV signatures, regardless of whether the planets are actually locked in a mean motion resonance. We distinguish these systems from the secular systems in our explorations. Additionally, we find that computing the autocorrelation function of a TTV signal can provide a useful diagnostic for identifying possible orbits for additional planets and suggest that this method could aid integration of TTV signals in future studies of particular exosystems.

preprint2010arXiv

Transit Timing Variations for Inclined and Retrograde Exoplanetary Systems

We perform numerical calculations of the expected transit timing variations (TTVs) induced on a Hot-Jupiter by an Earth-mass perturber. Motivated by the recent discoveries of retrograde transiting planets, we concentrate on an investigation of the effect of varying relative planetary inclinations, up to and including completely retrograde systems. We find that planets in low order (E.g. 2:1) mean-motion resonances (MMRs) retain approximately constant TTV amplitudes for $0<\,^{\circ}i<170\,^{\circ}$, only reducing in amplitude for $i > 170\,^{\circ}$. Systems in higher order MMRs (E.g. 5:1) increase in TTV amplitude as inclinations increase towards $45\,^{\circ}$, becoming approximately constant for $45 < i < 135$, and then declining for $i > 135\,^{\circ}$. Planets away from resonance slowly decrease in TTV amplitude as inclinations increase from 0 to 180, where-as planets adjacent to resonances can exhibit a huge range of variability in TTV amplitude as a function of both eccentricity and inclination. For highly retrograde systems ($135\,^{\circ} < i \leq 180\,^{\circ}$), TTV signals will be undetectable across almost the entirety of parameter space, with the exceptions occurring when the perturber has high eccentricity or is very close to a MMR. This high inclination decrease in TTV amplitude (on and away from resonance) is important for the analysis of the known retrograde and multi-planet transiting systems, as inclination effects need to be considered if TTVs are to be used to exclude the presence of any putative planetary companions: absence of evidence is not evidence of absence.

preprint2009arXiv

Outward Migration of Terrestrial Embryos in Binary Systems

We consider the formation and migration of protoplanetary embryos in disks around the stars in tight binary systems (separations ~ 20 AU. In such systems, the initial stages of runaway embryo formation are expected to only take place within some critical disk radius a_{crit}, due to the perturbing effect of the binary companions (Thebault et al. 2009). We perform n-body simulations of the evolution of such a population of inner-disk embryos surrounded by an outer-disk of smaller planetesimals. Taking Alpha Centauri-B as our fiducial reference example in which a_{crit} ~ 0.7 AU, and using a Minimum Mass Nebular Model with $Σ\propto a^{-3/2}$, we find that within 10^6 yrs (10^7 yrs), systems will on average contain embryos which have migrated out to 0.9 AU (1.2 AU), with the average outer-most body having a mass of 0.2 M_{earth} 0.4 M_{earth}. Changes to increase the surface density of solids or to use a flatter profile both produce increased embryo migration and growth. At a given time, the relative change in semi-major axis of the outer-most embryo in these simulations is found to be essentially independent of a_{crit}, and we note that little further embryo migration takes place beyond 10^7 years. We conclude that the suppression of runaway growth outside a_{crit} does not mean that the habitable zones in such tight binary systems will be devoid of detectable, terrestrial mass planets, even if a_{crit} lies significantly interior to the inner edge of the habitable zone.