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

J. C. B. Papaloizou

J. C. B. Papaloizou contributes to research discovery and scholarly infrastructure.

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

Trust snapshot

Quick read

Trust 21 - Emerging
25works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

25 published item(s)

preprint2016arXiv

Consequences of tidal interaction between disks and orbiting protoplanets for the evolution of multi-planet systems with architecture resembling that of Kepler 444

We study orbital evolution of multi-planet systems with masses in the terrestrial planet regime induced through tidal interaction with a protoplanetary disk assuming that this is the dominant mechanism for producing orbital migration and circularization. We develop a simple analytic model for a system that maintains consecutive pairs in resonance while undergoing orbital circularization and migration. Migration times for each planet may be estimated once planet masses, circularization times and the migration time for the innermost planet are given. We applied it to a model system with the current architecture of Kepler 444 interacting with a protoplanetary disk, the evolution time for the system as a whole being comparable to current protoplanetary disk lifetimes. In addition we performed numerical simulations with input data obtained from this model. These indicate that although the analytic model is inexact, relatively small corrections to estimated migration rates yield systems for which period ratios vary by a minimal extent. Because of relatively large deviations from exact resonance in the observed system of up to $2\%,$ the migration times obtained in this way indicate only weak convergent migration such that a system for which the planets did not interact would contract by only $\sim 1\%$ although undergoing significant inward migration as a whole. We performed additional simulations to investigate how the system could undergo significant convergent migration before reaching its final state. These indicate migration times have to be significantly shorter and resonances significantly closer. Relative migration rates would then have to decrease allowing period ratios to increase to become more distant from resonances as the system approached its final state in the inner regions of the protoplanetary disk (abridged).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the formation of a quasi-stationary twisted disc after a tidal disruption event

We investigate misaligned accretion discs formed after tidal disruption events that occur when a star encounters a supermassive black hole. We employ the linear theory of warped accretion discs to find the shape of a disc for which the stream arising from the disrupted star provides a source of angular momentum that is misaligned with that of the black hole. For quasi-steady configurations we find that when the warp diffusion or propagation time is large compared to the local mass accretion time and/or the natural disc alignment radius is small, misalignment is favoured. These results have been verified using SPH simulations. We also simulated 1D model discs including gas and radiation pressure. As accretion rates initially exceed the Eddington limit the disc is initially advection dominated. Assuming the $α$ model for the disc, where it can be thermally unstable it subsequently undergoes cyclic transitions between high and low states. During these transitions the aspect ratio varies from $\sim 1$ to $\sim 10^{-3}$ which is reflected in changes in the degree of disc misalignment at the stream impact location. For maximal black hole rotation and sufficiently large values of viscosity parameter $α> \sim 0.01-0.1$ the ratio of the disc inclination to that of the initial stellar orbit is estimated to be $0.1-0.2$ in the advection dominated state, while reaching of order unity in the low state. Misalignment descreases with decrease of $α$, but increases as the black hole rotation parameter decreases. Thus, it is always significant when the latter is small.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the orbital evolution of a pair of giant planets in mean motion resonance

Pairs of extrasolar giant planets in a mean motion commensurability are common with 2:1 resonance occurring most frequently. Disc-planet interaction provides a mechanism for their origin. However, the time scale on which this could operate in particular cases is unclear. We perform 2D and 3D numerical simulations of pairs of giant planets in a protoplanetary disc as they form and maintain a mean motion commensurability. We consider systems with current parameters similar to those of HD 155358, 24 Sextantis and HD 60532, and disc models of varying mass, decreasing mass corresponding to increasing age. For the lowest mass discs, systems with planets in the Jovian mass range migrate inwards maintaining a 2:1 commensurability. Systems with the inner planet currently at around 1 au from the central star could have originated at a few au and migrated inwards on a time scale comparable to protoplanetary disc lifetimes. Systems of larger mass planets such as HD 60532 attain 3:1 resonance as observed. For a given mass accretion rate, results are insensitive to the disc model for the range of viscosity prescriptions adopted, there being good agreement between 2D and 3D simulations. However, in a higher mass disc a pair of Jovian mass planets passes through 2:1 resonance before attaining a temporary phase lasting a few thousand orbits in an unstable 5:3 resonance prior to undergoing a scattering. Thus finding systems in this commensurability is unlikely.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A new general normal mode approach to dynamic tides in rotating stars with realistic structure and its applications

