BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Man Hoi Lee

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

ResearcherUnclaimed source record
astro-ph.EPastro-phastro-ph.SRcond-mat.stat-mech

Catalog footprint

What is connected

13works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

13 published item(s)

preprint2022arXiv

Dynamical architecture of the HD 107148 system

We present an independent Doppler validation and dynamical orbital analysis of the two-planet system HD 107148, which was recently announced in Rosenthal et al. (2021). Our detailed analyses are based on literature HIRES data and newly obtained HARPS and CARMENES radial velocity (RV) measurements as part of our survey in search for additional planets around single planet systems. We perform a periodogram analysis of the available HIRES and HARPS precise RVs and stellar activity indicators. We do not find any apparent correlation between the RV measurements and the stellar activity indicators, thus linking the two strong periodicities to a moderately compact multiple-planet system. We carry out orbital fitting analysis by testing various one- and two-planet orbital configurations and studying the posterior probability distribution of the fitted parameters. Our results solidify the existence of a Saturn-mass planet (HD 107148 b, discovered first) with a period $P_b\sim77.2$ d, and a second, eccentric ($e_c \sim$ 0.4), Neptune-mass exoplanet (HD 107148 c), with an orbital period of $P_c\sim18.3$ d. Finally, we investigate the two-planet system's long-term stability and overall orbital dynamics with the posterior distribution of our preferred orbital configuration. Our N-body stability simulations show that the system is long-term stable and exhibits large secular osculations in eccentricity but in no particular mean-motion resonance configuration. The HD 107148 system, consisting of a Solar-type main sequence star with two giant planets in a rare configuration, features a common proper motion white dwarf companion and is, therefore, a valuable target for understanding the formation and evolution of planetary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Numerical Method for Determining the Elements of Circumbinary Orbits and Its Application to Circumbinary Planets and the Satellites of Pluto-Charon

Planets and satellites orbiting a binary system exist in the solar system and extrasolar planetary systems. Their orbits can be significantly different from Keplerian orbits, if they are close to the binary and the secondary-to-primary mass ratio is high. A proper description of a circumbinary orbit is in terms of the free eccentricity $e_{\rm free}$ at the epicyclic frequency $κ_0$, forced eccentricity $e_{\rm forced}$ at the mean motion $n_0$, and oscillations at higher frequencies forced by the non-axisymmetric components of the binary's potential. We show that accurate numerical values for the amplitudes and frequencies of these terms can be extracted from numerical orbit integrations by applying fast Fourier transformation (FFT) to the cylindrical distance between the circumbinary object and the center of mass of the binary as a function of time. We apply this method to three Kepler circumbinary planets and the satellites of Pluto-Charon. For the satellite Styx of Pluto-Charon, the FFT results for $κ_0$ and $e_{\rm free}$ differ significantly from the first-order analytic value and the value reported by Showalter & Hamilton (2015), respectively. We show that the deviation in $κ_0$ is likely due to the effect of the 3:1 mean-motion resonance and discuss the implications of the lower value for $e_{\rm free}$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Precise radial velocities of giant stars IX. HD 59686 Ab: a massive circumstellar planet orbiting a giant star in a ~13.6 au eccentric binary system

Context: For over 12 yr, we have carried out a precise radial velocity survey of a sample of 373 G and K giant stars using the Hamilton Échelle Spectrograph at Lick Observatory. There are, among others, a number of multiple planetary systems in our sample as well as several planetary candidates in stellar binaries. Aims: We aim at detecting and characterizing substellar+stellar companions to the giant star HD 59686 A (HR 2877, HIP 36616). Methods: We obtained high precision radial velocity (RV) measurements of the star HD 59686 A. By fitting a Keplerian model to the periodic changes in the RVs, we can assess the nature of companions in the system. In order to discriminate between RV variations due to non-radial pulsation or stellar spots we used infrared RVs taken with the CRIRES spectrograph at the Very Large Telescope. Additionally, to further characterize the system, we obtain high-resolution images with LMIRCam at the Large Binocular Telescope. Results: We report the likely discovery of a giant planet with a mass of $m_{p}~\sin i=6.92_{-0.24}^{+0.18}~M_{Jup}$ orbiting at $a_{p}=1.0860_{-0.0007}^{+0.0006}$ au from the giant star HD 59686 A. Besides the planetary signal, we discover an eccentric ($e_{B}=0.729_{-0.003}^{+0.004}$) binary companion with a mass of $m_{B}~\sin i=0.5296_{-0.0008}^{+0.0011}~M_{Sun}$ orbiting at a semi-major axis of just $a_{B}=13.56_{-0.14}^{+0.18}$ au. Conclusions: The existence of the planet HD 59686 Ab in a tight eccentric binary system severely challenges standard giant planet formation theories and requires substantial improvements to such theories in tight binaries. Otherwise, alternative planet formation scenarios such as second generation planets or dynamical interactions in an early phase of the system's lifetime should be seriously considered in order to better understand the origin of this enigmatic planet.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Tilting Saturn without tilting Jupiter: Constraints on giant planet migration

