Search

We report the discovery of an ultrahot Jupiter with an extremely short orbital period of $0.67247414\,\pm\,0.00000028$ days ($\sim$16 hr). The $1.347 \pm 0.047$ $R_{\rm Jup}$ planet, initially identified by the Transiting Exoplanet Survey Satellite (TESS) mission, orbits TOI-2109 (TIC 392476080): a $T_{\rm eff} \sim 6500$ K F-type star with a mass of $1.447 \pm 0.077$ $M_{\rm Sun}$, a radius of $1.698 \pm 0.060$ $R_{\rm Sun}$, and a rotational velocity of $v\sin i_* = 81.9 \pm 1.7$ km s$^{-1}$. The planetary nature of TOI-2109b was confirmed through radial velocity measurements, which yielded a planet mass of $5.02 \pm 0.75$ $M_{\rm Jup}$. Analysis of the Doppler shadow in spectroscopic transit observations indicates a well-aligned system, with a sky-projected obliquity of $\lambda = 1\overset{\circ}{.}7 \pm 1\overset{\circ}{.}7$. From the TESS full-orbit light curve, we measured a secondary eclipse depth of $731 \pm 46$ ppm, as well as phase-curve variations from the planet's longitudinal brightness modulation and ellipsoidal distortion of the host star. Combining the TESS-band occultation measurement with a $K_s$-band secondary eclipse depth ($2012 \pm 80$ ppm) derived from ground-based observations, we find that the dayside emission of TOI-2109b is consistent with a brightness temperature of $3631 \pm 69$ K, making it the second hottest exoplanet hitherto discovered. By virtue of its extreme irradiation and strong planet-star gravitational interaction, TOI-2109b is an exceptionally promising target for intensive follow-up studies using current and near-future telescope facilities to probe for orbital decay, detect tidally driven atmospheric escape, and assess the impacts of H$_2$ dissociation and recombination on the global heat transport., Comment: 30 pages, 17 figures, published in AJ

A key process in planet formation is the exchange of angular momentum between a growing planet and the protoplanetary disc, which makes the planet migrate through the disc. Several works show that in general low-mass and intermediate-mass planets migrate towards the central star, unless corotation torques become dominant. Recently, a new kind of torque, called the thermal torque, was proposed as a new source that can generate outward migration of low-mass planets. While the Lindblad and corotation torques depend mostly on the properties of the protoplanetary disc and on the planet mass, the thermal torque depends also on the luminosity of the planet, arising mainly from the accretion of solids. Thus, the accretion of solids plays an important role not only in the formation of the planet but also in its migration process. In a previous work, we evaluated the thermal torque effects on planetary growth and migration mainly in the planetesimal accretion paradigm. In this new work, we study the role of the thermal torque within the pebble accretion paradigm. Computations are carried out consistently in the framework of a global model of planet formation that includes disc evolution, dust growth and evolution, and pebble formation. We also incorporate updated prescriptions of the thermal torque derived from high resolution hydrodynamical simulations. Our simulations show that the thermal torque generates extended regions of outward migration in low viscosity discs. This has a significant impact in the formation of the planets., Comment: Accepted for publication in MNRAS

Super-Earths and sub-Neptunes are commonly thought to have accreted hydrogen/helium envelopes, consisting of a few to ten percent of their total mass, from the primordial gas disk. Subsequently, hydrodynamic escape driven by core-powered mass-loss and/or photo-evaporation likely stripped much of these primordial envelopes from the lower-mass and closer-in planets to form the super-Earth population. In this work we show that after undergoing core-powered mass-loss, some super-Earths can retain small residual H/He envelopes. This retention is possible because, for significantly depleted atmospheres, the density at the radiative-convective boundary drops sufficiently such that thhe cooling time-scale becomes less than the mass-loss time-scale. The residual envelope is therefore able to contract, terminating further mass loss. Using analytic calculations and numerical simulations, we show that the mass of primordial H/He envelope retained as a fraction of the planet's total mass, $f_{ret}$, increases with increasing planet mass, $M_{c}$, and decreases with increasing equilibrium temperature, $T_{eq}$, scaling as $f_{ret} \propto M_{c}^{3/2} T_{eq}^{-1/2} \exp{[M_{c}^{3/4} T_{eq}^{-1}]}$. $f_{ret}$ varies from $< 10^{-8}$ to about $10^{-3}$ for typical super-Earth parameters. To first order, the exact amount of left-over H/He depends on the initial envelope mass, the planet mass, its equilibrium temperature, and the envelope's opacity. These residual hydrogen envelopes reduce the atmosphere's mean molecular weight compared to a purely secondary atmosphere, a signature observable by current and future facilities. These remnant atmospheres may, however, in many cases be vulnerable to long-term erosion by photo-evaporation. Any residual hydrogen envelope likely plays an important role in the long-term physical evolution of super-Earths, including their geology and geochemistry., 17 pages, 11 figures. Resubmitted to MNRAS after revision

Previous work concerning planet formation around low-mass stars has often been limited to large planets and individual systems. As current surveys routinely detect planets down to terrestrial size in these systems, a more holistic approach that reflects their diverse architectures is timely. Here, we investigate planet formation around low-mass stars and identify differences in the statistical distribution of planets. We compare the synthetic planet populations to observed exoplanets. We used the Generation III Bern model of planet formation and evolution to calculate synthetic populations varying the central star from solar-like stars to ultra-late M dwarfs. This model includes planetary migration, N-body interactions between embryos, accretion of planetesimals and gas, and long-term contraction and loss of the gaseous atmospheres. We find that temperate, Earth-sized planets are most frequent around early M dwarfs and more rare for solar-type stars and late M dwarfs. The planetary mass distribution does not linearly scale with the disk mass. The reason is the emergence of giant planets for M*>0.5 Msol, which leads to the ejection of smaller planets. For M*>0.3 Msol there is sufficient mass in the majority of systems to form Earth-like planets, leading to a similar amount of Exo-Earths going from M to G dwarfs. In contrast, the number of super-Earths and larger planets increases monotonically with stellar mass. We further identify a regime of disk parameters that reproduces observed M-dwarf systems such as TRAPPIST-1. However, giant planets around late M dwarfs such as GJ 3512b only form when type I migration is substantially reduced. We quantify the stellar mass dependence of multi-planet systems using global simulations of planet formation and evolution. The results compare well to current observational data and predicts trends that can be tested with future observations., Updated references, corrected inverted colorbar in Fig. 1. Data available on dace.unige.ch -> Formation & evolution

Context.Kepler-444 is one of the oldest planetary systems known thus far. Its peculiar configuration consisting of five sub-Earth-sized planets orbiting the companion to a binary stellar system makes its early history puzzling. Moreover, observations of HI-Lyαvariations raise many questions about the potential presence of escaping atmospheres today.Aims.We aim to study the orbital evolution of Kepler-444-d and Kepler-444-e and the impact of atmospheric evaporation on Kepler-444-e.Methods.Rotating stellar models of Kepler-444-A were computed with the Geneva stellar evolution code and coupled to an orbital evolution code, accounting for the effects of dynamical, equilibrium tides and atmospheric evaporation. The impacts of multiple stellar rotational histories and X-ray and extreme ultraviolet (XUV) luminosity evolutionary tracks are explored.Results.Using detailed rotating stellar models able to reproduce the rotation rate of Kepler-444-A, we find that its observed rotation rate is perfectly in line with what is expected for this old K0-type star, indicating that there is no reason for it to be exceptionally active as would be required to explain the observed HI-Lyαvariations from a stellar origin. We show that given the low planetary mass (~0.03 M⊕) and relatively large orbital distance (~0.06 AU) of Kepler-444-d and e, dynamical tides negligibly affect their orbits, regardless of the stellar rotational history considered. We point out instead how remarkable the impact is of the stellar rotational history on the estimation of the lifetime mass loss for Kepler-444-e. We show that, even in the case of an extremely slow rotating star, it seems unlikely that such a planet could retain a fraction of the initial water-ice content if we assume that it formed with a Ganymede-like composition.

