Giant Planet Formation Research Articles

CANDIDATE WATER VAPOR LINES TO LOCATE THE H2O SNOWLINE THROUGH HIGH-DISPERSION SPECTROSCOPIC OBSERVATIONS. I. THE CASE OF A T TAURI STAR

ABSTRACT Inside the H2O snowline of protoplanetary disks, water evaporates from the dust-grain surface into the gas phase, whereas it is frozen out onto the dust in the cold region beyond the snowline. H2O ice enhances the solid material in the cold outer part of a disk, which promotes the formation of gas-giant planet cores. We can regard the H2O snowline as the surface that divides the regions between rocky and gaseous giant planet formation. Thus observationally measuring the location of the H2O snowline is crucial for understanding the planetesimal and planet formation processes, and the origin of water on Earth. In this paper, we find candidate water lines to locate the H2O snowline through future high-dispersion spectroscopic observations. First, we calculate the chemical composition of the disk and investigate the abundance distributions of H2O gas and ice, and the position of the H2O snowline. We confirm that the abundance of H2O gas is high not only in the hot midplane region inside the H2O snowline but also in the hot surface layer of the outer disk. Second, we calculate the H2O line profiles and identify those H2O lines that are promising for locating the H2O snowline: the identified lines are those that have small Einstein A coefficients and high upper state energies. The wavelengths of the candidate H2O lines range from mid-infrared to sub-millimeter, and they overlap with the regions accessible to the Atacama Large Millimeter/sub-millimeter Array and future mid-infrared high-dispersion spectrographs (e.g., TMT/MICHI, SPICA).

FRIENDS OF HOT JUPITERS. IV. STELLAR COMPANIONS BEYOND 50 au MIGHT FACILITATE GIANT PLANET FORMATION, BUT MOST ARE UNLIKELY TO CAUSE KOZAI–LIDOV MIGRATION

ABSTRACT Stellar companions can influence the formation and evolution of planetary systems, but there are currently few observational constraints on the properties of planet-hosting binary star systems. We search for stellar companions around 77 transiting hot Jupiter systems to explore the statistical properties of this population of companions as compared to field stars of similar spectral type. After correcting for survey incompleteness, we find that of hot Jupiter systems have stellar companions with semimajor axes between 50 and 2000 au. This is 2.9 times larger than the field star companion fraction in this separation range, with a significance of . In the 1–50 au range, only of hot Jupiters host stellar companions, compared to the field star value of , which is a difference. We find that the distribution of mass ratios for stellar companions to hot Jupiter systems peaks at small values and therefore differs from that of field star binaries which tend to be uniformly distributed across all mass ratios. We conclude that either wide separation stellar binaries are more favorable sites for gas giant planet formation at all separations, or that the presence of stellar companions preferentially causes the inward migration of gas giant planets that formed farther out in the disk via dynamical processes such as Kozai–Lidov oscillations. We determine that less than 20% of hot Jupiters have stellar companions capable of inducing Kozai–Lidov oscillations assuming initial semimajor axes between 1 and 5 au, implying that the enhanced companion occurrence is likely correlated with environments where gas giants can form efficiently.

GAS GIANT PLANET FORMATION IN THE PHOTOEVAPORATING DISK. I. GAP FORMATION

ABSTRACT Planet formation and photoevaporation have both been considered as gap opening mechanisms in protoplanetary disks. We have studied giant planet formation in a photoevaporating disk with long-term evolution. Our calculations suggest that the core accretion rate of a protoplanet declines and the trigger of the runaway gas accretion for a giant planet is delayed under the action of photoevaporation. We find that the final mass of a giant planet characterized by the “gap-limiting” case is not influenced by photoevaporation but the final mass of a giant planet characterized by the “diffusion-limiting” case is greatly influenced by photoevaporation. Considering the formation process of giant planets, we suggest that the locations of the gaps opened by giant planets are within 30–40 au and the gap width in the “gap-limiting” case is wider than that in the “diffusion-limiting” case. We also find that gaps in photoevaporating disks are wider than those in non-photoevaporating disks. Our calculations suggest that the origins of multiple gaps in a disk can be diverse depending on their formation locations. In the formation region of giant planets, gaps are opened by giant planets. The outer gap beyond the giant planet formation region may be opened under the action of photoevaporation. A gap may also be opened at 1–3 au under the actions of photoevaporating dissipation and gas accretion of the outer giant planets.

