Observations of pressure-gravity mixed modes, combined with a theoretical framework for understanding mode formation, can yield a wealth of information about deep stellar interiors. In this paper, we seek to develop a formalism for treating the effects of deeply buried core magnetic fields on mixed modes in evolved stars, where the fields are moderate, i.e. not strong enough to disrupt wave propagation, but where they may be too strong for non-degenerate first-order perturbation theory to be applied. The magnetic field is incorporated in a way that avoids having to use this. Inclusion of the Lorentz force term is shown to yield a system of differential equations that allows for the magnetically-affected eigenfunctions to be computed from scratch, rather than following the approach of first-order perturbation theory. For sufficiently weak fields, coupling between different spherical harmonics can be neglected, allowing for reduction to a second-order system of ordinary differential equations akin to the usual oscillation equations that can be solved analogously. We derive expressions for (i) the mixed-mode quantisation condition in the presence of a field and (ii) the frequency shift associated with the magnetic field. In addition, for modes of low degree we uncover an extra offset term in the quantisation condition that is sensitive to properties of the evanescent zone. These expressions may be inverted to extract information about the stellar structure and magnetic field from observational data.

We calculate the dynamical tides raised by a close planetary companion on non-rotating stars of $1$ $\text{M}_{\odot}$ and $1.4$ $\text{M}_{\odot}$. Using the Henyey method, we solve the fully non-adiabatic equations throughout the star. The horizontal Lagrangian displacement is found to be 10 to 100 times larger than the equilibrium tide value in a thin region near the surface of the star. This is because non--adiabatic effects dominate in a region that extends from below the outer edge of the convection zone up to the stellar surface, and the equilibrium tide approximation is inconsistent with non--adiabaticity. Although this approximation generally applies in the low frequency limit, it also fails in the parts of the convection zone where the forcing frequency is small but larger than the Brunt-V\"ais\"al\"a frequency. We derive analytical estimates which give a good approximation to the numerical values of the magnitude of the ratio of the horizontal and radial displacements at the surface. The relative surface flux perturbation is also significant, on the order of 0.1 % for a system modelled on 51 Pegasi b. Observations affected by the horizontal displacement may therefore be more achievable than previously thought, and brightness perturbations may be the result of flux perturbations rather than due to the radial displacement. We discuss the implication of this on the possibility of detecting such tidally excited oscillations, including the prospect of utilising the large horizontal motion for observations of systems such as 51 Pegasi.

The discovery of ring systems around objects of the outer Solar System provides a strong motivation to apply theoretical models in order to better estimate their physical and orbital parameters, which can constrain scenarios for their origin. We review the criterion for maintaining apse-alignment across a ring and the balance between the energy input rate provided by a close by satellite and the internal dissipation rate occurring through ring particle collisions that is required to maintain ring eccentricity, as derived from the equations of motion governing the Lagrangian-displacements of the ring-particle orbits. We use the case of the epsilon-ring of Uranus, to calibrate our theoretical discussion and illustrate the basic dynamics governing these types of ring. In the case of the ring system of (10199) Chariklo, where the evidence that the rings are eccentric is not conclusive, we apply the theory of apse-alignment to derive information about the most plausible combination of values of the surface density and eccentricity-gradient, as well as the masses and locations of their postulated but -presently undetected-shepherd-satellites. When the balance conditions that we predict are applied to the ring system of (10199) Chariklo, we are able to estimate the minimum mass of a shepherd satellite required to prevent eccentricity decay, as a function of its orbital location, for two different models of dissipation. We conclude that the satellite mass required to maintain the m = 1 eccentric mode in the ring, would be similar or smaller than that needed to confine the rings radially.

We study the formation of the 9:7 mean motion resonance in a system of two low-mass planets ($m_{1}=m_{2}=3M_{\oplus}$) embedded in a gaseous protoplanetary disk employing a full 2D hydrodynamic treatment of the disk-planet interactions. Our aim is to determine the disk properties that favor a capture of two equal-mass super-Earths into this second-order resonance. For this purpose, we have performed a series of numerical hydrodynamic simulations of the system of two super-Earths migrating in disks with a variety of different initial parameters and found conditions for the permanent or temporary locking in the 9:7 resonance. We observe that capture occurs during the convergent migration of planets if their resonance angle at the moment of arrival at the resonance assumes values in a certain range (inside a window of capture). The width of such a window depends on the relative migration and circularization rates that are determined by the disk parameters. The window is wide if the relative migration rate is slow, and it becomes narrower as the relative migration rate increases. The window will be closed if the migration rate is sufficiently high, and the capture will not take place. We illustrate also how the 9:7 resonance window of capture is affected by the initial eccentricities and the initial orbits of the planets.

This paper focuses on two-planet systems in a first-order $(q+1):q$ mean motion resonance and undergoing type-I migration in a disc. We present a detailed analysis of the resonance valid for any value of $q$. Expressions for the equilibrium eccentricities, mean motions and departure from exact resonance are derived in the case of smooth convergent migration. We show that this departure, not assumed to be small, is such that period ratio normally exceeds, but can also be less than, $ (q+1)/q.$ Departure from exact resonance as a function of time for systems starting in resonance and undergoing divergent migration is also calculated. We discuss observed systems in which two low mass planets are close to a first-order resonance. We argue that the data are consistent with only a small fraction of the systems having been captured in resonance. Furthermore, when capture does happen, it is not in general during smooth convergent migration through the disc but after the planets reach the disc inner parts. We show that although resonances may be disrupted when the inner planet enters a central cavity, this alone cannot explain the spread of observed separations. Disruption is found to result in either the system moving interior to the resonance by a few percent, or attaining another resonance. We postulate two populations of low mass planets: a small one for which extensive smooth migration has occurred, and a larger one that formed approximately in-situ with very limited migration.

We consider misaligned accretion discs formed after tidal disruption events occurring when a star encounters a supermassive rotating black hole. We use the linear theory of warped accretion discs to find the disc shape when the stream produced by the disrupted star provides a source of mass and angular momentum that is misaligned with the black hole. The evolution of the surface density and aspect ratio is found from a one dimensional vertically averaged model. We extend previous work which assumed a quasi-stationary disc to allow unrestricted dynamical propagation of disc tilt and twist through time dependent backgrounds. We consider a smaller value of the viscosity parameter, $\alpha =0.01,$ finding the dynamics varies significantly. At early times the disc inclination is found to be nearly uniform at small radii where the aspect ratio is large. However, since torques arise from the Lense-Thirring effect and the stream there is non uniform precession. We propose a simple model for this requiring only the background surface density and aspect ratio. At these times the $\alpha \sim 0.01$ disc exhibits a new feature. An inclined hot inner region joins an outer low inclination cool region via a thin transition front propagating outwards with a speed exceeding that of bending waves in the cool region. These waves accumulate where the propagation speeds match producing an inclination spike separating inner and outer discs. At late times a sequence of quasi-stationary configurations approximates disc shapes at small radii. We discuss observational implications of our results.

The spectrum of oscillation modes of a star provides information not only about its material properties (e.g. mean density), but also its symmetries. Spherical symmetry can be broken by rotation and/or magnetic fields. It has been postulated that strong magnetic fields in the cores of some red giants are responsible for their anomalously weak dipole mode amplitudes (the "dipole dichotomy" problem), but a detailed understanding of how gravity waves interact with strong fields is thus far lacking. In this work, we attack the problem through a variety of analytical and numerical techniques, applied to a localised region centred on a null line of a confined axisymmetric magnetic field which is approximated as being cylindrically symmetric. We uncover a rich variety of phenomena that manifest when the field strength exceeds a critical value, beyond which the symmetry is drastically broken by the Lorentz force. When this threshold is reached, the spatial structure of the g-modes becomes heavily altered. The dynamics of wave packet propagation transitions from regular to chaotic, which is expected to fundamentally change the organisation of the mode spectrum. In addition, depending on their frequency and the orientation of field lines with respect to the stratification, waves impinging on different parts of the magnetised region are found to undergo either reflection or trapping. Trapping regions provide an avenue for energy loss through Alfven wave phase mixing. Our results may find application in various astrophysical contexts, including the dipole dichotomy problem, the solar interior, and compact star oscillations.

We study the dynamical evolution of the TRAPPIST-1 system under the influence of orbital circularization through tidal interaction with the central star. We find that systems with parameters close to the observed one evolve into a state where consecutive planets are linked by first order resonances and consecutive triples, apart from planets c, d and e, by connected three body Laplace resonances. The system expands with period ratios increasing and mean eccentricities decreasing with time. This evolution is largely driven by tides acting on the innermost planets which then influence the outer ones. In order that deviations from commensurability become significant only on $Gy$ time scales or longer, we require that the tidal parameter associated with the planets has to be such that $Q' > \sim 10^{2-3}.$ At the same time, if we start with two subsystems, with the inner three planets comprising the inner one, $Q'$ associated with the planets has to be on the order (and not significantly exceeding) $10^{2-3}$ for the two subsystems to interact and end up in the observed configuration. This scenario is also supported by modelling of the evolution through disk migration which indicates that the whole system cannot have migrated inwards together. Also in order to avoid large departures from commensurabilities, the system cannot have stalled at a disk inner edge for significant time periods. We discuss the habitability consequences of the tidal dissipation implied by our modelling, concluding that planets d, e and f are potentially in habitable zones.

If sufficiently irradiated by its central star, a protoplanetary disks falls into an equilibrium state exhibiting vertical shear. This state may be subject to a hydrodynamical instability, the `vertical shear instability' (VSI), whose breakdown into turbulence transports a moderate amount of angular momentum while also facilitating planet formation, possibly via the production of small-scale vortices. In this paper, we show that VSI modes (a) exhibit arbitrary spatial profiles and (b) remain nonlinear solutions to the incompressible local equations, no matter their amplitude. The modes are themselves subject to parasitic Kelvin-Helmholtz instability, though the disk rotation significantly impedes the parasites and permits the VSI to attain large amplitudes (fluid velocities < 10% the sound speed). This `delay' in saturation probably explains the prominence of the VSI linear modes in global simulations. More generally, the parasites may set the amplitude of VSI turbulence in strongly irradiated disks. They are also important in breaking the axisymmetry of the flow, via the unavoidable formation of vortices. The vortices, however, are not aligned with the orbital plane and thus express a pronounced $z$-dependence. We also briefly demonstrate that the vertical shear has little effect on the magnetorotational instability, whereas magnetic fields easily quench the VSI, a potential issue in the ionised surface regions of the disk and also at larger radii.

Local models of gaseous accretion discs have been successfully employed for decades to describe an assortment of small-scale phenomena, from instabilities and turbulence, to dust dynamics and planet formation. For the most part, they have been derived in a physically motivated but essentially ad hoc fashion, with some of the mathematical assumptions never made explicit nor checked for consistency. This approach is susceptible to error, and it is easy to derive local models that support spurious instabilities or fail to conserve key quantities. In this paper we present rigorous derivations, based on an asympototic ordering, and formulate a hierarchy of local models (incompressible, Boussinesq, and compressible), making clear which is best suited for a particular flow or phenomenon while spelling out explicitly the assumptions and approximations of each. We also discuss the merits of the anelastic approximation, emphasising that anelastic systems struggle to conserve energy unless strong restrictions are imposed on the flow. The problems encountered by the anelastic approximation are exacerbated by the disk's differential rotation, but also attend non-rotating systems such as stellar interiors. We conclude with a defence of local models and their continued utility in astrophysical research.

The core accretion hypothesis posits that planets with significant gaseous envelopes accreted them from their protoplanetary discs after the formation of rocky/icy cores. Observations indicate that such exoplanets exist at a broad range of orbital radii, but it is not known whether they accreted their envelopes in situ, or originated elsewhere and migrated to their current locations. We consider the evolution of solid cores embedded in evolving viscous discs that undergo gaseous envelope accretion in situ with orbital radii in the range $0.1-10\rm au$. Additionally, we determine the long-term evolution of the planets that had no runaway gas accretion phase after disc dispersal. We find: (i) Planets with $5 \rm M_{\oplus}$ cores never undergo runaway accretion. The most massive envelope contained $2.8 \rm M_{\oplus}$ with the planet orbiting at $10 \rm au$. (ii) Accretion is more efficient onto $10 \rm M_{\oplus}$ and $15 \rm M_{\oplus}$ cores. For orbital radii $a_{\rm p} \ge 0.5 \rm au$, $15 \rm M_{\oplus}$ cores always experienced runaway gas accretion. For $a_{\rm p} \ge 5 \rm au$, all but one of the $10 \rm M_{\oplus}$ cores experienced runaway gas accretion. No planets experienced runaway growth at $a_{\rm p} = 0.1 \rm au$. (iii) We find that, after disc dispersal, planets with significant gaseous envelopes cool and contract on Gyr time-scales, the contraction time being sensitive to the opacity assumed. Our results indicate that Hot Jupiters with core masses $\lesssim 15 \rm M_{\oplus}$ at $\lesssim 0.1 \rm au$ likely accreted their gaseous envelopes at larger distances and migrated inwards. Consistently with the known exoplanet population, Super-Earths and mini-Neptunes at small radii during the disc lifetime, accrete only modest gaseous envelopes.

