BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

W. K. M. Rice

W. K. M. Rice contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.EPastro-ph.SRastro-ph.GAastro-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

12 published item(s)

preprint2014arXiv

Convergence of simulations of self-gravitating accretion discs II: Sensitivity to the implementation of radiative cooling and artificial viscosity

Recently it has been suggested that the fragmentation boundary in Smoothed Particle Hydrodynamic (SPH) and FARGO simulations of self-gravitating accretion discs with beta-cooling do not converge as resolution is increased. Furthermore, this recent work suggests that by carefully optimising the artificial viscosity parameters in these codes it can be shown that fragmentation may occur for much longer cooling times than earlier work suggests. If correct, this result is intriguing as it suggests that gas giant planets could form, via direct gravitational collapse, reasonably close to their parent stars. This result is, however, slightly surprising and there have been a number of recent studies suggesting that the result is likely an indication of a numerical problem with the simulations. One suggestion, in particular, is that the SPH results are influenced by the manner in which the cooling is implemented. We extend this work here and show that if the cooling is implemented in a manner that removes a known numerical artefact in the shock regions, the fragmentation boundary converges to a value consistent with earlier work and that fragmentation is unlikely for the long cooling times suggested by this recent work. We also investigate the optimisation of the artificial viscosity parameters and show that the values that appear optimal are likely introducing numerical problems in both the SPH and FARGO simulations. We therefore conclude that earlier predictions for the cooling times required for fragmentation are likely correct and that, as suggested by this earlier work, fragmentation cannot occur in the inner parts (r < 50 au) of typical protostellar discs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Planetesimal formation in self-gravitating discs: the effects of particle self-gravity and back-reaction

We study particle dynamics in self-gravitating gaseous discs with a simple cooling law prescription via two-dimensional simulations in the shearing sheet approximation. It is well known that structures arising in the gaseous component of the disc due to a gravitational instability can have a significant effect on the evolution of dust particles. Previous results have shown that spiral density waves can be highly efficient at collecting dust particles, creating significant local over-densities of particles. The degree of such concentrations has been shown to be dependent on two parameters: the size of the dust particles and the rate of gas cooling. We expand on these findings, including the self-gravity of dust particles, to see how these particle over-densities evolve. We use the PENCIL CODE to solve the local shearing sheet equations for gas on a fixed grid together with the equations of motion for solids coupled to the gas through an aerodynamic drag force. We find that the enhancements in the surface density of particles in spiral density wave crests can reach levels high enough to allow the solid component of the disc to collapse under its own self-gravity. This produces many gravitationally bound collections of particles within the spiral structure. The total mass contained in bound structures appears nearly independent of the cooling time, suggesting that the formation of planetesimals through dust particle trapping by self-gravitating density waves may be possible at a larger range of radii within a disc than previously thought. So, density waves due to gravitational instabilities in the early stages of star formation may provide excellent sites for the rapid formation of many large, planetesimal-sized objects.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Planetesimal Formation In Self-Gravitating Discs

We study particle dynamics in local two-dimensional simulations of self-gravitating accretion discs with a simple cooling law. It is well known that the structure which arises in the gaseous component of the disc due to a gravitational instability can have a significant effect on the evolution of dust particles. Previous results using global simulations indicate that spiral density waves are highly efficient at collecting dust particles, creating significant local over-densities which may be able to undergo gravitational collapse. We expand on these findings, using a range of cooling times to mimic the conditions at a large range of radii within the disc. Here we use the Pencil Code to solve the 2D local shearing sheet equations for gas on a fixed grid together with the equations of motion for solids coupled to the gas solely through aerodynamic drag force. We find that spiral density waves can create significant enhancements in the surface density of solids, equivalent to 1-10cm sized particles in a disc following the profiles of Clarke (2009) around a solar mass star, causing it to reach concentrations several orders of magnitude larger than the particles mean surface density. We also study the velocity dispersion of the particles, finding that the spiral structure can result in the particle velocities becoming highly ordered, having a narrow velocity dispersion. This implies low relative velocities between particles, which in turn suggests that collisions are typically low energy, lessening the likelihood of grain destruction. Both these findings suggest that the density waves that arise due to gravitational instabilities in the early stages of star formation provide excellent sites for the formation of large, planetesimal-sized objects.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Tidal evolution of close-in giant planets : Evidence of Type II migration?