We constrain the angular momentum architecture of HD 106906, a 13 $\pm$ 2 Myr old system in the ScoCen complex composed of a compact central binary, a widely separated planetary-mass tertiary HD 106906 b, and a debris disk nested between the binary and tertiary orbital planes. We measure the orientations of three vectors: the companion spin axis, companion orbit normal, and disk normal. Using near-IR high-resolution spectra from Gemini/IGRINS, we obtain a projected rotational velocity of $v\sin{i_p}$ = 9.5 $\pm$ 0.2 km/s for HD 106906 b. This measurement together with a published photometric rotation period implies the companion is viewed nearly pole-on, with a line-of-sight spin axis inclination of $i_p$ = 14 $\pm$ 4 degrees or 166 $\pm$ 4 degrees. By contrast, the debris disk is known to be viewed nearly edge-on. The likely misalignment of all three vectors suggests HD 106906 b formed by gravitational instability in a turbulent environment, either in a disk or cloud setting., Comment: accepted to AJ. 14 pages, 7 figures

A crucial phase during planetary growth is the migration when the planetary core has been assembled, but the planet did not open a deep gap yet. During this phase the planet is subject to fast type-I migration, which is mostly directed inwards, and the planet can lose a significant fraction of its semi-major axis. The duration of this phase is set by how long the planetary envelope needs to contract until it reaches a mass similar to the mass of the planetary core, which is when runaway gas accretion can set in and the planet can open a deeper gap in the disc, transitioning into the slower type-II migration. This envelope contraction phase depends crucially on the planetary mass and on the opacity inside the planetary envelope. Here we study how different opacity prescriptions influence the envelope contraction time and how this in turn influences how far the planet migrates through the disc. We find within our simulations that the size distribution of the grains as well as the chemical composition of the grains crucially influences how far the planet migrates before it can reach the runaway gas accretion phase. Grain size distributions with larger grain sizes result in less inward migration of the growing planet, due to faster gas accretion enabled by more efficient cooling. In addition, we find that planets forming in water poor environments can contract their envelope faster and thus migrate less, implying that gas giants forming in water poor environments might be located further away from their central star compared to gas giants forming in water rich environments. Future studies of planet formation that aim to investigate the chemical composition of formed gas giants need to take these effects self consistently into account., Comment: Accepted by A&A, 11 pages

Planets are thought to form at the early stage of stellar evolution when the mass accretion is still ongoing. RY Tau is a T Tauri type star at the age of a few Myr, with accretion disc seen at high inclination, so that line of sight crosses both the wind and the accretion gas flows. In a long series of spectroscopic monitoring of the star in 2013-2020, we detected variations in H-alpha and NaI D absorptions at radial velocities of infall (accretion) and outflow (wind) with a period of about 22 days. The absorptions in the infalling and the outflowing gas streams vary in anti-phase: an increase of infall is accompanied by a decrease of outflow, and vice versa. These flip-flop oscillations retain phase over several years of observations. We suggest that this may result from the MHD processes at the disk-magnetosphere boundary in the propeller mode. Another possibility is that a massive planet modulates some processes in the disc and provides the observed effects. The period, if Keplerian, corresponds to a distance of 0.2 AU, which is close to the dust sublimation radius in this star. The presence of the putative planet may be confirmed by radial velocity measurements: expected amplitude is > 90 m/s if a planet mass is > 2 Mj., 7 pages, 12 figures, accepted for publication in MNRAS

We present TOI-1259Ab, a 1.0 Rjup gas giant planet transiting a 0.71 Rsun K-dwarf on a 3.48 day orbit. The system also contains a bound white dwarf companion TOI-1259B with a projected distance of approximately 1600 AU from the planet host. Transits are observed in nine TESS sector and are 2.7 per cent deep - among the deepest known - making TOI-1259Ab a promising target for atmospheric characterization. Our follow-up radial velocity measurements indicate a variability of semiamplitude K = 71 m/s, implying a planet mass of 0.44 Mjup. By fitting the spectral energy distribution of the white dwarf we derive a total age of 4.08 (+1.21 -0.53) Gyr for the system. The K-dwarf's light curve reveals a rotational variability with a period of 28 days, which implies a gyrochronology age broadly consistent with the white dwarf's total age., Comment: Accepted to MNRAS. Some structural changes from first arXiv version but no significant changes to results. One figure got a bit prettier

Inferring planetary parameters from transit timing variations (TTVs) is challenging for small exoplanets because their transits may be so weak that determination of individual transit timing is difficult or impossible. We implement a useful combination of tools that together provide a numerically fast global photodynamical model. This is used to fit the TTV-bearing light curve, in order to constrain the masses of transiting exoplanets in low-eccentricity, multiplanet systems—and small planets in particular. We present inferred dynamical masses and orbital eccentricities in four multi-planet systems from Kepler's complete long-cadence data set. We test our model against Kepler-36/KOI-277, a system with some of the most precisely determined planetary masses through TTV inversion methods, and find masses of 5.56+0.41 −0.45 and 9.76+0.79 −0.89 for Kepler-36 b and c, respectively—consistent with literature in both value and error. We then improve the mass determination of the four planets in Kepler-79/KOI-152, where literature values were physically problematic to 12.5+4.5 −3.6, 9.5+2.3 −2.1, 11.3+2.2 −2.2 and 6.3+1.0 −1.0 for Kepler-79 b, c, d, and e, respectively. We provide new mass constraints where none existed before for two systems. These are 12.5+3.2 −2.6 for Kepler-450 c, and 3.3+1.7 −1.0 and 17.4+7.1 −3.8 for Kepler-595 c (previously KOI-547.03) and b, respectively. The photodynamical code used here, called PyDynamicaLC, is made publicly available.

The habitable zone is the main tool that mission architectures utilize to select potentially habitable planets for follow up spectroscopic observation. Given its importance, the precise size and location of the habitable zone remains a hot topic, as many studies, using a hierarchy of models, have assessed various factors including: atmospheric composition, time, and planetary mass. However, little work has assessed how the habitable zone changes with variations in background nitrogen pressure, which is directly connected to the habitability and life bearing potential of planets. Here, I use an advanced energy balance model with clouds to show that our solar system habitable zone is about 0.9 to 1.7 AU, assuming a 5 bar nitrogen background pressure and a maximum 100 percent cloud cover at the inner edge. This width is about 20 percent wider than the conservative habitable zone estimate. Similar extensions are calculated for A to M stars. I also show that cooling clouds and hazes and high background pressures can decrease the runaway greenhouse threshold temperature to approximately 300 K (or less) for planets orbiting any star type. This is because the associated increase in planetary albedo enables stable climates closer to the star, where rapid destabilization can be triggered from a lower mean surface temperature. Enhanced longwave emission for planets with very high stratospheric temperatures also permits stable climates at smaller orbital distances. The model predicts a runaway greenhouse above approximately 330 K for planets orbiting the Sun, which is consistent with previous work. However, moist greenhouses only occur for planets orbiting A-stars., Comment: Published in The Monthly Notices of the Royal Astronomical Society (27 pages, 7 figures)

The most accurate method for modelling planetary migration and hence the formation of resonant systems is using hydrodynamical simulations. Usually, the force (torque) acting on a planet is calculated using the forces from the gas disc and the star, while the gas accelerations are computed using the pressure gradient, the star, and the planet's gravity, ignoring its own gravity. For the non-migrating the neglect of the disc gravity results in a consistent torque calculation while for the migrating case it is inconsistent. We aim to study how much this inconsistent torque calculation can affect the final configuration of a two-planet system. Our focus will be on low-mass planets because most of the multi-planetary systems, discovered by the Kepler survey, have masses around 10 Earth masses. Performing hydrodynamical simulations of planet-disc interaction, we measure the torques on non-migrating and migrating planets for various disc masses as well as density and temperature slopes with and without considering the disc self-gravity. Using this data, we find a relation that quantifies the inconsistency, use it in an N-body code, and perform an extended parameter study modelling the migration of a planetary system with different planet mass ratios and disc surface densities, in order to investigate the impact of the torque inconsistency on the architecture of the planetary system. Not considering disc self-gravity produces an artificially larger torque on the migrating planet that can result in tighter planetary systems. The deviation of this torque from the correct value is larger in discs with steeper surface density profiles. In hydrodynamical modelling of multi-planetary systems, it is crucial to account for the torque correction, otherwise the results favour more packed systems., Comment: 7 pages, 7 figures, accepted for publication in A&A

We have used infrared images obtained with the Wide Field Camera 3 on board the Hubble Space Telescope to search for planetary-mass brown dwarfs in the star-forming cluster IC 348. In those images, we have identified 12 objects that have colors indicative of spectral types later than M8, corresponding to masses of, Astronomical Journal, in press