We study the effect of dynamical tides associated with the excitation of gravity waves in an interior radiative region of the central star on orbital evolution in observed systems containing Hot Jupiters. We consider WASP-43, Ogle-tr-113, WASP-12, and WASP-18 which contain stars on the main sequence (MS). For these systems there are observational estimates regarding the rate of change of the orbital period. We also investigate Kepler-91 which contains an evolved giant star. We adopt the formalism of Ivanov et al. for calculating the orbital evolution. For the MS stars we determine expected rates of orbital evolution under different assumptions about the amount of dissipation acting on the tides, estimate the effect of stellar rotation for the two most rapidly rotating stars and compare results with observations. All cases apart from possibly WASP-43 are consistent with a regime in which gravity waves are damped during their propagation over the star. However, at present this is not definitive as observational errors are large. We find that although it is expected to apply to Kepler-91, linear radiative damping cannot explain this dis- sipation regime applying to MS stars. Thus, a nonlinear mechanism may be needed. Kepler-91 is found to be such that the time scale for evolution of the star is comparable to that for the orbit. This implies that significant orbital circularisation may have occurred through tides acting on the star. Quasi-static tides, stellar winds, hydrodynamic drag and tides acting on the planet have likely played a minor role.

Stars are self-gravitating fluids in which pressure, buoyancy, rotation and magnetic fields provide the restoring forces for global modes of oscillation. Pressure and buoyancy energetically dominate, while rotation and magnetism are generally assumed to be weak perturbations and often ignored. However, observations of anomalously weak dipole mode amplitudes in red giant stars suggest that a substantial fraction of these are subject to an additional source of damping localised to their core region, with indirect evidence pointing to the role of a deeply buried magnetic field. It is also known that in many instances the gravity-mode character of affected modes is preserved, but so far no effective damping mechanism has been proposed that accommodates this aspect. Here we present such a mechanism, which damps the oscillations of stars harbouring magnetised cores via resonant interactions with standing Alfv\'en modes of high harmonic index. The damping rates produced by this mechanism are quantitatively on par with those associated with turbulent convection, and in the range required to explain observations, for realistic stellar models and magnetic field strengths. Our results suggest that magnetic fields can provide an efficient means of damping stellar oscillations without needing to disrupt the internal structure of the modes, and lay the groundwork for an extension of the theory of global stellar oscillations that incorporates these effects.

We study, by means of numerical simulations and analysis, the details of the accretion process from a disc onto a binary system. We show that energy is dissipated at the edge of a circumbinary disc and this is associated with the tidal torque that maintains the cavity: angular momentum is transferred from the binary to the disc through the action of compressional shocks and viscous friction. These shocks can be viewed as being produced by fluid elements which drift into the cavity and, before being accreted, are accelerated onto trajectories that send them back to impact the disc. The rate of energy dissipation is approximately equal to the product of potential energy per unit mass at the disc's inner edge and the accretion rate, estimated from the disc parameters just beyond the cavity edge, that would occur without the binary. For very thin discs, the actual accretion rate onto the binary may be significantly less. We calculate the energy emitted by a circumbinary disc taking into account energy dissipation at the inner edge and also irradiation arising there from reprocessing of light from the stars. We find that, for tight PMS binaries, the SED is dominated by emission from the inner edge at wavelengths between 1-4 and 10 $\mu$m. This may apply to systems like CoRoT 223992193 and V1481 Ori.

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 $\alpha$ 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 $\alpha > \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 $\alpha$, but increases as the black hole rotation parameter decreases. Thus, it is always significant when the latter is small.

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.

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

In this note we discuss the main results of a study of a massive binary with unequal mass ratio, q, embedded in an accretion disk, with its orbital rotation being opposed to that of the disk. When the mass ratio is sufficiently large, a gap opens in the disk, but the mechanism of gap formation is very different from the prograde case. Inward migration occurs on a timescale of t_ev ~ M_p/(dot M), where M_p is the mass of the less massive component (the perturber), and dot M is the accretion rate. When q<< 1, the accretion takes place mostly onto the more massive component, with the accretion rate onto the perturber being smaller than, or of order of, q^(1/3)M. However, this rate increases when supermassive binary black holes are considered and gravitational wave emission is important. We estimate a typical duration of time for which the accretion onto the perturber and gravitational waves could be detected.

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.

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.

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

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\Omega_P$, where $\Omega_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.

We study the orbital evolution of a three planet system with masses in the super-Earth regime resulting from the action of tides on the planets induced by the central star which cause orbital circularization. We consider systems either in or near to a three body commensurability for which adjacent pairs of planets are in a first order commensurability. We develop a simple analytic solution, derived from a time averaged set of equations, that describes the expansion of the system away from strict commensurability as a function of time, once a state where relevant resonant angles undergo small amplitude librations has been attained. We perform numerical simulations that show the attainment of such resonant states focusing on the Kepler 60 system. The results of the simulations confirm many of the scalings predicted by the appropriate analytic solution. We go on to indicate how the results can be applied to put constraints on the amount of tidal dissipation that has occurred in the system. For example, if the system has been in a librating state since its formation, we find that its present period ratios imply an upper limit on the time average of 1/Q', with Q' being the tidal dissipation parameter. On the other hand if a librating state has not been attained, a lower upper bound applies.

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

The great diversity of extrasolar planetary systems has challenged our understanding of how planets form, and how their orbits evolve as they form. Among the various processes that may account for this diversity, the gravitational interaction between planets and their parent protoplanetary disc plays a prominent role in shaping young planetary systems. Planet-disc forces are large, and the characteristic times for the evolution of planets orbital elements are much shorter than the lifetime of protoplanetary discs. The determination of such forces is challenging, because it involves many physical mechanisms and it requires a detailed knowledge of the disc structure. Yet, the intense research of the past few years, with the exploration of many new avenues, represents a very significant improvement on the state of the discipline. This chapter reviews current understanding of planet-disc interactions, and highlights their role in setting the properties and architecture of observed planetary systems.

The Kepler mission is dramatically increasing the number of planets known in multi-planetary systems. Many adjacent planets have orbital period ratios near resonant values, with a tendency to be larger than required for exact first-order mean-motion resonances. This intriguing feature has been shown to be a natural outcome of orbital circularization of resonant planetary pairs due to star-planet tidal interactions. However, this feature holds in multi-planetary systems with periods longer than ten days, for which tidal circularization is unlikely to provide efficient divergent evolution of the planets orbits. Gravitational interactions between planets and their parent protoplanetary disk may instead provide efficient divergent evolution. For a planet pair embedded in a disk, we show that interactions between a planet and the wake of its companion can reverse convergent migration, and significantly increase the period ratio from a near-resonant value. Divergent evolution due to wake-planet interactions is particularly efficient when at least one of the planets opens a partial gap around its orbit. This mechanism could help account for the diversity of period ratios in Kepler's multiple systems comprising super-Earth to sub-jovian planets with periods greater than about ten days. Diversity is also expected for planet pairs massive enough to merge their gap. The efficiency of wake-planet interactions is then much reduced, but convergent migration may stall with a variety of period ratios depending on the density structure in the common gap. This is illustrated for the Kepler-46 system, for which we reproduce the period ratio of Kepler-46b and c.

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 'Hot Jupiters'.

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.

We study the interaction between massive planets and a gas disc with a mass in the range expected for protoplanetary discs. We use SPH simulations to study the orbital evolution of a massive planet as well as the dynamical response of the disc for planet masses between 1 and $6\ \rmn{M_J}$ and the full range of initial relative orbital inclinations. Gap formation can occur for planets in inclined orbits. For given planet mass, a threshold relative orbital inclination exists under which a gap forms. At high relative inclinations, the inclination decay rate increases for increasing planet mass and decreasing initial relative inclination. For an initial semi-major axis of 5 AU and relative inclination of $i_0=80^\circ,$ the times required for the inclination to decay by $10^\circ$ is $\sim10^{6}\ \rmn{yr}$ and $\sim10^{5}\ \rmn{yr}$ for $1\ \rmn{M_J}$ and $6\ \rmn{M_J}$. Planets on inclined orbits warp the disc by an extent that is negligible for $1\ \rmn{M_J}$ but increases with increasing mass becoming quite significant for a planet of mass $6\ \rmn{M_J}$. We also find a solid body precession of both the total disc angular momentum vector and the planet orbital momentum vector about the total angular momentum vector. Our results illustrate that the influence of an inclined massive planet on a protoplanetary disc can lead to significant changes of the disc structure and orientation which can in turn affect the orbital evolution of the planet significantly.

[Abridged] The migration of low mass planets has been studied in hydrodynamical disc models for more than three decades, but the impact of a magnetic field in the protoplanetary disc is less known. When the disc's magnetic field is strong enough to prevent horseshoe motion, the corotation torque is replaced by a torque arising from magnetic resonances. For weak enough magnetic fields, horseshoe motion and a corotation torque exist, and recent turbulent MHD simulations have reported the existence of a new component of the corotation torque in the presence of a mean toroidal field. The aim of this paper is to investigate the physical origin and the properties of this new corotation torque. We performed MHD simulations of a low mass planet embedded in a 2D laminar disc threaded by a weak toroidal magnetic field, with the effects of turbulence modelled by a viscosity and a resistivity. We confirm that the interaction between the magnetic field and the horseshoe motion results in an additional corotation torque on the planet, which we dub the MHD torque excess. It is caused by the accumulation of the magnetic field along the downstream separatrices of the planet's horseshoe region, which gives rise to an azimuthally asymmetric underdense region at that location. The properties of the MHD torque excess are characterised by varying the slope of the density, temperature and magnetic field profiles, as well as the diffusion coefficients and the strength of the magnetic field. The sign of the MHD torque excess depends on the density and temperature gradients only, and is positive for profiles expected in protoplanetary discs. Its magnitude is in turn mainly determined by the strength of the magnetic field and the turbulent resistivity. The MHD torque excess can be strong enough to reverse migration, even when the magnetic pressure is less than one percent of the thermal pressure.

We study the dynamics of a planet on an orbit inclined with respect to a disc. If the initial inclination of the orbit is larger than some critical value, the gravitational force exerted by the disc on the planet leads to a Kozai cycle in which the eccentricity of the orbit is pumped up to large values and oscillates with time in antiphase with the inclination. On the other hand, both the inclination and the eccentricity are damped by the frictional force that the planet is subject to when it crosses the disc. We show that, by maintaining either the inclination or the eccentricity at large values, the Kozai effect provides a way of delaying alignment with the disc and circularization of the orbit. We find the critical value to be characteristically as small as about 20 degrees. Typically, Neptune or lower mass planets would remain on inclined and eccentric orbits over the disc lifetime, whereas orbits of Jupiter or higher mass planets would align and circularize. This could play a significant role in planet formation scenarios.