The migration and encounter histories of the giant planets in our Solar System can be constrained by the obliquities of Jupiter and Saturn. We have performed secular simulations with imposed migration and N-body simulations with planetesimals to study the expected obliquity distribution of migrating planets with initial conditions resembling those of the smooth migration model, the resonant Nice model and two models with five giant planets initially in resonance (one compact and one loose configuration). For smooth migration, the secular spin-orbit resonance mechanism can tilt Saturn's spin axis to the current obliquity if the product of the migration time scale and the orbital inclinations is sufficiently large (exceeding 30 Myr deg). For the resonant Nice model with imposed migration, it is difficult to reproduce today's obliquity values, because the compactness of the initial system raises the frequency that tilts Saturn above the spin precession frequency of Jupiter, causing a Jupiter spin-orbit resonance crossing. Migration time scales sufficiently long to tilt Saturn generally suffice to tilt Jupiter more than is observed. The full N-body simulations tell a somewhat different story, with Jupiter generally being tilted as often as Saturn, but on average having a higher obliquity. The main obstacle is the final orbital spacing of the giant planets, coupled with the tail of Neptune's migration. The resonant Nice case is barely able to simultaneously reproduce the {orbital and spin} properties of the giant planets, with a probability ~0.15%. The loose five planet model is unable to match all our constraints (probability <0.08%). The compact five planet model has the highest chance of matching the orbital and obliquity constraints simultaneously (probability ~0.3%).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Complete Tidal Evolution of Pluto-Charon

Both Pluto and its satellite Charon have rotation rates synchronous with their orbital mean motion. This is the theoretical end point of tidal evolution where transfer of angular momentum has ceased. Here we follow Pluto's tidal evolution from an initial state having the current total angular momentum of the system but with Charon in an eccentric orbit with semimajor axis $a \approx 4R_P$ (where $R_P$ is the radius of Pluto), consistent with its impact origin. Two tidal models are used, where the tidal dissipation function $Q \propto$ 1/frequency and $Q=$ constant, where details of the evolution are strongly model dependent. The inclusion of the gravitational harmonic coefficient $C_{22}$ of both bodies in the analysis allows smooth, self consistent evolution to the dual synchronous state, whereas its omission frustrates successful evolution in some cases. The zonal harmonic $J_2$ can also be included, but does not cause a significant effect on the overall evolution. The ratio of dissipation in Charon to that in Pluto controls the behavior of the orbital eccentricity, where a judicious choice leads to a nearly constant eccentricity until the final approach to dual synchronous rotation. The tidal models are complete in the sense that every nuance of tidal evolution is realized while conserving total angular momentum - including temporary capture into spin-orbit resonances as Charon's spin decreases and damped librations about the same.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the Origin of Pluto's Small Satellites by Resonant Transport