It is well accepted that &#39;hot Jupiters&#39; did not form in situ, as the temperature in the protoplanetary disc at the radius at which they now orbit would have been too high for planet formation to have occurred. These planets, instead, form at larger radii and then move into the region in which they now orbit. The exact process that leads to the formation of these close-in planets is, however, unclear and it seems that there may be more than one mechanism that can produce these short-period systems. Dynamical interactions in multiple-planet systems can scatter planets into highly eccentric orbits which, if the pericentre is sufficiently close to the parent star, can be tidally circularised by tidal interactions between the planet and star. Furthermore, systems with distant planetary or stellar companions can undergo Kozai cycles which can result in a planet orbiting very close to its parent star. However, the most developed model for the origin of short period planets is one in which the planet exchanges angular momentum with the surrounding protoplanetary disc and spirals in towards the central star. In the case of &#39;hot Jupiters&#39;, the planet is expected to open a gap in the disc and migrate through Type II .migration. If this is the dominant mechanism for producing `hot Jupiters&#39; then we would expect the currect properties of observed close-in giant planets to be consistent with an initial population resulting from Type II migration followed by evolution due to tidal interactions with the central star. We consider initial distributions that are consistent with Type II migration and find that after tidal evolution, the final distributions can be consistent with that observed. Our results suggest that a modest initial pile-up at a ~ 0.05 au is required and that the initial eccentricity distribution must peak at e \sim 0.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Convergence of SPH simulations of self-gravitating accretion discs: Sensitivity to the implementation of radiative cooling

Recent simulations of self-gravitating accretion discs, carried out using a three-dimensional Smoothed Particle Hydrodynamics (SPH) code by Meru and Bate, have been interpreted as implying that three-dimensional global discs fragment much more easily than would be expected from a two-dimensional local model. Subsequently, global and local two-dimensional models have been shown to display similar fragmentation properties, leaving it unclear whether the three-dimensional results reflect a physical effect or a numerical problem associated with the treatment of cooling or artificial viscosity in SPH. Here, we study how fragmentation of self-gravitating disc flows in SPH depends upon the implementation of cooling. We run disc simulations that compare a simple cooling scheme, in which each particle loses energy based upon its internal energy per unit mass, with a method in which the cooling is derived from a smoothed internal energy density field. For the simple per particle cooling scheme, we find a significant increase in the minimum cooling time scale for fragmentation with increasing resolution, matching previous results. Switching to smoothed cooling, however, results in lower critical cooling time scales, and tentative evidence for convergence at the highest spatial resolution tested. We conclude that precision studies of fragmentation using SPH require careful consideration of how cooling (and, probably, artificial viscosity) is implemented, and that the apparent non-convergence of the fragmentation boundary seen in prior simulations is likely a numerical effect. In real discs, where cooling is physically smoothed by radiative transfer effects, the fragmentation boundary is probably displaced from the two-dimensional value by a factor that is only of the order of unity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Excitation of spiral density waves by convection in accretion discs

Motivated by the recent results of \citet{Lesur_Ogilvie10} on the transport properties of incompressible convection in protoplanetary discs, in this paper we study the role of compressibility and hence of another basic mode -- spiral density waves -- in convective instability in discs. We analyse the linear dynamics of non-axisymmetric convection and spiral density waves in a Keplerian disc with superadiabatic vertical stratification using the local shearing box approach. It is demonstrated that the shear associated with Keplerian differential rotation introduces a novel phenomenon, it causes these two perturbation modes to become coupled: during evolution the convective mode generates (trailing) spiral density waves and can therefore be regarded as a new source of spiral density waves in discs. The wave generation process studied here owes its existence solely to shear of the disc&#39;s differential rotation, and is a special manifestation of a more general linear mode coupling phenomena universally taking place in flows with an inhomogeneous velocity profile. We quantify the efficiency of spiral density wave generation by convection as a function of azimuthal and vertical wavenumbers of these modes and find that it is maximal and most powerful when both these length-scales are comparable to the disc scale height. We also show that unlike the convective mode, which tends to transport angular momentum inwards in the linear regime, the spiral density waves transport angular momentum outwards. Based on these findings, we suggest that in the non-linear regime spiral density waves generated by convection may play a role in enhancing the transport of angular momentum due the convective mode alone, which is actually being changed to outward by non-linearity, as indicated by above-mentioned recent developments.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Stability of self-gravitating discs under irradiation