Search for giant planets in M67

Since 2008 we used high-precision radial velocity (RV) measurements obtained with different telescopes to detect signatures of massive planets around main-sequence and evolved stars of the open cluster (OC) M67. We aimed to perform a long-term study on giant planet formation in open clusters and determine how this formation depends on stellar mass and chemical composition. A new hot Jupiter (HJ) around the main-sequence star YBP401 is reported in this work. An update of the RV measurements for the two HJ host-stars YBP1194 and YBP1514 is also discussed. Our sample of 66 main-sequence and turnoff stars includes 3 HJs, which indicates a high rate of HJs in this cluster (~5.6% for single stars and ~4.5% for the full sample ). This rate is much higher than what has been discovered in the field, either with RV surveys or by transits. High metallicity is not a cause for the excess of HJs in M67, nor can the excess be attributed to high stellar masses. When combining this rate with the non-zero eccentricity of the orbits, our results are qualitatively consistent with a HJ formation scenario dominated by strong encounters with other stars or binary companions and subsequent planet-planet scattering, as predicted by N-body simulations.

PROTOPLANETARY DISKS IN THE HOSTILE ENVIRONMENT OF CARINA

ABSTRACT We report the first direct imaging of protoplanetary disks in the star-forming region of Carina, the most distant massive cluster in which disks have been imaged. Using the Atacama Large Millimeter/sub-millimeter Array (ALMA), the disks are observed around two young stellar objects (YSOs) that are embedded inside evaporating gaseous globules and exhibit jet activity. The disks have an average radius of 60 au and total masses of 30 and 50 . Given the measured masses, the minimum timescale required for planet formation (∼1–2 Myr) and the average age of the Carina population (∼1–4 Myr), it is plausible that young planets are present or their formation is currently ongoing in these disks. The non-detection of millimeter emission above the 4σ threshold ( ) in the core of the massive cluster Trumpler 14, an area containing previously identified proplyd candidates, suggests evidence for rapid photo-evaporative disk destruction in the cluster’s harsh radiation field. This would prevent the formation of giant gas planets in disks located in the cores of Carina’s dense subclusters, whereas the majority of YSO disks in the wider Carina region remain unaffected by external photoevaporation.

DIRECT DETECTION OF PRECURSORS OF GAS GIANTS FORMED BY GRAVITATIONAL INSTABILITY WITH THE ATACAMA LARGE MILLIMETER/SUBMILLIMETER ARRAY

ABSTRACT Phases of gravitational instability are expected in the early phases of disk evolution, when the disk mass is still a substantial fraction of the mass of the star. Disk fragmentation into sub-stellar objects could occur in the cold exterior part of the disk. Direct detection of massive gaseous clumps on their way to collapse into gas giant planets would offer an unprecedented test of the disk instability model. Here we use state-of-the-art 3D radiation-hydro simulations of disks undergoing fragmentation into massive gas giants, post-processed with RADMC-3D to produce dust continuum emission maps. These are then fed into the Common Astronomy Software Applications (CASA) ALMA simulator. The synthetic maps show that both overdense spiral arms and actual clumps at different stages of collapse can be detected with the Atacama Large Millimeter/submillimeter Array (ALMA) in the full configuration at the distance of the Ophiuchus star forming region (125 pc). The detection of clumps is particularly effective at shorter wavelengths (690 GHz) combining two resolutions with multi-scale clean. Furthermore, we show that a flux-based estimate of the mass of a protoplanetary clump can be comparable to a factor of three higher than the gravitationally bound clump mass. The estimated mass depends on the assumed opacity, and on the gas temperature, which should be set using the input of radiation-hydro simulations. We conclude that ALMA has the capability to detect “smoking gun” systems that are a signpost of the disk instability model for gas giant planet formation.

