Continuing our study of the effects of secular resonances on the formation of terrestrial planets in moderately close binary stars, we present here the results of an extensive numerical simulations of the formation of these objects. Considering a binary with two giant planets and a protoplanetary disk around its primary star, we have simulated the late stage of terrestrial planet formation for different types of the secondary, and different orbital elements of the binary and giant planets. Results demonstrate that terrestrial planet formation can indeed proceed constructively in such systems; however, as predicted by the general theory, secular resonances are suppressed and do not contribute to the formation process. Simulations show that it is in fact the mean-motion resonances of the inner giant planet that drive the dynamics of the protoplanetary disk and the mass and orbital architecture of the final bodies. Simulations also show that in the majority of the cases, the final systems contain only one terrestrial planet with a mass of 0.6-1.7 Earth masses. Multiple planets appear on rare occasions in the form of Earth-Mars analogs with the smaller planet in an exterior orbit. When
With four giant planets ($m\sim5-10~M_{\rm Jup}$, $T_\rm{eff}\sim900-1200$ K) orbiting between 15-70 au, HR 8799 provides an unparalleled testbed for studying giant planet formation and probing compositional trends across the protoplanetary disk. We present new JWST/NIRSpec IFU observations ($2.85-5.3~μ$m, $R\approx2700$) that now include the spectrum of HR 8799 b, and higher S/N spectra for HR 8799 c, d, and e compared to that in Ruffio & Xuan et al. We detect CO, CH$_4$, H$_2$O, H$_2$S, CO$_2$, and for planet b, NH$_3$. We combine the NIRSpec spectra with $1-5 μ$m photometry to perform atmospheric retrievals that account for disequilibrium chemistry and clouds, and allow C/H, O/H, N/H, and S/H to scale independently. While the four planets are similarly enriched in carbon and oxygen, with C/H and O/H between $3-5\times$ stellar, we observe a tentative trend of increasing S/H - a tracer of refractory solids - from $2-5 \times$ stellar with increasing orbital distance. From HR 8799 b's NH$_3$ abundance, we estimate $\rm N/H=21.2^{+16.2}_{-8.8}\times$ stellar, suggesting the outer planet accreted significant amounts of N-rich gas. Overall, the elemental abundance patterns we obs
Motivated by the diversity of circumstellar planets in binary stars and the strong effects of the secular resonances of Jupiter and Saturn on the formation and architecture of the inner solar system, we have launched an expansive project on studying the effects of secular resonances on the formation of terrestrial planets around a star of a moderately close binary. As the first phase of our project, we present here the general theory of secular resonances in dual-star systems where the primary hosts two giant planets. Using the concept of generalized disturbing function, we derive the formula for the locations of secular resonances and show that in systems where the perturbation of the secondary star is stronger, the locations of secular resonances are farther way from the primary and closer to the giant planets. The latter implies that in such systems, terrestrial planet formation has a larger area to proceed with more of the protoplanetary disk being available to it. To demonstrate the validity of our theoretical results, we simulated the evolution of a protoplanetary disk interior to the inner giant planet. Results, in addition to confirming our theoretical predictions, pointed
Arid terrestrial exoplanets are potentially abundant and are thus interesting targets in the search for life. In particular, M-dwarf planets such as those in the TRAPPIST-1 system may possess limited surface water, whereas early solar system terrestrials may have had small surface water inventories postmagma ocean solidification. On modern Earth, there is enough surface water for a balanced geologic carbon cycle, meaning silicate weathering balances the volcanic outgassing of CO2. However, on arid planets, there may not be enough surface water for this silicate weathering thermostat to maintain habitable conditions. Here, we show that arid planets enter a regime where weathering cannot keep up with volcanic degassing of CO2. Using a coupled model of the geologic carbon cycle, we find that terrestrial Earth-like planets require an initial surface water inventory of at least ~20-50% of Earth's ocean mass to maintain a balanced geologic carbon cycle and temperate surface temperature over 4.5 Gyr of evolution. Arid planets with less than ~20-50% of Earth's oceans cannot maintain high silicate weathering fluxes, potentially causing a runaway increase in atmospheric CO2. In addition, we
GJ 536 is a low-mass star, located 10 pc away from the Sun, that hosts a low-mass planet orbiting with a period of 8.71 days. Based on an analysis of the radial velocity (RV) time series obtained from the available data of the spectrographs HARPS, HARPS-N, CARMENES and HIRES, we announce the discovery of a second low-mass planet orbiting the star. We performed a RV global analysis on RV, spectroscopic activity indicators, and ASAS photometry, within the multidimensional Gaussian process framework, updated the parameters of GJ 536 b, and found significant evidence of the presence of a second planet. GJ 536 c is a low-mass planet ($m_{p} \sin i$ = 5.89 $\pm$ 0.70 M$_{\oplus}$), orbiting with a period of 32.761 $\pm$ 0.015 days, at a distance of 0.1617 $\pm$ 0.0028 au from its parent star. It induces an RV semi-amplitude of 1.80 $\pm$ 0.20 m$\cdot$s$^{-1}$. Given its distance to the star, it receives a flux of 1.692 $\pm$ 0.069 F$_{\oplus}$, for an equilibrium temperature of 290.5 $\pm$ 9.5 K. We update the mass of the planet GJ 536 b to $m_{p} \sin i$ = 6.37 $\pm$ 0.38 M$_{\oplus}$. The orbits of both planets are consistent with circular. We explored the use of statistical Doppler im
TOI-6255~b (GJ 4256) is an Earth-sized planet (1.079$\pm0.065$ $R_\oplus$) with an orbital period of only 5.7 hours. With the newly commissioned Keck Planet Finder (KPF) and CARMENES spectrographs, we determined the planet's mass to be 1.44$\pm$0.14 $M_{\oplus}$. The planet is just outside the Roche limit, with $P_{\rm orb}/P_{\rm Roche}$ = 1.13 $\pm0.10$. The strong tidal force likely deforms the planet into a triaxial ellipsoid with a long axis that is $\sim$10\% longer than the short axis. Assuming a reduced stellar tidal quality factor $Q_\star^\prime \approx10^7$, we predict that tidal orbital decay will cause TOI-6255 to reach the Roche limit in roughly 400 Myr. Such tidal disruptions may produce the possible signatures of planet engulfment that have been on stars with anomalously high refractory elemental abundances compared to its conatal binary companion. TOI-6255 b is also a favorable target for searching for star-planet magnetic interactions, which might cause interior melting and hasten orbital decay. TOI-6255 b is a top target (Emission Spectroscopy Metric of about 24) for phase curve observations with the James Webb Space Telescope.
We use the optical and near-infrared photometry from the Kepler Input Catalog to provide improved estimates of the stellar characteristics of the smallest stars in the Kepler target list. We find 3897 dwarfs with temperatures below 4000K, including 64 planet candidate host stars orbited by 95 transiting planet candidates. We refit the transit events in the Kepler light curves for these planet candidates and combine the revised planet/star radius ratios with our improved stellar radii to revise the radii of the planet candidates orbiting the cool target stars. We then compare the number of observed planet candidates to the number of stars around which such planets could have been detected in order to estimate the planet occurrence rate around cool stars. We find that the occurrence rate of 0.5-4 Earth radius planets with orbital periods shorter than 50 days is 0.90 (+0.04/-0.03) planets per star. The occurrence rate of Earth-size (0.5-1.4 Earth radius) planets is constant across the temperature range of our sample at 0.51 (+0.06/-0.05) Earth-size planets per star, but the occurrence of 1.4-4 Earth radius planets decreases significantly at cooler temperatures. Our sample includes 2 E
Atmospheric metal enrichment (i.e., elements heavier than helium, also called "metallicity") is a key diagnostic of the formation of giant planets. The giant planets of the solar system exhibit an inverse relationship between mass and both their bulk metallicities and atmospheric metallicities. Extrasolar giant planets also display an inverse relationship between mass and bulk metallicity. However, there is significant scatter in the relationship and it is not known how atmospheric metallicity correlates with either planet mass or bulk metallicity. Here we show that the Saturn-mass exoplanet HD 149026b has an atmospheric metallicity 59 - 276 times solar (at 1 $σ$), which is greater than Saturn's atmospheric metallicity of ~7.5 times solar at >4 $σ$ confidence. This result is based on modeling CO$_2$ and H$_2$O absorption features in the thermal emission spectrum of the planet measured by JWST. HD 149026b is the most metal-rich giant planet known, with an estimated bulk heavy element abundance of 66 $\pm$ 2% by mass. We find that the atmospheric metallicities of both HD 149026b and the solar system giant planets are more correlated with bulk metallicity than planet mass.
We use a multiannulus planetesimal accretion code to investigate the growth of icy planets in the outer regions of a planetesimal disk. In a quiescent minimum mass solar nebula, icy planets grow to sizes of 1000--3000 km on a timescale t = 15-20 Myr (a/30 AU)^3 where a is the distance from the central star. Planets form faster in more massive nebulae. Newly-formed planets stir up leftover planetesimals along their orbits and produce a collisional cascade where icy planetesimals are slowly ground to dust. The dusty debris of planet formation has physical characteristics similar to those observed in beta Pic, HR 4796A, and other debris disks. We derive dust masses for small particles, 1 mm and smaller, and large particles, 1 mm and larger, as a function of the initial conditions in the planetesimal disk. The dust luminosities derived from these masses are similar to those observed in Vega, HR 4796A, and other debris disks. The calculations produce bright rings and dark gaps. Bright rings occur where 1000 km and larger planets have recently formed. Dark gaps are regions where planets have cleared out dust or shadows where planets have yet to form.
We present the discovery of a Neptune-mass planet OGLE-2007-BLG-368Lb with a planet-star mass ratio of q=[9.5 +/- 2.1] x 10^{-5} via gravitational microlensing. The planetary deviation was detected in real-time thanks to the high cadence of the MOA survey, real-time light curve monitoring and intensive follow-up observations. A Bayesian analysis returns the stellar mass and distance at M_l = 0.64_{-0.26}^{+0.21} M_\sun and D_l = 5.9_{-1.4}^{+0.9} kpc, respectively, so the mass and separation of the planet are M_p = 20_{-8}^{+7} M_\oplus and a = 3.3_{-0.8}^{+1.4} AU, respectively. This discovery adds another cold Neptune-mass planet to the planetary sample discovered by microlensing, which now comprise four cold Neptune/Super-Earths, five gas giant planets, and another sub-Saturn mass planet whose nature is unclear. The discovery of these ten cold exoplanets by the microlensing method implies that the mass ratio function of cold exoplanets scales as dN_{\rm pl}/d\log q \propto q^{-0.7 +/- 0.2} with a 95% confidence level upper limit of n < -0.35 (where dN_{\rm pl}/d\log q \propto q^n). As microlensing is most sensitive to planets beyond the snow-line, this implies that Neptune-ma
The exoplanet sub-Neptune population currently poses a conundrum, as to whether small-size planets are volatile-rich cores without an atmosphere, or rocky cores surrounded by a H-He envelope. To test the different hypotheses from an observational point of view, a large sample of small-size planets with precise mass and radius measurements is the first step. On top of that, much more information will likely be needed, including atmospheric characterisation and a demographic perspective on their bulk properties. We present here the concept and strategy of the THIRSTEE project, which aims to shed light on the composition of the sub-Neptune population across stellar types by increasing their number and improving the accuracy of bulk density measurements, as well as investigating their atmospheres and performing statistical, demographic analysis. We report the first results of the program, characterising a new two-planet system around the M-dwarf TOI-406. We analyse TESS and ground-based photometry, together with ESPRESSO and NIRPS/HARPS RVs to derive the orbital parameters and investigate the internal composition of the 2 planets orbiting TOI-406, which have radii and masses of $R_c =
$K2$ greatly extended $Kepler$'s ability to find new planets, but it was typically limited to identifying transiting planets with orbital periods below 40 days. While analyzing $K2$ data through the Exoplanet Explorers project, citizen scientists helped discover one super-Earth and four sub-Neptune sized planets in the relatively bright ($V=12.21$, $K=10.3$) K2-138 system, all which orbit near 3:2 mean motion resonances. The $K2$ light curve showed two additional transit events consistent with a sixth planet. Using $Spitzer$ photometry, we validate the sixth planet's orbital period of $41.966\pm0.006$ days and measure a radius of $3.44^{+0.32}_{-0.31}\,R_{\oplus}$, solidifying K2-138 as the $K2$ system with the most currently known planets. There is a sizeable gap between the outer two planets, since the fifth planet in the system, K2-138 f, orbits at 12.76 days. We explore the possibility of additional non-transiting planets in the gap between f and g. Due to the relative brightness of the K2-138 host star, and the near resonance of the inner planets, K2-138 could be a key benchmark system for both radial velocity and transit timing variation mass measurements, and indeed radial v
For minor planet observations to be archived and used by the scientific community, observations of individual objects must be linked together into groups called tracklets. This linking is nontrivial, as linking software must find real tracklets in noisy data within a reasonable amount of time. We describe FindPOTATOs, a linking software written in Python that can assemble minor planet tracklets. With the appropriate parameters, FindPOTATOs assembles tracklets for a variety of objects, including close-approaching near-Earth objects and trans-Neptunian objects. FindPOTATOs is ideal for processing data sets taken at small observatories, processing archival data, or finding tracklets of minor planets on unusual orbits. This paper describes the code structure, usage, and validation.
Systems of two gravitationally bound exoplanets orbiting a common barycenter outside their physical radii ("binary planets") may result from tidal capture during planet-planet scattering. These objects are expected to form in tight orbits of just a few times their summed radii due to dynamical tides. As a result of their close proximity, their transits overlap heavily, leading to the deceptive illusion of a single planet of larger effective size, an effect compounded in the presence of noisy data and/or long integration times. We show that these illusory single-component transits, dubbed "chimera transits", exhibit large-amplitude Transit Duration Variation (TDV) effects on the order of hours, as well as smaller Transit Timing Variations (TTVs). We compute an analytic approximation for the transit duration upper bound, assuming binary planets with low impact parameter and orbits coplanar with the stellarcentric orbit. We verify the accuracy of our expressions against dynamical simulations of binary Jupiters using the luna algorithm, and provide a Python code for numerical calculations of the TDV signal in binary planet systems (github.com/joheenc/binary-planet-transits). Additional
This paper describes the Habitable Energy balance model for eXoplaneT ObseRvations (HEXTOR), which is a model for calculating latitudinal temperature profiles on Earth and other rapidly rotating planets. HEXTOR includes a lookup table method for calculating the outgoing infrared radiative flux and planetary albedo, which provides improvements over other approaches at parameterizing radiative transfer in an energy balance model. Validation cases are presented for present-day Earth and other Earth-sized planets with aquaplanet and land planet conditions from 0 to 45 degrees obliquity. A tidally locked coordinate system is also implemented in the energy balance model, which enables calculation of the horizontal temperature profile for planets in synchronous rotation around low mass stars. This coordinate transformed model is applied to cases for TRAPPIST-1e as defined by the TRAPPIST Habitable Atmosphere Intercomparison protocol, which demonstrates better agreement with general circulation models compared to the latitudinal energy balance model. Advances in applying energy balance models to exoplanets can be made by using general circulation models as a benchmark for tuning as well as
Magnetic interactions between a planet and its environment are known to lead to aurorae and shocks in the solar system. The large number of close-in exoplanets that have been discovered so far triggered a renewed interest in understanding magnetic interactions in other star-planet systems. Multiple magnetic effects were then unveiled, such as planet inflation or heating, planet migration, planetary material escape, and even some modifications of the host star apparent activity. Our goal here is to lay out the basic physical principles underlying star-planet magnetic interactions. We first briefly review the hot exoplanets' population as we know it. We then move to a general description of star-planet magnetic interactions, and finally focus on the fundamental concept of Alfvén wings and its implication for exosystems.
This review is based on lectures given at the 45th Saas-Fee Advanced Course 'From Protoplanetary Disks to Planet Formation' held in March 2015 in Les Diablerets, Switzerland. Starting with an overview of the main characterictics of the Solar System and extrasolar planets, we describe the planet formation process in terms of the sequential accretion scenario. First the growth processes of dust particles to planetesimals and subsequently to terrestrial planets or planetary cores are presented. This is followed by the formation process of the giant planets either by core accretion or gravitational instability. Finally, the dynamical evolution of the orbital elements as driven by disk-planet interaction and the overall evolution of multi-object systems is presented.
We investigate gas accretion onto a protoplanet, by considering the thermal effect of gas in three-dimensional hydrodynamical simulations, in which the wide region from a protoplanetary gas disk to a Jovian radius planet is resolved using the nested-grid method. We estimate the mass accretion rate and growth timescale of gas giant planets. The mass accretion rate increases with protoplanet mass for M_p < M_cri, while it becomes saturated or decreases for M_p > M_cri, where M_cri = 0.036 M_Jup (a_p/1AU)^0.75, and M_Jup and a_p are the Jovian mass and the orbital radius, respectively. The growth timescale of a gas giant planet or the timescale of the gas accretion onto the protoplanet is about 10^5 yr, that is two orders of magnitude shorter than the growth timescale of the solid core. The thermal effects barely affect the mass accretion rate because the gravitational energy dominates the thermal energy around the protoplanet. The mass accretion rate obtained in our local simulations agrees quantitatively well with those obtained in global simulations with coarser spatial resolution. The mass accretion rate is mainly determined by the protoplanet mass and the property of the pr
Thirteen exo-planets have been discovered using the gravitational microlensing technique (out of which 7 have been published). These planets already demonstrate that super-Earths (with mass up to ~10 Earth masses) beyond the snow line are common and multiple planet systems are not rare. In this White Paper we introduce the basic concepts of the gravitational microlensing technique, summarise the current mode of discovery and outline future steps towards a complete census of planets including Earth-mass planets. In the near-term (over the next 5 years) we advocate a strategy of automated follow-up with existing and upgraded telescopes which will significantly increase the current planet detection efficiency. In the medium 5-10 year term, we envision an international network of wide-field 2m class telescopes to discover Earth-mass and free-floating exo-planets. In the long (10-15 year) term, we strongly advocate a space microlensing telescope which, when combined with Kepler, will provide a complete census of planets down to Earth mass at almost all separations. Such a survey could be undertaken as a science programme on Euclid, a dark energy probe with a wide-field imager which has
We report the discovery of a planet with a high planet-to-star mass ratio in the microlensing event MOA-2009-BLG-387, which exhibited pronounced deviations over a 12-day interval, one of the longest for any planetary event. The host is an M dwarf, with a mass in the range 0.07 M_sun < M_host < 0.49M_sun at 90% confidence. The planet-star mass ratio q = 0.0132 +- 0.003 has been measured extremely well, so at the best-estimated host mass, the planet mass is m_p = 2.6 Jupiter masses for the median host mass, M = 0.19 M_sun. The host mass is determined from two "higher order" microlensing parameters. One of these, the angular Einstein radius θ_E = 0.31 +- 0.03 mas, is very well measured, but the other (the microlens parallax π_E, which is due to the Earth's orbital motion) is highly degenate with the orbital motion of the planet. We statistically resolve the degeneracy between Earth and planet orbital effects by imposing priors from a Galactic model that specifies the positions and velocities of lenses and sources and a Kepler model of orbits. The 90% confidence intervals for the distance, semi-major axis, and period of the planet are 3.5 kpc < D_L < 7.9 kpc, 1.1 AU < a