BZPEERPapers, people, evidence and research paths
Menu

Discover

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

Research tools

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

Network

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

Opportunities

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

Account

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

Source author record

Ken Rice

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

ResearcherUnclaimed source record
astro-ph.EPastro-ph.SRastro-ph.GAastro-ph.IMMachine Learning

Catalog footprint

What is connected

40works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

40 published item(s)

preprint2022arXiv

Binary companions triggering fragmentation in self-gravitating discs

Observations of systems hosting close in ($<1$ AU) giant planets and brown dwarfs ($M\gtrsim7$ M$_{\rm Jup}$) find an excess of binary star companions, indicating that stellar multiplicity may play an important role in their formation. There is now increasing evidence that some of these objects may have formed via fragmentation in gravitationally unstable discs. We present a suite of 3D smoothed particle hydrodynamics (SPH) simulations of binary star systems with circumprimary self-gravitating discs, which include a realistic approximation to radiation transport, and extensively explore the companion's orbital parameter space for configurations which may trigger fragmentation. We identify a "sweet spot" where intermediate separation binary companions ($100$ AU $\lesssim a\lesssim400$ AU) can cause a marginally stable disc to fragment. The exact range of ideal binary separations is a function of the companion's eccentricity, inclination and mass. Heating is balanced by efficient cooling, and fragmentation occurs inside a spiral mode driven by the companion. Short separation, disc penetrating binary encounters ($a\lesssim100$ AU) are prohibitive to fragmentation, as mass stripping and disc heating quench any instability. This is also true of binary companions with high orbital eccentricities ($e\gtrsim0.75$). Wide separation companions ($a\gtrsim500$ AU) have little effect on the disc properties for the setup parameters considered here. The sweet spot found is consistent with the range of binary separations which display an excess of close in giant planets and brown dwarfs. Hence we suggest that fragmentation triggered by a binary companion may contribute to the formation of these substellar objects.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Identifying Exoplanets with Deep Learning. IV. Removing Stellar Activity Signals from Radial Velocity Measurements Using Neural Networks

Exoplanet detection with precise radial velocity (RV) observations is currently limited by spurious RV signals introduced by stellar activity. We show that machine learning techniques such as linear regression and neural networks can effectively remove the activity signals (due to starspots/faculae) from RV observations. Previous efforts focused on carefully filtering out activity signals in time using modeling techniques like Gaussian Process regression (e.g. Haywood et al. 2014). Instead, we systematically remove activity signals using only changes to the average shape of spectral lines, and no information about when the observations were collected. We trained our machine learning models on both simulated data (generated with the SOAP 2.0 software; Dumusque et al. 2014) and observations of the Sun from the HARPS-N Solar Telescope (Dumusque et al. 2015; Phillips et al. 2016; Collier Cameron et al. 2019). We find that these techniques can predict and remove stellar activity from both simulated data (improving RV scatter from 82 cm/s to 3 cm/s) and from more than 600 real observations taken nearly daily over three years with the HARPS-N Solar Telescope (improving the RV scatter from 1.753 m/s to 1.039 m/s, a factor of ~ 1.7 improvement). In the future, these or similar techniques could remove activity signals from observations of stars outside our solar system and eventually help detect habitable-zone Earth-mass exoplanets around Sun-like stars.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Independent validation of the temperate Super-Earth HD79211 b using HARPS-N

We present high-precision radial velocities (RVs) from the HARPS-N spectrograph for HD79210 and HD79211, two M0V members of a gravitationally-bound binary system. We detect a planet candidate with a period of $24.421^{+0.016}_{-0.017}$ days around HD79211 in these HARPS-N RVs, validating the planet candidate originally identified in CARMENES RV data alone. Using HARPS-N, CARMENES and HIRES RVs spanning a total of 25 years, we further refine the planet candidate parameters to $P=24.422\pm0.014$ days, $K=3.19\pm0.27$ m/s, $M$ sin $i = 10.6 \pm 1.2 M_\oplus$, and $a = 0.142 \pm0.005$ au. We do not find any additional planet candidate signals in the data of HD79211 nor do we find any planet candidate signals in HD79210. This system adds to the number of exoplanets detected in binaries with M dwarf members, and serves as a case study for planet formation in stellar binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

K2-79b and K2-222b: Mass measurements of two small exoplanets with periods beyond 10 days that overlap with periodic magnetic activity signals