We review our recent results on a unified normal mode approach to dynamic tides proposed in Ivanov, Papaloizou $\&$ Chernov (2013) and Chernov, Papaloizou $\&$ Ivanov (2013). Our formalism can be used whenever the tidal interactions are mainly determined by normal modes of a star with identifiable regular spectrum of low frequency modes. We provide in the text basic expressions for tidal energy and angular momentum transfer valid both for periodic and parabolic orbits, and different assumptions about efficiency of normal mode damping due to viscosity and/or non-linear effects and discuss applications to binary stars and close orbiting extrasolar planets.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Evolutionary outcomes for pairs of planets undergoing orbital migration and circularization: second order resonances and observed period ratios in Kepler's planetary systems

In order to study the origin of the architectures of low mass planetary systems, we perform numerical surveys of the evolution of pairs of coplanar planets in the mass range $(1-4)\ \rmn{M}_{\oplus}.$ These evolve for up to $2\times10^7 \rmn{yr}$ under a range of orbital migration torques and circularization rates assumed to arise through interaction with a protoplanetary disc. Near the inner disc boundary, significant variations of viscosity, interaction with density waves or with the stellar magnetic field could occur and halt migration, but allow ircularization to continue. This was modelled by modifying the migration and circularization rates. Runs terminated without an extended period of circularization in the absence of migration torques gave rise to either a collision, or a system close to a resonance. These were mostly first order with a few $\%$ terminating in second order resonances. Both planetary eccentricities were small $< 0.1$ and all resonant angles liberated. This type of survey produced only a limited range of period ratios and cannot reproduce Kepler observations. When circularization alone operates in the final stages, divergent migration occurs causing period ratios to increase. Depending on its strength the whole period ratio range between $1$ and $2$ can be obtained. A few systems close to second order commensurabilities also occur. In contrast to when arising through convergent migration, resonant trapping does not occur and resonant angles circulate. Thus the behaviour of the resonant angles may indicate the form of migration that led to near resonance.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Evolution of a disc-planet system with a binary companion on an inclined orbit

We study orbital inclination changes associated with the precession of a disc-planet system that occurs through gravitational interaction with a binary companion on an inclined orbit. We investigate whether this scenario can account for giant planets on close orbits highly inclined to the stellar equatorial plane. We obtain conditions for maintaining approximate coplanarity and test them with SPH-simulations. For parameters of interest, the system undergoes approximate rigid body precession with modest warping while the planets migrate inwards. Because of pressure forces, disc self-gravity is not needed to maintain the configuration. We consider a disc and single planet for different initial inclinations of the binary orbit to the midplane of the combined system and a system of three planets for which migration leads to dynamical instability that reorders the planets. As the interaction is dominated by the time averaged quadrupole component of the binary&#39;s perturbing potential, results for a circular orbit can be scaled to apply to eccentric orbits. The system responded adiabatically when changes to binary orbital parameters occurred on time scales exceeding the orbital period. Accordingly inclination changes are maintained under its slow removal. Thus the scenario for generating high inclination planetary orbits studied here, is promising.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the local stability of vortices in differentially rotating discs

In order to circumvent the loss of solid material through radial drift towards the central star, the trapping of dust inside persistent vortices in protoplanetary discs has often been suggested as a process that can eventually lead to planetesimal formation. Although a few special cases have been discussed, exhaustive studies of possible quasi-steady configurations available for dust-laden vortices and their stability have yet to be undertaken, thus their viability or otherwise as locations for the gravitational instability to take hold and seed planet formation is unclear. In this paper we generalise and extend the well known Kida solution to obtain a series of steady state solutions with varying vorticity and dust density distributions in their cores, in the limit of perfectly coupled dust and gas. We then present a local stability analysis of these configurations, considering perturbations localised on streamlines. Typical parametric instabilities found have growthrates of $~0.05Ω_P$, where $Ω_P$ is the angular velocity at the centre of the vortex. Models with density excess can exhibit many narrow parametric instability bands while those with a concentrated vorticity source display internal shear which significantly affects their stability. However, the existence of these parametric instabilities may not necessarily prevent the possibility of dust accumulation in vortices.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The evolution of a binary in a retrograde circular orbit embedded in an accretion disk