It is widely known that giant planets have the capacity to open deep gaps in their natal gaseous protoplanetary discs. It is unclear, however, how gas accretion onto growing planets influences the shape and depth of their growing gaps. We performed isothermal hydrodynamical simulations with the Fargo-2D1D code, which assumes planets accreting gas within full discs that range from 0.1 to 260 AU. The gas accretion routine uses a sink cell approach, in which different accretion rates are used to cope with the broad range of gas accretion rates cited in the literature. We find that the planetary gas accretion rate increases for larger disc aspect ratios and greater viscosities. Our main results show that gas accretion has an important impact on the gap-opening mass: we find that when the disc responds slowly to a change in planetary mass (i.e., at low viscosity), the gap-opening mass scales with the planetary accretion rate, with a higher gas accretion rate resulting in a larger gap-opening mass. On the other hand, if the disc response time is short (i.e., at high viscosity), then gas accretion helps the planet carve a deep gap. As a consequence, higher planetary gas accretion rates result in smaller gap-opening masses. Our results have important implications for the derivation of planet masses from disc observations: depending on the planetary gas accretion rate, the derived masses from ALMA observations might be off by up to a factor of two. We discuss the consequences of the change in the gap-opening mass on the evolution of planetary systems based on the example of the grand tack scenario. Planetary gas accretion also impacts stellar gas accretion, where the influence is minimal due to the presence of a gas-accreting planet.

We predict magnitudes for young planets embedded in transition discs, still affected by extinction due to material in the disc. We focus on Jupiter-size planets at a late stage of their formation, when the planet has carved a deep gap in the gas and dust distributions and the disc starts being transparent to the planet flux in the infrared (IR). Column densities are estimated by means of three-dimensional hydrodynamical models, performed for several planet masses. Expected magnitudes are obtained by using typical extinction properties of the disc material and evolutionary models of giant planets. For the simulated cases located at $5.2$ AU in a disc with local unperturbed surface density of $127$ $\mathrm{g} \cdot \mathrm{cm}^{-2}$, a $1$ $M_J$ planet is highly extincted in J-, H- and K-bands, with predicted absolute magnitudes $\ge 50$ mag. In L- and M-bands extinction decreases, with planet magnitudes between $25$ and $35$ mag. In the N-band, due to the silicate feature on the dust opacities, the expected magnitude increases to $40$ mag. For a $2$ $M_J$ planet, the magnitudes in J-, H- and K-bands are above $22$ mag, while for L-, M- and N-bands the planet magnitudes are between $15$ and $20$ mag. For the $5$ $M_J$ planet, extinction does not play a role in any IR band, due to its ability to open deep gaps. Contrast curves are derived for the transition discs in CQ Tau, PDS70, HL Tau, TW Hya and HD163296. Planet mass upper-limits are estimated for the known gaps in the last two systems., Comment: Accepted for publication on January 8, 2020 in MNRAS. 15 pages of main text with 14 figures, and 5 pages of appendices A and B with 4 figures

Recent exoplanet statistics indicate that photo-evaporation has a great impact on the mass and bulk composition of close-in low-mass planets. While there are many studies addressing photo-evaporation of hydrogen-rich or water-rich atmospheres, no detailed investigation regarding rocky vapor atmospheres (or mineral atmospheres) has been conducted. Here, we develop a new 1-D hydrodynamic model of the UV-irradiated mineral atmosphere composed of Na, Mg, O, Si, their ions and electrons, includin molecular diffusion, thermal conduction, photo-/thermo-chemistry, X--ray and UV heating, and radiative line cooling (i.e., the effects of the optical thickness and non-LTE). The focus of this paper is on describing our methodology but presents some new findings. Our hydrodynamic simulations demonstrate that almost all of the incident X-ray and UV energy from the host-star is converted into and lost by the radiative emission of the coolant gas species such as Na, Mg, Mg$^+$, Si$^{2+}$, Na$^{3+}$ and Si$^{3+}$. For an Earth-size planet orbiting 0.02~AU around a young solar-type star, we find that the X-ray and UV heating efficiency is as small as $1 \times 10^{-3}$, which corresponds to 0.3~$\Mearth$/Gyr of the mass loss rate simply integrated over all the directions. Because of such efficient cooling, the photo-evaporation of the mineral atmosphere on hot rocky exoplanets with masses of $1\Mearth$ is not massive enough to exert a great influence on the planetary mass and bulk composition. This suggests that close-in high-density exoplanets with sizes larger than the Earth radius survive in the high-UV environments., Comment: 23 pages, accepted for publication in the Monthly Notices of the Royal Astronomical Society

Comparing chemical abundances of a planet and the host star reveals the origin and formation path. Stellar abundance is measured with high-resolution spectroscopy. Planet abundance, on the other hand, is usually inferred from low-resolution data. For directly imaged exoplanets, the data are available from a slew of high-contrast imaging/spectroscopy instruments. Here, we study the chemical abundance of HR 8799 and its planet c. We measure stellar abundance using LBT/PEPSI (R=120,000) and archival HARPS data: stellar [C/H], [O/H], and C/O are 0.11$\pm$0.12, 0.12$\pm$0.14, and 0.54$^{+0.12}_{-0.09}$, all consistent with solar values. We conduct atmospheric retrieval using newly obtained Subaru/CHARIS data together with archival Gemini/GPI and Keck/OSIRIS data. We model the planet spectrum with petitRADTRANS and conduct retrieval using PyMultiNest. Retrieved planetary abundance can vary by $\sim$0.5 dex, from sub-stellar to stellar C and O abundances. The variation depends on whether strong priors are chosen to ensure a reasonable planet mass. Moreover, comparison with previous works also reveals inconsistency in abundance measurements. We discuss potential issues that can cause the inconsistency, e.g., systematics in individual data sets and different assumptions in the physics and chemistry in retrieval. We conclude that no robust retrieval can be obtained unless the issues are fully resolved., Comment: Published on AJ

Observations of bright protoplanetary disks often show annular gaps in their dust emission. One interpretation of these gaps is disk-planet interaction. If so, fitting models of planetary gaps to observed protoplanetary disk gaps can reveal the presence of hidden planets. However, future surveys are expected to produce an ever-increasing number of protoplanetary disks with gaps. In this case, performing a customized fitting for each target becomes impractical owing to the complexity of disk-planet interaction. To this end, we introduce DPNNet (Disk Planet Neural Network), an efficient model of planetary gaps by exploiting the power of machine learning. We train a deep neural network with a large number of dusty disk-planet hydrodynamic simulations across a range of planet masses, disk temperatures, disk viscosities, disk surface density profiles, particle Stokes numbers, and dust abundances. The network can then be deployed to extract the planet mass for a given gap morphology. In this work, first in a series, we focus on the basic concepts of our machine learning framework. We demonstrate its utility by applying it to the dust gaps observed in the protoplanetary disk around HL Tau at $10$ au, $30$ au, and $80$ au. Our network predict planet masses of $80 \, M_{\rm Earth}$, $63 \, M_{\rm Earth}$, and $70 \, M_{\rm Earth}$, respectively, which are comparable to other studies based on specialized simulations. We discuss the key advantages of our DPNNet in its flexibility to incorporate new physics, any number of parameters and predictions, and its potential to ultimately replace hydrodynamical simulations for disk observers and modelers., Comment: 16 Pages, 8 Figures, to appear in ApJ

Gravitational waves can probe the existence of planetary-mass primordial black holes. Considering a mass range of [10−7−10−2]M⊙, inspiraling primordial black holes could emit either continuous gravitational waves, quasi-monochromatic signals that last for many years, or transient continuous waves, signals whose frequency evolution follows a power law and last for O(hours-months). We show that primordial black hole binaries in our galaxy may produce detectable gravitational waves for different mass functions and formation mechanisms. In order to detect these inspirals, we adapt methods originally designed to search for gravitational waves from asymmetrically rotating neutron stars. The first method, the Frequency-Hough, exploits the continuous, quasi-monochromatic nature of inspiraling black holes that are sufficiently light and far apart such that their orbital frequencies can be approximated as linear with a small spin-up. The second method, the Generalized Frequency-Hough, drops the assumption of linearity and allows the signal frequency to follow a power-law evolution. We explore the parameter space to which each method is sensitive, derive a theoretical sensitivity estimate, determine optimal search parameters and calculate the computational cost of all-sky and directed searches. We forecast limits on the abundance of primordial black holes within our galaxy, showing that we can constrain the fraction of dark matter that primordial black holes compose, fPBH, to be fPBH≲1 for chirp masses between [4×10−5−10−3]M⊙ for current detectors. For the Einstein Telescope, we expect the constraints to improve to fPBH≲10−2 for chirp masses between [10−4−10−3]M⊙., SCOPUS: ar.j, info:eu-repo/semantics/published