Astrophysical fluid dynamics

These lecture notes and example problems are based on a course given at the University of Cambridge in Part III of the Mathematical Tripos. Fluid dynamics is involved in a very wide range of astrophysical phenomena, such as the formation and internal dynamics of stars and giant planets, the workings of jets and accretion discs around stars and black holes and the dynamics of the expanding Universe. Effects that can be important in astrophysical fluids include compressibility, self-gravitation and the dynamical influence of the magnetic field that is ‘frozen in’ to a highly conducting plasma. The basic models introduced and applied in this course are Newtonian gas dynamics and magnetohydrodynamics (MHD) for an ideal compressible fluid. The mathematical structure of the governing equations and the associated conservation laws are explored in some detail because of their importance for both analytical and numerical methods of solution, as well as for physical interpretation. Linear and nonlinear waves, including shocks and other discontinuities, are discussed. The spherical blast wave resulting from a supernova, and involving a strong shock, is a classic problem that can be solved analytically. Steady solutions with spherical or axial symmetry reveal the physics of winds and jets from stars and discs. The linearized equations determine the oscillation modes of astrophysical bodies, as well as their stability and their response to tidal forcing.

Giant planet formation in radially structured protoplanetary discs

Our recent N-body simulations of planetary system formation, incorporating models for the main physical processes thought to be important during the building of planets (i.e. gas disc evolution, migration, planetesimal/boulder accretion, gas accretion onto cores, etc.), have been successful in reproducing some of the broad features of the observed exoplanet population (e.g. compact systems of low mass planets, hot Jupiters), but fail completely to form any surviving cold Jupiters. The primary reason for this failure is rapid inward migration of growing protoplanets during the gas accretion phase, resulting in the delivery of these bodies onto orbits close to the star. Here, we present the results of simulations that examine the formation of gas giant planets in protoplanetary discs that are radially structured due to spatial and temporal variations in the effective viscous stresses, and show that such a model results in the formation of a population of cold gas giants. Furthermore, when combined with models for disc photoevaporation and a central magnetospheric cavity, the simulations reproduce the well-known hot-Jupiter/cold-Jupiter dichotomy in the observed period distribution of giant exoplanets, with a period valley between 10-100 days.

The HARPS search for southern extra-solar planets

Stellar metallicity -- as a probe of the metallicity of proto-planetary disks -- is an important ingredient for giant planet formation, likely through its effect on the timescales in which rocky/icy planet cores can form. Giant planets have been found to be more frequent around metal-rich stars, in agreement with predictions based on the core-accretion theory. In the metal-poor regime, however, the frequency of planets, especially low-mass planets, and how it depends on metallicity are still largely unknown. As part of a planet search programme focused on metal-poor stars, we study the targets from this survey that were observed with HARPS on more than 75 nights. The main goals are to assess the presence of low-mass planets and provide a first estimate of the frequency of Neptunes and super-Earths around metal-poor stars. We perform a systematic search for planetary companions, both by analysing the periodograms of the radial-velocities and by comparing, in a statistically-meaningful way, models with an increasing number of Keplerians. A first constraint on the frequency of planets in our metal-poor sample is calculated considering the previous detection (in our sample) of a Neptune-sized planet around HD175607 and one candidate planet (with an orbital period of 68.42d and minimum mass $M_p \sin i = 11.14 \pm 2.47 M_{\oplus}$) for HD87838, announced in the present study. This frequency is determined to be close to 13% and is compared with results for solar-metallicity stars.

PEERING INTO THE GIANT-PLANET-FORMING REGION OF THE TW HYDRAE DISK WITH THE GEMINI PLANET IMAGER