We present mass and radius measurements of K2-79b and K2-222b, two transiting exoplanets orbiting active G-type stars. Their respective 10.99d and 15.39d orbital periods fall near periods of signals induced by stellar magnetic activity. The two signals might therefore interfere and lead to an inaccurate estimate of exoplanet mass. We present a method to mitigate these effects when radial velocity and activity indicator observations are available over multiple observing seasons and the orbital period of the exoplanet is known. We perform correlation and periodogram analyses on sub-sets composed of each target's two observing seasons, in addition to the full data sets. For both targets, these analyses reveal an optimal season with little to no interference at the orbital period of the known exoplanet. We make a confident mass detection of each exoplanet by confirming agreement between fits to the full radial velocity set and the optimal season. For K2-79b, we measure a mass of 11.8 $\pm$ 3.6 $M_{Earth}$ and a radius of 4.09 $\pm$ 0.17 $R_{Earth}$. For K2-222b, we measure a mass of 8.0 $\pm$ 1.8 $M_{Earth}$ and a radius of 2.35 $\pm$ 0.08 $R_{Earth}$. According to model predictions, K2-79b is a highly irradiated Uranus-analog and K2-222b hosts significant amounts of water ice. We also present an RV solution for a candidate second companion orbiting K2-222 at 147.5d.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A pair of TESS planets spanning the radius valley around the nearby mid-M dwarf LTT 3780

We present the confirmation of two new planets transiting the nearby mid-M dwarf LTT 3780 (TIC 36724087, TOI-732, $V=13.07$, $K_s=8.204$, $R_s$=0.374 R$_{\odot}$, $M_s$=0.401 M$_{\odot}$, d=22 pc). The two planet candidates are identified in a single TESS sector and are validated with reconnaissance spectroscopy, ground-based photometric follow-up, and high-resolution imaging. With measured orbital periods of $P_b=0.77$ days, $P_c=12.25$ days and sizes $r_{p,b}=1.33\pm 0.07$ R$_{\oplus}$, $r_{p,c}=2.30\pm 0.16$ R$_{\oplus}$, the two planets span the radius valley in period-radius space around low mass stars thus making the system a laboratory to test competing theories of the emergence of the radius valley in that stellar mass regime. By combining 63 precise radial-velocity measurements from HARPS and HARPS-N, we measure planet masses of $m_{p,b}=2.62^{+0.48}_{-0.46}$ M$_{\oplus}$ and $m_{p,c}=8.6^{+1.6}_{-1.3}$ M$_{\oplus}$, which indicates that LTT 3780b has a bulk composition consistent with being Earth-like, while LTT 3780c likely hosts an extended H/He envelope. We show that the recovered planetary masses are consistent with predictions from both photoevaporation and from core-powered mass loss models. The brightness and small size of LTT 3780, along with the measured planetary parameters, render LTT 3780b and c as accessible targets for atmospheric characterization of planets within the same planetary system and spanning the radius valley.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fragmentation favoured in discs around higher mass stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Massive discs around low-mass stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The observational impact of dust trapping in self-gravitating discs

We present a 3D semi-analytic model of self-gravitating discs, and include a prescription for dust trapping in the disc spiral arms. Using Monte-Carlo radiative transfer we produce synthetic ALMA observations of these discs. In doing so we demonstrate that our model is capable of producing observational predictions, and able to model real image data of potentially self-gravitating discs. For a disc to generate spiral structure that would be observable with ALMA requires that the disc's dust mass budget is dominated by millimetre and centimetre-sized grains. Discs in which grains have grown to the grain fragmentation threshold may satisfy this criterion, thus we predict that signatures of gravitational instability may be detectable in discs of lower mass than has previously been suggested. For example, we find that discs with disc-to-star mass ratios as low as $0.10$ are capable of driving observable spiral arms. Substructure becomes challenging to detect in discs where no grain growth has occurred or in which grain growth has proceeded well beyond the grain fragmentation threshold. We demonstrate how we can use our model to retrieve information about dust trapping and grain growth through multi-wavelength observations of discs, and using estimates of the opacity spectral index. Applying our disc model to the Elias 27, WaOph 6 and IM Lup systems we find gravitational instability to be a plausible explanation for the observed substructure in all 3 discs, if sufficient grain growth has indeed occurred.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

TOI-1235 b: a keystone super-Earth for testing radius valley emergence models around early M dwarfs