A number of protoplanetary disks observed with ALMA potentially provide direct examples of initial conditions for planetary systems. In particular, the HL Tau disk has been intensively studied, and its rings/gaps are conventionally interpreted to be a result of unseen massive planets embedded in the gaps. Based on this interpretation, we carried out N-body simulations to investigate orbital evolution of planets within the protoplanetary disk and after the disk dispersal. Before the disk dispersal, our N-body simulations include both migration and mass-growth of the planet coupled with evolution of the disk. By varying the disk parameters, we produce a variety of widely-separated planetary systems consisting of three super-Jupiters at the end of disk dispersal. We found the outer planet is more massive than the inner one, and the migration of the innermost planet is inefficient due to the accretion of outer planet(s). We also showed how the final configuration and the final planetary mass depend on disk parameters. The migration is found to be convergent and no planet-pair has a period ratio less than 2. After the disk dispersal, we switch to pure gravitational N-body simulations and integrate the orbits up to 10 Gyr. Most simulated systems remain stable for at least 10 Gyr. We discuss implications of our result in terms of the observed widely-separated planetary systems HR 8799 and PDS 70., Comment: 16 pages, 9 figures. Accepted for publication in the ApJ

Some Jupiter-mass exoplanets contain $\sim$$100\, M_\oplus$ of metals, well above the $\sim$$10\, M_\oplus$ typically needed in a solid core to trigger giant planet formation by runaway gas accretion. We demonstrate that such `heavy-metal Jupiters' can result from planetary mergers near $\sim$10 au. Multiple cores accreting gas at runaway rates gravitationally perturb one another onto crossing orbits such that the average merger rate equals the gas accretion rate. Concurrent mergers and gas accretion implies the core mass scales with the total planet mass as $M_{\rm core} \propto M^{1/5}$ - heavier planets harbour heavier cores, in agreement with the observed mass-metallicity relation. While the average gas giant merges about once to double its core, others may merge multiple times, as merger trees grow chaotically. We show that the dispersion of outcomes inherent in mergers can reproduce the large scatter in observed planet metallicities, assuming $3-30\, M_\oplus$ pre-runaway cores. Mergers potentially correlate metallicity, eccentricity, and spin., Comment: Submitted to MNRAS

We report the discovery of TOI 837b and its validation as a transiting planet. We characterize the system using data from the NASA TESS mission, the ESA Gaia mission, ground-based photometry from El Sauce and ASTEP400, and spectroscopy from CHIRON, FEROS, and Veloce. We find that TOI 837 is a $T=9.9$ mag G0/F9 dwarf in the southern open cluster IC 2602. The star and planet are therefore $35^{+11}_{-5}$ million years old. Combining the transit photometry with a prior on the stellar parameters derived from the cluster color-magnitude diagram, we find that the planet has an orbital period of $8.3\,{\rm d}$ and is slightly smaller than Jupiter ($R_{\rm p} = 0.77^{+0.09}_{-0.07} \,R_{\rm Jup}$). From radial velocity monitoring, we limit $M_{\rm p}\sin i$ to less than 1.20 $M_{\rm Jup}$ (3-$\sigma$). The transits either graze or nearly graze the stellar limb. Grazing transits are a cause for concern, as they are often indicative of astrophysical false positive scenarios. Our follow-up data show that such scenarios are unlikely. Our combined multi-color photometry, high-resolution imaging, and radial velocities rule out hierarchical eclipsing binary scenarios. Background eclipsing binary scenarios, though limited by speckle imaging, remain a 0.2% possibility. TOI 837b is therefore a validated adolescent exoplanet. The planetary nature of the system can be confirmed or refuted through observations of the stellar obliquity and the planetary mass. Such observations may also improve our understanding of how the physical and orbital properties of exoplanets change in time., Comment: AJ accepted. Figure 11 is my favorite. Comments welcome!

We present a velocimetric and spectropolarimetric analysis of 27 observations of the 22-Myr M1 star AU Microscopii (Au Mic) collected with the high-resolution $YJHK$ (0.98-2.35 $\mu$m) spectropolarimeter SPIRou from 2019 September 18 to November 14. Our radial velocity (RV) time-series exhibits activity-induced fluctuations of 45 m/s RMS, about three times smaller than those measured in the optical domain, that we filter using Gaussian Process Regression. We report a 3.9$\sigma$-detection of the recently-discovered 8.46-d transiting planet AU Mic b, with an estimated mass of $17.1^{+4.7}_{-4.5}$ M$_{\odot}$ and a bulk density of $1.3 \pm 0.4$ g/cm$^{-3}$, inducing a RV signature of semi-amplitude $K=8.5^{+2.3}_{-2.2}$ m/s in the spectrum of its host star. A consistent detection is independently obtained when we simultaneously image stellar surface inhomogeneities and estimate the planet parameters with Zeeman-Doppler Imaging (ZDI). Using ZDI, we invert the time series of unpolarized and circularly-polarized spectra into surface brightness and large-scale magnetic maps. We find a mainly poloidal and axisymmetric field of 475 G, featuring, in particular, a dipole of 450 G tilted at 19{\deg} to the rotation axis. Moreover, we detect a strong differential rotation of d$\Omega = 0.167 \pm 0.009$ rad/d shearing the large-scale field, about twice stronger than that shearing the brightness distribution, suggesting that both observables probe different layers of the convective zone. Even though we caution that more RV measurements are needed to accurately pin down the planet mass, AU Mic b already appears as a prime target for constraining planet formation models, studying the interactions with the surrounding debris disk, and characterizing its atmosphere with upcoming space- and ground-based missions., Comment: 19 pages, 20 figures, MNRAS, in press

The origin of the observed diversity of planetary system architectures is one of the main topic of the exoplanetary research. The detection of a statistically significant sample of planets around young stars allows us to study the early stages of planet formation and evolution, but only a handful of them is known so far. In this regard, a considerable contribution is expected from the NASA TESS satellite, which is now performing a survey of $\sim 85 \%$ of the sky to search for short-period transiting planets In its first month of operations, TESS found a planet candidate with an orbital period of 8.14 days around a member of the Tuc-Hor young association ($\sim$ 40 Myr), the G6V main component of the binary system DS\,Tuc. If confirmed, it would be the first transiting planet around a young star suitable for radial velocity and/or atmospheric characterization. We aim to validate the planetary nature of this companion and to measure its orbital and physical parameters. We obtain accurate planet parameters by coupling an independent reprocessing of the TESS light curve with improved stellar parameters and the dilution caused by the binary companion; we analyse high precision archival radial velocities to impose an upper limit of about 0.1 M$_{\rm Jup}$ on the planet mass; we finally rule out the presence of external companions beyond 40 au with adaptive optics images. We confirm the presence of a young, giant (${\rm R} = 0.50$ R$_{\rm Jup}$) planet having a not negligible possibility to be inflated (theoretical mass $\lesssim 20$ M$_{\oplus}$) around DS\,Tuc~A. We discuss the feasibility of mass determination, Rossiter-McLaughlin analysis and atmosphere characterization, allowed by the brightness of the star., Comment: 12 pages, 11 figures. Accepted for publication in A&A

We explore the probable chemical signature of planet formation in the remarkable binary system HD 106515. The A star hosts a massive long-period planet with 9 MJup detected by radial velocity. We also refine stellar and planetary parameters by using non-solar-scaled opacities when modeling the stars. Methods. We carried out a simultaneous determination of stellar parameters and abundances, by applying for the first time non-solar-scaled opacities in this binary system, in order to reach the highest possible precision. Results. The stars A and B in the binary system HD 106515 do not seem to be depleted in refractory elements, which is different when comparing the Sun with solar-twins. Then, the terrestrial planet formation would have been less efficient in the stars of this binary system. Together with HD 80606/7, this is the second binary system which does not seem to present a (terrestrial) signature of planet formation, and hosting both systems an eccentric giant planet. This is in agreement with numerical simulations, where the early dynamical evolution of eccentric giant planets clear out most of the possible terrestrial planets in the inner zone. We refined the stellar mass, radius and age for both stars and found a notable difference of 78% in R compared to previous works. We also refined the planet mass to mp sini = 9.08 +/- 0.20 MJup, which differs by 6% compared with literature. In addition, we showed that the non-solar-scaled solution is not compatible with the classical solar-scaled method, and some abundance differences are comparable to NLTE or GCE effects specially when using the Sun as reference. Then, we encourage the use of non-solar-scaled opacities in high-precision studies such as the detection of Tc trends.[abridged], Comment: 9 pages, 10 figures, A&A accepted. arXiv admin note: text overlap with arXiv:1507.08125

As of today ten circumbinary planets orbiting solar type main sequence stars have been discovered. Nearly all of them orbit around the central binary very closely to the region of instability where it is difficult to form them in situ. It is assumed that they formed further out and migrated to their observed position. We extend previous studies to a more realistic thermal disc structure and determine what parameter influence the final parking location of a planet around a binary star. We perform two-dimensional numerical simulations of viscous accretion discs around a central binary that include viscous heating and radiative cooling from the disc surfaces. We vary the binary eccentricity as well as disc viscosity and mass. Concerning the disc evolution we find that it can take well over 100000 binary orbits until an equilibrium state is reached. As seen previously, we find that the central cavity opened by the binary becomes eccentric and precesses slowly in a prograde sense. Embedded planets migrate to the inner edge of the disc. In cases of lower disc viscosity they migrate further in maintaining a circular orbit, while for high viscosity they are parked further out on an eccentric orbit. The final location of an embedded planet is linked to its ability to open a gap in the disc. Gap opening planets separate inner from outer disc, preventing eccentricity excitation in the latter and making it more circular. This allows embedded planets to migrate closer to the binary, in agreement with the observations. The necessary condition for gap opening and the final planet position depend on the planet mass and disc viscosity., Comment: 15 pages, 16 figures, accepted by Astronomy and Astrophysics

We observe the dust continuum at 225 GHz and CO isotopologue (12CO, 13CO, and C18O) J=2-1 emission lines toward the CR Cha protoplanetary disk using the Atacama Large Millimeter/Submillimeter Array (ALMA). The dust continuum image shows a dust gap-ring structure in the outer region of the dust disk. A faint dust ring is also detected around 120 au beyond the dust gap. The CO isotopologue lines indicate that the gas disk is more extended than the dust disk. The peak brightness temperature of the 13CO line shows a small bump around 130 au while 12CO and C18O lines do not show. We investigate two possible mechanisms for reproducing the observed dust gap-ring structure and a gas temperature bump. First, the observed gap structure can be opened by a Jupiter mass planet using the relation between the planet mass and the gap depth and width. Meanwhile, the radiative transfer calculations based on the observed dust surface density profile show that the observed dust ring could be formed by dust accumulation at the gas temperature bump, that is, the gas pressure bump produced beyond the outer edge of the dust disk., Comment: 21 pages, figures

Several detections of wide-orbit planet-mass/sub-stellar companions around young solar-like stars were reported in the last decade. The origin of those possible planets is still unclear but accretion tracers and VLT/SPHERE observations indicate that they are surrounded by circumplanetary material or even a circumplanetary disk. We want to investigate if the gas component of disks around wide-orbit companions is detectable with current and future (sub)mm telescopes and what constraints such gas observations can provide on the nature of the circumplanetary material and on the mass of the companion. We applied the radiation thermo-chemical disk code ProDiMo to model the dust and gas component of passive circumplanetary disks and produced realistic synthetic observables. We considered different companion properties, disk parameters and radiative environments and compared the resulting synthetic observables to telescope sensitivities and to existing dust observations. The main criterion for a successful detection is the size of the circumplanetary disk. At a distance of about 150 pc, a circumplanetary disk with an outer radius of about 10 au is detectable with ALMA in about 6 hours in optically thick CO lines. Other aspects such as the companion's luminosity, disk inclination and background radiation fields are also relevant and should be considered to optimize the observing strategy for detection experiments. For most of the known wide-orbit planet-mass companions, their maximum theoretical disk size of one third of the Hill radius would be sufficient to allow detection of CO lines. It is therefore feasible to detect their gas disks and constrain the mass of the companion through the kinematic signature. Even in the case of non-detections such observations will provide stringent constraints on disk size and gas mass, information crucial for formation theories., Comment: 13 pages, 11 figures, accepted for publication in A&A

While thousands of exoplanets have been confirmed, the known properties about individual discoveries remain sparse and depend on detection technique. To utilize more than a small section of the exoplanet dataset, tools need to be developed to estimate missing values based on the known measurements. Here, we demonstrate the use of a neural network that models the density of planets in a space of six properties that is then used to impute a probability distribution for missing values. Our results focus on planetary mass which neither the radial velocity nor transit techniques for planet identification can provide alone. The neural network can impute mass across the four orders of magnitude in the exoplanet archive, and return a distribution of masses for each planet that can inform about trends in the underlying dataset. The average error on this mass estimate from a radial velocity detection is a factor of 1.5 of the observed value, and 2.7 for a transit observation. The mass of Proxima Centauri b found by this method is $1.6^{\rm +0.46}_{\rm -0.36}$ M$_\oplus$, where the upper and lower bounds are derived from the root mean square deviation from the log mass probability distribution. The network can similarly impute the other potentially missing properties, and we use this to predict planet radius for radial velocity measurements, with an average error of a factor 1.4 of the observed value. The ability of neural networks to search for patterns in multidimensional data means that such techniques have the potential to greatly expand the use of the exoplanet catalogue., Comment: 17 page, 15 figures, accepted for publication in AJ

We present adaptive optics imaging from the NIRC2 instrument on the Keck-2 telescope that resolves the exoplanet host (and lens) star as it separates from the brighter source star. These observations yield the $K$-band brightness of the lens and planetary host star, as well as the lens-source relative proper motion, $��_{\rm rel,H}$. in the heliocentric reference frame. The $��_{\rm rel,H}$ measurement allows determination of the microlensing parallax vector, $��_E$, which had only a single component determined by the microlensing light curve. The combined measurements of $��_{\rm rel,H}$ and $K_L$ provide the masses of the host stat, $M_{\rm host} = 0.426\pm 0.037 M_\odot$, and planet, $m_p = 3.27 \pm 0.32 M_{\rm Jup}$ with a projected separation of $3.4\pm 0.5\,$AU. This confirms the tentative conclusion of a previous paper (Dong et al. 2009) that this super-Jupiter mass planet, OGLE-2005-BLG-071Lb, orbits an M-dwarf. Such planets are predicted to be rare by the core accretion theory and have been difficult to find with other methods, but there are two such planets with firm mass measurements from microlensing, and an additional 11 planetary microlens events with host mass estimates $< 0.5M_\odot$ and planet mass estimates $> 2$ Jupiter masses that could be confirmed by high angular follow-up observations. We also point out that OGLE-2005-BLG-071L has separated far enough from its host star that it should be possible to measure the host star metallicity withspectra from a high angular resolution telescope such as Keck, the VLT, the Hubble Space Telescope or the James Webb Space Telescope., accepted by AJ

Gas giants are thought to form by runaway accretion: an instability driven by the self-gravity of growing atmospheres that causes accretion rates to rise super-linearly with planet mass. Why runaway should stop at a Jupiter or any other mass is unknown. We consider the proposal that final masses are controlled by circumstellar disc gaps (cavities) opened by planetary gravitational torques. We develop a fully time-dependent theory of gap formation and couple it self-consistently to planetary growth rates. When gaps first open, planetary torques overwhelm viscous torques, and gas depletes as if it were inviscid. In low-viscosity discs, of the kind motivated by recent observations and theory, gaps stay predominantly in this inviscid phase and planet masses finalize at $M_{\rm final}/M_\star\sim(\Omega t_{\rm disc})^{0.07}(H/a)^{2.73}(G\rho_0/\Omega^2)^{1/3}$, with $M_\star$ the host stellar mass, $\Omega$ the planet's orbital angular velocity, $t_{\rm disc}$ the gas disc's lifetime, $H/a$ its aspect ratio, and $\rho_0$ its unperturbed density. This final mass is independent of the dimensionless viscosity $\alpha$ and applies to large orbital distances, typically beyond $\sim$10 AU, where disc scale heights exceed planet radii. It evaluates to a few Jupiter masses at 10-100 AU, increasing gradually with distance as gaps become harder to open., Comment: Accepted to MNRAS

We present our observational results of the 1.1 mm continuum and the HCO+ (3-2) line in HL Tau at angular resolutions of 0.1" obtained with ALMA and our data analysis of the 2.9 mm and 1.1 mm continuum and the HCO+ (3-2) and (1-0) lines of the HL Tau disk. The Keplerian rotation of the HL Tau disk is well resolved in the HCO+ (3-2) emission, and the stellar mass is estimated to be 2.1+/-0.2 Msun with a disk inclination angle of 47 deg. The radial profiles of the HCO+ column density and excitation temperature are measured with the LTE analysis of the two transitions of the HCO+ emission. An HCO+ gas gap at a radius of 30 au, where the column density drops by a factor of 4-8, is found in the HCO+ column density profile, coincident with the dust gap traced by the continuum emission. No other clear HCO+ gas gap is seen. This HCO+ gas gap can be opened by a planet with mass of 0.5-0.8 Mj, which is comparable to the planet mass adopted in numerical simulations to form the dust gap at the same radius in the HL Tau disk. In addition to the disk component, a one-arm spiral with a length of ~3" (520 au) stretching out from the inner disk is observed in the HCO+ (3-2) emission. The observed velocity structures along the spiral suggest an infalling and rotating gas stream toward the inner disk., Comment: Accepted by ApJ

As of today, hundreds of hot Jupiters have been found, yet the inflated radii of a large fraction of them remain unexplained. It is still unclear whether a single inflation mechanism is enough to explain the entire distribution of radii, or whether a combination of them is needed. We seek to understand the relationship of hot Jupiter radii with stellar irradiation and mass. We also aim to find the intrinsic physical scatter in their radii caused by unobservable parameters. By constructing a hierarchical Bayesian model, we infer the probabilistic relation between planet radius, mass and incident flux for a sample of 286 gas giants. We separately incorporate the observational uncertainties of the data and the intrinsic physical scatter in the population. This allows us to treat the intrinsic physical scatter in radii (due to latent parameters such as the heavy element fraction) as a parameter to be inferred. We find that the planetary mass plays a key role in the inflation extent, with planets in the range between 0.37 - 0.98 MJ showing the most inflated radii. At higher masses, the radius response to incident flux begins to decrease. Below a threshold of 0.37 +- 0.03 MJ we find that giant exoplanets as a population are unable to maintain inflated radii greater than 1.4 RJ but instead exhibit smaller sizes as the incident flux is increased beyond 10^6 Wm-2. We also find that below 1 MJ, there is a cutoff point at high incident flux beyond which we find no more inflated planets, and that this cutoff point decreases as the mass decreases. At incident fluxes higher than 1.6*10^6 Wm -2 and in a mass range 0.37 - 0.98 MJ, we find no evidence for a population of non-inflated hot Jupiters. Our study demonstrates the importance of population-level analysis to grasp the underlying properties of exoplanets., Comment: Accepted in A&A on 27/03/2018, 14 pages, 12 figures

Massive planets form within the lifetime of protoplanetary disks, and therefore, they are subject to orbital migration due to planet–disk interactions. When the first planet reaches the inner edge of the disk, its migration stops and consequently the second planet ends up locked in resonance with the first one. We detail how the resonant trapping works comparing semi-analytical formulae and numerical simulations. We restrict to the case of two equal-mass coplanar planets trapped in first-order resonances, but the method can be easily generalized. We first describe the family of resonant stable equilibrium points (zero-amplitude libration orbits) using series expansions up to different orders in eccentricity as well as a non-expanded Hamiltonian. Then we show that during convergent migration the planets evolve along these families of equilibrium points. Eccentricity damping from the disk leads to a final equilibrium configuration that we predict precisely analytically. The fact that observed multi-exoplanetary systems are rarely seen in resonances suggests that in most cases the resonant configurations achieved by migration become unstable after the removal of the protoplanetary disk. Here we probe the stability of the resonances as a function of planetary mass. For this purpose, we fictitiously increase the masses of resonant planets, adiabatically maintaining the low-amplitude libration regime until instability occurs. We discuss two hypotheses for the instability, that of a low-order secondary resonance of the libration frequency with a fast synodic frequency of the system, and that of minimal approach distance between planets. We show that secondary resonances do not seem to impact resonant systems at low amplitude of libration. Resonant systems are more stable than non-resonant ones for a given minimal distance at close encounters, but we show that the latter nevertheless play the decisive role in the destabilization of resonant pairs. We show evidence that as the planetary mass increases and the minimal distance between planets gets smaller in terms of mutual Hill radius, the region of stability around the resonance center shrinks, until the equilibrium point itself becomes unstable.

Space-based high contrast imaging mission concepts for studying rocky exoplanets in reflected light are currently under community study. We develop an inverse modeling framework to estimate the science return of such missions given different instrument design considerations. By combining an exoplanet albedo model, an instrument noise model, and an ensemble Markov chain Monte Carlo sampler, we explore retrievals of atmospheric and planetary properties for Earth twins as a function of signal-to-noise ratio (SNR) and resolution ($R$). Our forward model includes Rayleigh scattering, single-layer water clouds with patchy coverage, and pressure-dependent absorption due to water vapor, oxygen, and ozone. We simulate data at $R = 70$ and $R = 140$ from 0.4-1.0 $\mu$m with SNR $ = 5, 10, 15, 20$ at 550 nm (i.e., for HabEx/LUVOIR-type instruments). At these same SNR, we simulate data for WFIRST paired with a starshade, which includes two photometric points between 0.48-0.6 $\mu$m and $R = 50$ spectroscopy from 0.6-0.97 $\mu$m. Given our noise model for WFIRST-type detectors, we find that weak detections of water vapor, ozone, and oxygen can be achieved with observations with at least $R = 70$ / SNR$\ = 15$, or $R = 140$ / SNR$\ = 10$ for improved detections. Meaningful constraints are only achieved with $R = 140$ / SNR$\ = 20$ data. The WFIRST data offer limited diagnostic information, needing at least SNR = 20 to weakly detect gases. Most scenarios place limits on planetary radius, but cannot constrain surface gravity and, thus, planetary mass., Comment: Resubmitted to AAS Journals after incorporating reviewer feedback. 26 pages, 18 figure, 9 tables

We present simulations of multi-year RV follow up campaigns of the {\it TESS} small exoplanet yield on the Automated Planet Finder telescope, using four different schemes to sample the transiting planets' RV phase curves. For planets below roughly 10 M$_\oplus$ we see a systematic bias of measured masses that are higher than the true planet mass, regardless of the observing scheme used. This produces a statistically significant difference in the mass-radius relation we recover, where planet masses are predicted to be too high and too similar across the entire super-Earth to Neptune radius range. This bias is due in part to only reporting masses that are measured with high statistical significance. Incorporating all mass measurements, even those that are essentially only upper limits, significantly mitigates this bias. We also find statistically significant differences between the mean number of planets measured at the 1-, 3-, and 5$\sigma_{\rm K}$ level by the different prioritization schemes. Our results show that prioritization schemes that more evenly sample the RV phase curves produce a larger number of significant mass detections. The scheme that aims to most uniformly sample the phase curve performs best, followed closely by the scheme which randomly samples, and then an in quadrature sampling approach. The fourth scheme, out of quadrature, performs noticeably worse. These results have important implications for determining accurate planet compositions and for designing effective RV follow up campaigns in the era of large planet detection surveys such as $K2$, $TESS$, and $PLATO$., Comment: 17 pages, 13 figures

We present the analyses of two short-timescale $(t_{\rm E} \sim 5~{\rm days})$ microlensing events, KMT-2016-BLG-1820 and KMT-2016-BLG-2142. In both light curves, the brief anomalies were clearly captured and densely covered by the Korea Microlensing Telescope Network survey. From these analyses, we find that both events have small Einstein radii of $\theta_{\rm E} = 0.12~{\rm mas}$, suggesting that the binary-lens systems are composed of very low-mass components and/or are located much closer to the lensed stars than to Earth. From Bayesian analyses, we find that these binaries have total system masses of $0.043_{-0.018}^{+0.043}~M_{\odot}$ and $0.088_{-0.041}^{+0.120}~M_{\odot}$, implying that they are well within the very-low-mass regime. The estimated lens-component masses indicate that the binary lenses consist of a giant-planet/brown-dwarf pair (KMT-2016-BLG-1820), and a dark/faint object pair (KMT-2016-BLG-2140) that are located near the deuterium-burning and hydrogen-burning mass limits, respectively. Both lens systems are likely to be in the Galactic disk with estimated distances of about $6$ kpc and $7$ kpc. The projected lens-components separations are $1.1$ AU and $0.8$ AU, and the mass ratios are $0.11$ and $0.20$. These prove that the microlensing method is effective to identify these closely-separated very-low-mass binaries having low mass-ratios., Comment: 20 pages, 7 figures, 3 tables. Accepted for publication in AJ

One of the main scenarios of planet formation is the core accretion model where a massive core forms first and then accretes a gaseous envelope. This core forms by accreting solids, either planetesimals, or pebbles. A key constraint in this model is that the accretion of gas must proceed before the dissipation of the gas disc. Classical planetesimal accretion scenario predicts that the time needed to form a giant planets core is much longer than the time needed to dissipate the disc. This difficulty led to the development of another accretion scenario, in which cores grow by accretion of pebbles, which are much smaller and thus more easily accreted, leading to a more rapid formation. The aim of this paper is to compare our updated pebble-based planet formation model with observations, in particular the well studied metallicity effect. We adopt the Bitsch et al. 2015a disc model and the Bitsch et al. 2015b pebble model and use a population synthesis approach to compare the formed planets with observations. We find that keeping the same parameters as in Bitsch et al. 2015b leads to no planet growth due to a computation mistake in the pebble flux (Bitsch et al. 2017). Indeed a large fraction of the heavy elements should be put into pebbles (Zpeb/Ztot = 0.9) in order to form massive planets using this approach. The resulting mass functions show a huge amount of giants and a lack of Neptune mass planets, which are abundant according to observations. To overcome this issue we include the computation of the internal structure for the planetary atmosphere to our model. This leads to the formation of Neptune mass planets but no observable giants. Reducing the opacity of the planetary envelope finally matches observations better. We conclude that modeling the internal structure for the planetary atmosphere is necessary to reproduce observations., Comment: 13 pages, 22 figures

Kepler-419 is a planetary system discovered by the Kepler photometry which is known to harbour two massive giant planets: [...] Here we present new radial velocity (RV) measurements secured over more than two years with the SOPHIE spectrograph, where both planets are clearly detected. The RV data is modelled together with the Kepler photometry using a photodynamical model. The inclusion of velocity information breaks the $MR^{-3}$ degeneracy inherent in timing data alone, allowing us to measure the absolute stellar and planetary radii and masses. With uncertainties of 12% and 13% for the stellar and inner planet radii, and 35%, 24%, and 35% for the masses of the star, planet b, and planet c respectively, these measurements are the most precise to date for a single host star system using this technique. The transiting planet mass is determined at better precision than the star mass. This shows that modelling the radial velocities and the light curve together in systems of dynamically interacting planets provides a way of characterising both the star and the planets without being limited by knowledge of the star. On the other hand, the period ratio and eccentricities place the Kepler-419 system in a sweet spot; had around twice as many transits been observed, the mass of the transiting planet could have been measured using its own TTVs. Finally, the origin of the Kepler-419 system is discussed. We show that the system is near a coplanar high-eccentricity secular fixed point, related to the alignment of the orbits, which has prevented the inner orbit from circularising. For most other relative apsidal orientations, planet b's orbit would be circular with a semi-major axis of 0.03 au. This suggests a mechanism for forming hot Jupiters in multiplanetary systems without the need of high mutual inclinations., Comment: 16 pages, 17 figures, accepted for publication in A&A

Atmospheric characterization through spectroscopic analysis, an essential tool of modern exoplanet science, can benefit significantly from the context provided by interior structure models. In particular, the planet's bulk metallicity $Z_p$ places an upper limit on potential atmospheric metallicity. Here we construct interior structure models to derive $Z_p$ and atmospheric metallicity upper limits for 403 known transiting giant exoplanets. These limits are low enough that they can usefully inform atmosphere models. Additionally, we argue that comparing $Z_p$ to the observed atmospheric metallicity gives a useful measure of how well-mixed metals are within the planet. This represents a new avenue for learning about planetary interiors. To aid in the future characterization of new planet discoveries we derive analytic prior predictions of atmosphere metallicity as a function of planet mass, and evaluate the effectiveness of our approach on Jupiter and Saturn. We we include log-linear fits for approximating the metallicities of planets not in our catalog., Comment: Submitted to ApJ. Revision includes more planets and clarity improvements, scientific results unchanged

We run numerical simulations to study the accretion of gas and dust grains onto gas giant planets embedded into massive protoplanetary discs. The outcome is found to depend on the disc cooling rate, planet mass, grain size and irradiative feedback from the planet. If radiative cooling is efficient, planets accrete both gas and pebbles rapidly, open a gap and usually become massive brown dwarfs. In the inefficient cooling case, gas is too hot to accrete onto the planet but pebble accretion continues and the planets migrate inward rapidly. Radiative feedback from the planet tends to suppress gas accretion. Our simulations predict that metal enrichment of planets by dust grain accretion inversely correlates with the final planet mass, in accordance with the observed trend in the inferred bulk composition of Solar System and exosolar giant planets. To account for observations, however, as much as ~30-50% of the dust mass should be in the form of large grains., Comment: Accepted to MNRAS

We aim to see if the difference between equilibrium and disequilibrium chemistry is observable in the atmospheres of transiting planets by the James Webb Space Telescope (JWST). We perform a case study comparing the dayside emission spectra of three planets like HD 189733b, WASP-80b, and GJ436b, in and out of chemical equilibrium at two metallicities each. These three planets were chosen because they span a large range of planetary masses and equilibrium temperatures, from hot and Jupiter-sized to warm and Neptune-sized. We link the one-dimensional disequilibrium chemistry model from Venot et al. (2012) in which thermochemical kinetics, vertical transport, and photochemistry are taken into account, to the one-dimensional, pseudo line-by-line radiative transfer model, Pyrat Bay, developed especially for hot Jupiters, and then simulate JWST spectra using PandExo for comparing the effects of temperature, metallicity, and radius. We find the most significant differences from 4 to 5 $\mu$m due to disequilibrium from CO and CO$_{2}$ abundances, and also H$_{2}$O for select cases. Our case study shows a certain "sweet spot" of planetary mass, temperature, and metallicity where the difference between equilibrium and disequilibrium is observable. For a planet similar to WASP-80b, JWST's NIRSpec G395M can detect differences due to disequilibrium chemistry with one eclipse event. For a planet similar to GJ 436b, the observability of differences due to disequilibrium chemistry is possible at low metallicity given five eclipse events, but not possible at the higher metallicity., Comment: Accepted for publication to ApJ on November 27, 2017

Searches for stellar companions to hot Jupiters (HJs) have revealed that planetary systems hosting a HJ are approximately three times more likely to have a stellar companion with a semimajor axis between 50 and 2000 AU, compared to field stars. This correlation suggests that HJ formation is affected by the stellar binary companion. A potential model is high-eccentricity migration, in which the binary companion induces high-eccentricity Lidov-Kozai (LK) oscillations in the proto-HJ orbit, triggering orbital migration driven by tides. A pitfall of this `binary-LK' model is that the observed stellar binaries hosting HJs are typically too wide to produce HJs in sufficient numbers, because of suppression by short-range forces. We propose a modification to the binary-LK model in which there is a second giant planet orbiting the proto-HJ at a semimajor axis of several tens of AU. Such companions are currently hidden to observations, but their presence could be manifested by a propagation of the perturbation of the stellar binary companion inwards to the proto-HJ, thereby overcoming the barrier imposed by short-range forces. Our model does not require the planetary companion orbit to be eccentric and/or inclined with respect to the proto-HJ, but its semimajor axis should lie in a specific range given the planetary mass and binary semimajor axis, and the inclination with respect to the binary should be near $40^\circ$ or $140^\circ$. Our prediction for planetary companions to HJs in stellar binaries should be testable by future observations., Comment: Accepted for publication in ApJL. 6 pages, 4 figures

