sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Rice, Ken

Normalized to: Rice, K.

72 article(s) in total. 399 co-authors, from 1 to 29 common article(s). Median position in authors list is 5,0.

[1]  [pdf] - 2040756
Massive discs around low-mass stars
Comments: 21 pages, submitted to MNRAS
Submitted: 2020-01-17, last modified: 2020-01-31
We use a suite of SPH simulations to investigate the susceptibility of protoplanetary discs to the effects of self-gravity as a function of star-disc properties. We also include passive irradiation from the host star using different models for the stellar luminosities. The critical disc-to-star mass ratio for axisymmetry (for which we produce criteria) increases significantly for low-mass stars. This could have important consequences for increasing the potential mass reservoir in a proto Trappist-1 system, since even the efficient Ormel et al. (2017) formation model will be influenced by processes like external photoevaporation, which can rapidly and dramatically deplete the dust reservoir. The aforementioned scaling of the critical $M_d/M_*$ for axisymmetry occurs in part because the Toomre $Q$ parameter has a linear dependence on surface density (which promotes instability) and only an $M_*^{1/2}$ dependence on shear (which reduces instability), but also occurs because, for a given $M_d/M_*$, the thermal evolution depends on the host star mass. The early phase stellar irradiation of the disc (for which the luminosity is much higher than at the zero age main sequence, particularly at low stellar masses) can also play a key role in significantly reducing the role of self-gravity, meaning that even Solar mass stars could support axisymmetric discs a factor two higher in mass than usually considered possible. We apply our criteria to the DSHARP discs with spirals, finding that self-gravity can explain the observed spirals so long as the discs are optically thick to the host star irradiation.
[2]  [pdf] - 2038459
Fragmentation favoured in discs around higher mass stars
Comments: Accepted for publication in MNRAS
Submitted: 2020-01-17
We investigate how a protoplanetary disc's susceptibility to gravitational instabilities and fragmentation depends on the mass of its host star. We use 1D disc models in conjunction with 3D SPH simulations to determine the critical disc-to-star mass ratios at which discs become unstable against fragmentation, finding that discs become increasingly prone to the effects of self-gravity as we increase the host star mass. The actual limit for stability is sensitive to the disc temperature, so if the disc is optically thin stellar irradiation can dramatically stabilise discs against gravitational instability. However, even when this is the case we find that discs around $2$M$_{\odot}$ stars are prone to fragmentation, which will act to produce wide-orbit giant planets and brown dwarfs. The consequences of this work are two-fold: that low mass stars could in principle support high disc-to-star mass ratios, and that higher mass stars have discs that are more prone to fragmentation, which is qualitatively consistent with observations that favour high-mass wide-orbit planets around higher mass stars. We also find that the initial masses of these planets depends on the temperature in the disc at large radii, which itself depends on the level of stellar irradiation.
[3]  [pdf] - 2042365
An ultra-short period rocky super-Earth orbiting the G2-star HD 80653
Comments: 15 pages, 12 figures; accepted by A&A
Submitted: 2020-01-07
Ultra-short period (USP) planets are a class of exoplanets with periods shorter than one day. The origin of this sub-population of planets is still unclear, with different formation scenarios highly dependent on the composition of the USP planets. A better understanding of this class of exoplanets will, therefore, require an increase in the sample of such planets that have accurate and precise masses and radii, which also includes estimates of the level of irradiation and information about possible companions. Here we report a detailed characterization of a USP planet around the solar-type star HD 80653 $\equiv$ EP 251279430 using the K2 light curve and 108 precise radial velocities obtained with the HARPS-N spectrograph, installed on the Telescopio Nazionale Galileo. From the K2 C16 data, we found one super-Earth planet ($R_{b}=1.613\pm0.071 R_{\oplus}$) transiting the star on a short-period orbit ($P_{\rm b}=0.719573\pm0.000021$ d). From our radial velocity measurements, we constrained the mass of HD 80653 b to $M_{b}=5.60\pm0.43 M_{\oplus}$. We also detected a clear long-term trend in the radial velocity data. We derived the fundamental stellar parameters and determined a radius of $R_{\star}=1.22\pm0.01 R_{\odot}$ and mass of $M_{\star}=1.18\pm0.04 M_{\odot}$, suggesting that HD 80653, has an age of $2.7\pm1.2$ Gyr. The bulk density ($\rho_{b} = 7.4 \pm 1.1$ g cm$^{-3}$) of the planet is consistent with an Earth-like composition of rock and iron with no thick atmosphere. Our analysis of the K2 photometry also suggests hints of a shallow secondary eclipse with a depth of 8.1$\pm$3.7 ppm. Flux variations along the orbital phase are consistent with zero. The most important contribution might come from the day-side thermal emission from the surface of the planet at $T\sim3480$ K.
[4]  [pdf] - 2013826
An 11 Earth-Mass, Long-Period Sub-Neptune Orbiting a Sun-like Star
Comments: Published in AJ, 17 pages, 8 figures, 3 tables
Submitted: 2019-08-22, last modified: 2019-12-11
Although several thousands of exoplanets have now been detected and characterized, observational biases have led to a paucity of long-period, low-mass exoplanets with measured masses and a corresponding lag in our understanding of such planets. In this paper we report the mass estimation and characterization of the long-period exoplanet Kepler-538b. This planet orbits a Sun-like star (V = 11.27) with M_* = 0.892 +/- (0.051, 0.035) M_sun and R_* = 0.8717 +/- (0.0064, 0.0061) R_sun. Kepler-538b is a 2.215 +/- (0.040, 0.034) R_earth sub-Neptune with a period of P = 81.73778 +/- 0.00013 d. It is the only known planet in the system. We collected radial velocity (RV) observations with HIRES on Keck I and HARPS-N on the TNG. We characterized stellar activity by a Gaussian process with a quasi-periodic kernel applied to our RV and cross correlation function full width at half maximum (FWHM) observations. By simultaneously modeling Kepler photometry, RV, and FWHM observations, we found a semi-amplitude of K = 1.68 +/- (0.39, 0.38) m s^-1 and a planet mass of M_p = 10.6 +/- (2.5, 2.4) M_earth. Kepler-538b is the smallest planet beyond P = 50 d with an RV mass measurement. The planet likely consists of a significant fraction of ices (dominated by water ice), in addition to rocks/metals, and a small amount of gas. Sophisticated modeling techniques such as those used in this paper, combined with future spectrographs with ultra high-precision and stability will be vital for yielding more mass measurements in this poorly understood exoplanet regime. This in turn will improve our understanding of the relationship between planet composition and insolation flux and how the rocky to gaseous transition depends on planetary equilibrium temperature.
[5]  [pdf] - 1999746
An extremely low-density and temperate giant exoplanet
Comments: Preprint submitted to Nature Astronomy. The results have not been peer-reviewed yet. Supplementary informations available as ancillary file
Submitted: 2019-11-17
Transiting extrasolar planets are key objects in the study of the formation, migration, and evolution of planetary systems. In particular, the exploration of the atmospheres of giant planets, through transmission spectroscopy or direct imaging, has revealed a large diversity in their chemical composition and physical properties. Studying these giant planets allows one to test the global climate models that are used for the Earth and other solar system planets. However, these studies are mostly limited either to highly-irradiated transiting giant planets or directly-imaged giant planets at large separations. Here we report the physical characterisation of the planets in a bright multi-planetary system (HIP41378) in which the outer planet, HIP41378 f is a Saturn-sized planet (9.2 $\pm$ 0.1 R$_\oplus$) with an anomalously low density of 0.09 $\pm$ 0.02 g cm$^{-3}$ that is not yet understood. Its equilibrium temperature is about 300 K. Therefore, it represents a planet with a mild temperature, in between the hot Jupiters and the colder giant planets of the Solar System. It opens a new window for atmospheric characterisation of giant exoplanets with a moderate irradiation, with the next-generation space telescopes such as JWST and ARIEL as well as the extremely-large ground-based telescopes. HIP41378 f is thus an important laboratory to understand the effect of the irradiation on the physical properties and chemical composition of the atmosphere of planets.
[6]  [pdf] - 1984005
Using HARPS-N to characterise the long-period planets in the PH-2 and Kepler-103 systems
Comments: Accepted for publication in MNRAS
Submitted: 2019-10-08
We present confirmation of the planetary nature of PH-2b, as well as the first mass estimates for the two planets in the Kepler-103 system. PH-2b and Kepler-103c are both long-period and transiting, a sparsely-populated category of exoplanet. We use {\it Kepler} light-curve data to estimate a radius, and then use HARPS-N radial velocities to determine the semi-amplitude of the stellar reflex motion and, hence, the planet mass. For PH-2b we recover a 3.5-$\sigma$ mass estimate of $M_p = 109^{+30}_{-32}$ M$_\oplus$ and a radius of $R_p = 9.49\pm0.16$ R$_\oplus$. This means that PH-2b has a Saturn-like bulk density and is the only planet of this type with an orbital period $P > 200$ days that orbits a single star. We find that Kepler-103b has a mass of $M_{\text{p,b}} = 11.7^{+4.31}_{-4.72}$ M$_{\oplus}$ and Kepler-103c has a mass of $M_{\text{p,c}} = 58.5^{+11.2}_{-11.4}$ M$_{\oplus}$. These are 2.5$\sigma$ and 5$\sigma$ results, respectively. With radii of $R_{\text{p,b}} = 3.49^{+0.06}_{-0.05}$ R$_\oplus$, and $R_{\text{p,c}} = 5.45^{+0.18}_{-0.17}$ R$_\oplus$, these results suggest that Kepler-103b has a Neptune-like density, while Kepler-103c is one of the highest density planets with a period $P > 100$ days. By providing high-precision estimates for the masses of the long-period, intermediate-mass planets PH-2b and Kepler-103c, we increase the sample of long-period planets with known masses and radii, which will improve our understanding of the mass-radius relation across the full range of exoplanet masses and radii.
[7]  [pdf] - 1897015
The temporal requirements of directly observing self-gravitating spiral waves in protoplanetary discs with ALMA
Comments: 13 pages, 6 figures
Submitted: 2019-01-08, last modified: 2019-06-08
We investigate how the detectability of signatures of self-gravity in a protoplanetary disc depends on its temporal evolution. We run a one-dimensional model for secular timescales to follow the disc mass as a function of time. We then combine this with three-dimensional global hydrodynamics simulations that employ a hybrid radiative transfer method to approximate realistic heating and cooling. We simulate ALMA continuum observations of these systems, and find that structures induced by the gravitational instability (GI) are readily detectable when $q=M_\mathrm{disc}/M_*\gtrsim 0.25$ and $R_\mathrm{outer}\lesssim 100$ au. The high accretion rate generated by gravito-turbulence in such a massive disc drains its mass to below the detection threshold in $\sim10^4$ years, or approximately 1 % of the typical disc lifetime. Therefore, discs with spiral arms detected in ALMA dust observations, if generated by self-gravity, must either be still receiving infall to maintain a high $q$ value, or have just emerged from their natal envelope. Detection of substructure in systems with lower $q$ is possible, but would require a specialist integration with the most extended configuration over several days. This disfavours the possibility of GI-caused spiral structure in systems with $q<0.25$ being detected in relatively short integration times, such as those found in the DSHARP ALMA survey (Andrews et al. 2018; Huang et al. 2018). We find no temporal dependence of detectability on dynamical timescales
[8]  [pdf] - 1896913
How formation time-scales affect the period dependence of the transition between rocky super-Earths and gaseous sub-Neptunes and implications for $\eta_\oplus$
Comments: 8 pages, 4 figures, accepted for publication in the Monthly Notices of the Royal Astronomical Society
Submitted: 2016-10-28, last modified: 2019-06-08
One of the most significant advances by NASA's ${\mathit Kepler}$ Mission was the discovery of an abundant new population of highly irradiated planets with sizes between the Earth and Neptune. Subsequent analysis showed that at ~1.5 Earth radii there is a transition from a population of predominantly rocky super-Earths to non-rocky sub-Neptunes, which must have substantial volatile envelopes. Determining the origin of these highly irradiated rocky planets will be critical to our understanding of low-mass planet formation and the frequency of potentially habitable Earth-like planets. These short-period rocky super-Earths could simply be the stripped cores of sub-Neptunes, which have lost their envelopes due to atmospheric photo-evaporation or other processes, or they might instead be a separate population of inherently rocky planets, which never had significant envelopes. Using models of atmospheric photo-evaporation, we show that if most bare rocky planets are the evaporated cores of sub-Neptunes then the transition radius should decrease as surveys push to longer orbital periods, since on wider orbits only planets with smaller less massive cores can be stripped. On the other hand, if most rocky planets formed after their disks dissipate then these planets will have formed without initial gaseous envelopes. In this case, we use N-body simulations of planet formation to show that the transition radius should increase with orbital period, due to the increasing solid mass available in their disks. Moreover, we show that distinguishing between these two scenarios should be possible in coming years with radial velocity follow-up of planets found by TESS. Finally, we discuss the broader implications of this work for current efforts to measure $\eta_{\mathrm{\oplus}}$, which may yield significant overestimates if most rocky planets form as evaporated cores.
[9]  [pdf] - 1912913
Temporal evolution and correlations of optical activity indicators measured in Sun-as-a-star observations
Comments: Accepted by Astronomy & Astrophysics
Submitted: 2019-06-07
(Abridged) We perform a detailed study of the main optical activity indicators (Ca II H & K, Balmer lines, Na I D$_{\rm 1}$ D$_{\rm 2}$, and He I D$_{\rm 3}$) measured for the Sun using the data provided by the HARPS-N solar-telescope feed at the Telescopio Nazionale Galileo. The value of the solar rotation period is found in all the activity indicators, with the only exception being H$\delta$. The derived values vary from 26.29 days (H$\gamma$ line) to 31.23 days (He I). From an analysis of sliding periodograms we find that in most of the activity indicators the spectral power is split into several "bands" of periods around 26 and 30 days, that might be explained by the migration of active regions between the equator and a latitude of $\sim$ 30$^{\circ}$, spot evolution or a combination of both effects. In agreement with previous works a typical lifetime of active regions of $\sim$ ten rotation periods is inferred from the pooled variance diagrams. We find that H$\alpha$, H$\beta$, H$\gamma$, H$\epsilon$, and He I show a significant correlation with the S index. Significant correlations between the contrast, bisector span, and the heliocentric radial velocity with the activity indexes are also found. We show that the full width at half maximum, the bisector, and the disc-integrated magnetic field correlate with the radial velocity variations. The correlation of the S index and H$\alpha$ changes with time, increasing with larger sun spot numbers and solar irradiance. A similar tendency with the S index - radial velocity correlation is also present in the data. Our results are consistent with a scenario in which higher activity favours the correlation between the S index and the H$\alpha$ activity indicators and between the S index and radial velocity variations.
[10]  [pdf] - 1894339
Three years of Sun-as-a-star radial-velocity observations on the approach to solar minimum
Comments: 19 pages, 15 figures, accepted for publication in MNRAS
Submitted: 2019-04-27
The time-variable velocity fields of solar-type stars limit the precision of radial-velocity determinations of their planets' masses, obstructing detection of Earth twins. Since 2015 July we have been monitoring disc-integrated sunlight in daytime using a purpose-built solar telescope and fibre feed to the HARPS-N stellar radial-velocity spectrometer. We present and analyse the solar radial-velocity measurements and cross-correlation function (CCF) parameters obtained in the first 3 years of observation, interpreting them in the context of spatially-resolved solar observations. We describe a Bayesian mixture-model approach to automated data-quality monitoring. We provide dynamical and daily differential-extinction corrections to place the radial velocities in the heliocentric reference frame, and the CCF shape parameters in the sidereal frame. We achieve a photon-noise limited radial-velocity precision better than 0.43 m s$^{-1}$ per 5-minute observation. The day-to-day precision is limited by zero-point calibration uncertainty with an RMS scatter of about 0.4 m s$^{-1}$. We find significant signals from granulation and solar activity. Within a day, granulation noise dominates, with an amplitude of about 0.4 m s$^{-1}$ and an autocorrelation half-life of 15 minutes. On longer timescales, activity dominates. Sunspot groups broaden the CCF as they cross the solar disc. Facular regions temporarily reduce the intrinsic asymmetry of the CCF. The radial-velocity increase that accompanies an active-region passage has a typical amplitude of 5 m s$^{-1}$ and is correlated with the line asymmetry, but leads it by 3 days. Spectral line-shape variability thus shows promise as a proxy for recovering the true radial velocity.
[11]  [pdf] - 1857740
A high binary fraction for the most massive close-in giant planets and brown dwarf desert members
Comments: Accepted for publication in MNRAS. 30 pages, 20 figures. Updated to include proof corrections
Submitted: 2019-03-06, last modified: 2019-03-28
Stellar multiplicity is believed to influence planetary formation and evolution, although the precise nature and extent of this role remain ambiguous. We present a study aimed at testing the role of stellar multiplicity in the formation and/or evolution of the most massive, close-in planetary and substellar companions. Using direct imaging observations, as well as the Gaia DR2 catalogue, we searched for wide binary companions to 38 stars hosting massive giant planets or brown dwarfs (M > 7 MJup) on orbits shorter than ~1 AU. We report the discovery of a new component in the WASP-14 system, and present an independent confirmation of a comoving companion to WASP-18. From a robust Bayesian statistical analysis, we derived a binary fraction of 79.0+13.2-14.7% between 20-10,000 AU for our sample, twice as high as for field stars with a 3-{\sigma} significance. This binary frequency was found to be larger than for lower-mass planets on similar orbits, and we observed a marginally higher binary rate for inner companions with periods shorter than 10 days. These results demonstrate that stellar companions greatly influence the formation and/or evolution of these systems, suggesting that the role played by binary companions becomes more important for higher-mass planets, and that this trend may be enhanced for systems with tighter orbits. Our analysis also revealed a peak in binary separation at 250 AU, highlighting a shortfall of close binaries among our sample. This indicates that the mechanisms affecting planet and brown dwarf formation or evolution in binaries must operate from wide separations, although we found that the Kozai-Lidov mechanism is unlikely to be the dominant underlying process. We conclude that binarity plays a crucial role in the existence of very massive short-period giant planets and brown dwarf desert inhabitants, which are almost exclusively observed in multiple systems.
[12]  [pdf] - 1834236
K2-291 b: A rocky super-Earth in a 2.2 day orbit
Comments: Accepted to AJ, 15 pages, 8 figures
Submitted: 2019-01-14, last modified: 2019-02-14
K2-291 (EPIC 247418783) is a solar-type star with a radius of R_star = 0.899 $\pm$ 0.034 R_sun and mass of M_star=0.934 $\pm$ 0.038 M_sun. From K2 C13 data, we found one super-Earth planet (R_p = 1.589+0.095-0.072 R_Earth) transiting this star on a short period orbit (P = 2.225177 +6.6e-5 -6.8e-5 days). We followed this system up with adaptive-optic imaging and spectroscopy to derive stellar parameters, search for stellar companions, and determine a planet mass. From our 75 radial velocity measurements using HIRES on Keck I and HARPS-N on Telescopio Nazionale Galileo, we constrained the mass of EPIC 247418783b to M_p = 6.49 $\pm$ 1.16 M_Earth. We found it necessary to model correlated stellar activity radial velocity signals with a Gaussian process in order to more accurately model the effect of stellar noise on our data; the addition of the Gaussian process also improved the precision of this mass measurement. With a bulk density of 8.84+2.50-2.03 g cm-3, the planet is consistent with an Earth-like rock/iron composition and no substantial gaseous envelope. Such an envelope, if it existed in the past, was likely eroded away by photo-evaporation during the first billion years of the star's lifetime.
[13]  [pdf] - 1860009
HARPS-N Solar Radial-Velocity Variations Are Dominated By Large, Bright Magnetic Regions
Comments: 16 pages, 7 figures, accepted by ApJ
Submitted: 2019-02-11
State of the art radial-velocity (RV) exoplanet searches are currently limited by RV signals arising from stellar magnetic activity. We analyze solar observations acquired over a 3-year period during the decline of Carrington Cycle 24 to test models of RV variation of Sun-like stars. A purpose-built solar telescope at the High Accuracy Radial velocity Planet Searcher for the Northern hemisphere (HARPS-N) provides disk-integrated solar spectra, from which we extract RVs and $\log{R'_{\rm HK}}$. The Solar Dynamics Observatory (SDO) provides disk-resolved images of magnetic activity. The Solar Radiation and Climate Experiment (SORCE) provides near-continuous solar photometry, analogous to a Kepler light curve. We verify that the SORCE photometry and HARPS-N $\log{R'_{\rm HK}}$ correlate strongly with the SDO-derived magnetic filling factor, while the HARPS-N RV variations do not. To explain this discrepancy, we test existing models of RV variations. We estimate the contributions of the suppression of convective blueshift and the rotational imbalance due to brightness inhomogeneities to the observed HARPS-N RVs. We investigate the time variation of these contributions over several rotation periods, and how these contributions depend on the area of active regions. We find that magnetic active regions smaller than $60 \ \rm Mm^2$ do not significantly suppress convective blueshift. Our area-dependent model reduces the amplitude of activity-induced RV variations by a factor of two. The present study highlights the need to identify a proxy that correlates specifically with large, bright magnetic regions on the surfaces of exoplanet-hosting stars.
[14]  [pdf] - 1863974
So close, so different: characterization of the K2-36 planetary system with HARPS-N
Comments: Accepted for publication on Astronomy $\&$ Astrophysics
Submitted: 2019-02-05
K2-36 is a K dwarf orbited by two small ($R_{\rm b}=1.43\pm0.08$ $R_\oplus$ and $R_{\rm c}=3.2\pm0.3$ $R_\oplus$), close-in ($a_{\rm b}$=0.022 AU and $a_{\rm c}$=0.054 AU) transiting planets discovered by Kepler/K2. They are representatives of two families of small planets ($R_{\rm p}$<4 $R_\oplus$) recently emerged from the analysis of Kepler data, with likely a different structure, composition and evolutionary pathways. We revise the fundamental stellar parameters and the sizes of the planets, and provide the first measurement of their masses and bulk densities, which we use to infer their structure and composition. We observed K2-36 with the HARPS-N spectrograph over $\sim$3.5 years, collecting 81 useful radial velocity measurements. The star is active, with evidence for increasing levels of magnetic activity during the observing time span. The radial velocity scatter is $\sim$17 \ms due to the stellar activity contribution, which is much larger that the semi-amplitudes of the planetary signals. We tested different methods for mitigating the stellar activity contribution to the radial velocity time variations and measuring the planet masses with good precision. We found that K2-36 is likely a $\sim$1 Gyr old system, and by treating the stellar activity through a Gaussian process regression, we measured the planet masses $m_{\rm b}$=3.9$\pm$1.1 $M_\oplus$ and $m_{\rm c}$=7.8$\pm$2.3 $M_\oplus$. The derived planet bulk densities $\rho_{\rm b}$=7.2$^{+2.5}_{-2.1}$ $g/cm^{3}$ and $\rho_{\rm c}$=1.3$^{+0.7}_{-0.5}$ $g/cm^{3}$ point out that K2-36\,b has a rocky, Earth-like composition, and K2-36\,c is a low-density sub-Neptune. Composed of two planets with similar orbital separations but different densities, K2-36 represents an optimal laboratory for testing the role of the atmospheric escape in driving the evolution of close-in, low-mass planets after $\sim$1 Gyr from their formation.
[15]  [pdf] - 1826285
A giant impact as the likely origin of different twins in the Kepler-107 exoplanet system
Comments: Published in Nature Astronomy on 4 February 2019, 35 pages including Supplementary Information material
Submitted: 2019-02-04
Measures of exoplanet bulk densities indicate that small exoplanets with radius less than 3 Earth radii ($R_\oplus$) range from low-density sub-Neptunes containing volatile elements to higher density rocky planets with Earth-like or iron-rich (Mercury-like) compositions. Such astonishing diversity in observed small exoplanet compositions may be the product of different initial conditions of the planet-formation process and/or different evolutionary paths that altered the planetary properties after formation. Planet evolution may be especially affected by either photoevaporative mass loss induced by high stellar X-ray and extreme ultraviolet (XUV) flux or giant impacts. Although there is some evidence for the former, there are no unambiguous findings so far about the occurrence of giant impacts in an exoplanet system. Here, we characterize the two innermost planets of the compact and near-resonant system Kepler-107. We show that they have nearly identical radii (about $1.5-1.6~R_\oplus$), but the outer planet Kepler-107c is more than twice as dense (about $12.6~\rm g\,cm^{-3}$) as the innermost Kepler-107b (about $5.3~\rm g\,cm^{-3}$). In consequence, Kepler-107c must have a larger iron core fraction than Kepler-107b. This imbalance cannot be explained by the stellar XUV irradiation, which would conversely make the more-irradiated and less-massive planet Kepler-107b denser than Kepler-107c. Instead, the dissimilar densities are consistent with a giant impact event on Kepler-107c that would have stripped off part of its silicate mantle. This hypothesis is supported by theoretical predictions from collisional mantle stripping, which match the mass and radius of Kepler-107c.
[16]  [pdf] - 1875304
Giant planets and brown dwarfs on wide orbits: a code comparison project
Comments: submitted to MNRAS. This is version 2 of the paper after considering referees comments and changes
Submitted: 2019-01-23
Gas clumps formed within massive gravitationally unstable circumstellar discs are potential seeds of gas giant planets, brown dwarfs and companion stars. Simulations show that competition between three processes -- migration, gas accretion and tidal disruption -- establishes what grows from a given seed. Here we investigate the robustness of numerical modelling of clump migration and accretion with the codes PHANTOM, GADGET, SPHINX, SEREN, GIZMO-MFM, SPHNG and FARGO. The test problem comprises a clump embedded in a massive disc at an initial separation of 120 AU. There is a general qualitative agreement between the codes, but the quantitative agreement in the planet migration rate ranges from $\sim 10$% to $\sim 50$%, depending on the numerical setup. We find that the artificial viscosity treatment and the sink particle prescription may account for much of the differences between the codes. In order to understand the wider implications of our work, we also attempt to reproduce the planet evolution tracks from our hydrodynamical simulations with prescriptions from three previous population synthesis studies. We find that the disagreement amongst the population synthesis models is far greater than that between our hydrodynamical simulations. The results of our code comparison project are therefore encouraging in that uncertainties in the given problem are probably dominated by the physics not yet included in the codes rather than by how hydrodynamics is modelled in them.
[17]  [pdf] - 1817204
HARPS-N radial velocities confirm the low densities of the Kepler-9 planets
Comments: 11 pages, 9 figures, 5 tables. Accepted for publication in MNRAS
Submitted: 2019-01-16
We investigated the discrepancy between planetary mass determination using the transit timing variations (TTVs) and radial velocities (RVs), by analysing the multi-planet system Kepler-9. Despite being the first system characterised with TTVs, there are several discrepant solutions in the literature, with those reporting lower planetary densities being apparently in disagreement with high-precision RV observations. To resolve this, we gathered HARPS-N RVs at epochs that maximised the difference between the predicted RV curves from discrepant solutions in the literature. We also re-analysed the full Kepler data-set and performed a dynamical fit, within a Bayesian framework, using the newly derived central and duration times of the transits. We compared these results with the RV data and found that our solution better describes the RV observations, despite the masses of the planets being nearly half that presented in the discovery paper. We therefore confirm that the TTV method can provide mass determinations that agree with those determined using high-precision RVs. The low densities of the planets place them in the scarcely populated region of the super-Neptunes / inflated sub-Saturns in the mass-radius diagram.
[18]  [pdf] - 1818822
Masses and radii for the three super-Earths orbiting GJ 9827, and implications for the composition of small exoplanets
Comments: 16 pages, 9 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society
Submitted: 2018-12-18, last modified: 2019-01-11
Super-Earths belong to a class of planet not found in the Solar System, but which appear common in the Galaxy. Given that some super-Earths are rocky, while others retain substantial atmospheres, their study can provide clues as to the formation of both rocky planets and gaseous planets, and - in particular - they can help to constrain the role of photo-evaporation in sculpting the exoplanet population. GJ 9827 is a system already known to host 3 super-Earths with orbital periods of 1.2, 3.6 and 6.2 days. Here we use new HARPS-N radial velocity measurements, together with previously published radial velocities, to better constrain the properties of the GJ 9827 planets. Our analysis can't place a strong constraint on the mass of GJ 9827 c, but does indicate that GJ 9827 b is rocky with a composition that is probably similar to that of the Earth, while GJ 9827 d almost certainly retains a volatile envelope. Therefore, GJ 9827 hosts planets on either side of the radius gap that appears to divide super-Earths into pre-dominantly rocky ones that have radii below $\sim 1.5 R_\oplus$, and ones that still retain a substantial atmosphere and/or volatile components, and have radii above $\sim 2 R_\oplus$. That the less heavily irradiated of the 3 planets still retains an atmosphere, may indicate that photoevaporation has played a key role in the evolution of the planets in this system.
[19]  [pdf] - 1791071
The fate of formamide in a fragmenting protoplanetary disc
Comments: 10 pages, 5 figures, accepted for publication in ApJ
Submitted: 2018-09-26
Recent high-sensitivity observations carried out with ALMA have revealed the presence of complex organic molecules (COMs) such as methyl cyanide (CH$_{\rm 3}$CN) and methanol (CH$_{\rm 3}$OH) in relatively evolved protoplanetary discs. The behaviour and abundance of COMs in earlier phases of disc evolution remains unclear. Here we combine a smoothed particle hydrodynamics simulation of a fragmenting, gravitationally unstable disc with a gas-grain chemical code. We use this to investigate the evolution of formamide (NH$_{\rm 2}$CHO), a pre-biotic species, in both the disc and in the fragments that form within it. Our results show that formamide remains frozen onto grains in the majority of the disc where the temperatures are $<$100 K, with a predicted solid-phase abundance that matches those observed in comets. Formamide is present in the gas-phase in three fragments as a result of the high temperatures ($\geq$200\,K), but remains in the solid-phase in one colder ($\leq$150 K) fragment. The timescale over which this occurs is comparable to the dust sedimentation timescales, suggesting that any rocky core which is formed would inherit their formamide content directly from the protosolar nebula.
[20]  [pdf] - 1757275
K2-263b: A 50-day period sub-Neptune with a mass measurement using HARPS-N
Comments: 10 pages, 6 figures, accepted for publication in MNRAS. v2: EPIC ID changed to assigned K2 name
Submitted: 2018-08-24, last modified: 2018-09-26
This paper reports on the validation and mass measurement of K2-263b, a sub-Neptune orbiting a quiet G9V star. Using K2 data from campaigns C5 and C16, we find this planet to have a period of $50.818947\pm 0.000094$ days and a radius of $2.41\pm0.12$ R$_{\oplus}$. We followed this system with HARPS-N to obtain 67 precise radial velocities. A combined fit of the transit and radial velocity data reveals that K2-263b has a mass of $14.8\pm3.1$ M$_{\oplus}$. Its bulk density ($5.7_{-1.4}^{+1.6}$ g cm$^{-3}$) implies that this planet has a significant envelope of water or other volatiles around a rocky core. EPIC211682544b likely formed in a similar way as the cores of the four giant planets in our own Solar System, but for some reason, did not accrete much gas. The planetary mass was confirmed by an independent Gaussian process-based fit to both the radial velocities and the spectroscopic activity indicators. K2-263b belongs to only a handful of confirmed K2 exoplanets with periods longer than 40 days. It is among the longest periods for a small planet with a precisely determined mass using radial velocities.
[21]  [pdf] - 1674979
An accurate mass determination for Kepler-1655b, a moderately-irradiated world with a significant volatile envelope
Comments: 21 pages, 14 figures, 4 tables; accepted for publication in the Astronomical Journal
Submitted: 2018-03-23
We present the confirmation of a small, moderately-irradiated (F = 155 +/- 7 Fearth) Neptune with a substantial gas envelope in a P=11.8728787+/-0.0000085-day orbit about a quiet, Sun-like G0V star Kepler-1655. Based on our analysis of the Kepler light curve, we determined Kepler-1655b's radius to be 2.213+/-0.082 Rearth. We acquired 95 high-resolution spectra with TNG/HARPS-N, enabling us to characterize the host star and determine an accurate mass for Kepler-1655b of 5.0+3.1/-2.8 Mearth via Gaussian-process regression. Our mass determination excludes an Earth-like composition with 98\% confidence. Kepler-1655b falls on the upper edge of the evaporation valley, in the relatively sparsely occupied transition region between rocky and gas-rich planets. It is therefore part of a population of planets that we should actively seek to characterize further.
[22]  [pdf] - 1648939
Is the spiral morphology of the Elias 2-27 circumstellar disc due to gravitational instability?
Comments: 12 pages, 5 figures
Submitted: 2018-02-26
A recent ALMA observation of the Elias 2-27 system revealed a two-armed structure extending out to ~300 au in radius. The protostellar disc surrounding the central star is unusually massive, raising the possibility that the system is gravitationally unstable. Recent work has shown that the observed morphology of the system can be explained by disc self-gravity, so we examine the physical properties of the disc necessary to detect self-gravitating spiral waves. Using three-dimensional Smoothed Particle Hydrodynamics, coupled with radiative transfer and synthetic ALMA imaging, we find that observable spiral structure can only be explained by self-gravity if the disc has a low opacity (and therefore efficient cooling), and is minimally supported by external irradiation. This corresponds to a very narrow region of parameter space, suggesting that, although it is possible for the spiral structure to be due to disc self-gravity, other explanations, such as an external perturbation, may be preferred.
[23]  [pdf] - 1717019
Eyes on K2-3: A system of three likely sub-Neptunes characterized with HARPS-N and HARPS
Comments: 14 pages, 15 figures, accepted for publication in A&A
Submitted: 2018-02-22
M-dwarf stars are promising targets for identifying and characterizing potentially habitable planets. K2-3 is a nearby (45 pc), early-type M dwarf hosting three small transiting planets, the outermost of which orbits close to the inner edge of the stellar (optimistic) habitable zone. The K2-3 system is well suited for follow-up characterization studies aimed at determining accurate masses and bulk densities of the three planets. Using a total of 329 radial velocity measurements collected over 2.5 years with the HARPS-N and HARPS spectrographs and a proper treatment of the stellar activity signal, we aim to improve measurements of the masses and bulk densities of the K2-3 planets. We use our results to investigate the physical structure of the planets. We analyse radial velocity time series extracted with two independent pipelines by using Gaussian process regression. We adopt a quasi-periodic kernel to model the stellar magnetic activity jointly with the planetary signals. We use Monte Carlo simulations to investigate the robustness of our mass measurements of K2-3\,c and K2-3\,d, and to explore how additional high-cadence radial velocity observations might improve them. Despite the stellar activity component being the strongest signal present in the radial velocity time series, we are able to derive masses for both planet b ($M_{\rm b}=6.6\pm1.1$ $M_{\rm \oplus}$) and planet c ($M_{\rm c}=3.1^{+1.3}_{-1.2}$ $M_{\rm \oplus}$). The Doppler signal due to K2-3\,d remains undetected, likely because of its low amplitude compared to the radial velocity signal induced by the stellar activity. The closeness of the orbital period of K2-3\,d to the stellar rotation period could also make the detection of the planetary signal complicated. [...]
[24]  [pdf] - 1637648
An ultra-short period rocky super-Earth with a secondary eclipse and a Neptune-like companion around K2-141
Comments: 16 pages, 10 figures., accepted for publication in AJ
Submitted: 2018-01-10
Ultra-short period (USP) planets are a class of low mass planets with periods shorter than one day. Their origin is still unknown, with photo-evaporation of mini-Neptunes and in-situ formation being the most credited hypotheses. Formation scenarios differ radically in the predicted composition of USP planets, it is therefore extremely important to increase the still limited sample of USP planets with precise and accurate mass and density measurements. We report here the characterization of an USP planet with a period of 0.28 days around K2-141 (EPIC 246393474), and the validation of an outer planet with a period of 7.7 days in a grazing transit configuration. We derived the radii of the planets from the K2 light curve and used high-precision radial velocities gathered with the HARPS-N spectrograph for mass measurements. For K2-141b we thus inferred a radius of $1.51\pm0.05~R_\oplus$ and a mass of $5.08\pm0.41~M_\oplus$, consistent with a rocky composition and lack of a thick atmosphere. K2-141c is likely a Neptune-like planet, although due to the grazing transits and the non-detection in the RV dataset, we were not able to put a strong constraint on its density. We also report the detection of secondary eclipses and phase curve variations for K2-141b. The phase variation can be modeled either by a planet with a geometric albedo of $0.30 \pm 0.06$ in the Kepler bandpass, or by thermal emission from the surface of the planet at $\sim$3000K. Only follow-up observations at longer wavelengths will allow us to distinguish between these two scenarios.
[25]  [pdf] - 1626459
On fragmentation of turbulent self-gravitating discs in the long cooling time regime
Comments: 12 pages, 12 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society
Submitted: 2017-12-15
It has recently been suggested that in the presence of driven turbulence discs may be much less stable against gravitational collapse than their non turbulent analogs, due to stochastic density fluctuations in turbulent flows. This mode of fragmentation would be especially important for gas giant planet formation. Here we argue, however, that stochastic density fluctuations due to turbulence do not enhance gravitational instability and disc fragmentation in the long cooling time limit appropriate for planet forming discs. These fluctuations evolve adiabatically and dissipate away by decompression faster than they could collapse. We investigate these issues numerically in 2D via shearing box simulations with driven turbulence and also in 3D with a model of instantaneously applied turbulent velocity kicks. In the former setting turbulent driving leads to additional disc heating that tends to make discs more, rather than less, stable to gravitational instability. In the latter setting, the formation of high density regions due to convergent velocity kicks is found to be quickly followed by decompression, as expected. We therefore conclude that driven turbulence does not promote disc fragmentation in protoplanetary discs and instead tends to make the discs more stable. We also argue that sustaining supersonic turbulence is very difficult in discs that cool slowly.
[26]  [pdf] - 1593660
Precise Masses in the WASP-47 System
Comments: 15 pages, 3 figures, 3 tables. Accepted in AJ
Submitted: 2017-09-29
We present precise radial velocity observations of WASP-47, a star known to host a hot Jupiter, a distant Jovian companion, and, uniquely, two additional transiting planets in short-period orbits: a super-Earth in a ~19 hour orbit, and a Neptune in a ~9 day orbit. We analyze our observations from the HARPS-N spectrograph along with previously published data to measure the most precise planet masses yet for this system. When combined with new stellar parameters and reanalyzed transit photometry, our mass measurements place strong constraints on the compositions of the two small planets. We find unlike most other ultra-short-period planets, the inner planet, WASP-47 e, has a mass (6.83 +/- 0.66 Me) and radius (1.810 +/- 0.027 Re) inconsistent with an Earth-like composition. Instead, WASP-47 e likely has a volatile-rich envelope surrounding an Earth-like core and mantle. We also perform a dynamical analysis to constrain the orbital inclination of WASP-47 c, the outer Jovian planet. This planet likely orbits close to the plane of the inner three planets, suggesting a quiet dynamical history for the system. Our dynamical constraints also imply that WASP-47 c is much more likely to transit than a geometric calculation would suggest. We calculate a transit probability for WASP-47 c of about 10%, more than an order of magnitude larger than the geometric transit probability of 0.6%.
[27]  [pdf] - 1584322
Three's Company: An additional non-transiting super-Earth in the bright HD 3167 system, and masses for all three planets
Comments: 22 pages, 14 figures, 5 tables. Submitted to AJ March 3rd, 2017. Accepted April 28th, 2017. In press
Submitted: 2017-06-06
HD 3167 is a bright (V = 8.9), nearby K0 star observed by the NASA K2 mission (EPIC 220383386), hosting two small, short-period transiting planets. Here we present the results of a multi-site, multi-instrument radial velocity campaign to characterize the HD 3167 system. The masses of the transiting planets are 5.02+/-0.38 MEarth for HD 3167 b, a hot super-Earth with a likely rocky composition (rho_b = 5.60+2.15-1.43 g/cm^3), and 9.80+1.30-1.24 MEarth for HD 3167 c, a warm sub-Neptune with a likely substantial volatile complement (rho_c = 1.97+0.94-0.59 g/cm^3). We explore the possibility of atmospheric composition analysis and determine that planet c is amenable to transmission spectroscopy measurements, and planet b is a potential thermal emission target. We detect a third, non-transiting planet, HD 3167 d, with a period of 8.509+/-0.045 d (between planets b and c) and a minimum mass of 6.90+/-0.71 MEarth. We are able to constrain the mutual inclination of planet d with planets b and c: we rule out mutual inclinations below 1.3 degrees as we do not observe transits of planet d. From 1.3-40 degrees, there are viewing geometries invoking special nodal configurations which result in planet d not transiting some fraction of the time. From 40-60 degrees, Kozai-Lidov oscillations increase the system's instability, but it can remain stable for up to 100Myr. Above 60 degrees, the system is unstable. HD 3167 promises to be a fruitful system for further study and a preview of the many exciting systems expected from the upcoming NASA TESS mission.
[28]  [pdf] - 1583598
Identifying and Analysing Protostellar Disc Fragments in Smoothed Particle Hydrodynamics Simulations
Comments: 22 pages, 22 figures
Submitted: 2017-05-18
We present a new method of identifying protostellar disc fragments in a simulation based on density derivatives, and analyse our data using this and the existing CLUMPFIND method, which is based on an ordered search over all particles in gravitational potential energy. Using smoothed particle hydrodynamics, we carry out 9 simulations of a $0.25$ M$_{\odot}$ disc around a 1 M$_{\odot}$ star, all of which fragment to form at least 2 bound objects. We find that when using all particles ordered in gravitational potential space, only fragments that survive the duration of the simulation are detected. When we use the density derivative method, all fragments are detected, so the two methods are complementary, as using the two methods together allows us to identify all fragments, and to then determine those that are likely to be destroyed. We find a tentative empirical relationship between the dominant azimuthal wavenumber in the disc $m$ and the maximum semi-major axis a fragment may achieve in a simulation, such that $a_{\mathrm{max}}\propto\frac{1}{m}$. We find the fragment destruction rate to be around half that predicted from population synthesis models. This is due to fragment-fragment interactions in the early gas phase of the disc, which can cause scattering and eccentricity pumping on short timescales, and affects the fragment's internal structure. We therefore caution that measurements of eccentricity as a function of semi-major axis may not necessarily constrain the formation mechanism of giant planets and brown dwarfs.
[29]  [pdf] - 1571060
The Kepler-19 system: a thick-envelope super-Earth with two Neptune-mass companions characterized using Radial Velocities and Transit Timing Variations
Comments: Accepted by AJ, 15 pages, 10 figures, 6 tables. Mass-Radius diagram for Exoplanets and code available at and
Submitted: 2017-03-20
We report a detailed characterization of the Kepler-19 system. This star was previously known to host a transiting planet with a period of 9.29 days, a radius of 2.2 R$_\oplus$ and an upper limit on the mass of 20 M$_\oplus$. The presence of a second, non-transiting planet was inferred from the transit time variations (TTVs) of Kepler-19b, over 8 quarters of Kepler photometry, although neither mass nor period could be determined. By combining new TTVs measurements from all the Kepler quarters and 91 high-precision radial velocities obtained with the HARPS-N spectrograph, we measured through dynamical simulations a mass of $8.4 \pm 1.6$ M$_\oplus$ for Kepler-19b. From the same data, assuming system coplanarity, we determined an orbital period of 28.7 days and a mass of $13.1 \pm 2.7$ M$_\oplus$ for Kepler-19c and discovered a Neptune-like planet with a mass of $20.3 \pm 3.4$ M$_\oplus$ on a 63 days orbit. By comparing dynamical simulations with non-interacting Keplerian orbits, we concluded that neglecting interactions between planets may lead to systematic errors that could hamper the precision in the orbital parameters when the dataset spans several years. With a density of $4.32 \pm 0.87$ g cm$^{-3}$ ($0.78 \pm 0.16$ $\rho_\oplus$) Kepler-19b belongs to the group of planets with a rocky core and a significant fraction of volatiles, in opposition to low-density planets characterized by transit-time variations only and the increasing number of rocky planets with Earth-like density. Kepler-19 joins the small number of systems that reconcile transit timing variation and radial velocity measurements.
[30]  [pdf] - 1581953
The VLT/NaCo large program to probe the occurrence of exoplanets and brown dwarfs at wide orbits. IV. Gravitational instability rarely forms wide, giant planets
Comments: 20 pages, 11 figures. Accepted for publication in A&A
Submitted: 2017-03-15
Understanding the formation and evolution of giant planets ($\ge$1 $M_{Jup}$) at wide orbital separation ($\ge$5 AU) is one of the goals of direct imaging. Over the past 15 years, many surveys have placed strong constraints on the occurrence rate of wide-orbit giants, mostly based on non-detections, but very few have tried to make a direct link with planet formation theories. In the present work, we combine the results of our previously published VLT/NaCo large program with the results of 12 past imaging surveys to constitute a statistical sample of 199 FGK stars within 100 pc, including 3 stars with sub-stellar companions. Using Monte Carlo simulations and assuming linear flat distributions for the mass and semi-major axis of planets, we estimate the sub-stellar companion frequency to be within 0.75-5.7% at the 68% confidence level (CL) within 20-300 AU and 0.5-75 $M_{Jup}$, which is compatible with previously published results. We also compare our results with the predictions of state-of-the-art population synthesis models based on the gravitational instability (GI) formation scenario by Forgan & Rice (2013), with and without scattering. We estimate that in both the scattered and non-scattered populations, we would be able to detect more than 30% of companions in the 1-75 $M_{Jup}$ range (95% CL). With the 3 sub-stellar detections in our sample, we estimate the fraction of stars that host a planetary system formed by GI to be within 1.0-8.6% (95% CL). We also conclude that even though GI is not common, it predicts a mass distribution of wide-orbit massive companions that is much closer to what is observed than what the core accretion scenario predicts. Finally, we associate the present paper with the release of the Direct Imaging Virtual Archive (DIVA,, a public database that aims at gathering the results of past, present, and future direct imaging surveys.
[31]  [pdf] - 1531727
Kepler-21b: A rocky planet around a V = 8.25 magnitude star
Comments: 52 pages, 13 figures, 1 long table. Accepted for publication in ApJ
Submitted: 2016-09-24
HD 179070, aka Kepler-21, is a V = 8.25 F6IV star and the brightest exoplanet host discovered by Kepler. An early detailed analysis by Howell et al. (2012) of the first thirteen months (Q0 - Q5) of Kepler light curves revealed transits of a planetary companion, Kepler-21b, with a radius of about 1.60 +/- 0.04 R_earth and an orbital period of about 2.7857 days. However, they could not determine the mass of the planet from the initial radial velocity observations with Keck-HIRES, and were only able to impose a 2-sigma upper limit of 10 M_earth. Here we present results from the analysis of 82 new radial velocity observations of this system obtained with HARPS-N, together with the existing 14 HIRES data points. We detect the Doppler signal of Kepler-21b with a radial velocity semi-amplitude K = 2.00 +/- 0.65 m/s, which corresponds to a planetary mass of 5.1 +/- 1.7 M_earth. We also measure an improved radius for the planet of 1.639 (+0.019, -0.015) R_earth, in agreement with the radius reported by Howell et al. (2012). We conclude that Kepler-21b, with a density of 6.4 +/- 2.1 g/cm^3, belongs to the population of terrestrial planets with iron, magnesium silicate interiors, which have lost the majority of their envelope volatiles via stellar winds or gravitational escape. The radial velocity analysis presented in this paper serves as example of the type of analysis that will be necessary to confirm the masses of TESS small planet candidates.
[32]  [pdf] - 1521296
The Orbit and Mass of the Third Planet in the Kepler-56 System
Comments: 7 pages, 1 figure, 2 tables; accepted for publication in AJ. Minor edits made after referee report
Submitted: 2016-08-11, last modified: 2016-09-22
While the vast majority of multiple-planet systems have their orbital angular momentum axes aligned with the spin axis of their host star, Kepler-56 is an exception: its two transiting planets are coplanar yet misaligned by at least 40 degrees with respect to their host star. Additional follow-up observations of Kepler-56 suggest the presence of a massive, non-transiting companion that may help explain this misalignment. We model the transit data along with Keck/HIRES and HARPS-N radial velocity data to update the masses of the two transiting planets and infer the physical properties of the third, non-transiting planet. We employ a Markov Chain Monte Carlo sampler to calculate the best-fitting orbital parameters and their uncertainties for each planet. We find the outer planet has a period of 1002 $\pm$ 5 days and minimum mass of 5.61 $\pm$ 0.38 Jupiter masses. We also place a 95% upper limit of 0.80 m/s/yr on long-term trends caused by additional, more distant companions.
[33]  [pdf] - 1502043
Grand challenges in protoplanetary disc modelling
Comments: Accepted for publication in PASA. 25 pages
Submitted: 2016-08-03, last modified: 2016-09-20
The Protoplanetary Discussions conference --- held in Edinburgh, UK, from 7th --11th March 2016 --- included several open sessions led by participants. This paper reports on the discussions collectively concerned with the multiphysics modelling of protoplanetary discs, including the self-consistent calculation of gas and dust dynamics, radiative transfer and chemistry. After a short introduction to each of these disciplines in isolation, we identify a series of burning questions and grand challenges associated with their continuing development and integration. We then discuss potential pathways towards solving these challenges, grouped by strategical, technical and collaborative developments. This paper is not intended to be a review, but rather to motivate and direct future research and collaboration across typically distinct fields based on \textit{community driven input}, to encourage further progress in our understanding of circumstellar and protoplanetary discs.
[34]  [pdf] - 1521311
A 1.9 Earth radius rocky planet and the discovery of a non-transiting planet in the Kepler-20 system
Comments: 33 pages, 8 figures. Accepted for publication in AJ
Submitted: 2016-08-24
Kepler-20 is a solar-type star (V = 12.5) hosting a compact system of five transiting planets, all packed within the orbital distance of Mercury in our own Solar System. A transition from rocky to gaseous planets with a planetary transition radius of ~1.6 REarth has recently been proposed by several publications in the literature (Rogers 2015;Weiss & Marcy 2014). Kepler-20b (Rp ~ 1.9 REarth) has a size beyond this transition radius, however previous mass measurements were not sufficiently precise to allow definite conclusions to be drawn regarding its composition. We present new mass measurements of three of the planets in the Kepler-20 system facilitated by 104 radial velocity measurements from the HARPS-N spectrograph and 30 archival Keck/HIRES observations, as well as an updated photometric analysis of the Kepler data and an asteroseismic analysis of the host star (MStar = 0.948+-0.051 Msun and Rstar = 0.964+-0.018 Rsun). Kepler-20b is a 1.868+0.066-0.034 REarth planet in a 3.7 day period with a mass of 9.70+1.41-1.44 MEarth resulting in a mean density of 8.2+1.5-1.3 g/cc indicating a rocky composition with an iron to silicate ratio consistent with that of the Earth. This makes Kepler-20b the most massive planet with a rocky composition found to date. Furthermore, we report the discovery of an additional non-transiting planet with a minimum mass of 19.96+3.08-3.61 MEarth and an orbital period of ~34 days in the gap between Kepler-20f (P ~ 11 days) and Kepler-20d (P ~ 78 days).
[35]  [pdf] - 1388960
The evolution of self-gravitating accretion discs
Comments: 12 pages, 9 figures, accepted for publication in PASA as part of the special issue on "Disc dynamics and planet formation"
Submitted: 2016-02-26
It is quite likely that self-gravity will play an important role in the evolution of accretion discs, in particular those around young stars, and those around supermassive black holes. We summarise, here, our current understanding of the evolution of such discs, focussing more on discs in young stellar system, than on discs in active galactic nuclei. We consider the conditions under which such discs may fragment to form bound objects, and when they might, instead, be expected to settle into a quasi-steady, self-regulated state. We also discuss how this understanding may depend on the mass of the disc relative to the mass of the central object, and how it might depend on the presence of external irradiation. Additionally, we consider whether or not fragmentation might be stochastic, where we might expect it to occur in an actual protostellar disc, and if there is any evidence for fragmentation actually playing a role in the formation of planetary-mass bodies. Although there are still a number of outstanding issue, such as the convergence of simulations of self-gravitating discs, whether or not there is more than one mode of fragmentation, and quite what role self-gravitating discs may play in the planet formation process, our general understanding of these systems seems quite robust.
[36]  [pdf] - 1359284
Directly observing continuum emission from self-gravitating spiral waves
Comments: 13 pages, 9 figures
Submitted: 2016-02-03
We use a simple, self-consistent, self-gravitating semi-analytic disc model to conduct an examination of the parameter space in which self-gravitating discs may exist. We then use Monte-Carlo radiative transfer to generate synthetic ALMA images of these self-gravitating discs to determine the subset of this parameter space in which they generate non-axisymmetric structure that is potentially detectable by ALMA. Recently, several transition discs have been observed to have non-axisymmetric structure that extends out to large radii. It has been suggested that one possible origin of these asymmetries could be spiral density waves induced by disc self-gravity. We use our simple model to see if these discs exist in the region of parameter space where self-gravity could feasibly explain these spiral features. We find that for self-gravity to play a role in these systems typically requires a disc mass around an order of magnitude higher than the observed disc masses for the systems. The spiral amplitudes produced by self-gravity in the local approximation are relatively weak when compared to amplitudes produced by tidal interactions, or spirals launched at Lindblad resonances due to embedded planets in the disc. As such, we ultimately caution against diagnosing spiral features as being due to self-gravity, unless the disc exists in the very narrow region of parameter space where the spiral wave amplitudes are large enough to produce detectable features, but not so large as to cause the disc to fragment.
[37]  [pdf] - 1359194
Tensor classification of structure in smoothed particle hydrodynamics density fields
Comments: 14 pages, 28 figures, accepted for publication in MNRAS
Submitted: 2016-01-11
As hydrodynamic simulations increase in scale and resolution, identifying structures with non-trivial geometries or regions of general interest becomes increasingly challenging. There is a growing need for algorithms that identify a variety of different features in a simulation without requiring a "by-eye" search. We present tensor classification as such a technique for smoothed particle hydrodynamics (SPH). These methods have already been used to great effect in N-Body cosmological simulations, which require smoothing defined as an input free parameter. We show that tensor classification successfully identifies a wide range of structures in SPH density fields using its native smoothing, removing a free parameter from the analysis and preventing the need for tesselation of the density field, as required by some classification algorithms. As examples, we show that tensor classification using the tidal tensor and the velocity shear tensor successfully identifies filaments, shells and sheet structures in giant molecular cloud simulations, as well as spiral arms in discs. The relationship between structures identified using different tensors illustrates how different forces compete and co-operate to produce the observed density field. We therefore advocate the use of multiple tensors to classify structure in SPH simulations, to shed light on the interplay of multiple physical processes.
[38]  [pdf] - 1347708
The Kepler-454 System: A Small, Not-rocky Inner Planet, a Jovian World, and a Distant Companion
Comments: 40 pages, 13 figures, 6 pages; ApJ in press
Submitted: 2015-11-29
Kepler-454 (KOI-273) is a relatively bright (V = 11.69 mag), Sun-like starthat hosts a transiting planet candidate in a 10.6 d orbit. From spectroscopy, we estimate the stellar temperature to be 5687 +/- 50 K, its metallicity to be [m/H] = 0.32 +/- 0.08, and the projected rotational velocity to be v sin i <2.4 km s-1. We combine these values with a study of the asteroseismic frequencies from short cadence Kepler data to estimate the stellar mass to be 1.028+0:04-0:03 M_Sun, the radius to be 1.066 +/- 0.012 R_Sun and the age to be 5.25+1:41-1:39 Gyr. We estimate the radius of the 10.6 d planet as 2.37 +/- 0.13 R_Earth. Using 63 radial velocity observations obtained with the HARPS-N spectrograph on the Telescopio Nazionale Galileo and 36 observations made with the HIRES spectrograph at Keck Observatory, we measure the mass of this planet to be 6.8 +/- 1.4M_Earth. We also detect two additional non-transiting companions, a planet with a minimum mass of 4.46 +/- 0.12 M_J in a nearly circular 524 d orbit and a massive companion with a period >10 years and mass >12.1M_J . The twelve exoplanets with radii <2.7 R_Earth and precise mass measurements appear to fall into two populations, with those <1.6 R_Earth following an Earth-like composition curve and larger planets requiring a significant fraction of volatiles. With a density of 2.76 +/- 0.73 g cm-3, Kepler-454b lies near the mass transition between these two populations and requires the presence of volatiles and/or H/He gas.
[39]  [pdf] - 1300377
Spiral Arms in Gravitationally Unstable Protoplanetary Disks as Imaged in Scattered Light
Comments: 18 pages, 4 figures, 1 table, ApJ Letter accepted
Submitted: 2015-10-01
Combining 3D smoothed-particle hydrodynamics and Monte Carlo radiative transfer calculations, we examine the morphology of spiral density waves induced by gravitational instability (GI) in protoplanetary disks, as they would appear in direct images at near-infrared (NIR) wavelengths. We find that systems with disk-to-star-mass ratios q=M_disk/M_star that are ~0.25 or more may produce prominent spiral arms in NIR imaging, remarkably resembling features observed in the MWC 758 and SAO 206462 systems. The contrast of GI-induced arms at NIR wavelengths can reach a factor of ~3, and their pitch angles are about 10-15 degree. The dominant azimuthal wavenumber of GI-induced spiral arms roughly obeys m~1/q in the range 2<~1/q<~8. In particular, a massive disk with q~0.5 can exhibit grand-design m=2 spirals. GI-induced arms are in approximate corotation with the local disk, and may therefore trap dust particles by pressure drag. Although GI can produce NIR spiral arms with morphologies, contrasts, and pitch angles similar to those reported in recent observations, it also makes other demands that may or may not be satisfied in any given system. A GI origin requires that the spirals be relatively compact, on scales <~100 AU; that the disk be massive, q>~0.25; and that the accretion rate Mdot be high, on the order of 1e-6 solar mass per year.
[40]  [pdf] - 1296262
Disc fragmentation rarely forms planetary-mass objects
Comments: 8 pages, 5 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society
Submitted: 2015-08-26
It is now reasonably clear that disc fragmentation can only operate in the outer parts of protostellar discs ($r > 50$ au). It is also expected that any object that forms via disc fragmentation will have an initial mass greater than that of Jupiter. However, whether or not such a process actually operates, or can play a significant role in the formation of planetary-mass objects, is still unclear. We do have a few examples of directly imaged objects that may have formed in this way, but we have yet to constrain how often disc fragmentation may actually form such objects. What we want to consider here is whether or not we can constrain the likely population of planetary-mass objects formed via disc fragmentation by considering how a population of objects at large radii ($a > 50$) au - if they do exist - would evolve under perturbations from more distant stellar companions. We find that there is a specific region of parameter space to which such objects would be scattered and show that the known exoplanets in that region have properties more consistent with that of the bulk exoplanet population, than with having been formed via disc fragmentation at large radii. Along with the scarcity of directly-imaged objects at large radii, our results provide a similar, but independent, constraint on the frequency of objects formed via disc fragmentation.
[41]  [pdf] - 1317011
The HARPS-N Rocky Planet Search I. HD219134b: A transiting rocky planet in a multi-planet system at 6.5 pc from the Sun
Comments: Accepted for publication in A&A
Submitted: 2015-07-30
We present here the detection of a system of four low-mass planets around the bright (V=5.5) and close-by (6.5 pc) star HD219134. This is the first result of the Rocky Planet Search program with HARPS-N on the TNG in La Palma. The inner planet orbits the star in 3.0937 +/-0.0004 days, on a quasi-circular orbit with a semi-major axis of 0.0382 +/- 0.0003 AU. Spitzer observations allowed us to detect the transit of the planet in front of the star making HD219134b the nearest known transiting planet to date. From the amplitude of the radial-velocity variation (2.33 +/- 0.24 m/s) and observed depth of the transit (359 +/- 38 ppm), the planet mass and radius are estimated to be 4.46 +/- 0.47 M_{\oplus} and 1.606 +/- 0.086 R_{\oplus} leading to a mean density of 5.89 +/- 1.17 g/cc, suggesting a rocky composition. One additional planet with minimum mass of 2.67 +/- 0.59 M_{\oplus} moves on a close-in, quasi-circular orbit with a period of 6.765 +/- 0.005 days. The third planet in the system has a period of 46.78 +/- 0.16 days and a minimum mass of 8.7 +/- 1.1 M{\oplus}, at 0.234 +/- 0.002 AU from the star. Its eccentricity is 0.32 +/- 0.14. The period of this planet is close to the rotational period of the star estimated from variations of activity indicators (42.3 +/- 0.1 days). The planetary origin of the signal is, however, the preferred solution as no indication of variation at the corresponding frequency is observed for activity-sensitive parameters. Finally, a fourth additional longer-period planet of mass of 62 +/- 6 M_{\oplus} orbits the star in 1190 days, on an eccentric orbit (e=0.27 +/- 0.11) at a distance of 2.14 +/- 0.27 AU.
[42]  [pdf] - 1273216
The quest for cradles of life: using the fundamental metallicity relation to hunt for the most habitable type of galaxy
Comments: Accepted to APJ letters. Comments welcome
Submitted: 2015-07-15
The field of astrobiology has made huge strides in understanding the habitable zones around stars (Stellar Habitable Zones) where life can begin, sustain its existence and evolve into complex forms. A few studies have extended this idea by modelling galactic-scale habitable zones (Galactic Habitable Zones) for our Milky Way and specific elliptical galaxies. However, estimating the habitability for galaxies spanning a wide range of physical properties has so far remained an outstanding issue. Here, we present a "cosmobiological" framework that allows us to sift through the entire galaxy population in the local Universe and answer the question "Which type of galaxy is most likely to host complex life in the cosmos"? Interestingly, the three key astrophysical criteria governing habitability (total mass in stars, total metal mass and ongoing star formation rate) are found to be intricately linked through the "fundamental metallicity relation" as shown by SDSS (Sloan Digital Sky Survey) observations of more than a hundred thousand galaxies in the local Universe. Using this relation we show that metal-rich, shapeless giant elliptical galaxies at least twice as massive as the Milky Way (with a tenth of its star formation rate) can potentially host ten thousand times as many habitable (earth-like) planets, making them the most probable "cradles of life" in the Universe.
[43]  [pdf] - 1269301
The Gemini NICI Planet-Finding Campaign: Asymmetries in the HD 141569 disc
Comments: 14 pages, 13 figures, accepted to MNRAS
Submitted: 2015-04-20, last modified: 2015-05-05
We report here the highest resolution near-IR imaging to date of the HD 141569A disc taken as part of the NICI Science Campaign. We recover 4 main features in the NICI images of the HD 141569 disc discovered in previous HST imaging: 1) an inner ring / spiral feature. Once deprojected, this feature does not appear circular. 2) an outer ring which is considerably brighter on the western side compared to the eastern side, but looks fairly circular in the deprojected image. 3) an additional arc-like feature between the inner and outer ring only evident on the east side. In the deprojected image, this feature appears to complete the circle of the west side inner ring and 4) an evacuated cavity from 175 AU inwards. Compared to the previous HST imaging with relatively large coronagraphic inner working angles (IWA), the NICI coronagraph allows imaging down to an IWA of 0.3". Thus, the inner edge of the inner ring/spiral feature is well resolved and we do not find any additional disc structures within 175 AU. We note some additional asymmetries in this system. Specifically, while the outer ring structure looks circular in this deprojection, the inner bright ring looks rather elliptical. This suggests that a single deprojection angle is not appropriate for this system and that there may be an offset in inclination between the two ring / spiral features. We find an offset of 4+-2 AU between the inner ring and the star center, potentially pointing to unseen inner companions.
[44]  [pdf] - 1224016
Can Kozai-Lidov cycles explain Kepler-78b?
Comments: 10 pages, 9 figures, accepted for publication in MNRAS
Submitted: 2015-01-14, last modified: 2015-02-24
Kepler-78b is one of a growing sample of planets similar, in composition and size, to the Earth. It was first detected with NASA's \emph{Kepler} spacecraft and then characterised in more detail using radial velocity follow-up observations. Not only is its size very similar to that of the Earth ($1.2 R_\oplus$), it also has a very similar density ($5.6$ g cm$^{-2}$). What makes this planet particularly interesting is that it orbits its host star every $8.5$ hours, giving it an orbital distance of only $0.0089$ au. What we investigate here is whether or not such a planet could have been perturbed into this orbit by an outer companion on an inclined orbit. In this scenario, the outer perturber causes the inner orbit to undergo Kozai-Lidov cycles which, if the periapse comes sufficiently close to the host star, can then lead to the planet being tidally circularised into a close orbit. We find that this process can indeed produce such very-close-in planets within the age of the host star ($\sim 600 - 900$ Myr), but it is more likely to find such ultra-short-period planets around slightly older stars ($> 1$ Gyr). However, given the size of the Kepler sample and the likely binarity, our results suggest that Kepler-78b may indeed have been perturbed into its current orbit by an outer stellar companion. The likelihood of this happening, however, is low enough that other processes - such as planet-planet scattering - could also be responsible.
[45]  [pdf] - 1223821
The Mass of Kepler-93b and The Composition of Terrestrial Planets
Comments: 8 pages, 4 figures. Accepted for publication in ApJ
Submitted: 2014-12-30
Kepler-93b is a 1.478 +/- 0.019 Earth radius planet with a 4.7 day period around a bright (V=10.2), astroseismically-characterized host star with a mass of 0.911+/-0.033 solar masses and a radius of 0.919+/-0.011 solar radii. Based on 86 radial velocity observations obtained with the HARPS-N spectrograph on the Telescopio Nazionale Galileo and 32 archival Keck/HIRES observations, we present a precise mass estimate of 4.02+/-0.68 Earth masses. The corresponding high density of 6.88+/-1.18 g/cc is consistent with a rocky composition of primarily iron and magnesium silicate. We compare Kepler-93b to other dense planets with well-constrained parameters and find that between 1-6 Earth masses, all dense planets including the Earth and Venus are well-described by the same fixed ratio of iron to magnesium silicate. There are as of yet no examples of such planets with masses > 6 Earth masses: All known planets in this mass regime have lower densities requiring significant fractions of volatiles or H/He gas. We also constrain the mass and period of the outer companion in the Kepler-93 system from the long-term radial velocity trend and archival adaptive optics images. As the sample of dense planets with well-constrained masses and radii continues to grow, we will be able to test whether the fixed compositional model found for the seven dense planets considered in this paper extends to the full population of 1-6 Earth mass planets.
[46]  [pdf] - 1223702
Characterizing K2 Planet Discoveries: A super-Earth transiting the bright K-dwarf HIP 116454
Comments: 16 pages, 8 figures. Accepted by ApJ
Submitted: 2014-12-17
We report the first planet discovery from the two-wheeled Kepler (K2) mission: HIP 116454 b. The host star HIP 116454 is a bright (V = 10.1, K = 8.0) K1-dwarf with high proper motion, and a parallax-based distance of 55.2 +/- 5.4 pc. Based on high-resolution optical spectroscopy, we find that the host star is metal-poor with [Fe/H] = -.16 +/- .18, and has a radius R = 0.716 +/- .0024 R_sun and mass M = .775 +/- .027 Msun. The star was observed by the Kepler spacecraft during its Two-Wheeled Concept Engineering Test in February 2014. During the 9 days of observations, K2 observed a single transit event. Using a new K2 photometric analysis technique we are able to correct small telescope drifts and recover the observed transit at high confidence, corresponding to a planetary radius of Rp = 2.53 +/- 0.18 Rearth. Radial velocity observations with the HARPS-N spectrograph reveal a 11.82 +/- 1.33 Mearth planet in a 9.1 day orbit, consistent with the transit depth, duration, and ephemeris. Follow-up photometric measurements from the MOST satellite confirm the transit observed in the K2 photometry and provide a refined ephemeris, making HIP 116454 b amenable for future follow-up observations of this latest addition to the growing population of transiting super-Earths around nearby, bright stars.
[47]  [pdf] - 1223357
The Dynamical Fate of Self-Gravitating Disc Fragments After Tidal Downsizing
Comments: 11 pages, 16 figures, accepted for publication in MNRAS
Submitted: 2014-11-26
The gravitational instability model of planet/brown dwarf formation proposes that protostellar discs can fragment into objects with masses above a few Jupiter masses at large semimajor axis. Tidal downsizing may reduce both the object mass and semimajor axis. However, most studies of tidal downsizing end when the protostellar disc disperses, while the system is embedded in its parent star-forming region. To compare disc fragment descendants with exoplanet and brown dwarf observations, the subsequent dynamical evolution must be explored. We carry out N-Body integrations of fragment-fragment scattering in multi-object star systems, and star systems embedded in substructured clusters. In both cases, we use initial conditions generated by population synthesis models of tidal downsizing. The scattering simulations produce a wide range of eccentricities. The ejection rate is around 25%. The ejecta mass distribution is similar to that for all objects, with a velocity dispersion consistent with those produced by full hydrodynamic simulations. The semimajor axis distribution after scattering extends to parsec scales. In the cluster simulations, 13% of objects are ejected from their planetary system, and around 10% experience significant orbit modification. A small number of objects are recaptured on high eccentricity, high inclination orbits. The velocity distribution of ejecta is similar to that produced by fragment-fragment scattering. If fragment-fragment scattering and cluster stripping act together, then disc fragmentation should be efficient at producing free-floating substellar objects, and hence characterising the free-floating planet population will provide strong constraints on the frequency of disc fragmentation.
[48]  [pdf] - 903740
Characterization of the Kepler-101 planetary system with HARPS-N. A hot super-Neptune with an Earth-sized low-mass companion
Comments: 7 pages, 3 figures, accepted in A&A
Submitted: 2014-09-16
We report on the characterization of the Kepler-101 planetary system, thanks to a combined DE-MCMC analysis of Kepler data and forty radial velocities obtained with the HARPS-N spectrograph. This system was previously validated by Rowe et al. (2014) and is composed of a hot super-Neptune, Kepler-101b, and an Earth-sized planet, Kepler-101c. These two planets orbit the slightly evolved and metal-rich G-type star in 3.49 and 6.03 days, respectively. With mass $M_{\rm p}=51.1_{-4.7}^{+5.1}~M_{\oplus}$, radius $R_{\rm p}=5.77_{-0.79}^{+0.85}~R_{\oplus}$, and density $\rho_{\rm p}=1.45_{-0.48}^{+0.83} \rm g\;cm^{-3}$, Kepler-101b is the first fully-characterized super-Neptune, and its density suggests that heavy elements make up a significant fraction of its interior; more than $60\%$ of its total mass. Kepler-101c has a radius of $1.25_{-0.17}^{+0.19}~R_{\oplus}$, which implies the absence of any H/He envelope, but its mass could not be determined due to the relative faintness of the parent star for highly precise radial-velocity measurements ($K_{\rm p}=13.8$) and the limited number of radial velocities. The $1~\sigma$ upper limit, $M_{\rm p} < 3.8~M_{\oplus}$, excludes a pure iron composition with a $68.3\%$ probability. The architecture of the Kepler-101 planetary system - containing a close-in giant planet and an outer Earth-sized planet with a period ratio slightly larger than the 3:2 resonance - is certainly of interest for planet formation and evolution scenarios. This system does not follow the trend, seen by Ciardi et al. (2013), that in the majority of Kepler systems of planet pairs with at least one Neptune-size or larger planet, the larger planet has the longer period.
[49]  [pdf] - 1209846
The Kepler-10 planetary system revisited by HARPS-N: A hot rocky world and a solid Neptune-mass planet
Comments: 44 pages, 8 figures, accepted for publication in ApJ
Submitted: 2014-05-30
Kepler-10b was the first rocky planet detected by the Kepler satellite and con- firmed with radial velocity follow-up observations from Keck-HIRES. The mass of the planet was measured with a precision of around 30%, which was insufficient to constrain models of its internal structure and composition in detail. In addition to Kepler-10b, a second planet transiting the same star with a period of 45 days was sta- tistically validated, but the radial velocities were only good enough to set an upper limit of 20 Mearth for the mass of Kepler-10c. To improve the precision on the mass for planet b, the HARPS-N Collaboration decided to observe Kepler-10 intensively with the HARPS-N spectrograph on the Telescopio Nazionale Galileo on La Palma. In to- tal, 148 high-quality radial-velocity measurements were obtained over two observing seasons. These new data allow us to improve the precision of the mass determina- tion for Kepler-10b to 15%. With a mass of 3.33 +/- 0.49 Mearth and an updated radius of 1.47 +0.03 -0.02 Rearth, Kepler-10b has a density of 5.8 +/- 0.8 g cm-3, very close to the value -0.02 predicted by models with the same internal structure and composition as the Earth. We were also able to determine a mass for the 45-day period planet Kepler-10c, with an even better precision of 11%. With a mass of 17.2 +/- 1.9 Mearth and radius of 2.35 +0.09 -0.04 Rearth, -0.04 Kepler-10c has a density of 7.1 +/- 1.0 g cm-3. Kepler-10c appears to be the first strong evidence of a class of more massive solid planets with longer orbital periods.
[50]  [pdf] - 1202880
Forming misaligned stellar discs around a massive black hole: Cloud infall in the Galactic Centre
Comments: To appear in the Proceedings of the IAU Symposium 303: The Galactic Center: Feeding and Feedback in a Normal Galactic Nucleus
Submitted: 2014-01-20
The innermost parsec around Sgr A* has been found to play host to two discs or streamers of O and W-R stars. They are misaligned by an angle approaching 90 degrees. That the stars are approximately coeval indicates that they formed in the same event rather than independently. We have performed SPH simulations of the infall of a single prolate cloud towards a massive black hole. As the cloud is disrupted, the large spread in angular momentum can, if conditions allow, lead to the creation of misaligned gas discs. In turn, stars may form within those discs. We are now investigating the origins of these clouds in the Galactic Centre (GC) region.
[51]  [pdf] - 740151
An Earth-sized planet with an Earth-like density
Comments: 10 Pages, 8 figures, 5 tables, data. Nature, 2013 (published online Oct. 30)
Submitted: 2013-10-29
Kepler-78 (KIC 8435766) was identified by Sanchis-Ojeda et al. (2013) as harbouring a transiting planet of 1.16 times the size of the Earth and an orbital period of only 8.5 hours. While the exquisite Kepler photometry was able to determine its radius and period, the mass of the planet (and thus its mean density) remained unconstrained in the absence of precise radial-velocity measurements. Here we present an accurate mass measurement of Kepler-78b using the HARPS-N spectrograph, recently installed on the Telescopio Nazionale Galileo (INAF) at the Roque de los Muchachos Observatory, La Palma, Spain. These new data yield a mass of 1.86 Earth masses. The resulting mean density of the planet is 5.57 grams per cubic centimetre, which is similar to that of the Earth and implies a composition of iron and rock. Kepler-78b, which orbits a Sun-like star called Kepler 78 located in the Cygnus constellation at a distance of about 400 light years from us, is now the smallest exoplanet for which both the mass and radius are known accurately.
[52]  [pdf] - 1166330
Misaligned streamers around a galactic centre black hole from a single cloud's infall
Comments: 14 pages, 17 figures, 1 table, accepted for publication in MNRAS
Submitted: 2013-04-30
We follow the near radial infall of a prolate cloud onto a 4 x 10^6 Msun supermassive black hole in the Galactic Centre using smoothed particle hydrodynamics (SPH). We show that a prolate cloud oriented perpendicular to its orbital plane naturally produces a spread in angular momenta in the gas which can translate into misaligned discs as is seen in the young stars orbiting Sagittarius A*. A turbulent or otherwise highly structured cloud is necessary to avoid cancelling too much angular momentum through shocks at closest approach. Our standard model of a 2 x 10^4 Msun gas cloud brought about the formation of a disc within 0.3 pc from the black hole and a larger, misaligned streamer at 0.5 pc. A total of 1.5 x 10^4 Msun of gas formed these structures. Our exploration of the simulation parameter space showed that when star formation occurred, it resulted in top-heavy IMFs with stars on eccentric orbits with semi-major axes 0.02 to 0.3 pc and inclinations following the gas discs and streamers. We suggest that the single event of an infalling prolate cloud can explain the occurrence of multiple misaligned discs of young stars.
[53]  [pdf] - 1166264
A Self-Gravitating Disc Around L1527 IRS?
Comments: 6 pages, 10 figures, accepted for publication in MNRAS
Submitted: 2013-04-27
Recent observations of the Class 0 protostar L1527 IRS have revealed a rotationally supported disc with an outer radius of at least 100 au. Measurements of the integrated flux at 870 microns suggest a disc mass that is too low for gravitational instability to govern angular momentum transport. However, if parts of the disc are optically thick at sub-mm wavelengths, the sub-mm fluxes will underestimate the disc mass, and the disc's actual mass may be substantially larger, potentially sufficient to be self-gravitating. We investigate this possibility using simple self-gravitating disc models. To match the observed mass accretion rates requires a disc-to-star mass ratio of at least ~0.5, which produces sub-mm fluxes that are similar to those observed for L1527 IRS in the absence of irradiation from the envelope or central star. If irradiation is significant, then the predicted fluxes exceed the observed fluxes by around an order of magnitude. Our model also indicates that the stresses produced by the gravitational instability are low enough to prevent disc fragmentation. As such, we conclude that observations do not rule out the possibility that the disc around L1527 IRS is self-gravitating, but it is more likely that despite being a very young system, this disc may already have left the self-gravitating phase.
[54]  [pdf] - 1166065
Towards a Population Synthesis Model of Objects formed by Self-Gravitating Disc Fragmentation and Tidal Downsizing
Comments: 19 pages, 28 figures, accepted for publication in MNRAS
Submitted: 2013-04-17
Recently, the gravitational instability (GI) model of giant planet and brown dwarf formation has been revisited and recast into what is often referred to as the "tidal downsizing" hypothesis. The fragmentation of self-gravitating protostellar discs into gravitationally bound embryos - with masses of a few to tens of Jupiter masses, at semi major axes above 30 - 40 au - is followed by grain sedimentation inside the embryo, radial migration towards the central star and tidal disruption of the embryo's upper layers. The properties of the resultant object depends sensitively on the timescales upon which each process occurs. Therefore, GI followed by tidal downsizing can theoretically produce objects spanning a large mass range, from terrestrial planets to giant planets and brown dwarfs. Whether such objects can be formed in practice, and what proportions of the observed population they would represent, requires a more involved statistical analysis. We present a simple population synthesis model of star and planet formation via GI and tidal downsizing. We couple a semi-analytic model of protostellar disc evolution to analytic calculations of fragmentation, initial embryo mass, grain growth and sedimentation, embryo migration and tidal disruption. While there are key pieces of physics yet to be incorporated, it represents a first step towards a mature statistical model of GI and tidal downsizing as a mode of star and planet formation. We show results from four runs of the population synthesis model, varying the opacity law and the strength of migration, as well as investigating the effect of disc truncation during the fragmentation process. Our early results suggest that GI plus tidal downsizing is not the principal mode of planet formation, but remains an excellent means of forming gas giant planets, brown dwarfs and low mass stars at large semimajor axes. (Abridged)
[55]  [pdf] - 1158852
The Effect of Irradiation on the Jeans Mass in Fragmenting Self-Gravitating Protostellar Discs
Comments: 8 pages, 14 figures, accepted for publication in MNRAS
Submitted: 2013-01-07
When a self-gravitating disc is subject to irradiation, its propensity to fragmentation will be affected. The strength of self-gravitating disc stresses is expected to dictate disc fragmentation: as the strength of these torques typically decrease with increasing sound speed, it is reasonable to assume, to first-order, that disc fragmentation is suppressed when compared to the non-irradiated case, although previous work has shown that the details are complicated by the source of the irradiation. We expand on previous analysis of the Jeans mass inside spiral structures in self-gravitating discs, incorporating the effects of stellar irradiation and background irradiation. If irradiation is present, fragmentation is suppressed for marginally unstable discs at low accretion rates (compared to the no-irradiation case), but these lower accretion rates correspond to higher mass discs. Fragmentation can still occur for high accretion rates, but is consequently suppressed at lower disc surface densities, and the subsequent Jeans mass is boosted. These results further bolster the consensus that, without subsequent fragment disruption or mass loss, the gravitational instability is more likely to form brown dwarfs and low-mass stars than gas giant planets.
[56]  [pdf] - 1151662
How fast do Jupiters grow? Signatures of the snowline and growth rate in the distribution of gas giant planets
Comments: 7 pages, 6 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society
Submitted: 2012-09-26
We present here observational evidence that the snowline plays a significant role in the formation and evolution of gas giant planets. When considering the population of observed exoplanets, we find a boundary in mass-semimajor axis space that suggests planets are preferentially found beyond the snowline prior to undergoing gap-opening inward migration and associated gas accretion. This is consistent with theoretical models suggesting that sudden changes in opacity -- as would occur at the snowline -- can influence core migration. Furthermore, population synthesis modelling suggests that this boundary implies that gas giant planets accrete ~ 70 % of the inward flowing gas, allowing ~ 30$ % through to the inner disc. This is qualitatively consistent with observations of transition discs suggesting the presence of inner holes, despite there being ongoing gas accretion.
[57]  [pdf] - 1084969
A Lower Angular Momentum Limit for Self-Gravitating Protostellar Disc Fragmentation
Comments: 10 pages, 19 figures, 1 table, accepted for publication in MNRAS
Submitted: 2011-10-19
We attempt to verify recent claims (made using semi-analytic models) that for the collapse of spherical homogeneous molecular clouds, fragmentation of the self-gravitating disc that subsequently forms can be predicted using the cloud's initial angular momentum alone. In effect, this condition is equivalent to requiring the resulting disc be sufficiently extended in order to fragment, in line with studies of isolated discs. We use smoothed particle hydrodynamics with hybrid radiative transfer to investigate this claim, confirming that in general, homogeneous spherical molecular clouds will produce fragmenting self-gravitating discs if the ratio of rotational kinetic energy to gravitational potential energy is greater than ~ 5e-3, where this result is relatively insensitive to the initial thermal energy. This condition begins to fail at higher cloud masses, suggesting that sufficient mass at large radii governs fragmentation. While these results are based on highly idealised initial conditions, and may not hold if the disc's accretion from the surrounding envelope is sufficiently asymmetric, or if the density structure is perturbed, they provide a sensible lower limit for the minimum angular momentum required to fragment a disc in the absence of significant external turbulence.
[58]  [pdf] - 1077763
The Jeans Mass as a Fundamental Measure of Self-Gravitating Disc Fragmentation and Initial Fragment Mass
Comments: 12 pages, 4 figures, accepted for publication in MNRAS
Submitted: 2011-07-05
As a formation route for objects such as giant planets and low-mass stars in protostellar discs (as well as stars in AGN discs), theories of self-gravitating disc fragmentation need to be able to predict the initial masses of fragments. We describe a means by which the local Jeans mass inside the spiral structure of a self-gravitating disc can be estimated. If such a self-gravitating disc satisfies the criteria for disc fragmentation, this estimate provides a lower limit for the initial mass of any fragments formed. We apply this approach to a series of self-gravitating protostellar disc models, to map out the typical masses of fragments produced by this formation mode. We find a minimum fragment mass of around 3 Jupiter masses, which is insensitive to the stellar mass, and that - within the parameter space surveyed - fragments with masses between 10 and 20 Jupiter masses are the most common. We also describe how the Jeans mass allows us to derive a more general criterion for disc fragmentation, which accounts for the processes of viscous heating, radiative cooling, accretion and the disc's thermal history. We demonstrate how such a criterion can be determined, and show that in limiting cases it recovers several fragmentation criteria that have been posited in the past, including the minimum cooling time/maximum stress criterion.
[59]  [pdf] - 1034184
The Nature of Angular Momentum Transport in Radiative Self-Gravitating Protostellar Discs
Comments: 13 pages, 44 figures, accepted for publication in MNRAS
Submitted: 2010-08-09
Semi-analytic models of self-gravitating discs often approximate the angular momentum transport generated by the gravitational instability using the phenomenology of viscosity. This allows the employment of the standard viscous evolution equations, and gives promising results. It is, however, still not clear when such an approximation is appropriate. This paper tests this approximation using high resolution 3D smoothed particle hydrodynamics (SPH) simulations of self-gravitating protostellar discs with radiative transfer. The nature of angular momentum transport associated with the gravitational instability is characterised as a function of both the stellar mass and the disc-to-star mass ratio. The effective viscosity is calculated from the Reynolds and gravitational stresses in the disc. This is then compared to what would be expected if the effective viscosity were determined by assuming local thermodynamic equilibrium or, equivalently, that the local dissipation rate matches the local cooling rate. In general, all the discs considered here settle into a self-regulated state where the heating generated by the gravitational instability is modulated by the local radiative cooling. It is found that low-mass discs can indeed be represented by a local "alpha-parametrisation", provided that the disc aspect ratio is small (H/R < 0.1) which is generally the case when the disc-to-star mass ratio q <0.5. However, this result does not extend to discs with masses approaching that of the central object. These are subject to transient burst events and global wave transport, and the effective viscosity is not well modelled by assuming local thermodynamic equilibrium. In spite of these effects, it is shown that massive (compact) discs can remain stable and not fragment, evolving rapidly to reduce their disc-to-star mass ratios through stellar accretion and radial spreading.
[60]  [pdf] - 1032656
Herschel-PACS observation of the 10 Myr old T Tauri disk TW Hya: Constraining the disk gas mass
Comments: 6 pages, 6 figures, accepted for publication in the A&A Herschel special issue
Submitted: 2010-05-20
Planets are formed in disks around young stars. With an age of ~10 Myr, TW Hya is one of the nearest T Tauri stars that is still surrounded by a relatively massive disk. In addition a large number of molecules has been found in the TW Hya disk, making TW Hya the perfect test case in a large survey of disks with Herschel-PACS to directly study their gaseous component. We aim to constrain the gas and dust mass of the circumstellar disk around TW Hya. We observed the fine-structure lines of [OI] and [CII] as part of the Open-time large program GASPS. We complement this with continuum data and ground-based 12CO 3-2 and 13CO 3-2 observations. We simultaneously model the continuum and the line fluxes with the 3D Monte-Carlo code MCFOST and the thermo-chemical code ProDiMo to derive the gas and dust masses. We detect the [OI] line at 63 micron. The other lines that were observed, [OI] at 145 micron and [CII] at 157 micron, are not detected. No extended emission has been found. Preliminary modeling of the photometric and line data assuming [12CO]/[13CO]=69 suggests a dust mass for grains with radius < 1 mm of ~1.9 times 10^-4 Msun (total solid mass of 3 times 10^-3 Msun) and a gas mass of (0.5--5) times 10^-3 Msun. The gas-to-dust mass may be lower than the standard interstellar value of 100.
[61]  [pdf] - 1032652
The Herschel view of Gas in Protoplanetary Systems (GASPS). First comparisons with a large grid of models
Comments: 5 pages, 4 figures, accepted for publication in the A&A Herschel special issue
Submitted: 2010-05-20
The Herschel GASPS Key Program is a survey of the gas phase of protoplanetary discs, targeting 240 objects which cover a large range of ages, spectral types, and disc properties. To interpret this large quantity of data and initiate self-consistent analyses of the gas and dust properties of protoplanetary discs, we have combined the capabilities of the radiative transfer code MCFOST with the gas thermal balance and chemistry code ProDiMo to compute a grid of 300 000 disc models (DENT). We present a comparison of the first Herschel/GASPS line and continuum data with the predictions from the DENT grid of models. Our objective is to test some of the main trends already identified in the DENT grid, as well as to define better empirical diagnostics to estimate the total gas mass of protoplanetary discs. Photospheric UV radiation appears to be the dominant gas-heating mechanism for Herbig stars, whereas UV excess and/or X-rays emission dominates for T Tauri stars. The DENT grid reveals the complexity in the analysis of far-IR lines and the difficulty to invert these observations into physical quantities. The combination of Herschel line observations with continuum data and/or with rotational lines in the (sub-)millimetre regime, in particular CO lines, is required for a detailed characterisation of the physical and chemical properties of circumstellar discs.
[62]  [pdf] - 1032654
Gas in the protoplanetary disc of HD 169142: Herschel's view
Comments: 5 pages, 3 Postscript figures, accepted for publication in the A&A Herschel special issue
Submitted: 2010-05-20
In an effort to simultaneously study the gas and dust components of the disc surrounding the young Herbig Ae star HD 169142, we present far-IR observations obtained with the PACS instrument onboard the Herschel Space Observatory. This work is part of the Open Time Key Project GASPS, which is aimed at studying the evolution of protoplanetary discs. To constrain the gas properties in the outer disc, we observed the star at several key gas-lines, including [OI] 63.2 and 145.5 micron, [CII] 157.7 micron, CO 72.8 and 90.2 micron, and o-H2O 78.7 and 179.5 micron. We only detect the [OI] 63.2 micron line in our spectra, and derive upper limits for the other lines. We complement our data set with PACS photometry and 12/13CO data obtained with the Submillimeter Array. Furthermore, we derive accurate stellar parameters from optical spectra and UV to mm photometry. We model the dust continuum with the 3D radiative transfer code MCFOST and use this model as an input to analyse the gas lines with the thermo-chemical code ProDiMo. Our dataset is consistent with a simple model in which the gas and dust are well-mixed in a disc with a continuous structure between 20 and 200 AU, but this is not a unique solution. Our modelling effort allows us to constrain the gas-to-dust mass ratio as well as the relative abundance of the PAHs in the disc by simultaneously fitting the lines of several species that originate in different regions. Our results are inconsistent with a gas-poor disc with a large UV excess; a gas mass of 5.0 +/- 2.0 times 10^(-3) Msun is still present in this disc, in agreement with earlier CO observations.
[63]  [pdf] - 1032661
Gas in Protoplanetary Discs (GASPS) 1. First results
Comments: Comments: 5 pages, 2 Postscript figures, accepted for publication in the A&A Herschel special issue
Submitted: 2010-05-20
Context - Circumstellar discs are ubiquitous around young stars, but rapidly dissipate their gas and dust on timescales of a few Myr. The Herschel space observatory allows for the study of the warm disc atmosphere, using far-infrared spectroscopy to measure gas content and excitation conditions, and far-IR photometry to constrain the dust distribution. Aims - We aim to detect and characterize the gas content of circumstellar discs in four targets as part of the Herschel science demonstration phase. Methods - We carried out sensitive medium resolution spectroscopy and high sensitivity photometry at lambda ~60-190 micron using the Photodetector Array Camera and Spectrometer instrument on the Herschel space observatory. Results - We detect [OI] 63 micron emission from the young stars HD 169142, TW Hydrae, and RECX 15, but not HD 181327. No other lines, including [CII] 158 and [OI] 145, are significantly detected. All four stars are detected in photometry at 70 and 160 micron. Extensive models are presented in associated papers.
[64]  [pdf] - 1026360
Native Synthetic Imaging of Smoothed Particle Hydrodynamics density fields using gridless Monte Carlo Radiative Transfer
Comments: 10 pages, 10 figures, accepted for publication in MNRAS
Submitted: 2010-04-16
An algorithm for creating synthetic telescope images of Smoothed Particle Hydrodynamics (SPH) density fields is presented, which utilises the adaptive nature of the SPH formalism in full. The imaging process uses Monte Carlo Radiative Transfer (MCRT) methods to model the scattering and absorption of photon packets in the density field, which then exit the system and are captured on a pixelated image plane, creating a 2D image (or a 3D datacube, if the photons are also binned by their wavelength). The algorithm is implemented on the density field directly: no gridding of the field is required, allowing the density field to be described to an identical level of accuracy as the simulations that generated it. Some applications of the method to star and planet formation simulations are presented to illustrate the advantages of this new technique, and suggestions as to how this framework could support a Radiative Equilibrium algorithm are also given as an indication for future development.
[65]  [pdf] - 149551
Numerical Testing of The Rare Earth Hypothesis using Monte Carlo Realisation Techniques
Comments: 13 pages, 10 figures, accepted for publication in the International Journal of Astrobiology
Submitted: 2010-01-11
The Search for Extraterrestrial Intelligence (SETI) has thus far failed to provide a convincing detection of intelligent life. In the wake of this null signal, many "contact pessimistic" hypotheses have been formulated, the most famous of which is the Rare Earth Hypothesis. It postulates that although terrestrial planets may be common, the exact environmental conditions that Earth enjoys are rare, perhaps unique. As a result, simple microbial life may be common, but complex metazoans (and hence intelligence) will be rare. This paper uses Monte Carlo Realisation Techniques to investigate the Rare Earth Hypothesis, in particular the environmental criteria considered imperative to the existence of intelligence on Earth. By comparing with a less restrictive, more optimistic hypothesis, the data indicates that if the Rare Earth hypothesis is correct, intelligent civilisation will indeed be relatively rare. Studying the separations of pairs of civilisations shows that most intelligent civilisation pairs (ICPs) are unconnected: that is, they will not be able to exchange signals at lightspeed in the limited time that both are extant. However, the few ICPs that are connected are strongly connected, being able to participate in numerous exchanges of signals. This may provide encouragement for SETI researchers: although the Rare Earth Hypothesis is in general a contact-pessimistic hypothesis, it may be a "soft" or "exclusive" hypothesis, i.e. it may contain facets that are latently contact-optimistic.
[66]  [pdf] - 1018193
Stellar Encounters in the Context of Outburst Phenomena
Comments: 8 pages, 11 figures, 1 table, accepted for publication in MNRAS
Submitted: 2009-11-03
Young stellar systems are known to undergo outbursts, where the star experiences an increased accretion rate, and the system's luminosity increases accordingly. The archetype is the FU Orionis (FU Ori) outburst, where the accretion rate can increase by three orders of magnitude (and the brightness of the system by five magnitudes). The cause appears to be instability in the circumstellar disc, but there is currently some debate as to the nature of this instability (e.g. thermal, gravitational, magneto-rotational). This paper details high resolution Smoothed Particle Hydrodynamics (SPH) simulations that were carried out to investigate the influence of stellar encounters on disc dynamics. Star-star encounters (where the primary has a self-gravitating, marginally stable protostellar disc) were simulated with various orbital parameters to investigate the resulting disc structure and dynamics. Crucially, the simulations include the effects of radiative transfer to realistically model the resulting thermodynamics. Our results show that the accretion history and luminosity of the system during the encounter displays many of the features of outburst phenomena. In particular, the magnitudes and decay times seen are comparable to those of FU Ori. There are two caveats to this assertion: the first is that these events are not expected to occur frequently enough to explain all FU Ori or EX Lupi; the second is that the inner discs of these simulations are subject to numerical viscosity, which will act to reduce the accretion rate (although it has less of an effect on the total mass accreted). In short, these results cannot rule out binary interactions as a potential source of some FU Ori-esque outbursts.
[67]  [pdf] - 1017179
Stellar Encounters: A Stimulus for Disc Fragmentation?
Comments: 16 figures, 2 tables, accepted for publication in MNRAS
Submitted: 2009-08-23
An interaction between a star-disc system and another star will perturb the disc, possibly resulting in a significant modification of the disc structure and its properties. It is still unclear if such an encounter can trigger fragmentation of the disc to form brown dwarfs or gas giant planets. This paper details high resolution Smoothed Particle Hydrodynamics (SPH) simulations investigating the influence of stellar encounters on disc dynamics. Star-star encounters (where the primary has a self-gravitating, marginally stable protostellar disc, and the secondary has no disc) were simulated with various orbital parameters to investigate the resulting disc structure and dynamics. This work is the first of its kind to incorporate realistic radiative transfer techniques to realistically model the resulting thermodynamics. The results suggest that the effect of stellar encounters is to prohibit fragmentation - compressive and shock heating stabilises the disc, and the radiative cooling is insufficient to trigger gravitational instability. The encounter strips the outer regions of the disc (either through tidal tails or by capture of matter to form a disc around the secondary), which triggers a readjustment of the primary disc to a steeper surface density profile (and a flatter Toomre Q profile). The disc around the secondary plays a role in the potential capture of the secondary to form a binary. However, this applies only to orbits that are parabolic - hyperbolic encounters do not form a secondary disc, and are not captured.
[68]  [pdf] - 19044
Introducing a Hybrid Method of Radiative Transfer for Smoothed Particle Hydrodynamics
Comments: 11 pages, 22 figures, accepted for publication in MNRAS
Submitted: 2008-12-01
A new means of incorporating radiative transfer into smoothed particle hydrodynamics (SPH) is introduced, which builds on the success of two previous methods - the polytropic cooling approximation as devised by Stamatellos et al (2007), and flux limited diffusion (e.g. Mayer et al 2007). This hybrid method preserves the strengths of its individual components, while removing the need for atmosphere matching or other boundary conditions to marry optically thick and optically thin regions. The code uses a non-trivial equation of state to calculate temperatures and opacities of SPH particles, which captures the effects of molecular hydrogen dissociation, atomic hydrogen ionisation, He0 and He+ ionisation, ice evaporation, dust sublimation, molecular absorption, bound-free and free-free transitions and electron scattering. The method is tested in several scenarios, including: (1) the evolution of a 0.07 solar mass protoplanetary disc surrounding a 0.5 solar mass star; (2) the collapse of a 1 solar mass protostellar cloud, and (3) the thermal relaxation of temperature fluctuations in a static homogeneous sphere.
[69]  [pdf] - 16790
Introducing a Hybrid Method of Radiative Transfer in Smoothed Particle Hydrodynamics
Comments: 4 pages, 6 figures, to appear in the proceedings of the Cool Stars 15 conference
Submitted: 2008-09-26
We present a new method of incorporating radiative transfer into Smoothed Particle Hydrodynamics (SPH). There have been many recent attempts at radiative transfer in SPH (Stamatellos et al 2005, 2005, Mayer et al 2007, Whitehouse and Bate 2006), however these are becoming increasingly complex, with some methods requiring the photosphere to be mapped (which is often of non-trivial geometric shape), and extra conditions to be applied there (matching atmospheres as in Cai et al (2008), or specifying cooling at the photosphere as in Mayer et al (2007)). The method of identifying the photosphere is usually a significant addition to the total simulation runtime, and often requires extra free parameters, the changing of which will affect the final results. Our method is not affected by such concerns, as the photosphere is constructed implicitly by the algorithm without the need for extra free parameters. The algorithm used is a synergy of two current formalisms for radiative effects: a) the polytropic cooling formalism proposed by Stamatellos et al (2007), and b) flux-limited diffusion, used by many authors to simulate radiation transport in the optically thick regime (e.g. Mayer et al 2007). We present several tests of this method: (1) The evolution of a 0.07 solar mass protoplanetary disc around a 0.5 solarmass star (Pickett et al 2003, Mejia et al 2005, Boley et al 2006, Cai et al 2008); (2) The collapse of a non-rotating 1 solar mass molecular cloud (Masunaga & Inutsuka 2000, Stamatellos et al 2007); (3) The thermal relaxation of temperature fluctuations in an static homogeneous sphere (Masunaga et al 1998, Spiegel 1957, Stamatellos et al 2007)
[70]  [pdf] - 16674
Gas Evolution in Protoplanetary Disks
Comments: 9 pages, 6 figures
Submitted: 2008-09-24
This article summarizes a Splinter Session at the Cool Stars XV conference in St. Andrews with 3 review and 4 contributed talks. The speakers have discussed various approaches to understand the structure and evolution of the gas component in protoplanetary disks. These ranged from observational spectroscopy in the UV, infrared and millimeter, through to chemical and hydrodynamical models. The focus was on disks around low-mass stars, ranging from classical T Tauri stars to transitional disks and debris disks. Emphasis was put on water and organic molecules, the relation to planet formation, and the formation of holes and gaps in the inner regions.
[71]  [pdf] - 81925
A comparative study of disc-planet interaction
Comments: 32 pages, accepted for publication in MNRAS
Submitted: 2006-05-09, last modified: 2006-06-01
We perform numerical simulations of a disc-planet system using various grid-based and smoothed particle hydrodynamics (SPH) codes. The tests are run for a simple setup where Jupiter and Neptune mass planets on a circular orbit open a gap in a protoplanetary disc during a few hundred orbital periods. We compare the surface density contours, potential vorticity and smoothed radial profiles at several times. The disc mass and gravitational torque time evolution are analyzed with high temporal resolution. There is overall consistency between the codes. The density profiles agree within about 5% for the Eulerian simulations while the SPH results predict the correct shape of the gap although have less resolution in the low density regions and weaker planetary wakes. The disc masses after 200 orbital periods agree within 10%. The spread is larger in the tidal torques acting on the planet which agree within a factor 2 at the end of the simulation. In the Neptune case the dispersion in the torques is greater than for Jupiter, possibly owing to the contribution from the not completely cleared region close to the planet.
[72]  [pdf] - 80388
Gravitational Instabilities in Gaseous Protoplanetary Disks and Implications for Giant Planet Formation
Comments: 17 pages, 7 figures, Invited Review Chapter to appear on Protostars and Planets V. A version with high resolution figures can be found at
Submitted: 2006-03-08
Protoplanetary gas disks are likely to experience gravitational instabilites (GI's) during some phase of their evolution. Density perturbations in an unstable disk grow on a dynamic time scale into spiral arms that produce efficient outward transfer of angular momentum and inward transfer of mass through gravitational torques. In a cool disk with rapid enough cooling, the spiral arms in an unstable disk form self-gravitating clumps. Whether gas giant protoplanets can form by such a disk instability process is the primary question addressed by this review. We discuss the wide range of calculations undertaken by ourselves and others using various numerical techniques, and we report preliminary results from a large multi-code collaboration. Additional topics include -- triggering mechanisms for GI's, disk heating and cooling, orbital survival of dense clumps, interactions of solids with GI-driven waves and shocks, and hybrid scenarios where GI's facilitate core accretion. The review ends with a discussion of how well disk instability and core accretion fare in meeting observational constraints.