Small planets on close-in orbits tend to exhibit envelope mass fractions of either effectively zero or up to a few percent depending on their size and orbital period. Models of thermally-driven atmospheric mass loss and of terrestrial planet formation in a gas-poor environment make distinct predictions regarding the location of this rocky/non-rocky transition in period-radius space. Here we present the confirmation of TOI-1235 b ($P=3.44$ days, $r_p=1.738^{+0.087}_{-0.076}$ R$_{\oplus}$), a planet whose size and period are intermediate between the competing model predictions thus making the system an important test case for emergence models of the rocky/non-rocky transition around early M dwarfs ($R_s=0.630\pm 0.015$ R$_{\odot}$, $M_s=0.640\pm 0.016$ M$_{\odot}$). We confirm the TESS planet discovery using reconnaissance spectroscopy, ground-based photometry, high-resolution imaging, and a set of 38 precise radial-velocities from HARPS-N and HIRES. We measure a planet mass of $6.91^{+0.75}_{-0.85}$ M$_{\oplus}$, which implies an iron core mass fraction of $20^{+15}_{-12}$% in the absence of a gaseous envelope. The bulk composition of TOI-1235 b is therefore consistent with being Earth-like and we constrain a H/He envelope mass fraction to be $<0.5$% at 90% confidence. Our results are consistent with model predictions from thermally-driven atmospheric mass loss but not with gas-poor formation, suggesting that the former class of processes remain efficient at sculpting close-in planets around early M dwarfs. Our RV analysis also reveals a strong periodicity close to the first harmonic of the photometrically-determined stellar rotation period that we treat as stellar activity, despite other lines of evidence favoring a planetary origin ($P=21.8^{+0.9}_{-0.8}$ days, $m_p\sin{i}=13.0^{+3.8}_{-5.3}$ M$_{\oplus}$) that cannot be firmly ruled out by our data.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A 1.9 Earth radius rocky planet and the discovery of a non-transiting planet in the Kepler-20 system

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).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Directly observing continuum emission from self-gravitating spiral waves

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Grand challenges in protoplanetary disc modelling

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Kepler-21b: A rocky planet around a V = 8.25 magnitude star

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Tensor classification of structure in smoothed particle hydrodynamics density fields

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The evolution of self-gravitating accretion discs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Orbit and Mass of the Third Planet in the Kepler-56 System

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Can Kozai-Lidov cycles explain Kepler-78b?

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Disc fragmentation rarely forms planetary-mass objects

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Spiral Arms in Gravitationally Unstable Protoplanetary Disks as Imaged in Scattered Light

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Gemini NICI Planet-Finding Campaign: Asymmetries in the HD 141569 disc

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Kepler-454 System: A Small, Not-rocky Inner Planet, a Jovian World, and a Distant Companion

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The quest for cradles of life: using the fundamental metallicity relation to hunt for the most habitable type of galaxy

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Characterizing K2 Planet Discoveries: A super-Earth transiting the bright K-dwarf HIP 116454

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Forming misaligned stellar discs around a massive black hole: Cloud infall in the Galactic Centre

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Dynamical Fate of Self-Gravitating Disc Fragments After Tidal Downsizing

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Kepler-10 planetary system revisited by HARPS-N: A hot rocky world and a solid Neptune-mass planet

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Mass of Kepler-93b and The Composition of Terrestrial Planets

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Self-Gravitating Disc Around L1527 IRS?

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An Earth-sized planet with an Earth-like density

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Misaligned streamers around a galactic centre black hole from a single cloud's infall

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Effect of Irradiation on the Jeans Mass in Fragmenting Self-Gravitating Protostellar Discs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Towards a Population Synthesis Model of Objects formed by Self-Gravitating Disc Fragmentation and Tidal Downsizing

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)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

How fast do Jupiters grow? Signatures of the snowline and growth rate in the distribution of gas giant planets

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Lower Angular Momentum Limit for Self-Gravitating Protostellar Disc Fragmentation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Jeans Mass as a Fundamental Measure of Self-Gravitating Disc Fragmentation and Initial Fragment Mass

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Native Synthetic Imaging of Smoothed Particle Hydrodynamics density fields using gridless Monte Carlo Radiative Transfer

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerical Testing of The Rare Earth Hypothesis using Monte Carlo Realisation Techniques

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Nature of Angular Momentum Transport in Radiative Self-Gravitating Protostellar Discs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Stellar Encounters in the Context of Outburst Phenomena

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Stellar Encounters: A Stimulus for Disc Fragmentation?

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.

0 citations0 reviewsNo code linkedOpen access