The anomalously large radii of strongly irradiated exoplanets have remained a major puzzle in astronomy. Based on a 2D steady state atmospheric circulation model, the validity of which is assessed by comparison to 3D calculations, we reveal a new mechanism, namely the advection of the potential temperature due to mass and longitudinal momentum conservation, a process occuring in the Earth's atmosphere or oceans. At depth, the vanishing heating flux forces the atmospheric structure to converge to a hotter adiabat than the one obtained with 1D calculations, implying a larger radius for the planet. Not only do the calculations reproduce the observed radius of HD209458b, but also the observed correlation between radius inflation and irradiation for transiting planets. Vertical advection of potential temperature induced by non uniform atmospheric heating thus provides a robust mechanism explaining the inflated radii of irradiated hot Jupiters., Comment: accepted in ApJ

In the presence of a magnetic field and weakly ionizing winds, ohmic dissipation is expected to take place in the envelopes of Jovian and lower-mass planets alike. While the process has been investigated on the former, there have been no studies done on mini-Neptunes so far. From structure and thermal evolution models, we determine that the required energy deposition for halting the contraction of mini-Neptunes increases with planetary mass and envelope fraction. Scaled to the insolation power, the ohmic heating needed is small $\sim10^{-5}$ -- orders of magnitude lower than for exo-Jupiters $\sim 10^{-2}$. Conversely, from solving the magnetic induction equation, we find that ohmic energy is dissipated more readily for lower-mass planets and those with larger envelope fractions. Combining these two trends, we find that ohmic dissipation in hot mini-Neptunes is strong enough to inflate their radii ($\sim 10^{15}$ W for $T_{eq}=1400K$). The implication is that the radii of hot mini-Neptunes may be attributed in part to ohmic heating. Thus, there is a trade-off between ohmic dissipation and H/He content for hot mini-Neptunes, adding a new degeneracy for the interpretation of the composition of such planets. In addition, ohmic dissipation would make mini-Neptunes more vulnerable to atmospheric evaporation., Comment: Published in ApJ