The recurrent outbursts that characterise low-mass binary systems reflect thermal state changes in their associated accretion discs. The observed outbursts are connected to the strong variation in disc opacity as hydrogen ionises near 5000 K. This physics leads to accretion disc models that exhibit bistability and thermal limit cycles, whereby the disc jumps between a family of cool and low accreting states and a family of hot and efficiently accreting states. Previous models have parametrised the disc turbulence via an alpha (or `eddy') viscosity. In this paper we treat the turbulence more realistically via a suite of numerical simulations of the magnetorotational instability (MRI) in local geometry. Radiative cooling is included via a simple but physically motivated prescription. We show the existence of bistable equilibria and thus the prospect of thermal limit cycles, and in so doing demonstrate that MRI-induced turbulence is compatible with the classical theory. Our simulations also show that the turbulent stress and pressure perturbations are only weakly dependent on each other on orbital times; as a consequence, thermal instability connected to variations in turbulent heating (as opposed to radiative cooling) is unlikely to operate, in agreement with previous numerical results. Our work presents a first step towards unifying simulations of full MHD turbulence with the correct thermal and radiative physics of the outbursting discs associated with dwarf novae, low-mass X-ray binaries, and possibly young stellar objects.

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.

We present numerical simulations of disc-planet interactions where the planet opens a gravitationally unstable gap in an otherwise gravitationally stable disc. In our disc models, where the outer gap edge can be unstable to global spiral modes, we find that as we increase the surface density scale the gap becomes more unstable and the planet migrates outwards more rapidly. We show that the positive torque is provided by material brought into the planet's coorbital region by the spiral arms. This material is expected to execute horseshoe turns upon approaching the planet and hence torque it. Our results suggest that standard type II migration, applicable to giant planets in non-self-gravitating viscous discs, is likely to be significantly modified in massive discs when gravitational instabilities associated with the gap occur.

We develop an analytic theory to describe spiral density waves propagating in a shearing disc in the weakly nonlinear regime. Such waves are generically found to be excited in simulations of turbulent accretion disks, in particular if said turbulence arises from the magneto-rotational instability (MRI). We derive a modified Burgers equation governing their dynamics, which includes the effects of nonlinear steepening, dispersion, and a bulk viscosity to support shocks. We solve this equation approximately to obtain nonlinear sawtooth solutions that are asymptotically valid at late times. In this limit, the presence of shocks is found to cause the wave amplitude to decrease with time as 1/t^2. The validity of the analytic description is confirmed by direct numerical solution of the full nonlinear equations of motion. The asymptotic forms of the wave profiles of the state variables are also found to occur in MRI simulations indicating that dissipation due to shocks plays a significant role apart from any effects arising from direct coupling to the turbulence.

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.

(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'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.

We study the effect of disc self-gravity on vortex-forming instabilities associated with gaps opened by a Saturn mass planet in a protoplanetary disc. It is shown analytically and confirmed through linear calculations that vortex modes with low azimuthal mode number $,m,$ are stabilised by increasing self-gravity if the basic state is fixed. However, linear calculations show that the combined effect of self-gravity through the background and through the linear response shifts the most unstable vortex mode to higher $m.$ Nonlinear hydrodynamic simulations of planetary gaps show more vortices develop with increasing strength of self-gravity. For sufficiently large disc mass the vortex instability is suppressed and replaced by a new global instability, consistent with analytical expectations. In the nonlinear regime, vortex merging is increasingly delayed as the disc mass increases and multiple vortices may persist until the end of simulations. With self-gravity, the post-merger vortex is localised in azimuth and has similar structure to a Kida-like vortex. This is unlike the case without self-gravity where vortices merge to form a single vortex extended in azimuth. We also performed a series of supplementary simulations of co-orbital Kida-like vortices and found that self-gravity enables such vortices to execute horseshoe turns upon encountering each other. As a result vortex merging is avoided on time-scales where it would occur without self-gravity. Thus we suggest that mutual repulsion of self-gravitating vortices in a rotating flow is responsible for the delayed vortex merging above. The effect of self-gravity on vortex-induced migration is briefly discussed. We found that when self-gravity is included, the vortex-induced type III migration of Lin & Papaloizou (2010) is delayed but the extent of migration is unchanged.

We study the stability of gaps opened by a giant planet in a self-gravitating protoplanetary disc. We find a linear instability associated with both the self-gravity of the disc and local vortensity maxima which coincide with gap edges. For our models, these edge modes develop and extend to twice the orbital radius of a Saturn mass planet in discs with disc-to-star mass ratio >0.06, corresponding to a Toomre Q < 1.5 at the outer disc boundary. Unlike the local vortex-forming instabilities associated with gap edges in weakly or non-self-gravitating low viscosity discs, the edge modes are global and exist only in sufficiently massive discs, but for the typical viscosity values adopted for protoplanetary discs. Analytic modelling and linear calculations show edge modes may be interpreted as a localised disturbance associated with a gap edge inducing activity in the extended disc, through the launching of density waves excited at Lindblad resonances. Nonlinear hydrodynamic simulations are performed to investigate the evolution of edge modes in disc-planet systems. The form and growth rates of unstable modes are consistent with linear theory. Their dependence on viscosity and gravitational softening is also explored. We also performed a first study of the effect of edge modes on planetary migration. We found that if edge modes develop, then the average disc-on-planet torque becomes more positive with increasing disc mass. In simulations where the planet was allowed to migrate, although a fast type III migration could be seen that was similar to that seen in non-self-gravitating discs, we found that it was possible for the planet to interact gravitationally with the spiral arms associated with an edge mode and that this could result in the planet being scattered outwards. Thus orbital migration is likely to be complex and non monotonic in massive discs of the type we consider.

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.

Hundred meter sized objects have been identified by the Cassini spacecraft in Saturn's A ring through the so-called "propeller" 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'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's orbit in Saturn'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't exceed 1 to 100 million years, depending on the dominant migration regime.

Many small moonlets, creating propeller structures, have been found in Saturn's rings by the Cassini spacecraft. We study the dynamical evolution of such 20-50m sized bodies which are embedded in Saturn's rings. We estimate the importance of various interaction processes with the ring particles on the moonlet's eccentricity and semi-major axis analytically. For low ring surface densities, the main effects on the evolution of the eccentricity and the semi-major axis are found to be due to collisions and the gravitational interaction with particles in the vicinity of the moonlet. For large surface densities, the gravitational interaction with self-gravitating wakes becomes important. We also perform realistic three dimensional, collisional N-body simulations with up to a quarter of a million particles. A new set of pseudo shear periodic boundary conditions is used which reduces the computational costs by an order of magnitude compared to previous studies. Our analytic estimates are confirmed to within a factor of two. On short timescales the evolution is always dominated by stochastic effects caused by collisions and gravitational interaction with self-gravitating ring particles. These result in a random walk of the moonlet's semi-major axis. The eccentricity of the moonlet quickly reaches an equilibrium value due to collisional damping. The average change in semi-major axis of the moonlet after 100 orbital periods is 10-100m. This translates to an offset in the azimuthal direction of several hundred kilometres. We expect that such a shift is easily observable.

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

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.

We study the type III migration of a Saturn mass planet in low viscosity discs. The planet is found to experience cyclic episodes of rapid decay in orbital radius, each amounting to a few Hill radii. We find this to be due to the scattering of large- scale vortices present in the disc. The origin and role of vortices in the context of type III migration is explored. It is shown through numerical simulations and semi- analytical modelling that spiral shocks induced by a sufficiently massive planet will extend close to the planet orbital radius. The production of vortensity across shock tips results in thin high vortensity rings with a characteristic width of the local scale height. For planets with masses equal to and above that of Saturn, the rings are co-orbital features extending the entire azimuth. Linear stability analysis show there exists unstable modes that are localised about local vortensity minima which coincide with gap edges. Simulations show that vortices are non-linear a outcome. We used hydrodynamic simulations to examine vortex-planet interactions. Their effect is present in discs with kinematic viscosity less than about an order of magnitude smaller than the typically adopted value of \nu = 10^{-5}\Omega_pr_p(0)^2, where r_p(0) and \Omega_p are the initial orbital radius and angular velocity of the planet respectively. We find that the magnitude of viscosity affects the nature of type III migration but not the extent of the orbital decay. The role of vortices as a function of initial disc mass is also explored and it is found that the amount of orbital decay during one episode of vortex-planet interaction is independent of initial disc mass. We incorporate the concept of the co-orbital mass deficit in the analysis of our results and link it to the presence of vortices at gap edges.

[Abridged] We consider the dynamics of a system of hot super-Earths or Neptunes such as HD 40307. We show that, as tidal interaction with the central star leads to small eccentricities, the planets in this system could be undergoing resonant coupling even though the period ratios depart significantly from very precise commensurability. In a three planet system, this is indicated by the fact that resonant angles librate or are associated with long term changes to the orbital elements. We propose that the planets in HD 40307 were in a strict Laplace resonance while they migrated through the disc. After entering the disc inner cavity, tidal interaction would cause the period ratios to increase from two but with the inner pair deviating less than the outer pair, counter to what occurs in HD 40307. However, the relationship between these pairs that occurs in HD 40307 might be produced if the resonance is impulsively modified by an event like a close encounter shortly after the planetary system decouples from the disc. We find this to be in principle possible for a small relative perturbation on the order of a few 1.d-3 but then we find that evolution to the present system in a reasonable time is possible only if the masses are significantly larger than the minimum masses and tidal dissipation is very effective. On the other hand we found that a system like HD 40307 with minimum masses and more realistic tidal dissipation could be produced if the eccentricity of the outermost planet was impulsively increased to about 0.15. The form of resonantly coupled tidal evolution we consider here is quite general and could be of greater significance for systems with inner planets on significantly shorter orbital periods characteristic of for example CoRoT 7 b.

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

The aim of this talk is to present the most recent advances in establishing plausible planetary system architectures determined by the gravitational tidal interactions between the planets and the disc in which they are embedded during the early epoch of planetary system formation. We concentrate on a very well defined and intensively studied process of the disc-planet interaction leading to the planet migration. We focus on the dynamics of the systems in which low-mass planets are present. Particular attention is devoted to investigation of the role of resonant configurations. Our studies, apart from being complementary to the fast progress occurring just now in observing the whole variety of planetary systems and uncovering their structure and origin, can also constitute a valuable contribution in support of the missions planned to enhance the number of detected multiple systems.

The dynamical interactions that occur in newly formed planetary systems may reflect the conditions occurring in the protoplanetary disk out of which they formed. With this in mind, we explore the attainment and maintenance of orbital resonances by migrating planets in the terrestrial mass range. Migration time scales varying between millions of years and thousands of years are considered. In the former case, for which the migration time is comparable to the lifetime of the protoplanetary gas disk, a 2:1 resonance may be formed. In the latter, relatively rapid migration regime commensurabilities of high degree such as 8:7 or 11:10 may be formed. However, in any one large-scale migration several different commensurabilities may be formed sequentially, each being associated with significant orbital evolution. We also use a simple analytic theory to develop conditions for first order commensurabilities to be formed. These depend on the degree of the commensurability, the imposed migration and circularization rates, and the planet mass ratios. These conditions are found to be consistent with the results of our simulations.

The recently discovered planetary system HD45364 which consists of a Jupiter and Saturn mass planet is very likely in a 3:2 mean motion resonance. The standard scenario to form planetary commensurabilities is convergent migration of two planets embedded in a protoplanetary disc. When the planets are initially separated by a period ratio larger than two, convergent migration will most likely lead to a very stable 2:1 resonance for moderate migration rates. To avoid this fate, formation of the planets close enough to prevent this resonance may be proposed. However, such a simultaneous formation of the planets within a small annulus, seems to be very unlikely. Rapid type III migration of the outer planet crossing the 2:1 resonance is one possible way around this problem. In this paper, we investigate this idea in detail. We present an estimate for the required convergent migration rate and confirm this with N-body and hydrodynamical simulations. If the dynamical history of the planetary system had a phase of rapid inward migration that forms a resonant configuration, we predict that the orbital parameters of the two planets are always very similar and hence should show evidence of that. We use the orbital parameters from our simulation to calculate a radial velocity curve and compare it to observations. Our model can explain the observational data as good as the previously reported fit. The eccentricities of both planets are considerably smaller and the libration pattern is different. Within a few years, it will be possible to observe the planet-planet interaction directly and thus distinguish between these different dynamical states.

Many multi-planet systems have been discovered in recent years. Some of them are in mean-motion resonances (MMR). Planet formation theory was successful in explaining the formation of 2:1, 3:1 and other low resonances as a result of convergent migration. However, higher order resonances require high initial orbital eccentricities in order to be formed by this process and these are in general unexpected in a dissipative disk. We present a way of generating large initial eccentricities using additional planets. This procedure allows us to form high order MMRs and predict new planets using a genetic N-body code.

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.

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.

We summarize the analytic model and numerical simulations of stochastically forced planets in a turbulent disk presented in a recent paper by Rein and Papaloizou. We identify two modes of libration in systems with planets in mean motion resonance which react differently to random forces. The slow mode, which mostly corresponds to motion of the angle between the apsidal lines of the two planets, is converted to circulation more readily than the fast mode which is associated with oscillations of the semi-major axes. We therefore conclude that stochastic forcing due to disk turbulence may have played an important role in shaping the configuration of observed systems in mean motion resonance. For example, it naturally provides a mechanism for accounting for the HD128311 system for which the fast mode librates and the slow mode does not.

(Abriged) The existence of large-scale and long-lived 2D vortices in accretion discs has been debated for more than a decade. They appear spontaneously in several 2D disc simulations and they are known to accelerate planetesimal formation through a dust trapping process. However, the issue of the stability of these structures to the imposition of 3D disturbances is still not fully understood, and it casts doubts on their long term survival. Aim: We present new results on the 3D stability of elliptical vortices embedded in accretion discs, based on a linear analysis and several non-linear simulations. Methods: We derive the linearised equations governing the 3D perturbations in the core of an elliptical vortex, and we show that they can be reduced to a Floquet problem. We solve this problem numerically in the astrophysical regime and we present several analytical limits for which the mechanism responsible for the instability can be explained. Finally, we compare the results of the linear analysis to some high resolution simulations. Results: We show that most anticyclonic vortices are unstable due to a resonance between the turnover time and the local epicyclic oscillation period. In addition, we demonstrate that a strong vertical stratification does not create any additional stable domain of aspect ratio, but it significantly reduces growth rates for relatively weak (and therefore elongated) vortices. Conclusions: Elliptical vortices are always unstable, whatever the horizontal or vertical aspect-ratio is. The instability can however be weak and is often found at small scales, making it difficult to detect in low-order finite-difference simulations.

Aims. We clarify the response of extrasolar planetary systems in a 2:1 mean motion commensurability with masses ranging from the super Jovian range to the terrestrial range to stochastic forcing that could result from protoplanetary disk turbulence. The behaviour of the different libration modes for a wide range of system parameters and stochastic forcing magnitudes is investigated. The growth of libration amplitudes is parameterized as a function of the relevant physical parameters. The results are applied to provide an explanation of the configuration of the HD128311 system. Methods. We first develop an analytic model from first principles without making the assumption that both eccentricities are small. We also perform numerical N-body simulations with additional stochastic forcing terms to represent the effects of putative disk turbulence. Results. Systems are quickly destabilized by large magnitudes of stochastic forcing but some stability is imparted should systems undergo a net orbital migration. The slow mode, which mostly corresponds to motion of the angle between the apsidal lines of the two planets, is converted to circulation more readily than the fast mode which is associated with oscillations of the semi-major axes. This mode is also vulnerable to the attainment of small eccentricities which causes oscillations between periods of libration and circulation. Conclusions. Stochastic forcing due to disk turbulence may have played a role in shaping the configurations of observed systems in mean motion resonance. It naturally provides a mechanism for accounting for the HD128311 system.

We study the coorbital flow for embedded, low mass planets. We provide a simple semi-analytic model for the corotation region, which is subsequently compared to high resolution numerical simulations. The model is used to derive an expression for the half-width of the horseshoe region, x_s, which in the limit of zero softening is given by x_s/r_p = 1.68(q/h)^(1/2), where q is the planet to central star mass ratio, h is the disc aspect ratio and r_p the orbital radius. This is in very good agreement with the same quantity measured from simulations. This result is used to show that horseshoe drag is about an order of magnitude larger than the linear corotation torque in the zero softening limit. Thus the horseshoe drag, the sign of which depends on the gradient of specific vorticity, is important for estimates of the total torque acting on the planet. We further show that phenomena, such as the Lindblad wakes, with a radial separation from corotation of ~ a pressure scale height H can affect x_s, even though for low-mass planets x_s << H. The effect is to distort streamlines and to reduce x_s through the action of a back pressure. This effect is reduced for smaller gravitational softening parameters and planets of higher mass, for which x_s becomes comparable to H.

We study the torque on low mass protoplanets on fixed circular orbits, embedded in a protoplanetary disc in the isothermal limit. For low mass protoplanets and large viscosity the corotation torque behaves as expected from linear theory. However, when the viscosity becomes small enough to enable horseshoe turns to occur, the linear corotation torque exists only temporarily after insertion of a planet into the disc, being replaced by the horseshoe drag first discussed by Ward. This happens after a time that is equal to the horseshoe libration period reduced by a factor amounting to about twice the disc aspect ratio. This torque scales with the radial gradient of specific vorticity, as does the linear torque, but we find it to be many times larger. If the viscosity is large enough for viscous diffusion across the coorbital region to occur within a libration period, we find that the horseshoe drag may be sustained. If not, the corotation torque saturates leaving only the linear Lindblad torques. As the magnitude of the non linear coorbital torque (horseshoe drag) is always found to be larger than the linear torque, we find that the sign of the total torque may change even for for mildly positive surface density gradients. In combination with a kinematic viscosity large enough to keep the torque from saturating, strong sustained deviations from linear theory and outward or stalled migration may occur in such cases (abridged).

We study the development of finite eccentricity in accretion disks in close binary systems using a two-dimensional grid-based numerical scheme. We perform detailed parameter studies to explore the dependence on viscosity, disk aspect ratio, the inclusion of a mass-transfer stream and the role of the boundary conditions. We consider mass ratios 0.05<q<0.3 appropriate to superoutbursting cataclysmic binary systems. Instability to the formation of a precessing eccentric disk that attains a quasi-steady state with mean eccentricity in the range 0.3-0.5 occurs readily. The shortest growth times are ~15 binary orbits for the largest viscosities and the instability mechanism is for the most part consistent with the mode-coupling mechanism associated with the 3:1 resonance proposed by Lubow. However, the results are sensitive to the treatment of the inner boundary and to the incorporation of the mass-transfer stream. In the presence of a stream we found a critical viscosity below which the disk remains circular. Incorporation of a mass-transfer stream tends to impart stability for small enough viscosity (or, equivalently, mass-transfer rate through the disk) and does assist in obtaining a prograde precession rate that is in agreement with observations. For the larger q the location of the 3:1 resonance is pushed outwards towards the Roche lobe where higher-order mode couplings and nonlinearity occur. It is likely that three-dimensional simulations that properly resolve the disk's vertical structure are required to make significant progress in this case.

We study the disc planet interactions of low-mass protoplanets embedded in a circumstellar disc. We extend the standard theory of planet migration from the usual locally isothermal assumption to include non-barotropic effects, focusing on the validity of linear theory. We compared solutions of the linear equations with results from non-linear hydrodynamic simulations, where in both cases we adopted a background entropy gradient and solved the energy equation. We show that the migration behavior of embedded planets depends critically on the background radial entropy gradient in the disc. The presence of such a gradient not only changes the corotation torque on the planet, but also always guarantees a departure from linear behavior, which gives a singular density response at corotation, in the absence of thermal or viscous diffusion. A negative entropy gradient tends to give rise to positive, non-linear corotation torques apparently produced as material executes horseshoe turns at approximately constant entropy. These torques have no counterpart in linear theory, but can be strong enough to push the planet outwards until saturation starts to occur after a horseshoe libration period. Increased thermal diffusion acts to reduce these non-linear torques, but, at the same time, it can help to prevent their saturation. In combination with a small kinematic viscosity that is able to maintain a smooth density profile the positive torque could be sustained.

(abbreviated) We consider how tight binaries consisting of a super-massive black hole of mass $M=10^{3}-10^{4}M_{\odot}$ and a white dwarf can be formed in a globular cluster. We point out that a major fraction of white dwarfs tidally captured by the black hole may be destroyed by tidal inflation during ongoing circularisation, and the formation of tight binaries is inhibited. However, some stars may survive being spun up to high rotation rates. Then the energy loss through gravitational wave emission induced by tidally excited pulsation modes and dissipation through non linear effects may compete with the increase of pulsation energy due to dynamic tides. The semi-major axes of these stars can be decreased below a 'critical' value where dynamic tides are not effective because pulsation modes retain phase coherence between successive pericentre passages. The rate of formation of such circularising stars is estimated assuming that they can be modelled as $n=1.5$ polytropes and that results of the tidal theory for slow rotators can be extrapolated to fast rotators. We estimate the total capture rate as $\sim \dot N\sim 2.5\cdot 10^{-8}M_{4}^{1.3}r_{0.1}^{-2.1}yr^{-1}$, where $M_{4}=M/10^4M_{\odot}$ and $r_{0.1}$ is the radius of influence of the black hole in units $0.1pc$. We find that the formation rate of tight pairs is approximately 10 times smaller than the total capture rate. It is used to estimate the probability of detection of gravitational waves coming from such tight binaries by LISA. We conclude that LISA may detect such binaries provided that the fraction of globular clusters with black holes in the mass range of interest is substantial and that the dispersion velocity of the cluster stars near the radius of influence of the black hole exceeds $\sim 20km/s$.

We study the dynamics of a twisted tilted disc under the influence of an external radiation field. Assuming the effect of absorption and reemission/scattering is that a pressure is applied to the disc surface where the local optical depth is of order unity, we determine the response of the vertical structure and the influence it has on the possibility of instability to warping. We derive a pair of equations describing the evolution of a small tilt as a function of radius in the small amplitude regime that applies to both the diffusive and bending wave regimes. We also study the non linear vertical response of the disc numerically using an analogous one dimensional slab model. For global warps, we find that in order for the disc vertical structure to respond as a quasi uniform shift or tilt, as has been assumed in previous work, the product of the ratio of the external radiation momentum flux to the local disc mid plane pressure, where it is absorbed, with the disc aspect ratio should be significantly less than unity. Namely, this quantity should be of the order of or smaller than the ratio of the disc gas density corresponding to the layer intercepting radiation to the mid plane density, $\lambda \ll 1$. When this condition is not satisfied the disc surface tends to adjust so that the local normal becomes perpendicular to the radiation propagation direction. In this case dynamical quantities determined by the disc twist and warp tend to oscillate with a large characteristic period $T_{*}\sim \lambda^{-1}T_{K}$, where $T_{K}$ is some 'typical' orbital period of a gas element in the disc. The possibility of warping instability then becomes significantly reduced. In addition, when the vertical response is non uniform, the possible production of shocks may lead to an important dissipation mechanism.

We study the properties of MHD turbulence driven by the magnetorotational instability (MRI) in accretion disks. We adopt the local shearing box model and focus on the special case for which the initial magnetic flux threading the disk vanishes. We employ the finite difference code ZEUS to evolve the ideal MHD equations. Performing a set of numerical simulations in a fixed computational domain with increasing resolution, we demonstrate that turbulent activity decreases as resolution increases. We quantify the turbulent activity by measuring the rate of angular momentum transport through evaluating the standard alpha parameter. We find alpha=0.004 when (N_x,N_y,N_z)=(64,100,64), alpha=0.002 when (N_x,N_y,N_z)=(128,200,128) and alpha=0.001 when (N_x,N_y,N_z)=(256,400,256). This steady decline is an indication that numerical dissipation, occurring at the grid scale is an important determinant of the saturated form of the MHD turbulence. Analysing the results in Fourier space, we demonstrate that this is due to the MRI forcing significant flow energy all the way down to the grid dissipation scale. We also use our results to study the properties of the numerical dissipation in ZEUS. Its amplitude is characterised by the magnitude of an effective magnetic Reynolds number Re_M which increases from 10^4 to 10^5 as the number of grid points is increased from 64 to 256 per scale height. The simulations we have carried out do not produce results that are independent of the numerical dissipation scale, even at the highest resolution studied. Thus it is important to use physical dissipation, both viscous and resistive, and to quantify contributions from numerical effects, when performing numerical simulations of MHD turbulence with zero net flux in accretion disks at the resolutions normally considered.

We study the influence of the choice of transport coefficients (viscosity and resistivity) on MHD turbulence driven by the magnetorotational instability (MRI) in accretion disks. We follow the methodology described in paper I: we adopt an unstratified shearing box model and focus on the case where the net vertical magnetic flux threading the box vanishes. For the most part we use the finite difference code ZEUS, including explicit transport coefficients in the calculations. However, we also compare our results with those obtained using other algorithms (NIRVANA, the PENCIL code and a spectral code) to demonstrate both the convergence of our results and their independence of the numerical scheme. We find that small scale dissipation affects the saturated state of MHD turbulence. In agreement with recent similar numerical simulations done in the presence of a net vertical magnetic flux, we find that turbulent activity (measured by the rate of angular momentum transport) is an increasing function of the magnetic Prandtl number Pm for all values of the Reynolds number Re that we investigated. We also found that turbulence disappears when the Prandtl number falls below a critical value Pm_c that is apparently a decreasing function of Re. For the limited region of parameter space that can be probed with current computational resources, we always obtained Pm_c>1. We conclude that the magnitudes of the transport coefficients are important in determining the properties of MHD turbulence in numerical simulations in the shearing box with zero net flux, at least for Reynolds numbers and magnetic Prandtl numbers that are such that transport is not dominated by numerical effects and thus can be probed using current computational resources.

In this paper, we study the propagation and stability of nonlinear sound waves in accretion disks. Using the shearing box approximation, we derive the form of these waves using a semi-analytic approach and go on to study their stability. The results are compared to those of numerical simulations performed using finite difference approaches such as employed by ZEUS as well as Godunov methods. When the wave frequency is between Omega and two Omega (where Omega is the disk orbital angular velocity), it can couple resonantly with a pair of linear inertial waves and thus undergo a parametric instability. Neglecting the disk vertical stratification, we derive an expression for the growth rate when the amplitude of the background wave is small. Good agreement is found with the results of numerical simulations performed both with finite difference and Godunov codes. During the nonlinear phase of the instability, the flow remains well organised if the amplitude of the background wave is small. However, strongly nonlinear waves break down into turbulence. In both cases, the background wave is damped and the disk eventually returns to a stationary state. Finally, we demonstrate that the instability also develops when density stratification is taken into account and so is robust. This destabilisation of freely propagating nonlinear sound waves may be important for understanding the complicated behaviour of density waves in disks that are unstable through the effects of self-gravity or magnetic fields and is likely to affect the propagation of waves that are tidally excited by objects such as a protoplanet or companion perturbing a protoplanetary disk. The nonlinear wave solutions described here as well as their stability properties were also found to be useful for testing and comparing the performance of different numerical codes.

In this paper, we study a simple model of an orbiting protoplanet in a central magnetospheric cavity, the entry into such a cavity having been proposed as a mechanism for halting inward orbital migration. We have calculated the gravitational interaction of the protoplanet with the magnetosphere using a local model and determined the rate of evolution of the orbit. The interaction is found to be determined by the outward flux of MHD waves and thus the possibility of the existence of such waves in the cavity is significant. The estimated orbital evolution rates due to gravitational and other interactions with the magnetosphere are unlikely to be significant during protoplanetary disk lifetimes.

We study the vertical settling of solid bodies in a turbulent protoplanetary disc. We consider the situation when the coupling to the gas is weak or equivalently when the particle stopping time tau_{st} due to friction with the gas is long compared to the orbital timescale Omega^{-1}. An analytical model, which takes into account the stochastic nature of the sedimentation process using a Fokker-Planck equation for the particle distribution function in phase space, is used to obtain the vertical scale height of the solid layer as a function of the vertical component of the turbulent gas velocity correlation function and the particle stopping time. This is found to be of the same form as the relation obtained for strongly coupled particles in previous work. We compare the predictions of this model with results obtained from local shearing box MHD simulations of solid particles embedded in a vertically stratified disc in which there is turbulence driven by the MRI. We find that the ratio of the dust disc thickness to the gas disc thickness satifies H_d/H=0.08 (Omega tau_{st})^{-1/2}, which is in very good agreement with the analytical model. By discussing the conditions for gravitational instability in the outer regions of protoplanetary discs in which there is a similar level of turbulence, we find that bodies in the size range 50 to 600 metres can aggregate to form Kuiper belt-like objects with characteristic radii ranging from tens to hundreds of kilometres.

The existence of extrasolar planets with short orbital periods suggests that planetary migration induced by tidal interaction with the protoplanetary disk is important. Cores and terrestrial planets may undergo migration as they form. In this paper we investigate the evolution of a population of cores with initial masses in the range 0.1-1 earth mass embedded in a disk. Mutual interactions lead to orbit crossing and mergers, so that the cores grow during their evolution. Interaction with the disk leads to orbital migration, which results in the cores capturing each other in mean motion resonances. As the cores migrate inside the disk inner edge, scatterings and mergers of planets on unstable orbits together with orbital circularization causes strict commensurability to be lost. Near commensurability however is usually maintained. All the simulations end with a population of typically between two and five planets, with masses depending on the initial mass. These results indicate that if hot super-Earths or Neptunes form by mergers of inwardly migrating cores, then such planets are most likely not isolated. We would expect to always find at least one, more likely a few, companions on close and often near-commensurable orbits. To test this hypothesis, it would be of interest to look for planets of a few to about 10 earth masses in systems where hot super-Earths or Neptunes have already been found.

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.

The discovery of close orbiting extrasolar giant planets led to extensive studies of disk planet interactions and the forms of migration that can result as a means of accounting for their location. Early work established the type I and type II migration regimes for low mass embedded planets and high mass gap forming planets respectively. While providing an attractive means of accounting for close orbiting planets intially formed at several AU, inward migration times for objects in the earth mass range were found to be disturbingly short, making the survival of giant planet cores an issue. Recent progress in this area has come from the application of modern numerical techniques which make use of up to date supercomputer resources. These have enabled higher resolution studies of the regions close to the planet and the initiation of studies of planets interacting with disks undergoing MHD turbulence. This work has led to indications of how the inward migration of low to intermediate mass planets could be slowed down or reversed. In addition, the possibility of a new very fast type III migration regime, that can be directed inwards or outwards, that is relevant to partial gap forming planets in massive disks has been investigated.

In this paper, we study the effect of MHD turbulence on the dynamics of dust particles in protoplanetary disks. We vary the size of the particles and relate the dust evolution to the turbulent velocity fluctuations. We performed numerical simulations using two Eulerian MHD codes, both based on finite difference techniques: ZEUS--3D and NIRVANA. These were local shearing box simulations incorporating vertical stratification. Both ideal and non ideal MHD simulations with midplane dead zones were carried out. The codes were extended to incorporate different models for the dust as an additional fluid component. Good agreement between results obtained using the different approaches was obtained. The simulations show that a thin layer of very small dust particles is diffusively spread over the full vertical extent of the disk. We show that a simple description obtained using the diffusion equation with a diffusion coefficient simply expressed in terms of the velocity correlations accurately matches the results. Dust settling starts to become apparent for particle sizes of the order of 1 to 10 centimeters for which the gas begins to decouple in a standard solar nebula model at 5.2 AU. However, for particles which are 10 centimeters in size, complete settling toward a very thin midplane layer is prevented by turbulent motions within the disk, even in the presence of a midplane dead zone of significant size. These results indicate that, when present, MHD turbulence affects dust dynamics in protoplanetary disks. We find that the evolution and settling of the dust can be accurately modelled using an advection diffusion equation that incorporates vertical settling. The value of the diffusion coefficient can be calculated from the turbulent velocity field when that is known for a time of several local orbits.

(abbreviated) We consider the problem of the tidal capture or circularisation from large eccentricity of a uniformly rotating object. We extend the self-adjoint formalism introduced in Papaloizou \& Ivanov 2005 (PI) to derive general expressions for the energy and angular momentum transfered when the planet or a star passes through periastron in a parabolic or highly eccentric orbit around a central mass, without making a low frequency approximation as was done in PI. We show how these can be adapted to the low frequency limit in which only inertial modes contribute to the energy and angular momentum transfer. We calculate the inertial mode eigenspectrum for planet models of one and five Jupiter masses $M_J,$ without a solid core, with different radii corresponding to different ages. We consider the multi-passage problem when there is no dissipation finding that stochastic instability resulting in the stochastic gain of inertial mode energy over many periastron passages occurs under similar conditions to those already found for the $f$ modes. We apply our calculations to the problem of the tidal circularisation of the orbits of the extra solar planets in a state of pseudo synchronisation, and find that inertial mode excitation dominates the tidal interaction for $1 M_J$ planets that start with semi- major axes less than $10 AU$ and end up on circular orbits with final period in the $4-6$ day range. It is potentially able to account for initial circularisation up to a final $6$ day period within a few $Gyr$ But in the case of $5M_J$ oscillation modes excited in the star are more important.

We review the observations of extrasolar planets, ongoing developments in theories of planet formation, orbital migration, and the evolution of multiplanet systems.

We consider the excitation of the inertial modes of a uniformly rotating fully convective body due to a close encounter with another object. This could lead to a tidal capture or orbital circularisation depending on whether the initial orbit is unbound or highly eccentric. We develop a general self-adjoint formalism for the response problem and thus solve it taking into account the inertial modes with $m=2$ for a full polytrope with $n=1.5.$ We are accordingly able to show in this case that the excitation of inertial modes dominates the response for large impact parameters and thus cannot be neglected in calculations of tidal energy and angular momentum exchange or orbital circularisation from large eccentricity.

We investigate orbital resonances expected to arise when a system of two planets, with masses in the range 1-4 Earth masses, undergoes convergent migration while embedded in a section of gaseous disc where the flow is laminar. We consider surface densities corresponding to 0.5-4 times that expected for a minimum mass solar nebula at 5.2 AU. Using hydrodynamic simulations we find that when the configuration is such that convergent migration occurs the planets can become locked in a first order commensurability for which the period ratio is (p+1)/p with p being an integer and migrate together maintaining it for many orbits. Relatively rapid convergent migration as tends to occur for disparate masses, results in commensurabilities with p larger than 2. However, in these cases the dynamics is found to have a stochastic character. When the convergent migration is slower, such as occurs in the equal mass case, lower p commensurabilities such as 3:2 are attained which show much greater stability. In one already known example of a system with nearly equal masses in the several Earth mass range (planets around pulsar PSR B1257+12) the two largest planets are intriguingly close to a 3:2 commensurability. A very similar behaviour is obtained when the systems are modeled using an N body code with simple prescriptions for the disc planet interaction. Using that, we found that an 8:7 resonance established in a hydrodynamic simulation run for 10-100 thousand orbits could be maintained for more than million orbits. Resonant capture leads to a rise in eccentricities that can be predicted using a simple analytic model constructed in this paper. We find that the system with the 8:7 commensurability is fully consistent with this prediction.

We study and review disk protoplanet interactions using local shearing box simulations. These suffer the disadvantage of having potential artefacts arising from periodic boundary conditions but the advantage, when compared to global simulations, of being able to capture much of the dynamics close to the protoplanet at high resolution for low computational cost. Cases with and without self sustained MHD turbulence are considered. The conditions for gap formation and the transition from type I migration are investigated and found to depend on whether the single parameter M_p R^3/(M_* H^3), with M_p, M_*, R and H being the protoplanet mass, the central mass, the orbital radius and the disk semi-thickness respectively exceeds a number of order unity. We also investigate the coorbital torques experienced by a moving protoplanet in an inviscid disk. This is done by demonstrating the equivalence of the problem for a moving protoplanet to one where the protoplanet is in a fixed orbit which the disk material flows through radially as a result of the action of an appropriate external torque. For sustainable coorbital torques to be realized a quasi steady state must be realized in which the planet migrates through the disk without accreting significant mass. In that case although there is sensitivity to computational parameters, in agreement with earlier work by Masset & Papaloizou (2003) based on global simulations, the coorbital torques are proportional to the migration speed and result in a positive feedback on the migration, enhancing it and potentially leading to a runaway. This could lead to a fast migration for protoplanets in the Saturn mass range in massive disks and may be relevant to the mass period correlation for extrasolar planets which gives a preponderance of sub Jovian masses at short orbital period.

We carry out a general study of the stability of astrophysical flows that appear steady in a uniformly rotating frame. Such a flow might correspond to a stellar pulsation mode or an accretion disk with a free global distortion giving it finite eccentricity. We consider perturbations arbitrarily localized in the neighbourhood of unperturbed fluid streamlines.When conditions do not vary around them, perturbations take the form of oscillatory inertial or gravity modes. However, when conditions do vary so that a circulating fluid element is subject to periodic variations, parametric instability may occur. For nearly circular streamlines, the dense spectra associated with inertial or gravity modes ensure that resonance conditions can always be satisfied when twice the period of circulation round a streamline falls within. We apply our formalism to a differentially rotating disk for which the streamlines are Keplerian ellipses, with free eccentricity up to 0.7, which do not precess in an inertial frame. We show that for small $e,$ the instability involves parametric excitation of two modes with azimuthal mode number differing by unity in magnitude which have a period of twice the period of variation as viewed from a circulating unperturbed fluid element. Instability persists over a widening range of wave numbers with increasing growth rates for larger eccentricities. The nonlinear outcome is studied in a follow up paper which indicates development of small scale subsonic turbulence.

(Abridged) We consider models of gas giant planets forming in protoplanetary disks consisting of solid cores with gaseous envelopes in contact with their critical Hill spheres while accreting gas from the surrounding disk.We suppose the luminosity derives from gas accretion alone.We label such models as type A and follow their evolution which may occur on a time scale similar to the protostellar disk lifetime until rapid gas accretion. We consider another set of models, we label type B, with a free surface, powered by gravitational contraction, while accreting through a disk.We find these models rapidly attain a radius <~ 2x10^(10)cm without subsequent expansion.We speculate that giant planet formation is initially described by models of type A, until at the onset of rapid gas accretion, there is a transition to models of type B. Protoplanet migration in standard models tends to be most effective near this transition where it also changes from type I to type II.If a mechanism prevents type I migration of low mass protoplanets, a rapid inward migration might occur near the transitional mass regime. Such protoplanets would end up in the inner disk regions undergoing type II migration and further accretion potentially becoming sub Jovian close orbiting planets. Noting that dustier more massive cores spend longer at a larger transitional mass where faster migration is expected, these may be more prone to end in close orbiters.We find the luminosity of the protoplanets during the later stages is dominated by the circumplanetary disk and protoplanet disk boundary layer.For one Jupiter mass the luminosity range is 10^-(1.5-4) L_sun$ depending on the evolutionary stage and external conditions.

We use a simple model of the dynamics of a narrow-eccentric ring, to put some constraints on some of the observable properties of the real systems.In this work we concentrate on the case of the `Titan ringlet of Saturn'.Our approach is fluid-like, since our description is based on normal modes of oscillation rather than in individual particle orbits. Thus, the rigid precession of the ring is described as a global $m=1$ mode, which originates from a standing wave superposed on an axisymmetric background. An integral balance condition for the maintenance of the $m=1$ standing-wave can be set up, in which the differential precession induced by the oblateness of the central planet must cancel the contributions of self-gravity, the resonant satellite forcing and collisional effects. We expect that in nearly-circular narrow rings dominated by self-gravity, the eccentricity varies linearly across the ring. Thus, we take a first order expansion and we derive two integral relationships from the rigid-precession condition. These relate the surface density of the ring with the eccentricity at the center, the eccentricity gradient and the location of the secular resonance. These relationships are applied to the Titan ringlet of Saturn, which has a secular resonance with the satellite Titan in which the ring precession period is close to Titan's orbital period. In this case, we estimate the mean surface density and the location of the secular resonance.

(Abridged) A simple and general description of the dynamics of a narrow eccentric ring is presented.We view an eccentric ring which precesses uniformly at a slow rate as exhibiting a global $m=1$ mode originating from a standing wave superposed on an axisymmetric background.We adopt a continuum description using the language of fluid dynamics which gives equivalent results for the secular dynamics of thin rings as the the well known description using discrete elliptical streamlines formulated by Goldreich and Tremaine (1979). We use this to discuss the non linear mode interactions that appear through the excitation of higher $m$ modes due to the coupling of the $m=1$ mode with an external satellite potential, showing that they can lead to the excitation of the $m=1$ mode through a feedback process.Two conditions for the maintainance of a steady $m=1$ mode are obtained. One,being the condition for the normal mode pattern to precess uniformly requires a balance between the differential precession induced by the oblateness of the central planet,self-gravity and collisional effects and is the continuum form of that obtained from the $N$ streamline model of Goldreich and Tremaine (1979).The other condition is for the steady maintenance of the non-zero radial action of the ring on account of the normal mode.This requires a balance between input due to eccentric resonances due to external satellites and additional collisional damping associated with the presence of the $m=1$ mode We estimate that such a balance can occur in the $\epsilon-$ring of Uranus,given its currently observed physical and orbital parameters.

We present the results of both global cylindrical disc simulations and local shearing box simulations of protoplanets interacting with a disc undergoing MHD turbulence with zero net flux magnetic fields. We investigate the nature of the disc response and conditions for gap formation. This issue is an important one for determining the type and nature of the migration of the protoplanet, with the presence of a deep gap being believed to enable slower migration. For both types of simulation we find a common pattern of behaviour for which the main parameter determining the nature of the response is $M_p R^3/(M_* H^3)$, with $M_p$, $M_*$, $R$, and $H$ being the protoplanet mass, the central mass, the orbital radius and the disc semi-thickness respectively. We find that as this parameter is increased towards 0.1, the presence of the protoplanet is first indicated by the appearance of the well known trailing wake which, although it may appear erratic on account of the turbulence, appears to be well defined. Once the above parameter exceeds a number around unity a gap starts to develop inside which the magnetic energy density tends to be concentrated in the high density wakes. This gap formation condition can be understood from simple dimensional considerations of the conditions for nonlinearity, and the balance of angular momentum transport due to Maxwell and Reynolds' stresses with that due to tidal torques. An important result is that the basic flow morphology in the vicinity of the protoplanet is very similar in both the local and global simulations. This indicates that local shearing box simulations, which are computationally less demanding, capture much of the physics of disc-planet interaction. Thus they may provide a useful tool for studying the local interaction between forming protoplanets and turbulent, protostellar discs.

We consider the tidal interaction of a fully convective primary star and a point mass. Using a normal mode decomposition we calculate the evolution of the primary angular velocity and orbit for arbitrary eccentricity e. The dissipation acting on the tidal perturbation is associated with convective turbulence. A novel feature of the Paper is that, to take into account of the fact that there is a relaxation time t_{c}, being the turn-over time of convective eddies, associated with the process, this is allowed to act non locally in time, producing a dependence of the dissipation on tidal forcing frequency. Results are expressed in terms of the Fourier coefficients of the tidal potential. We find analytical approximations for these valid for $e>0.2$. When the tidal response is frequency independent, our results are equivalent to those obtained in the standard constant time lag approximation. When there is the frequency dependence of the dissipative response, the evolution can differ drastically. In that case the system can evolve through a sequence of spin-orbit corotation resonances with Omega_{r}/Omega_{orb}=n/2, where Omega_{r} and Omega_{orb} are the rotation and orbital frequencies and n is an integer. We study this case analytically and numerically.

In this article we present results from three on-going projects related to the formation of protoplanets in protostellar discs. We present the results of simulations that model the interaction between embedded protoplanets and disc models undergoing MHD turbulence. We review the similarities and differences that arise when the disc is turbulent as opposed to laminar (but viscous), and present the first results of simulations that examine the tidal interaction between low mass protoplanets and turbulent discs. We describe the results of simulations of Jovian mass protoplanets forming in circumbinary discs, and discuss the range of possible outcomes that arise in hydrodynamic simulations. Finally, we report on some preliminary simulations of three protoplanets of Jovian mass that form approximately coevally within a protostellar disc. We describe the conditions under which such a system can form a stable three planet resonance.

Scattered light images of the optically thin dust disk around the 5 Myr old star HD141569 have revealed its complex asymmetric structure. We show in this paper that the surface density inferred from the observations presents similarities with that expected from a circumprimary disk within a highly eccentric binary system. We assume that either the two M stars in the close vicinity of HD141569 are bound companions or at least one of them is an isolated binary companion. We discuss the resulting interaction with an initially axisymmetric disk. This scenario accounts for the formation of a spiral structure, a wide gap in the disk and a broad faint extension outside the truncation radius of the disk after 10-15 orbital periods with no need for massive companion(s) in the midst of the disk resolved in scattered light. The simulations match the observations and the star age if the perturber is on an elliptic orbit with a periastron distance of 930 AU and an eccentricity from 0.7 to 0.9. We find that the numerical results can be reasonably well reproduced using an analytical approach proposed to explain the formation of a spiral structure by secular perturbation of a circumprimary disk by an external bound companion. We also interpret the redness of the disk in the visible reported by Clampin et al.(2003) and show that short-lived grains one order of magnitude smaller than the blow-out size limit are abundant in the disk. The most probable reason for this is that the disk sustains high collisional activity. Finally we conclude that additional processes are required to clear out the disk inside 150 AU and that interactions with planetary companions possibly coupled with the remnant gas disk are likely candidates.

In this paper we develop a theory of disturbances induced by the stellar tidal field in a fully convective slowly rotating planet orbiting on a highly eccentric orbit around a central star. We show that there are two contributions to the mode energy and angular momentum gain due to impulsive tidal interaction: a) 'the quasi-static' contribution which requires dissipative processes operating in the planet; b) the dynamical contribution associated with excitation of modes of oscillation. These contributions are obtained self-consistently from a single set of the governing equations. We calculate a critical 'equilibrium' value of angular velocity of the planet \Omega_{crit} determined by the condition that action of the dynamical tides does not alter the angular velocity at that rotation rate. We show that this can be much larger than the corresponding rate associated with quasi-static tides and that at this angular velocity, the rate of energy exchange is minimised. We also investigate the conditions for the stochastic increase in oscillation energy that may occur if many periastron passages are considered. We make some simple estimates of time scale of circularization of initially eccentric orbit due to tides, using a realistic model of the planet, for orbits withperiods after circularization typical of those observed for extra-solar planets P_{obs} > 3days. We find that dynamic tides could have produced a very large decrease of the semi-major axis of a planet with mass of the order of the Jupiter mass M_{J} and final periods P_{obs} < 4.5days on a time-scale < a few Gyrs. We also discuss several unresolved issues in the context of the scenario of the orbit circularization due to dynamic tides.