Self-gravity becomes competitive as an angular momentum transport process in accretion discs at large radii, where the temperature is low enough that external irradiation likely contributes to the thermal balance. Irradiation is known to weaken the strength of disc self-gravity, and can suppress it entirely if the disc is maintained above the threshold for linear instability. However, its impact on the susceptibility of the disc to fragmentation is less clear. We use two-dimensional numerical simulations to investigate the evolution of self-gravitating discs as a function of the local cooling time and strength of irradiation. In the regime where the disc does not fragment, we show that local thermal equilibrium continues to determine the stress - which can be represented as an effective viscous alpha - out to very long cooling times (at least 240 dynamical times). In this regime, the power spectrum of the perturbations is uniquely set by the effective viscous alpha and not by the cooling rate. Fragmentation occurs for cooling times tau < beta_crit / Omega, where beta_crit is a weak function of the level of irradiation. We find that beta_crit declines by approximately a factor of two, as irradiation is increased from zero up to the level where instability is almost quenched. The numerical results imply that irradiation cannot generally avert fragmentation of self-gravitating discs at large radii; if other angular momentum transport sources are weak mass will build up until self-gravity sets in, and fragmentation will ensue.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Axisymmetric modes in vertically stratified self-gravitating discs

We perform linear analysis of axisymmetric vertical normal modes in stratified compressible self-gravitating polytropic discs in the shearing box approximation. We study specific dynamics for subadiabatic, adiabatic and superadiabatic vertical stratifications. In the absence of self-gravity, four well-known principal modes can be identified in a stratified disc: acoustic p-, surface gravity f-, buoyancy g- and inertial r-modes. After characterizing modes in the non-self-gravitating case, we include self-gravity and investigate how it modifies the properties of these modes. We find that self-gravity, to a certain degree, reduces their frequencies and changes the structure of the dispersion curves and eigenfunctions at radial wavelengths comparable to the disc height. Its influence on the basic branch of the r-mode, in the case of subadiabatic and adiabatic stratifications, and on the basic branch of the g-mode, in the case of superadiabatic stratification (which in addition exhibits convective instability), does appear to be strongest. Reducing the three-dimensional Toomre&#39;s parameter Q_{3D} results in the latter modes becoming unstable due to self-gravity, so that they determine the onset criterion and nature of gravitational instability of a stratified disc. By contrast, the p-, f- and convectively stable g-modes, although their corresponding ω^2 are reduced by self-gravity, never become unstable however small the value of Q_{3D}. This is a consequence of the three-dimensionality of the disc. The eigenfunctions corresponding to the gravitationally unstable modes are intrinsically 3D. We also contrast the more exact instability criterion based on our 3D model with that of density waves in 2D (razor-thin) discs. Based on these findings, we comment on the origin of surface distortions seen in numerical simulations of self-gravitating discs.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Have proto-planetary discs formed planets?