Supermassive black hole binaries may form as a consequence of galaxy mergers. Both prograde and retrograde orbits have been proposed. We study a binary of a small mass ratio, q, in a retrograde orbit immersed in and interacting with a gaseous accretion disk in order to estimate time scales for inward migration leading to coalescence and the accretion rate to the secondary component. We employ both semi-analytic methods and two dimensional numerical simulations, focusing on the case where the binary mass ratio is small but large enough to significantly perturb the disk. We develop the theory of type I migration for this case and determine conditions for gap formation finding that then inward migration occurs on a time scale equal to the time required for one half of the secondary mass to be accreted through the unperturbed disk, with accretion onto the secondary playing only a minor role. The semi-analytic and fully numerical approaches are in good agreement, the former being applicable over long time scales. Inward migration induced by interaction with the disk alleviates the final parsec problem. Accretion onto the secondary does not significantly affect the orbital evolution, but may have observational consequences for high accretion efficiency. The binary may then appear as two sources of radiation rotating around each other. This study should be extended to consider orbits with significant eccentricity and the effects of gravitational radiation at small length scales. Note too that torques acting between a circumbinary disk and a retrograde binary orbit may cause the mutual inclination to increase on a timescale that can be similar to, or smaller than that for orbital evolution, depending on detailed parameters. This is also an aspect for future study (abridged).

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A unified normal mode approach to dynamic tides and its application to rotating Sun-like stars

We determine the response of a uniformly rotating star to tidal perturbations due to a companion. General periodic orbits and parabolic flybys are considered. We evaluate energy and angular momentum exchange rates as a sum of contributions from normal modes allowing for dissipative processes. We consider the case when the response is dominated by the contribution of an identifiable regular spectrum of low frequency modes, such as gravity modes and evaluate it in the limit of very weak dissipation. Our formalism may be applied both to Sun-like stars with radiative cores and convective envelopes and to more massive stars with convective cores and radiative envelopes. We provide general expressions for transfer of energy and angular momentum valid for an orbit with any eccentricity. Detailed calculations are made for Sun-like stars in the slow rotation regime where centrifugal distortion is neglected in the equilibrium and the traditional approximation is made for the normal modes. We use both a WKBJ procedure and direct numerical evaluation which are found to be in good agreement for regimes of interest. Finally we use our formalism to determine the evolution time scales for an object, in an orbit of small eccentricity, around a Sun-like star in which the tidal response is assumed to occur. Systems with either no rotation or synchronous rotation are considered. Only rotationally modified gravity modes are taken into account under the assumption that wave dissipation occurs close to the stellar centre.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamical tides excited in rotating stars of different masses and ages and the formation of close in orbits

We study the tidal response of rotating solar mass stars, as well as more massive rotating stars, of different ages in the context of tidal captures leading to either giant exoplanets on close in orbits, or the formation of binary systems in star clusters. To do this, we adopt approaches based on normal mode and associated overlap integral evaluation, developed in a companion paper by Ivanov et al., and direct numerical simulation, to evaluate energy and angular momentum exchanges between the orbit and normal modes. The two approaches are found to be in essential agreement apart from when encounters occur near to pseudosynchronization, where the stellar angular velocity and the orbital angular velocity at periastron are approximately matched. We find that the strength of tidal interaction being expressed in dimensionless natural units is significantly weaker for the more massive stars, as compared to the solar mass stars, because of the lack of significant convective regions in the former case. On the other hand the interaction is found to be stronger for retrograde as opposed to prograde orbits in all cases. In addition, for a given pericentre distance, tidal interactions also strengthen for more evolved stars on account of their radial expansion. In agreement with previous work based on simplified polytropic models, we find that energy transferred to their central stars could play a significant role in the early stages of the circularisation of potential &#39;Hot Jupiters&#39;.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Close encounters of a rotating star with planets in parabolic orbits of varying inclination and the formation of Hot Jupiters