We examine whether various characteristics of planet-driven spiral arms can be used to constrain the masses of unseen planets and their positions within their disks. By carrying out two-dimensional hydrodynamic simulations varying planet mass and disk gas temperature, we find that a larger number of spiral arms form with a smaller planet mass and a lower disk temperature. A planet excites two or more spiral arms interior to its orbit for a range of disk temperature characterized by the disk aspect ratio $0.04\leq(h/r)_p\leq0.15$, whereas exterior to a planet's orbit multiple spiral arms can form only in cold disks with $(h/r)_p \lesssim 0.06$. Constraining the planet mass with the pitch angle of spiral arms requires accurate disk temperature measurements that might be challenging even with ALMA. However, the property that the pitch angle of planet-driven spiral arms decreases away from the planet can be a powerful diagnostic to determine whether the planet is located interior or exterior to the observed spirals. The arm-to-arm separations increase as a function of planet mass, consistent with previous studies; however, the exact slope depends on disk temperature as well as the radial location where the arm-to-arm separations are measured. We apply these diagnostics to the spiral arms seen in MWC 758 and Elias 2-27. As shown in Bae et al. (2017), planet-driven spiral arms can create concentric rings and gaps, which can produce more dominant observable signature than spiral arms under certain circumstances. We discuss the observability of planet-driven spiral arms versus rings and gaps., Comment: 15 pages, 10 figures, Figure 2 size reduced to meet the requirement, accepted for publication in the ApJ