(Abridged) We present global disc and local shearing box simulations of planets interacting with a MHD turbulent disc. We examine the torque exerted by the disc on the embedded planets as a function of planet mass, and thus make a first study of orbital migration of planets due to interaction with turbulent discs. Global simulations were performed for a disc with H/R=0.07 and planet masses M_p=3,10,30 Earth masses, and 3 Jupiter masses. Shearing box runs were performed for values of (M_p/M_*)/(H/R)^3=0.1,0.3,1.0 and 2.0, M_* being the central mass. These allow embedded and gap forming planets to be examined. In all cases the instantaneous torque exerted on a planet showed strong fluctuations. In the embedded cases it oscillated between negative and positive values, and migration occurs as a random walk, unlike the usual type I migration. Running time averages for embedded planets over 20-25 orbital periods show that large fluctuations occur on longer time scales, preventing convergence of the average torque to well defined values, or even to a well defined sign. Fluctuations become relatively smaller for larger masses, giving better convergence, due to the planet's perturbation of the disc becoming larger than the turbulence in its vicinity. Eventually gap formation occurs, with a transition to type II migration. The existence of significant fluctuations occurring in turbulent discs on long time scales is important for lower mass embedded protoplanets. If significant fluctuations occur on the longest disc evolutionary time scales, convergence of torque running averages for practical purposes will not occur, and the migration behaviour of low mass protoplanets considered as an ensemble would be very different from predictions of type I theory for laminar discs.