The orbits of Pluto's four small satellites (Styx, Nix, Kerberos, and Hydra) are nearly circular and coplanar with the orbit of the large satellite Charon, with orbital periods nearly in the ratios 3:1, 4:1, 5:1, and 6:1 with Charon's orbital period. These properties suggest that the small satellites were created during the same impact event that placed Charon in orbit and had been pushed to their current positions by being locked in mean-motion resonances with Charon as Charon's orbit was expanded by tidal interactions with Pluto. Using the Pluto-Charon tidal evolution models developed by Cheng et al. (2014), we show that stable capture and transport of a test particle in multiple resonances at the same mean-motion commensurability is possible at the 5:1, 6:1, and 7:1 commensurabilities, if Pluto's zonal harmonic $J_{2P} = 0$. However, the test particle has significant orbital eccentricity at the end of the tidal evolution of Pluto-Charon in almost all cases, and there are no stable captures and transports at the 3:1 and 4:1 commensurabilities. Furthermore, a non-zero hydrostatic value of $J_{2P}$ destroys the conditions necessary for multiple resonance migration. Simulations with finite but minimal masses of Nix and Hydra also fail to yield any survivors. We conclude that the placing of the small satellites at their current orbital positions by resonant transport is extremely unlikely.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Precise radial velocities of giant stars VI. A possible 2:1 resonant planet pair around the K giant star $η$ Cet

We report the discovery of a new planetary system around the K giant $η$ Cet (HIP 5364, HD 6805) based on 118 high-precision optical radial velocities taken at Lick Observatory since July 2000. Since October 2011 an additional nine near-infrared Doppler measurements have been taken using the ESO CRIRES spectrograph (VLT, UT1). The visible data set shows two clear periodicities. Although we cannot completely rule out that the shorter period is due to rotational modulation of stellar features, the infrared data show the same variations as in the optical, which strongly supports that the variations are caused by two planets. Assuming the mass of $η$ Cet to be 1.7 $M_\odot$, the best edge-on coplanar dynamical fit to the data is consistent with two massive planets ($m_b\sin i$ = 2.6 $\pm$ 0.2 $M_{\mathrm{Jup}}$, $m_c\sin i$ = 3.3 $\pm$ 0.2 $M_{\mathrm{Jup}}$), with periods of $P_b$ = 407 $\pm$ 3 days and $P_c$ = 740 $\pm$ 5 days and eccentricities of $e_b$ = 0.12 $\pm$ 0.05 and $e_c$ = 0.08 $\pm$ 0.03. We tested a wide variety of edge-on coplanar and inclined planetary configurations for stability, which agree with the derived radial velocities. We find that in certain coplanar orbital configurations with moderate $e_b$ eccentricity, the planets can be effectively trapped in an anti-aligned 2:1 mean motion resonance. A much larger non-resonant stable region exists in low-eccentricity parameter space, although it appears to be much farther from the best fit than the 2:1 resonant region. In all other cases, the system is categorized as unstable or chaotic. Another conclusion from the coplanar inclined dynamical test is that the planets can be at most a factor of $\sim$ 1.4 more massive than their suggested minimum masses. This stability constraint on the inclination excludes the possibility of two brown dwarfs, and strongly favors a planetary system.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Detailed Analysis of the HD 73526 2:1 Resonant Planetary System

We present six years of new radial-velocity data from the Anglo-Australian and Magellan Telescopes on the HD 73526 2:1 resonant planetary system. We investigate both Keplerian and dynamical (interacting) fits to these data, yielding four possible configurations for the system. The new data now show that both resonance angles are librating, with amplitudes of 40 degrees and 60 degrees, respectively. We then perform long-term dynamical stability tests to differentiate these solutions, which only differ significantly in the masses of the planets. We show that while there is no clearly preferred system inclination, the dynamical fit with i=90 degrees provides the best combination of goodness-of-fit and long-term dynamical stability.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Are the Kepler Near-Resonance Planet Pairs due to Tidal Dissipation?

The multiple-planet systems discovered by the Kepler mission show an excess of planet pairs with period ratios just wide of exact commensurability for first-order resonances like 2:1 and 3:2. In principle, these planet pairs could have both resonance angles associated with the resonance librating if the orbital eccentricities are sufficiently small, because the width of first-order resonances diverges in the limit of vanishingly small eccentricity. We consider a widely-held scenario in which pairs of planets were captured into first-order resonances by migration due to planet-disk interactions, and subsequently became detached from the resonances, due to tidal dissipation in the planets. In the context of this scenario, we find a constraint on the ratio of the planet's tidal dissipation function and Love number that implies that some of the Kepler planets are likely solid. However, tides are not strong enough to move many of the planet pairs to the observed separations, suggesting that additional dissipative processes are at play.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Planetary Migration and Eccentricity and Inclination Resonances in Extrasolar Planetary Systems