We present Gemini Planet Imager (GPI) adaptive optics near-infrared images of the giant planet-forming regions of the protoplanetary disk orbiting the nearby (D = 54 pc), pre-main sequence (classical T Tauri) star TW Hydrae. The GPI images, which were obtained in coronagraphic/polarimetric mode, exploit starlight scattered off small dust grains to elucidate the surface density structure of the TW Hya disk from 80 AU to within 10 AU of the star at 1.5 AU resolution. The GPI polarized intensity images unambiguously confirm the presence of a gap in the radial surface brightness distribution of the inner disk. The gap is centered near 23 AU, with a width of 5 AU and a depth of 50%. In the context of recent simulations of giant planet formation in gaseous, dusty disks orbiting pre-main sequence stars, these results indicate that at least one young planet with a mass 0.2 M_J could be present in the TW Hya disk at an orbital semi-major axis similar to that of Uranus. If this (proto)planet is actively accreting gas from the disk, it may be readily detectable by GPI or a similarly sensitive, high-resolution infrared imaging system.

Quantification of stochastic fragmentation of self-gravitating discs

Using 2D smoothed particle hydrodynamics, we investigate the distribution of wait times between strong shocks in a turbulent, self-gravitating accretion disc. We show the resulting distributions do not depend strongly on the cooling time or resolution of the disc and that they are consistent with the predictions of earlier work (Young & Clarke 2015; Cossins et al. 2009, 2010). We use the distribution of wait times between shocks to estimate the likelihood of stochastic fragmentation by gradual contraction of shear-resistant clumps on the cooling time scale. We conclude that the stochastic fragmentation mechanism (Paardekooper 2012) cannot change the radius at which fragmentation is possible by more than ~20%, restricting direct gravitational collapse as a mechanism for giant planet formation to the outer regions of protoplanetary discs.

MAPPING THE SHORES OF THE BROWN DWARF DESERT. IV. OPHIUCHUS

We conduct a multiplicity survey of members of the rho Ophiuchus cloud complex with high resolution imaging to characterize the multiple star population of this nearby star forming region and investigate the relation between stellar multiplicity and star and planet formation. Our aperture masking survey reveals the presence of 5 new stellar companions beyond the reach of previous studies, but does not result in the detection of any new substellar companions. We find that 43+/-6% of the 114 stars in our survey have stellar mass companions between 1.3-780AU, while 7 (+8 -5)% host brown dwarf companions in the same interval. By combining this information with knowledge of disk-hosting stars, we show that the presence of a close binary companion (separation < 40 AU) significantly influences the lifetime of protoplanetary disks, a phenomenon previously seen in older star forming regions. At the ~1-2Myr age of our Ophiuchus members ~2/3 of close binary systems have lost their disks, compared to only ~30% of single stars and wide binaries. This has significant impact on the formation of giant planets, which are expected to require much longer than 1 Myr to form via core accretion and thus planets formed via this pathway should be rare in close binary systems.

On the history and future of cosmic planet formation

We combine constraints on galaxy formation histories with planet formation models, yielding the Earth-like and giant planet formation histories of the Milky Way and the Universe as a whole. In the Hubble Volume (10^13 Mpc^3), we expect there to be ~10^20 Earth-like and ~10^20 giant planets; our own galaxy is expected to host ~10^9 and ~10^10 Earth-like and giant planets, respectively. Proposed metallicity thresholds for planet formation do not significantly affect these numbers. However, the metallicity dependence for giant planets results in later typical formation times and larger host galaxies than for Earth-like planets. The Solar System formed at the median age for existing giant planets in the Milky Way, and consistent with past estimates, formed after 80% of Earth-like planets. However, if existing gas within virialised dark matter haloes continues to collapse and form stars and planets, the Universe will form over 10 times more planets than currently exist. We show that this would imply at least a 92% chance that we are not the only civilisation the Universe will ever have, independent of arguments involving the Drake Equation.