[Abridged] We evaluate the coorbital corotation torque on a migrating protoplanet. The coorbital torque is assumed to come from orbit crossing fluid elements which exchange angular momentum with the planet when they execute a U-turn at the end of horseshoe streamlines. When the planet migrates inward, the fluid elements of the inner disk undergo one such exchange as they pass to the outer disk. The angular momentum they gain is removed from the planet, and this corresponds to a negative contribution to the corotation torque, which scales with the drift rate. In addition, the material trapped in the coorbital region drifts radially with the planet giving a positive contribution to the corotation torque, which also scales with the drift rate. These two contributions do not cancel out if the coorbital region is depleted, in which case there is a net corotation torque which scales with the drift rate and the mass deficit in the coorbital region, and which has same sign as the drift rate. This leads to a positive feedback on the migrating planet. In particular, if the coorbital mass deficit is larger than the planet mass, the migration rate undergoes a runaway which can vary the protoplanet semi-major axis by 50% over a few tens of orbits. This can happen only if the planet mass is sufficient to create a dip or gap in its surrounding region, and if the surrounding disk mass is larger than the planet mass. This typically corresponds to planet masses in the sub-Saturnian to Jovian mass range embedded in massive protoplanetary disks. Runaway migration is a good candidate to account for the orbital characteristics of close orbiting giant planets, most of which have sub-Jovian masses. Further, we show that in the runaway regime, migration can be directed outwards.