(abbreviated) We extend the theory of close encounters of a planet on a parabolic orbit with a star to include the effects of tides induced on the central rotating star. Orbits with arbitrary inclination to the stellar rotation axis are considered. We obtain results both from an analytic treatment and numerical one that are in satisfactory agreement. These results are applied to the initial phase of the tidal circularisation problem. We find that both tides induced in the star and planet can lead to a significant decrease of the orbital semi-major axis for orbits having periastron distances smaller than 5-6 stellar radii (corresponding to periods $\sim 4-5$ days after the circularisation has been completed) with tides in the star being much stronger for retrograde orbits compared to prograde orbits. We use the simple Skumanich law for the stellar rotation with its rotational period equal to one month at the age of 5Gyr. The strength of tidal interactions is characterised by circularisation time scale, $t_{ev}$ defined as a time scale of evolution of the planet&#39;s semi-major axis due to tides considered as a function of orbital period $P_{obs}$ after the process of tidal circularisation has been completed. We find that the ratio of the initial circularisation time scales corresponding to prograde and retrograde orbits is of order 1.5-2 for a planet of one Jupiter mass and $P_{obs}\sim $ four days. It grows with the mass of the planet, being of order five for a five Jupiter mass planet with the same $P_{orb}$. Thus, the effect of stellar rotation may provide a bias in the formation of planetary systems having planets on close orbits around their host stars, as a consequence of planet-planet scattering, favouring systems with retrograde orbits. The results may also be applied to the problem of tidal capture of stars in young stellar clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Outward migration of a super-Earth in a disc with outward propagating density waves excited by a giant planet

In this paper we consider a new mechanism for stopping the inward migration of a low-mass planet embedded in a gaseous protoplanetary disc. It operates when a low-mass planet (for example a super-Earth), encounters outgoing density waves excited by another source in the disc. This source could be a gas giant in an orbit interior to that of the low-mass planet. As the super-Earth passes through the wave field, angular momentum is transferred to the disc material and then communicated to the planet through coorbital dynamics, with the consequence that its inward migration can be halted or even reversed. We illustrate how the mechanism we consider works in a variety of different physical conditions employing global two-dimensional hydrodynamical calculations. We confirm our results by performing local shearing box simulations in which the super-Earth interacts with density waves excited by an independent harmonically varying potential. Finally, we discuss the constraints arising from the process considered here, on formation scenarios for systems containing a giant planet and lower mass planet in an outer orbit with a 2:1 commensurability such as GJ876.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Tidal interactions in multi-planet systems

We study systems of close orbiting planets evolving under the influence of tidal circularization. It is supposed that a commensurability forms through the action of disk induced migration and orbital circularization. After the system enters an inner cavity or the disk disperses the evolution continues under the influence of tides due to the central star which induce orbital circularization. We derive approximate analytic models that describe the evolution away from a general first order resonance that results from tidal circularization in a two planet system and which can be shown to be a direct consequence of the conservation of energy and angular momentum. We consider the situation when the system is initially very close to resonance and also when the system is between resonances. We also perform numerical simulations which confirm these models and then apply them to two and four planet systems chosen to have parameters related to the GJ581 and HD10180 systems. We also estimate the tidal dissipation rates through effective quality factors that could result in evolution to observed period ratios within the lifetimes of the systems. Thus the survival of, or degree of departure from, close commensurabilities in observed systems may be indicative of the effectiveness of tidal disipation, a feature which in turn may be related to the internal structure of the planets involved.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dynamic tides in rotating objects: a numerical investigation of inertial waves in fully convective or barotropic stars and planets

We perform direct numerical simulations of the tidal encounter of a rotating planet on a highly eccentric or parabolic orbit about a central star formulated as an initial value problem. This approach enables us to extend previous work of Ivanov & Papaloizou to consider planet models with solid cores and to avoid making an anelastic approximation. We obtain a power spectrum of the tidal response of coreless models which enables global inertial modes to be identified. Their frequencies are found to be in good agreement with those obtained using either a WKBJ approach or the anelastic spectral approach adopted in previous work for small planet rotation rates. We also find that the dependence of the normal mode frequencies on the planet angular velocity in case of higher rotation rates can for the most part be understood by applying first order perturbation theory to the anelastic modes. We calculate the energy and angular momentum exchanged as a result of the tidal encounter and for coreless models again find good agreement with results obtained using either the anelastic spectral method. Models with a solid core showed evidence of the emission of shear layers at critical latitudes and possibly wave attractors after the encounter but the total energy exchanged during the encounter did not differ dramatically from the coreless case as long as the ratio of the core radius to the total radius was less than 50%, there being hardly any difference at all when this ratio was less than 25% of the total radius. We give a physical and mathematical interpretation of this result. Finally we are able to validate the use of the anelastic approximation for both the work presented here and our previous work which led to estimates of circularisation rates for planets in highly eccentric orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Inertial waves in rotating bodies: a WKBJ formalism for inertial modes and a comparison with numerical results