Forming terrestrial planets and delivering water

AbstractBuilding models capable of successfully matching the Terrestrial Planet's basic orbital and physical properties has proven difficult. Meanwhile, improved estimates of the nature of water-rich material accreted by the Earth, along with the timing of its delivery, have added even more constraints for models to match. While the outer Asteroid Belt seemingly provides a source for water-rich planetesimals, models that delivered enough of them to the still-forming Terrestrial Planets typically failed on other basic constraints - such as the mass of Mars.Recent models of Terrestrial Planet Formation have explored how the gas-driven migration of the Giant Planets can solve long-standing issues with the Earth/Mars size ratio. This model is forced to reproduce the orbital and taxonomic distribution of bodies in the Asteroid Belt from a much wider range of semimajor axis than previously considered. In doing so, it also provides a mechanism to feed planetesimals from between and beyond the Giant Planet formation region to the still-forming Terrestrial Planets.

The Hera Saturn entry probe mission

The Hera Saturn entry probe mission is proposed as an M-class mission led by ESA with a contribution from NASA. It consists of one atmospheric probe to be sent into the atmosphere of Saturn, and a Carrier-Relay spacecraft. In this concept, the Hera probe is composed of ESA and NASA elements, and the Carrier-Relay Spacecraft is delivered by ESA. The probe is powered by batteries, and the Carrier-Relay Spacecraft is powered by solar panels and batteries. We anticipate two major subsystems to be supplied by the United States, either by direct procurement by ESA or by contribution from NASA: the solar electric power system (including solar arrays and the power management and distribution system), and the probe entry system (including the thermal protection shield and aeroshell). Hera is designed to perform in situ measurements of the chemical and isotopic compositions as well as the dynamics of Saturn's atmosphere using a single probe, with the goal of improving our understanding of the origin, formation, and evolution of Saturn, the giant planets and their satellite systems, with extrapolation to extrasolar planets. Hera's aim is to probe well into the cloud-forming region of the troposphere, below the region accessible to remote sensing, to the locations where certain cosmogenically abundant species are expected to be well mixed. By leading to an improved understanding of the processes by which giant planets formed, including the composition and properties of the local solar nebula at the time and location of giant planet formation, Hera will extend the legacy of the Galileo and Cassini missions by further addressing the creation, formation, and chemical, dynamical, and thermal evolution of the giant planets, the entire solar system including Earth and the other terrestrial planets, and formation of other planetary systems.

IN SEARCH OF ALIEN JUPITERS.

The article discusses efforts by astronomers Bruce Macintosh and Jean-Luc Beuzit to capture images of planets in other solar systems which resemble the planet Jupiter. Details on the Gemini Planet Imager (GPI) used by Macintosh and the Spectro-Polarimetric High-contrast Expolanet Research Instrument (SPHERE) used by Beuzit are presented. The formation of giant planets and their role in shaping their planetary systems are also discussed.

THE DISSIMILAR CHEMICAL COMPOSITION OF THE PLANET-HOSTING STARS OF THE XO-2 BINARY SYSTEM

Using high-quality spectra of the twin stars in the XO-2 binary system, we have detected significant differences in the chemical composition of their photospheres. The differences correlate strongly with the elements' dust condensation temperature. In XO-2N, volatiles are enhanced by about 0.015 dex and refractories are overabundant by up to 0.090 dex. On average, our error bar in relative abundance is 0.012 dex. We present an early metal-depletion scenario in which the formation of the gas giant planets known to exist around these stars is responsible for a 0.015 dex offset in the abundances of all elements while 20 M_Earth of non-detected rocky objects that formed around XO-2S explain the additional refractory-element difference. An alternative explanation involves the late accretion of at least 20 M_Earth of planet-like material by XO-2N, allegedly as a result of the migration of the hot Jupiter detected around that star. Dust cleansing by a nearby hot star as well as age or Galactic birthplace effects can be ruled out as valid explanations for this phenomenon.