We present a global MHD simulation of a turbulent accretion disc interacting with a protoplanet of 5 Jupiter masses. The disc model had H/r=0.1,and a value of the Shakura & Sunyaev alpha ~ 0.005. The protoplanet opened a gap in the disc, with the interaction leading to inward migration on the expected time scale. Spiral waves were launched by the protoplanet and although they were diffused and dissipated through interaction with the turbulence, they produced an outward angular momentum flow which compensated for a reduced flux associated with the turbulence, so maintaining the gap. When compared with laminar disc models with the same estimated alpha, the gap was found to be deeper and wider indicating that the turbulent disc behaved as if it possessed a smaller alpha. This may arise for two reasons. First, the turbulence does not provide a source of constantly acting friction in the near vicinity of the planet that leads to steady mass flow into the gap region. Instead the turbulence is characterised by large fluctuations in the radial velocity, and time averaging over significant time scales is required to recover the underlying mass flow through the disc. Near the planet the disc material experiences high amplitude perturbations on time scales that are short relative to the time scale required for averaging. The disc response is thus likely to be altered relative to a Navier--Stokes model. Second, the simulation indicates that an ordered magnetic connection between the inner and outer disc can occur enabling angular momentum to flow out across the gap, helping to maintain it independently of the protoplanet's tide. This type of effect may assist gap formation for smaller mass protoplanets which otherwise would not be able to maintain them.

This is the first of a series of papers aimed at developing and interpreting simulations of protoplanets interacting with turbulent accretion discs. Here we study the disc models prior to the introduction of a protoplanet.We study models in which a Keplerian domain is unstable to the magnetorotational instability (MRI). Various models with B-fields having zero net flux are considered.We relate the properties of the models to classical viscous disc theory.All models attain a turbulent state with volume averaged stress parameter alpha ~ 0.005. At any particular time the vertically and azimuthally averaged value exhibited large fluctuations in radius. Time averaging over periods exceeding 3 orbital periods at the outer boundary of the disc resulted in a smoother quantity with radial variations within a factor of two or so. The vertically and azimuthally averaged radial velocity showed much larger spatial and temporal fluctuations, requiring additional time averaging for 7-8 orbital periods at the outer boundary to limit them. Comparison with the value derived from the averaged stress using viscous disc theory yielded schematic agreement for feasible averaging times but with some indication that the effects of residual fluctuations remained. The behaviour described above must be borne in mind when considering laminar disc simulations with anomalous Navier--Stokes viscosity. This is because the operation of a viscosity as in classical viscous disc theory with anomalous viscosity coefficient cannot apply to a turbulent disc undergoing rapid changes due to external perturbation. The classical theory can only be used to describe the time averaged behaviour of the parts of the disc that are in a statistically steady condition for long enough for appropriate averaging to be carried out.

We calculate global $m=1$ modes with low pattern speed corresponding to introducing a finite eccentricity into a protoplanetary disc. We consider disc models which are either isolated or contain one or two protoplanets orbiting in an inner cavity. Global modes that are strongly coupled to inner protoplanets are found to have disc orbits which tend to have apsidal lines antialigned with respect to those of the inner protoplanets. Other modes corresponding to free disc modes may be global over a large range of length scales and accordingly be long lived. We consider the motion of a protoplanet in the earth mass range embedded in an eccentric disc and determine the equilibrium orbits which maintain fixed apsidal alignment with respect to the disc gas orbits. Equilibrium eccentricities are found to be comparable or possibly exceed the disc eccentricity. We then approximately calculate the tidal interaction with the disc in order to estimate the orbital migration rate. Results are found to deviate from the case of axisymmetric disc with near circular protoplanet orbit once eccentricities of protoplanet and disc orbits become comparable to the disc aspect ratio in magnitude. Aligned protoplanet orbits with very similar eccentricity to that of the gas disc are found to undergo litle eccentricity change while undergoing inward migration in general. However, for significantly larger orbital eccentricities, migration may be significantly reduced or even reverse from inwards to outwards. Thus the existence of global non circular motions in discs with radial excursions comparable to the semi-thickness may have important consequences for the migration and survival of protoplanetary cores in the earth mass range.

We investigate the possible commensurabilities to be expected when two protoplanets in the Jovian mass range, gravitationally interacting with each other and an external protoplanetary disc, are driven by disc induced orbital migration of the outer protoplanet into a commensurability which is then maintained in subsequent evolution. We find that for a variety of protoplanet masses and typical protoplanetary disc properties, as well as the setting up of 2:1 commensurabilities of the type recently observed in GJ876, 3:1 commensurabilities are often formed, in addition to 4:1, 5:1, and 5:2 commensurabilities which occur less frequently. The higher order commensurabilites are favoured when either one of the planets is massive, or the inner planet begins with a significant orbital eccentricity. Detection of such commensurabilities could yield important information relating to the operation of protoplanet disc interactions during and shortly post formation.