The distribution of exoplanet masses is not primordial. After the initial stage of planet formation is complete, the gravitational interactions between planets can lead to the physical collision of two planets, or the ejection of one or more planets from the system. When this occurs, the remaining planets are typically left in more eccentric orbits. Here we use present-day eccentricities of the observed exoplanet population to reconstruct the initial mass function of exoplanets before the onset of dynamical instability. We developed a Bayesian framework that combines data from N-body simulations with present-day observations to compute a probability distribution for the planets that were ejected or collided in the past. Integrating across the exoplanet population, we obtained an estimate of the initial mass function of exoplanets. We find that the ejected planets are primarily sub-Saturn type planets. While the present-day distribution appears to be bimodal, with peaks around $\sim 1 M_{\rm J}$ and $\sim 20 M_\oplus$, this bimodality does not seem to be primordial. Instead, planets around $\sim 60 M_\oplus$ appear to be preferentially removed by dynamical instabilities. Attempts to reproduce exoplanet populations using population synthesis codes should be mindful of the fact that the present population has been depleted of intermediate-mass planets. Future work should explore how the system architecture and multiplicity might alter our results., Comment: 10 pages, 9 figures; submitted to MNRAS

We have investigated the probability of temporary capture of asteroids in eccentric orbits by a planet in a circular or an eccentric orbit through analytical and numerical calculations. We found that in the limit of the circular orbit, the capture probability is $\sim 0.1\%$ of encounters to the planet's Hill sphere, independent of planetary mass and semimajor axis. In general, the temporary capture becomes more difficult as the planet's eccentricity ($e_{\rm p}$) increases. We found that the capture probability is almost independent of $e_{\rm p}$ until a critical value ($e_{\rm p}^{\rm c}$) that is given by $\simeq$ 5 times the Hill radius scaled by the planet's semimajor axis. For $e_{\rm p} > e_{\rm p}^{\rm c}$, the probability decreases approximately in proportion to $e_{\rm p}^{-1}$. The current orbital eccentricity of Mars is several times larger than $e_{\rm p}^{\rm c}$. However, since the range of secular change in Martian eccentricity overlaps $e_{\rm p}^{\rm c}$, the capture of minor bodies by the past Mars is not ruled out., Accepted for publication in AJ, 14 pages and 6 figures