(abbreviated) In this paper we develop a consistent WKBJ formalism, together with a formal first order perturbation theory for calculating the properties of the inertial modes of a uniformly rotating coreless body (modelled as a polytrope and referred hereafter to as a planet) under the assumption of a spherically symmetric structure. The eigenfrequencies, spatial form of the associated eigenfunctions and other properties we obtained analytically using the WKBJ eigenfunctions are in good agreement with corresponding results obtained by numerical means for a variety of planet models even for global modes with a large scale distribution of perturbed quantities. This indicates that even though they are embedded in a dense spectrum, such modes can be identified and followed as model parameters changed and that first order perturbation theory can be applied. This is used to estimate corrections to the eigenfrequencies as a consequence of the anelastic approximation, which we argue here to be small when the rotation frequency is small. These are compared with simulation results in an accompanying paper with a good agreement between theoretical and numerical results. The results reported here may provide a basis of theoretical investigations of inertial waves in many astrophysical and other applications, where a rotating body can be modelled as a uniformly rotating barotropic object, for which the density has, close to its surface, an approximately power law dependence on distance from the surface.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Migration of a moonlet in a ring of solid particles : Theory and application to Saturn&#39;s propellers

Hundred meter sized objects have been identified by the Cassini spacecraft in Saturn&#39;s A ring through the so-called &#34;propeller&#34; features they create in the ring. These moonlets should migrate, due to their gravitational interaction with the ring ; in fact, some orbital variation have been detected. The standard theory of type I migration of planets in protoplanetary disks can&#39;t be applied to the ring system, as it is pressureless. Thus, we compute the differential torque felt by a moonlet embedded in a two-dimensional disk of solid particles, with flat surface density profile, both analytically and numerically. We find that the corresponding migration rate is too small to explain the observed variations of the propeller&#39;s orbit in Saturn&#39;s A-ring. However, local density fluctuations (due to gravity wakes in the marginally gravitationally stable A-ring) may exert a stochastic torque on a moonlet. Our simulations show that this torque can be large enough to account for the observations, depending on the parameters of the rings. We find that on time scales of several years the migration of propellers is likely to be dominated by stochastic effects (while the former, non-stochastic migration dominates after ~ 10^{4-5} years). In that case, the migration rates provided by observations so far suggests that the surface density of the A ring should be of the order of 700 kg/m^2. The age of the propellers shouldn&#39;t exceed 1 to 100 million years, depending on the dominant migration regime.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The subcritical baroclinic instability in local accretion disc models

(abridged) Aims: We present new results exhibiting a subcritical baroclinic instability (SBI) in local shearing box models. We describe the 2D and 3D behaviour of this instability using numerical simulations and we present a simple analytical model describing the underlying physical process. Results: A subcritical baroclinic instability is observed in flows stable for the Solberg-Hoiland criterion using local simulations. This instability is found to be a nonlinear (or subcritical) instability, which cannot be described by ordinary linear approaches. It requires a radial entropy gradient weakly unstable for the Schwartzchild criterion and a strong thermal diffusivity (or equivalently a short cooling time). In compressible simulations, the instability produces density waves which transport angular momentum outward with typically alpha<3e-3, the exact value depending on the background temperature profile. Finally, the instability survives in 3D, vortex cores becoming turbulent due to parametric instabilities. Conclusions: The subcritical baroclinic instability is a robust phenomenon, which can be captured using local simulations. The instability survives in 3D thanks to a balance between the 2D SBI and 3D parametric instabilities. Finally, this instability can lead to a weak outward transport of angular momentum, due to the generation of density waves by the vortices.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Vortex migration in protoplanetary disks