We consider the evolution of a system containing a population of massive planets formed rapidly through a fragmentation process occurring on a scale on the order of 100 au and a lower mass planet that assembles in a disc on a much longer timescale. During the formation phase, the inner planet is kept on a circular orbit due to tidal interaction with the disc, while the outer planets undergo dynamical relaxation. Interaction with the massive planets left in the system after the inner planet forms may increase the eccentricity of the inner orbit to high values, producing systems similar to those observed.

We present global 3D MHD simulations of geometrically thin but unstratified accretion disks in which a near Keplerian disk rotates between two bounding regions with initial rotation profiles that are stable to the MRI. The inner region models the boundary layer between the disk and an assumed more slowly rotating central, non magnetic star. We investigate the dynamical evolution of this system in response to initial vertical and toroidal fields imposed in a variety of domains contained within the near Keplerian disk. Cases with both non zero and zero net magnetic flux are considered and sustained dynamo activity found in runs for up to fifty orbital periods at the outer boundary of the near Keplerian disk. Simulations starting from fields with small radial scale and with zero net flux lead to the lowest levels of turbulence and smoothest variation of disk mean state variables. For our computational set up, average values of the Shakura & Sunyaev (1973) $\alpha$ parameter in the Keplerian disk are typically $0.004\pm 0.002.$ Magnetic field eventually always diffuses into the boundary layer resulting in the build up of toroidal field inward angular momentum transport and the accretion of disk material. The mean radial velocity, while exhibiting large temporal fluctuations is always subsonic. Simulations starting with net toroidal flux may yield an average $\alpha \sim 0.04.$ While being characterized by one order of magnitude larger average $\alpha$, simulations starting from vertical fields with large radial scale and net flux may lead to the formation of persistent non-homogeneous, non-axisymmetric magnetically dominated regions of very low density.

We consider two protoplanets gravitationally interacting with each other and a protoplanetary disc. The two planets orbit interior to a tidally maintained disc cavity while the disc interaction indices inward migration. When the migration is slow enough, the more rapidly migrating outer protoplanet approaches and becomes locked in a 2:1 commensurability with the inner one. This is maintained in subsequent evolution. We study this evolution using a simple anaytic model, full hydrodynamic 2D simulations of the disc planet system and longer time N body integrations incorporating simple prescriptions for the effect of the disc on the planet orbits. The eccentricity of the protoplanets are found to be determined by the migration rate induced in the outer planet orbit by the external disc. We apply our results to the recently discovered resonant planets around GJ876. Simulation shows that a disc with parameters expected for protoplanetary discs causes trapping in the 2:1 commensurability when the planets orbit in an inner cavity and that eccentricities in the observed range may be obtained.

We report a new and complete model of the beta Pictoris disk, which succeeds in accounting for both the surface brightness distribution, warp characteristics, the outer ``butterfly'' asymmetry as observed by HST/STIS in scattered light, as well as the infrared emission. Our model includes the presence of a disk of planetesimals extending out to 120-150 AU, perturbed gravitationally by a giant planet on an inclined orbit, following the approach of Mouillet et al. (1997b). At any time, the planetesimal disk is assumed to be the source of a distribution of grains produced through collisional evolution, with the same initial orbital parameter distribution. The steady state spatial grain distribution is found incorporating the effects of radiation pressure which can cause the distribution of the smallest particles to become very distended. With realistic assumptions about the grains' chemical properties, the modeling confirms the previously evident need for an additional population of hot grains close to the star, to account for the 12 microns fluxes at short distances from the star. It also indicates that this population cannot explain the outer 12 microns flux distribution when the effects of gravity and radiation pressure determine the distribution. Very small grains, produced by collisions among aggregates, are tentatively proposed to account for this 12 microns outer emission.

Following the suggestion of Black (1997) that some massive extrasolar planets may be associated with the tail of the distribution of stellar companions, we investigate a scenario in which 5 < N < 100 planetary mass objects are assumed to form rapidly through a fragmentation process occuring in a disc or protostellar envelope on a scale of 100 au. These are assumed to have formed rapidly enough through gravitational instability or fragmentation that their orbits can undergo dynamical relaxation on a timescale of about 100 orbits. Under a wide range of initial conditions and assumptions the relaxation process ends with either (i) one potential 'hot Jupiter' plus up to two 'external' companions, i.e. planets orbiting near the outer edge of the initial distribution; (ii) one or two 'external' planets or even none at all; (iii) one planet on an orbit with a semi--major axis 10 to a 100 times smaller than the outer boundary radius of the inital distribution together with an 'external' companion. Most of the other objects are ejected and could contribute to a population of free floating planets. Apart from the potential 'hot Jupiters', all the bound objects are on highly eccentric orbits. We found that, apart from the close orbiters, the probability of ending up with a planet orbiting at a given distance from the central star increases with the distance. This is because of the tendency of the relaxation process to lead to collisions with the central star. We discuss the application of these results to some of the more massive extrasolar planets.

We compute the warping of a disc induced by an inclined dipole. We consider a magnetised star surrounded by a thin Keplerian diamagnetic disc with an inner edge that corotates with the star. We suppose the stellar field is a dipole with an axis that is slightly misaligned with the stellar rotation axis. The rotation axes of the disc material orbiting at large distances from the star and that of the star are supposed to coincide. The misalignment of the magnetic and rotation axes results in the magnetic pressure not being the same on the upper and lower surfaces of the disc. The resultant net vertical force produces a warp which appears stationary in a frame corotating with the star. We find that, if viscosity is large enough (alpha on the order of 0.01-0.1) to damp bending waves as they propagate away, a smoothly varying warp of the inner region of the disc is produced. The amplitude of the warp can easily be on the order of ten percent of the disc inner radius for reasonably small misalignment angles (less than 30 degrees). Viscous damping also introduces a phase shift between the warp and the forcing torque, which results in the locations of maximum elevation above the disc forming a trailing spiral pattern. We apply these results to recent observations of AA Tau, and show that the variability of its light curve, which occurs with a period comparable to the expected stellar rotation period, could be due to obscuration produced by a warp configuration of the type we obtain.

We calculate the structure of a force-free magnetosphere which is assumed to corotate with a central star and which interacts with an embedded differentially rotating accretion disc. The magnetic and rotation axes are aligned and the stellar field is assumed to be a dipole. We concentrate on the case when the amount of field line twisting through the disc-magnetosphere interaction is large and consider different outer boundary conditions. In general the field line twisting produces field line inflation (eg. Bardou & Heyvaerts 1996) and in some cases with large twisting many field lines can become open. We calculate the spin-down torque acting between the star and the disc and we find that it decreases significantly for cases with large field line twisting. This suggests that the oscillating torques observed for some accreting neutron stars could be due to the magnetosphere varying between states with low and high field line inflation. Calculations of the spin evolution of T Tauri stars may also have to be revised in light of the significant effect that field line twisting has on the magnetic torque resulting from star-disc interactions.

We present SPH simulations of accretion discs in orbit about rotating compact objects such as black holes and neutron stars, and study the structure of warped discs produced by the Bardeen-Petterson effect. We calculate the transition radius out to which the disc specific angular momentum vector is aligned with that of the black hole. We focus on the parameter regime where the warp dynamics are controlled by bending wave propagation, but also consider models in which warps are subject to diffusion rather than wave transport, and are able to consider the fully nonlinear regime. Because of hydrodynamic or pressure effects, for the parameter range investigated, the transition radius is always found to be much smaller than that obtained by Bardeen & Petterson (1975). For discs with midplane Mach numbers of about 10, the transition occurs between 10 - 16 gravitational radii, whereas for a Mach number of about 30 it occurs at around 30 gravitational radii. A thicker disc with a Mach number of 5 is found to produce no discernible warped structure. The rate of black hole - disc alignment is found to be consistent with the ideas of Ress (1978), with the alignment torque behaving as if it arises from the accreted material transferring its misaligned component of angular momentum at the larger transition radius of Bardeen & Petterson (1975). The inclusion of Einstein precession in the calculations modified both the warped disc structure and, consistent with linear analysis, produced an increased alignment rate by up to a factor of 4 because of the effect that a non Keplerian potential has on the propagation of warps.

We present a new computation of the linear tidal interaction of a protoplanetary core with a thin gaseous disc in which it is fully embedded. For the first time a discussion of the orbital evolution of cores with eccentricity (e) significantly larger than the gas-disc scale height to radius ratio (H/r) is given. We find that the direction of orbital migration reverses for e>1.1H/r. This occurs as a result of the orbital crossing of resonances in the disc that do not overlap the orbit when the eccentricity is very small. Simple expressions giving approximate fits to the eccentricity damping rate and the orbital migration rate are presented. We go on to calculate the rate of increase of the mean eccentricity for a system of protoplanetary cores due to dynamical relaxation. By equating the eccentricity damping time-scale with the dynamical relaxation time-scale we deduce that an equilibrium between eccentricity damping and excitation through scattering is attained on a 10^3 to 10^4 yr time-scale, at 1au. The equilibrium thickness of the protoplanet distribution is such that it is generally well confined within the gas disc. By use of a three dimensional N-body code we simulate the evolution of a system of protoplanetary cores, incorporating our eccentricity damping and migration rates. Assuming that collisions lead to agglomeration, we find that the vertical confinement of the protoplanet distribution permits cores to build up from 0.1 to 1 earth mass in only ~10^4 yr, within 1au. The time-scale required to achieve this is comparable to the migration time-scale. We deduce that it is not possible to build up a massive enough core to form a gas giant planet before orbital migration ultimately results in the preferential delivery of all such bodies to the neighbourhood of the central star. [Abridged]

We investigate the gravitational interaction of a Jovian mass protoplanet with a gaseous disc with aspect ratio and kinematic viscosity expected for the protoplanetary disc from which it formed. Different disc surface density distributions have been investigated. We focus on the tidal interaction with the disc with the consequent gap formation and orbital migration of the protoplanet. Nonlinear hydrodynamic simulations are employed using three independent numerical codes. A principal result is that the direction of the orbital migration is always inwards and such that the protoplanet reaches the central star in a near circular orbit after a characteristic viscous time scale of approximately 10,000 initial orbital periods. This was found to be independent of whether the protoplanet was allowed to accrete mass or not. Inward migration is helped through the disappearance of the inner disc, and therefore the positive torque it would exert, because of accretion onto the central star.Our results indicate that a realistic upper limit for the masses of closely orbiting giant planets is approximately 5 Jupiter masses, because of the reduced accretion rates obtained for planets of increasing mass. Assuming some process such as termination of the inner disc through a magnetospheric cavity stops the migration, the range of masses estimated for a number of close orbiting giant planets (Marcy, Cochran, & Mayor 1999; Marcy & Butler 1998) as well as their inward orbital migration can be accounted for by consideration of disc--protoplanet interactions during the late stages of giant planet formation. Maximally accreting protoplanets reached about four Jovian masses on reaching the neighbourhood of the central star.