It has recently been noted that many discs around T Tauri stars appear to comprise only a few Jupiter-masses of gas and dust. Using millimetre surveys of discs within six local star-formation regions, we confirm this result, and find that only a few percent of young stars have enough circumstellar material to build gas giant planets, in standard core accretion models. Since the frequency of observed exo-planets is greater than this, there is a `missing mass&#39; problem. As alternatives to simply adjusting the conversion of dust-flux to disc mass, we investigate three other classes of solution. Migration of planets could hypothetically sweep up the disc mass reservoir more efficiently, but trends in multi-planet systems do not support such a model, and theoretical models suggest that the gas accretion timescale is too short for migration to sweep the disc. Enhanced inner-disc mass reservoirs are possible, agreeing with predictions of disc evolution through self-gravity, but not adding to millimetre dust-flux as the inner disc is optically thick. Finally, the incidence of massive discs is shown to be higher at the {\it proto}stellar stages, Classes 0 and I, where discs substantial enough to form planets via core accretion are abundant enough to match the frequency of exo-planets. Gravitational instability may also operate in the Class 0 epoch, where half the objects have potentially unstable discs of $\ga$30 % of the stellar mass. However, recent calculations indicate that forming gas giants inside 50 AU by instability is unlikely, even in such massive discs. Overall, the results presented suggest that the canonically &#39;proto-planetary&#39; discs of Class II T Tauri stars {\bf have globally low masses in dust observable at millimetre wavelengths, and conversion to larger bodies (anywhere from small rocks up to planetary cores) must already have occurred.}

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The role of disc self-gravity in the formation of protostars and protostellar discs

We use time-dependent, one-dimensional disc models to investigate the evolution of protostellar discs that form through the collapse of molecular cloud cores and in which the primary transport mechanism is self-gravity. We assume that these discs settle into a state of thermal equilibrium with Q = 2 and that the strength of the angular momentum transport is set by the cooling rate of the disc. The results suggest that these discs will attain a quasi-steady state that persists for a number of free-fall times and in which most of the mass within 100 au is located inside 10-20 au. This pile-up of mass in the inner disc could result in temperatures that are high enough for the growth of MHD turbulence which could rapidly drain the inner disc and lead to FU Orionis-like outbursts. In all our simulations, the inner regions of the discs (r < 40 au) were stable against fragmentation, while fragmentation was possible in the outer regions (r > 40 au) of discs that formed from cores that had enough initial angular momentum to deposit sufficient mass in these outer regions. The large amounts of mass in these outer regions, however, suggests that fragmentation will lead to the formation of sub-stellar and stellar mass companions, rather than planetary mass objects. Although mass accretion rates were largely consistent with observations, the large disc masses suggest that an additional transport mechanism (such as MRI occuring in the upper layers of the disc) must operate in order to drain the remaining disc material within observed disc lifetimes.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Time-dependent models of the structure and evolution of self-gravitating protoplanetary discs

Angular momentum transport within young massive protoplanetary discs may be dominated by self-gravity at radii where the disk is too weakly ionized to allow the development of the magneto-rotational instability. We use time-dependent one-dimensional disc models, based on a local cooling time calculation of the efficiency of transport, to study the radial structure and stability (against fragmentation) of protoplanetary discs in which self-gravity is the sole transport mechanism. We find that self-gravitating discs rapidly attain a quasi-steady state in which the surface density in the inner disc is high and the strength of turbulence very low (alpha ~ 10^{-3} or less inside 5 au). Temperatures high enough to form crystalline silicates may extend out to several au at early times within these discs. None of our discs spontaneously develop regions that would be unambiguously unstable to fragmentation into substellar objects, though the outer regions (beyond 20 au) of the most massive discs are close enough to the threshold that fragmentation cannot be ruled out. We discuss how the mass accretion rates through such discs may vary with disc mass and with mass of the central star, and note that a determination of the \dot{M}-M_* relation for very young systems may allow a test of the model.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Constraints on a planetary origin for the gap in GM Aurigae&#39;s protoplanetary disc

The unusual spectral energy distribution (SED) of the classical T Tauri star GM Aurigae provides evidence for the presence of an inner disc hole extending to several au. Using a combination of hydrodynamical simulations and Monte Carlo radiative transport, we investigate whether the observed SED is consistent with the inner hole being created and maintained by an orbiting planet. We show that an ~ 2 M_Jupiter planet, orbiting at 2.5 au in a disc with mass 0.047 M_sun and radius 300 au, provides a good match both to the SED and to CO observations which constrain the velocity field in the disc. A range of planet masses is allowed by current data, but could in principle be distinguished with further observations between 3 and ~ 20 microns. Future high precision astrometric instruments should also be able to detect the motion of the central star due to an orbiting Jupiter mass planet. We argue that the small number of T Tauri stars with SEDs resembling that of GM Aur is broadly consistent with the expected statistics of embedded migrating planets.

0 citations0 reviewsNo code linkedOpen access