We present the results of our spectroscopic observations of eight detached eclipsing binaries (DEBs), selected from the $Kepler$ Eclipsing Binary Catalog. Radial velocities (RVs) were calculated from high resolution spectra obtained with the HIDES pectrograph, attached to the 1.88-m telescope of the Okayama Astrophysical Observatory, and were used to characterize the targets in combination with the $Kepler$ light curves. For each binary we obtained a full set of orbital and physical parameters, reaching precision below 3% in masses and radii for 5 pairs. By comparing our results with theoretical models, we assess the distance, age and evolutionary status of the researched objects. We also study eclipse timing variations of selected objects, and identify a new system with a $\gamma$Dor pulsator. Two systems are triples, and show lines coming from three components. In one case the motion of the outer star and the perturbation in the RVs of the inner binary are clearly visible and periodical, which allows us to directly calculate the mass of the third star, and inclination of the outer orbit. In the second case we only see a clear motion of the tertiary, and investigate two scenarios: that it is a linear trend coming from the orbital motion around the inner binary, and that it is caused by a planetary mass companion. When possible, we also compare our results with the literature, and conclude that only by combining photometry with RVs it is possible to obtain correct physical parameters of both components of a DEB., Comment: Accepted by MNRAS, uncorrected proof. 23 pages, 19 figures, 9 tables (2 incomplete - will be published entirely with the electronic version in MNRAS)

With the discovery of rocky planets in the temperate habitable zone (HZ) of the close-by cool star TRAPPIST-1 the question of whether such planets could harbour life arises. Habitable planets around red dwarf stars can orbit in radiation environments that can be life-sterilizing. UV flares from these stars are more frequent and intense than solar flares. Additionally, their temperate HZs are closer to the star. Here we present UV surface environment models for TRAPPIST-1's HZ planets and explore the implications for life. TRAPPIST-1 has high X-ray/EUV activity, placing planetary atmospheres at risk from erosion. If a dense Earth-like atmosphere with a protective ozone layer exists on planets in the HZ of TRAPPIST-1, UV surface environments would be similar to present-day Earth. However an eroded or an anoxic atmosphere, would allow more UV to reach the surface, making surface environments hostile even to highly UV-tolerant terrestrial extremophiles. If future observations detect ozone in the atmospheres of any of the planets in the HZ of TRAPPIST-1, these would be interesting targets for the search for surface life. We anticipate our assay to be a starting point for in-depth exploration of stellar and atmospheric observations of the TRAPPIST-1 planets to constrain their UV-surface-habitability., Comment: 5 pages, 4 figures. Revised version contains extra information on the model atmospheres and minor changes to figures. Accepted for publication in MNRAS Letters