CONFIRMATION AND CHARACTERIZATION OF THE PROTOPLANET HD 100546 b—DIRECT EVIDENCE FOR GAS GIANT PLANET FORMATION AT 50 AU

We present the first multi-wavelength, high-contrast imaging study confirming the protoplanet embedded in the disk around the Herbig Ae/Be star HD 100546. The object is detected at L? () and M? (), but not at Ks (), and the emission consists of a point source component surrounded by spatially resolved emission. For the point source component we derive apparent magnitudes of mag, mag, and mag (3? limit), and a separation and position angle of and , and and in L? and M?, respectively. We demonstrate that the object is co-moving with HD 100546 and can reject any (sub-)stellar fore-/background object. Fitting a single-temperature blackbody to the observed fluxes of the point source component yields an effective temperature of K and a radius for the emitting area of . The best-fit luminosity is . We quantitatively compare our findings with predictions from evolutionary and atmospheric models for young, gas giant planets, discuss the possible existence of a warm, circumplanetary disk, and note that the deprojected physical separation from the host star of AU poses a challenge to standard planet formation theories. Considering the suspected existence of an additional planet orbiting at ?13?14 AU, HD 100546 appears to be an unprecedented laboratory to study the formation of multiple gas giant planets empirically.

INDIRECT DETECTION OF FORMING PROTOPLANETS VIA CHEMICAL ASYMMETRIES IN DISKS

We examine changes in the molecular abundances resulting from increased heating due to a self-luminous planetary companion embedded within a narrow circumstellar disk gap. Using 3D models that include stellar and planetary irradiation, we find that luminous young planets locally heat up the parent circumstellar disk by many tens of Kelvin, resulting in efficient thermal desorption of molecular species that are otherwise locally frozen out. Furthermore, the heating is deposited over large regions of the disk, ±5 AU radially and spanning azimuthally. From the 3D chemical models, we compute rotational line emission models and full Atacama Large Millimeter Array simulations, and find that the chemical signatures of the young planet are detectable as chemical asymmetries in observations. HCN and its isotopologues are particularly clear tracers of planetary heating for the models considered here, and emission from multiple transitions of the same species is detectable, which encodes temperature information in addition to possible velocity information from the spectra itself. We find submillimeter molecular emission will be a useful tool to study gas giant planet formation in situ, especially beyond AU.

Modeling gravitational instabilities in self-gravitating protoplanetary disks with adaptive mesh refinement techniques

The astonishing diversity in the observed planetary population requires theoretical efforts and advances in planet formation theories. Numerical approaches provide a method to tackle the weaknesses of current planet formation models and are an important tool to close gaps in poorly constrained areas. We present a global disk setup to model the first stages of giant planet formation via gravitational instabilities (GI) in 3D with the block-structured adaptive mesh refinement (AMR) hydrodynamics code ENZO. With this setup, we explore the impact of AMR techniques on the fragmentation and clumping due to large-scale instabilities using different AMR configurations. Additionally, we seek to derive general resolution criteria for global simulations of self-gravitating disks of variable extent. We run a grid of simulations with varying AMR settings, including runs with a static grid for comparison, and study the effects of varying the disk radius. Adopting a marginally stable disk profile (Q_init=1), we validate the numerical robustness of our model for different spatial extensions, from compact to larger, extended disks (R_disk = 10, 100 and 300 AU, M_disk ~ 0.05 M_Sun, M_star = 0.646 M_Sun). By combining our findings from the resolution and parameter studies we find a lower limit of the resolution to be able to resolve GI induced fragmentation features and distinct, turbulence inducing clumps. Irrespective of the physical extension of the disk, topologically disconnected clump features are only resolved if the fragmentation-active zone of the disk is resolved with at least 100 cells, which holds as a minimum requirement for all global disk setups. Our simulations illustrate the capabilities of AMR-based modeling techniques for planet formation simulations and underline the importance of balanced refinement settings to reproduce fragmenting structures.