We present the results of a study of propagating warp or bending waves in accretion discs. Three dimensional hydrodynamic simulations were performed using SPH, and the results of these are compared with calculations based on the linear theory of warped discs. We consider primarily the physical regime in which the dimensionless viscosity parameter `alpha' < H/r, the disc aspect ratio, so that bending waves are expected to propagate. We also present calculations in which `alpha' > H/r, where the warps are expected to behave diffusively. Small amplitude perturbations are studied in both Keplerian and slightly non Keplerian discs, and we find that the SPH results can be reasonably well fitted by those of the linear theory. The main results of these calculations are: (1) the warp in Keplerian discs when `alpha' < H/r propagates with little dispersion and damps at a rate expected from estimates of the code viscosity, (2) warps evolve diffusively when `alpha' > H/r, (3) the non Keplerian discs exhibit a substantially more dispersive behaviour of the warps. Initially imposed higher amplitude nonlinear warping disturbances were studied in Keplerian discs. The results indicate that nonlinear warps can lead to the formation of shocks, and that the evolution of the warp becomes less wave-like and more diffusive in character. This work is relevant to the study of the warped accretion discs that may occur around Kerr black holes or in misaligned binary systems. The results indicate that SPH can accurately model the hydrodynamics of warped discs, even when using rather modest numbers of particles.

We investigate in details the properties of the disk surrounding the 8 Myr old star HR 4796A, one of the few stars bringing precious clues to better understand the scenario which lead to planetary system formation. We propose a model able to reproduce all the available observations : the full spectral energy distribution from the mid-infrared to the millimeter wavelengths, resolved scattered light and thermal emission observations. We show that the circumstellar matter splits into two dust components : a cold annulus, peaked at 70 AU from the star, made of ISM-like grains (amorphous composition, porosity ~ 0.6) larger than 10 microns and a population of hot dust close to the star (at about 9 AU) made of comet-like grains (crystalline composition, porosity ~ 0.97). Both dust populations are highly collisional and the grain size distribution in the cold annulus is found to be cut-off by radiation pressure. At 70 AU, bodies as large as a few meters are required to fit the data leading to a minimum disk mass of a few Earth masses and to a gas to dust ratio less than 1. We discuss aftewards some implications on the disk structure and effects of larger bodies.

We review results about protoplanetary disk models, protoplanet migration and formation of giant planets with migrating cores. We first model the protoplanetary nebula as an \alpha-accretion disk and present steady state calculations for different values of \alpha and gas accretion rate through the disk. We then review the current theories of protoplanet migration in the context of these models, focusing on the gaseous disk-protoplanet tidal interaction. According to these theories, the migration timescale may be shorter than the planetary formation timescale. Therefore we investigate planet formation in the context of a migrating core, considering both the growth of the core and the build-up of the envelope in the course of the migration.

We study a solid protoplanetary core of 1-10 earth masses migrating through a disk. We suppose the core luminosity is generated as a result of planetesimal accretion and calculate the structure of the gaseous envelope assuming equilibrium. This is a good approximation when the core mass is less than the critical value, M_{crit}, above which rapid gas accretion begins. We model the structure of the protoplanetary nebula as an accretion disk with constant \alpha. We present analytic fits for the steady state relation between disk surface density and mass accretion rate as a function of radius r. We calculate M_{crit} as a function of r, gas accretion rate through the disk, and planetesimal accretion rate onto the core \dot{M}. For a fixed \dot{M}, M_{crit} increases inwards, and it decreases with \dot{M}. We find that \dot{M} onto cores migrating inwards in a time 10^3-10^5 yr at 1 AU is sufficient to prevent the attainment of M_{crit} during the migration process. Only at small radii where planetesimals no longer exist can M_{crit} be attained. At small radii, the runaway gas accretion phase may become longer than the disk lifetime if the core mass is too small. However, massive cores can be built-up through the merger of additional incoming cores on a timescale shorter than for in situ formation. Therefore, feeding zone depletion in the neighborhood of a fixed orbit may be avoided. Accordingly, we suggest that giant planets may begin to form early in the life of the protostellar disk at small radii, on a timescale that may be significantly shorter than for in situ formation. (abridged)

The dynamical foundations of $\alpha$ disk models are described. At the heart of the viscous formalism of accretion disk models are correlations in the fluctuating components of the disk velocity, magnetic field, and gravitational potential. We relate these correlations to the large scale mean flow dynamics used in phenomenological viscous disk models. MHD turbulence readily lends itself to the $\alpha$ formalsim, but transport by self-gravity does not. Nonlocal transport is an intrinsic property of turbulent self-gravitating disks, which in general cannot be captured by an $\alpha$ model. Local energy dissipation and $\alpha$-like behavior can be re-established if the pattern speeds associated with the amplitudes of an azimuthal Fourier decomposition of the turbulence are everywhere close to the local rotation frequency. In this situation, global wave transport must be absent. Shearing box simulations, which employ boundary conditions forcing local behavior, are probably not an adequate tool for modeling the behavior of self-gravitating disks. As a matter of principle, it is possible that disks which hover near the edge of gravitational stability may behave in accord with a local $\alpha$ model, but global simulations performed to date suggest matters are not this simple.

The interaction of a massive binary and a non-self-gravitating circumbinary accretion disc is considered. The shape of the stationary twisted disc produced by the binary is calculated. It is shown that the inner part of the disc must lie in the binary orbital plane for any value of viscosity. When the inner disc midplane is aligned with the binary orbital plane on the scales of interest and it rotates in the same sense as the binary, the modification of the disc structure and the rate of decay of the binary orbit, assumed circular, due to tidal exchange of angular momentum with the disc, are calculated. It is shown that the modified disc structure is well described by a self-similar solution of the non-linear diffusion equation governing the evolution of the disc surface density. The calculated time scale for decay of the binary orbit is always smaller than the "accretion" time $t_{acc}=m/{\dot M}$ ($m$ is the mass of the secondary component, and $\dot M$ is the disc accretion rate), and is determined by ratio of secondary mass $m$, assumed to be much smaller than the primary mass, the disc mass inside the initial binary orbit, and the form of viscosity in the disc.

We consider several protostellar systems where either a precessing jet or at least two misaligned jets have been observed. We assume that the precession of jets is caused by tidal interactions in noncoplanar binary systems. For Cep E, V1331 Cyg and RNO 15-FIR the inferred orbital separations and disk radii are in the range 4-160 AU and 1-80 AU, respectively, consistent with those expected for pre-main sequence stars. Furthermore, we assume or use the fact that the source of misaligned outflows is a binary, and evaluate the lengthscale over which the jets should precess as a result of tidal interactions. For T Tau, HH1 VLA 1/2 and HH 24 SVS63, it may be possible to detect a bending of the jets rather than 'wiggling'. In HH 111 IRS and L1551 IRS5, 'wiggling' may be detected on the current observed scale. Our results are consistent with the existence of noncoplanar binary systems in which tidal interactions induce jets to precess.

We review results on the dynamics of warped gaseous discs. We consider tidal perturbation of a Keplerian disc by a companion star orbiting in a plane inclined to the disc. The perturbation induces the precession of the disc, and thus of any jet it could drive. In some conditions the precession rate is uniform, and as a result the disc settles into a warp mode. The tidal torque also leads to the truncation of the disc, to the evolution of the inclination angle (not necessarily towards alignment of the disc and orbital planes) and to a transport of angular momentum in the disc. We note that the spectral energy distribution of such a warped disc is different from that of a flat disc. We conclude by listing observational effects of warps in protostellar discs.

The formation, structure and evolution of protoplanetary discs is considered. The formation of giant planets within the environment of these models is also discussed.

We consider several processes operating during the late stages of planet formation that can affect observed orbital elements. Disk-planet interactions, tidal interactions with the central star, long term orbital instability and the Kozai mechanism are discussed.

Potential-density pair basis sets can be used for highly efficient N-body simulation codes, but they suffer from a lack of versatility, i.e. a basis set has to be constructed for each different class of stellar system. We present numerical techniques for generating a biorthonormal potential-density pair basis set which has a general specified pair as its lowest order member. We go on to demonstrate how the set can be used to construct N-body equilibria, which we then evolve using an N-body code which calculates forces using the basis set.

We calculate the dynamical tides raised on a non-rotating solar-type star by a close stellar or planetary companion. Dissipation arising from a turbulent viscosity operating in the convection zone and radiative damping in the radiative core are considered. We compute the torque exerted on the star by a companion in circular orbit, and determine the potentially observable magnitude of the tidally induced velocity at the stellar photosphere. These calculations are compared with the results obtained by assuming that a very small frequency limit can be taken in order to calculate the tidal response (equilibrium tide). For a standard solar model, the latter is found to give a relatively poor approximation at the periods of interest of several days, even when the system is far from resonance with a normal mode. It is shown that although the companion may go through a succession of resonances as it spirals in under the action of the tides, for a fixed spectrum of normal modes its migration is controlled essentially by the non-resonant interaction. We find that the turbulent viscosity that is required to provide the observed circularization rates of main sequence solar-type binaries is about fifty times larger than that simply estimated from mixing length theory for non-rotating stars. We discuss the means by which this enhanced viscosity might be realized. These calculations are applied to 51 Pegasi. We show that the perturbed velocity induced by the tides at the stellar surface is too small to be observed.

We examine the effect of dynamical tides raised by a companion on a solar-type star. In these binaries, gravity or g mode oscillations are excited by the companion in the radiative region beneath the convective envelope of the star. They become evanescent in the convection zone. This is of particular interest in connection with the newly discovered planets, some of which are found to orbit around solar-type stars with a period comparable to that of the high order g modes of the star. One such example is 51 Pegasi. Here, we determine the magnitude of the perturbed velocity induced by the tides at the stellar surface. We show that, in the case of 51 Pegasi, this velocity is too small to be observed. This result is insensitive to the magnitude of the stellar turbulent viscosity assumed and is not affected by the possibility of resonance, which occurs when the frequency of the tidal disturbance is close to that of some normal mode of the star. We also discuss the orbital evolution and synchronization timescales associated with the tidal interaction.

On the global warping of a thin self-gravitating near Keplerian gaseous disk with application to the disk in NGC 4258

We study the global bending modes of a thin annular disc subject to both an internally generated magnetic field and a magnetic field due to a dipole embedded in the central star with axis aligned with the disc rotation axis. When there is a significant inner region of the disc corotating with the star, we find spectra of unstable bending modes. These may lead to elevation of the disc above the original symmetry plane facilitating accretion along the magnetospheric field lines. The resulting non-axisymmetric disc configuration may result in the creation of hot spots on the stellar surface and the periodic photometric variations observed in many classical T Tauri stars (CTTS). Time-dependent behaviour may occur including the shadowing of the central source in magnetic accretors even when the dipole and rotation axes are aligned.

We study the fully non-adiabatic tidal response of a uniformly rotating unevolved 20 Msun star to the dominant l=m=2 component of the companion's perturbing potential. This is done numerically with a 2D implicit finite difference scheme. We assume the star is rotating slowly, so that the centrifugal force can be neglected, but the Coriolis force is taken fully into account. We study the low frequency `inertial' regime where the tidal oscillation frequency in the frame corotating with the star is less than twice the stellar spin rate. In this frequency range inertial modes are excited in the convective core which can interact with rotationally modified g- or r-modes in the radiative envelope and cause significant strengthening of the tidal interaction. Resonant interaction with quasi-toroidal (r-)modes in slightly super-synchronous stars causes efficient spin down towards corotation.

We describe a generalization of the asymptotic calculation of the tidal torques experienced by a massive star as a result of a companion in circular orbit originally considered by Zahn (1975,1977) to the case of a rotating star when the forcing frequency is small and in the inertial regime, that is it is less than twice the rotation frequency in magnitude. The results confirm the presence of a strong toroidal mode resonance feature for retrograde forcing and also, with a simple description of the convective core, the presence of some core inertial mode features in the response. These were found numerically by Savonije and Papaloizou (astro-ph/9706186).

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.

We describe some recent nonlinear three dimensional hydrodynamic simulations of accretion discs in binary systems where the orbit is circular and not necessarily coplanar with the disc midplane. 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.

The structure of the boundary layer region between the disc and a comparatively slowly rotating star is studied using a causal prescription for viscosity. The vertically integrated viscous stress relaxes towards its equilibrium value on a relaxation timescale $\tau$, which naturally yields a finite speed of propagation for viscous information. For a standard alpha prescription with alpha in the range 0.1-0.01, and ratio of viscous speed to sound speed in the range 0.02-0.5, details in the boundary layer are strongly affected by the causality constraint. We study both steady state polytropic models and time dependent models, taking into account energy dissipation and transport. Steady state solutions are always subviscous with a variety of $\Omega$ profiles which may exhibit near discontinuities. For alpha =0.01 and small viscous speeds, the boundary layer adjusted to a near steady state. A long wavelength oscillation generated by viscous overstability could be seen at times near the outer boundary. Being confined there, the boundary layer remained almost stationary. However, for alpha =0.1 and large viscous speeds, short wavelength disturbances were seen throughout which could significantly affect the power output in the boundary layer. This could be potentially important in producing time dependent behaviour in accreting systems such as CVs and protostars.

In this paper we present an analytical and numerical study of the response of a circumbinary disc subject to the tidal-forcing of a binary with a fixed circular orbit. We consider fluid discs with a range of thicknesses and binaries with a range of mass ratios, orbital separations and inclination angles. Our numerical simulations are implemented using a SPH code. For our unperturbed disc models, we write a scaling relation for the shear viscosity and deduce that the disc thickness cannot be varied without affecting the viscosity in these kinds of SPH disc models. It is found that maintainance of an inner cavity owing to the tidal truncation of the disc is effective for non-zero orbital inclinations. Also we show that our model discs may precess approximately like rigid bodies, provided that the disc is able to communicate on a length scale comparable to the inner boundary radius by either sonic or viscous effects, in a sufficiently small fraction of the local precession period. Furthermore, the disc precession period may tend to infinity if the disc outer edge is allowed to become arbitrarily large, the disc suffering only a modest quasi-steady warp near the inner boundary. When the disc response is linear, or weakly non-linear, the precession periods and the forms of warping that we measure yield reasonable quantitative agreement with the analytical expressions that we derive from a linear response calculation. For a stronger disc response the results can agree poorly with our linear analysis, although some qualitative features of the response remain intact. We show that the response of a disc of non-interacting particles is qualitatively different from this. The work presented in this paper is of relevance to a number of astrophysical phenomena of current interest in star and planet formation.

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.

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.

We present a new formulation of the equations of motion used in smoothed particle hydrodynamics (SPH). The spatial resolution in SPH is determined by the smoothing length, $h$, and it has become common practice for each particle to be given its own adaptive smoothing length, $h_i$. Consequently, the dynamic range that may be spatially resolved is greatly increased, but additional ($\nabla h$) terms, which account for the variability of the smoothing lengths, should appear in the equations of motion in order to satisy conservation requirements. Previous implementations of SPH have neglected these additional terms, whereas we have included them. This is achieved by defining a functional form for the $h_i$s, that depends on inter-particle distances only, and then deriving the equations of motion using a Hamiltonian formalism. A number of test calculations are presented. We find that including the $\nabla h$ terms has no detrimental effect on the ability of SPH to model known problems with reasonable accuracy. For problems where the energy conservation was rather poor, including the $\nabla h$ terms results in a dramatic improvement. In particluar non-conservation of energy during a collision between two polytropes can occur at the level of 10% when the $\nabla h$ terms are neglected. When the $\nabla h$ terms are included, then this error reduces to 0.8%.

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 '$\alpha $' 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 $\alpha $ between 0.1 and 10$^{-4}$.