We consider the radial migration of vortices in two-dimensional isothermal gaseous disks. We find that a vortex core, orbiting at the local gas velocity, induces velocity perturbations that propagate away from the vortex as density waves. The resulting spiral wave pattern is reminiscent of an embedded planet. There are two main causes for asymmetries in these wakes: geometrical effects tend to favor the outer wave, while a radial vortensity gradient leads to an asymmetric vortex core, which favors the wave at the side that has the lowest density. In the case of asymmetric waves, which we always find except for a disk of constant pressure, there is a net exchange of angular momentum between the vortex and the surrounding disk, which leads to orbital migration of the vortex. Numerical hydrodynamical simulations show that this migration can be very rapid, on a time scale of a few thousand orbits, for vortices with a size comparable to the scale height of the disk. We discuss the possible effects of vortex migration on planet formation scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The excitation of spiral density waves through turbulent fluctuations in accretion discs I: WKBJ theory

We study and elucidate the mechanism of spiral density wave excitation in a differentially rotating turbulent flow. We formulate a set of wave equations with sources that are only non-zero in the presence of turbulent fluctuations. We solve these in a shearing box domain using a WKBJ method. It is found that, for a particular azimuthal wave length, the wave excitation occurs through a sequence of regularly spaced swings during which the wave changes from leading to trailing form. This is a generic process that is expected to occur in shearing discs with turbulence. Trailing waves of equal amplitude propagating in opposite directions are produced, both of which produce an outward angular momentum flux that we give expressions for as functions of the disc parameters and azimuthal wave length. By solving the wave amplitude equations numerically we justify the WKBJ approach for a Keplerian rotation law for all parameter regimes of interest. In order to quantify the wave excitation completely the important wave source terms need to be specified. Assuming conditions of weak nonlinearity, these can be identified and are associated with a quantity related to the potential vorticity, being the only survivors in the linear regime. Under the additional assumption that the source has a flat power spectrum at long azimuthal wave lengths, the optimal azimuthal wave length produced is found to be determined solely by the WKBJ response and is estimated to be 2 pi H, with H being the nominal disc scale height.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The excitation of spiral density waves through turbulent fluctuations in accretion discs II: Numerical Simulations with MRI driven turbulence

We present fully three-dimensional local MRI simulations with the object of studying the excitation of non-axisymmetric spiral density waves that are observed to always be present in such simulations. They are potentially important for affecting protoplanetary migration through the action of associated stochastic gravitational forces and producing residual transport in MHD inactive regions. The simulations we perform are with zero net flux and produce mean activity levels of alpha ~ 0.005. We reveal the nature of the mechanism responsible for the excitation of these waves by determining the time dependent evolution of the Fourier transforms of the participating state variables. The dominant waves are found to have no vertical structure and to be excited during periodically repeating swings in which they change from leading to trailing. The initial phase of the evolution of such a swing is found to be in excellent agreement with that expected from the WKBJ theory developed in a preceding paper by Heinemann & Papaloizou. We demonstrate that the important source terms causing excitation of the waves are related to a quantity that reduces to the potential vorticity for small perturbations from the background state with no vertical dependence. We find that the root mean square density fluctuations associated with the waves are positively correlated with both this quantity and the general level of hydromagnetic turbulence. The mean angular momentum transport associated with spiral density waves generated in our simulations is estimated to be a significant fraction of that associated with the turbulent Reynolds stress.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

One-armed oscillations in Be star discs

In this paper we study the effect of the quadrupole-term in the gravitational potential of a rotationally deformed central Be star on one armed density waves in the circumstellar disc. The aim is to explain the observed long-term violet over red (V/R) intensity variations of the double peaked Balmer emission-lines, not only in cool Be star systems, but also in the hot systems like gamma Cas. We have carried out semi-analytic and numerical studies of low-frequency one armed global oscillations in near Keplerian discs around Be stars. In these we have investigated surface density profiles for the circumstellar disc which have inner narrow low surface density or gap regions, just interior to global maxima close to the rapidly rotating star, as well as the mode inner boundary conditions. Our results indicate that it is not necessary to invoke extra forces such as caused by line absorption from the stellar flux in order to explain the long-term V/R variations in the discs around massive Be stars. When there exists a narrow gap between the star and its circumstellar disc, with the result that the radial velocity perturbation is non-zero at the inner disc boundary, we find oscillation (and V/R) periods in the observed range for plausible magnitudes for the rotational quadrupole term.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