The differential migration of two planets due to planet-disk interaction can result in capture into the 2:1 eccentricity-type mean-motion resonances. Both the sequence of 2:1 eccentricity resonances that the system is driven through by continued migration and the possibility of a subsequent capture into the 4:2 inclination resonances are sensitive to the migration rate within the range expected for type II migration. If the migration rate is fast, the resonant pair can evolve into a family of 2:1 eccentricity resonances different from those found by Lee (2004). This new family has outer orbital eccentricity e_2 > 0.4-0.5, asymmetric librations of both eccentricity resonance variables, and orbits that intersect if they are exactly coplanar. Although this family exists for an inner-to-outer planet mass ratio m_1/m_2 > 0.2, it is possible to evolve into this family by fast migration only for m_1/m_2 > 2. Thommes and Lissauer (2003) have found that a capture into the 4:2 inclination resonances is possible only for m_1/m_2 < 2. We show that this capture is also possible for m_1/m_2 > 2 if the migration rate is slightly slower than that adopted by Thommes and Lissauer. There is significant theoretical uncertainty in both the sign and the magnitude of the net effect of planet-disk interaction on the orbital eccentricity of a planet. If the eccentricity is damped on a timescale comparable to or shorter than the migration timescale, e_2 may not be able to reach the values needed to enter either the new 2:1 eccentricity resonances or the 4:2 inclination resonances. Thus, if future observations of extrasolar planetary systems were to reveal certain combinations of mass ratio and resonant configuration, they would place a constraint on the strength of eccentricity damping during migration, as well as on the migration rate itself.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

A Primordial Origin of the Laplace Relation Among the Galilean Satellites

Understanding the origin of the orbital resonances of the Galilean satellites of Jupiter will constrain the longevity of the extensive volcanism on Io, may explain a liquid ocean on Europa, and may guide studies of the dissipative properties of stars and Jupiter-like planets. The differential migration of the newly formed Galilean satellites due to interactions with a circumjovian disk can lead to the primordial formation of the Laplace relation n_1 - 3 n_2 + 2 n_3 = 0, where the n_i are the mean orbital angular velocities of Io, Europa, and Ganymede, respectively. This contrasts with the formation of the resonances by differential expansion of the orbits from tidal torques from Jupiter.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

A Survey of Numerical Solutions to the Coagulation Equation

We present the results of a systematic survey of numerical solutions to the coagulation equation for a rate coefficient of the form A_ij \propto (i^mu j^nu + i^nu j^mu) and monodisperse initial conditions. The results confirm that there are three classes of rate coefficients with qualitatively different solutions. For nu \leq 1 and lambda = mu + nu \leq 1, the numerical solution evolves in an orderly fashion and tends toward a self-similar solution at large time t. The properties of the numerical solution in the scaling limit agree with the analytic predictions of van Dongen and Ernst. In particular, for the subset with mu > 0 and lambda < 1, we disagree with Krivitsky and find that the scaling function approaches the analytically predicted power-law behavior at small mass, but in a damped oscillatory fashion that was not known previously. For nu \leq 1 and lambda > 1, the numerical solution tends toward a self-similar solution as t approaches a finite time t_0. The mass spectrum n_k develops at t_0 a power-law tail n_k \propto k^{-tau} at large mass that violates mass conservation, and runaway growth/gelation is expected to start at t_crit = t_0 in the limit the initial number of particles n_0 -> \infty. The exponent tau is in general less than the analytic prediction (lambda + 3)/2, and t_0 = K/[(lambda - 1) n_0 A_11] with K = 1--2 if lambda > 1.1. For nu > 1, the behaviors of the numerical solution are similar to those found in a previous paper by us. They strongly suggest that there are no self-consistent solutions at any time and that runaway growth is instantaneous in the limit n_0 -> \infty. They also indicate that the time t_crit for the onset of runaway growth decreases slowly toward zero with increasing n_0.

0 citations0 reviewsNo code linkedOpen access