A planet on an inclined orbit as an explanation of the warp in the Beta Pictoris disk

We consider the deformation that has recently been observed in the inner part of the circumstellar disk around Beta Pictoris with the HST. Our recent ground based adaptive optics coronographic observations confirm that the inner disk is warped. We investigate the hypothesis that a yet undetected planet is responsible for the observed warp, through simulations of the effect of the gravitational perturbation due to a massive companion on the disk. The physical processes assumed in the simulations are discussed: since the observed particles do not survive collisions, the apparent disk shape is driven by the underlying collisionless parent population. The resulting possible parameters for the planet that are consistent with the observed disk deformation are reviewed.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

The tidally induced warping, precession and truncation of accretion discs in binary systems: three-dimensional simulations

We present the results of non linear, hydrodynamic simulations, in three dimensions, of the tidal perturbation of accretion discs in binary systems where the orbit is circular and not necessarily coplanar with the disc mid-plane. The accretion discs are assumed to be geometrically thin, and of low mass relative to the stellar mass so that they are governed by thermal pressure and viscosity, but not self-gravity. The parameters that we consider in our models are the ratio of the orbital distance to the disc radius, D/R, the binary mass ratio, the initial inclination angle between the orbit and disc planes and the Mach number in the outer parts of the unperturbed disc. For binary mass ratios of around unity and D/R in the range 3 to 4, we find that the global evolution of the discs is governed primarily by the value of the Mach number. For relatively low Mach numbers (i.e. 10 to 20) we find that the discs develop a mildly warped structure, are tidally truncated, and undergo a near rigid body precession at a rate which is in close agreement with analytical arguments. For higher Mach numbers (approx. 30), the evolution is towards a considerably more warped structure, but the disc none the less maintains itself as a long-lived, coherent entity. Additionally, it is found that the inclination angle between the disc and orbital angular momentum vectors evolves on a longer timescale which is probably the viscous evolution timescale of the disc. The calculations are relevant to a number of observed astrophysical phenomena, including the precession of jets associated with young stars, the high spectral index of some T-Tauri stars, and the light curves of X-ray binaries such as Hercules X-1 which suggest the presence of precessing accretion discs.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Interchange instability in an accretion disc with a poloidal magnetic field

We investigate the stability to nonaxisymmetric perturbations of an accretion disc in which a poloidal magnetic field provides part of the radial support against gravity. Interchange instability due to radial gradients in the magnetic field are strongly stabilized by the shear flow in the disc. For smooth field distributions this instability is restricted to discs in which the magnetic energy is comparable to the gravitational energy. An incompressible model for the instability akin to the Boussinesq approximation for convection is given which predicts the behaviour of the instability accurately. Global axisymmetric disturbances are also considered and found to be stable for a certain class of models. The results indicate that accretion discs may be able to support poloidal fields which are strong enough to suppress other forms of magnetic instability. These strong and stable field distributions are likely to be well suited for the magnetic acceleration of jets and winds.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

A comparison of one and two dimensional models of transonic accretion discs around collapsed objects

We construct models of the inner part of a transonic adiabatic accretion disc assuming constant specific angular momentum taking the vertical structure fully into account.\\ For comparison purposes, we construct the corresponding one dimensional viscous disc models derived under vertical averaging assumptions. The conditions under which a unique location for the critical/sonic point is obtained, given an appropriate set of exterior boundary conditions for these models, is also discussed. This is not unique if the standard &#39;$α$&#39; prescription with viscous stress proportional to the angular velocity gradient is used.\\ We use a simple model to discuss the possible limitations on the form of the viscous stress arising from the requirement that viscous information must travel at a finite speed. Contrary to results in the existing literature, the viscous stress tends to be {\it increased} rather than reduced for the type of flows we consider in which the angular momentum and angular velocity gradients have opposite signs. However, finite propagation effects may result in a unique location for the sonic point.\\ We found good agreement between the radial flow and specific angular momentum profiles in the inner regions of the one dimensional models and those in the equatorial plane for corresponding two dimensional models which may be matched for a range of $α$ between 0.1 and 10$^{-4}$.

0 citations0 reviewsNo code linkedOpen access