The nearest exoplanets to the Sun are our best possibilities for detailed characterization. We report the discovery of a compact multi-planet system of super-Earths orbiting the nearby red dwarf GJ 887, using radial velocity measurements. The planets have orbital periods of 9.3 and 21.8~days. Assuming an Earth-like albedo, the equilibrium temperature of the 21.8 day planet is approx 350 K; which is interior, but close to the inner edge, of the liquid-water habitable zone. We also detect a further unconfirmed signal with a period of 50 days which could correspond to a third super-Earth in a more temperate orbit. GJ 887 is an unusually magnetically quiet red dwarf with a photometric variability below 500 parts-per-million, making its planets amenable to phase-resolved photometric characterization.

Spin off events and impacts can eject boulders from an asteroid surface and rubble pile asteroids can accumulate from debris following a collision between large asteroids. These processes produce a population of gravitational bound objects in orbit that can impact an asteroid surface at low velocity and with a distribution of impact angles. We present laboratory experiments of low velocity spherical projectiles into a fine granular medium, sand. We delineate velocity and impact angles giving ricochets, those giving projectiles that roll-out from the impact crater and those that stop within their impact crater. With high speed camera images and fluorescent markers on the projectiles we track spin and projectile trajectories during impact. We find that the projectile only reaches a rolling without slipping condition well after the marble has reached peak penetration depth. The required friction coefficient during the penetration phase of impact is 4-5 times lower than that of the sand suggesting that the sand is fluidized near the projectile surface during penetration. We find that the critical grazing impact critical angle dividing ricochets from roll-outs, increases with increasing impact velocity. The critical angles for ricochet and for roll-out as a function of velocity can be matched by an empirical model during the rebound phase that balances a lift force against gravity. We estimate constraints on projectile radius, velocity and impact angle that would allow projectiles on asteroids to ricochet or roll away from impact, finally coming to rest distant from their initial impact sites.

Eccentric gaseous discs are unstable to a parametric instability involving the resonant interaction between inertial-gravity waves and the eccentric mode in the disc. We present 3D global hydrodynamical simulations of inviscid circumbinary discs that form an inner cavity and become eccentric through interaction with the central binary. The parametric instability grows and generates turbulence that transports angular momentum with stress parameter $\alpha \sim 5 \times 10^{-3}$ at distances $\lesssim 7 \;a_{bin} $, where $a_{bin}$ is the binary semi-major axis. Vertical turbulent diffusion occurs at a rate corresponding to $\alpha_{diff}\sim 1-2\times 10^{-3}$. We examine the impact of turbulent diffusion on the vertical settling of pebbles, and on the rate of pebble accretion by embedded planets. In steady state, dust particles with Stokes numbers ${\it St} \lesssim 0.1$ form a layer of finite thickness $H_d \gtrsim 0.1 H$, where $H$ is the gas scale height. Pebble accretion efficiency is then reduced by a factor $r_{acc}/H_d$, where $r_{acc}$ is the accretion radius, compared to the rate in a laminar disc. For accreting core masses with $m_p \lesssim 0.1\; M_\oplus$, pebble accretion for particles with ${\it St} \gtrsim 0.5$ is also reduced because of velocity kicks induced by the turbulence. These effects combine to make the time needed by a Ceres-mass object to grow to the pebble isolation mass, when significant gas accretion can occur, longer than typical disc lifetimes. Hence, the origins of circumbinary planets orbiting close to their central binary systems, as discovered by the Kepler mission, are difficult to explain using an in situ model that invokes a combination of the streaming instability and pebble accretion.

Theoretical models of the ionization state in protoplanetary disks suggest the existence of large areas with low ionization and weak coupling between the gas and magnetic fields. In this regime hydrodynamical instabilities may become important. In this work we investigate the gas and dust structure and dynamics for a typical T Tauri system under the influence of the vertical shear instability (VSI). We use global 3D radiation hydrodynamics simulations covering all $360^\circ$ of azimuth with embedded particles of 0.1 and 1mm size, evolved for 400 orbits. Stellar irradiation heating is included with opacities for 0.1- to 10-$\mu$m-sized dust. Saturated VSI turbulence produces a stress-to-pressure ratio of $\alpha \simeq 10^{-4}$. The value of $\alpha$ is lowest within 30~au of the star, where thermal relaxation is slower relative to the orbital period and approaches the rate below which VSI is cut off. The rise in $\alpha$ from 20 to 30~au causes a dip in the surface density near 35~au, leading to Rossby wave instability and the generation of a stationary, long-lived vortex spanning about 4~au in radius and 40~au in azimuth. Our results confirm previous findings that mm size grains are strongly vertically mixed by the VSI. The scale height aspect ratio for 1mm grains is determined to be 0.037, much higher than the value $H/r=0.007$ obtained from millimeter-wave observations of the HL~Tau system. The measured aspect ratio is better fit by non-ideal MHD models. In our VSI turbulence model, the mm grains drift radially inwards and many are trapped and concentrated inside the vortex. The turbulence induces a velocity dispersion of $\sim 12$~m/s for the mm grains, indicating that grain-grain collisions could lead to fragmentation.

Outflows driven by large-scale magnetic fields likely play an important role in the evolution and dispersal of protoplanetary disks, and in setting the conditions for planet formation. We extend our 2-D axisymmetric non-ideal MHD model of these outflows by incorporating radiative transfer and simplified thermochemistry, with the twin aims of exploring how heating influences wind launching, and illustrating how such models can be tested through observations of diagnostic spectral lines. Our model disks launch magnetocentrifugal outflows primarily through magnetic tension forces, so the mass-loss rate increases only moderately when thermochemical effects are switched on. For typical field strengths, thermochemical and irradiation heating are more important than magnetic dissipation. We furthermore find that the entrained vertical magnetic flux diffuses out of the disk on secular timescales as a result of non-ideal MHD. Through post-processing line radiative transfer, we demonstrate that spectral line intensities and moment-1 maps of atomic oxygen, the HCN molecule, and other species show potentially observable differences between a model with a magnetically driven outflow and one with a weaker, photoevaporative outflow. In particular, the line shapes and velocity asymmetries in the moment-1 maps could enable the identification of outflows emanating from the disk surface.

We present simulations of low-mass planet--disc interactions in inviscid three-dimensional discs. We show that a wind-driven laminar accretion flow through the surface layers of the disc does not significantly modify the migration torque experienced by embedded planets. More importantly, we find that 3D effects lead to a dramatic change in the behaviour of the dynamical corotation torque compared to earlier 2D theory and simulations. Although it was previously shown that the dynamical corotation torque could act to slow and essentially stall the inward migration of a low-mass planet, our results in 3D show that the dynamical corotation torque has the complete opposite effect and speeds up inward migration. Our numerical experiments implicate buoyancy resonances as the cause. These have two effects: (i) they exert a direct torque on the planet, whose magnitude relative to the Lindblad torque is measured in our simulations to be small; (ii) they torque the gas librating on horseshoe orbits in the corotation region and drive evolution of its vortensity, leading to the negative dynamical corotation torque. This indicates that at low turbulent viscosity, the detailed vertical thermal structure of the protoplanetary disc plays an important role in determining the migration behaviour of embedded planets. If this result holds up under a more refined treatment of disc thermal evolution, then it has important implications for understanding the formation and early evolution of planetary systems.

We consider how an impact generated seismic pulse affects the surface of an asteroid distant from the impact site. With laboratory experiments on dry polydisperse gravel mixtures, we track the trajectories of particles ejected from the surface by a single strong upward propagating pressure pulse. High speed video images show that ejecta trajectories are independent of particle size, and collisions primarily take place upon landing. When they land particles are ballistically sorted, as proposed by Shinbrot et al. (2017), leaving larger particles on the surface and smaller particles more widely dispersed. A single strong pulse can leave previously buried boulders stranded on the surface. Boulder stranding due to an impact excited seismic pulse is an additional mechanism that could leave large boulders present on the surface of rubble asteroids such as 162173 Ryugu, 101955 Bennu and 25143 Itokawa.

NASA's Kepler mission discovered $\sim700$ planets in multi-planet systems containing 3 or more transiting bodies, many of which are super-Earths and mini-Neptunes in compact configurations. Using $N$-body simulations, we examine the in situ, final stage assembly of multi-planet systems via the collisional accretion of protoplanets. Our initial conditions are constructed using a subset of the Kepler 5-planet systems as templates. Two different prescriptions for treating planetary collisions are adopted. The simulations address numerous questions: do the results depend on the accretion prescription?; do the resulting systems resemble the Kepler systems, and do they reproduce the observed distribution of planetary multiplicities when synthetically observed?; do collisions lead to significant modification of protoplanet compositions, or to stripping of gaseous envelopes?; do the eccentricity distributions agree with those inferred for the Kepler planets? We find the accretion prescription is unimportant in determining the outcomes. The final planetary systems look broadly similar to the Kepler templates adopted, but the observed distributions of planetary multiplicities or eccentricities are not reproduced, because scattering does not excite the systems sufficiently. In addition, we find that $\sim 1$% of our final systems contain a co-orbital planet pair in horseshoe or tadpole orbits. Post-processing the collision outcomes suggests they would not significantly change the ice fractions of initially ice-rich protoplanets, but significant stripping of gaseous envelopes appears likely. Hence, it may be difficult to reconcile the observation that many low mass Kepler planets have H/He envelopes with an in situ formation scenario that involves giant impacts after dispersal of the gas disc.

The increasing number of newly detected exoplanets at short orbital periods raises questions about their formation and migration histories. A particular puzzle that requires explanation arises from one of the key results of the Kepler mission, namely the increase in the planetary occurrence rate with orbital period up to 10 days for F, G, K and M stars. We investigate the conditions for planet formation and migration near the dust sublimation front in protostellar disks around young Sun-like stars. For this analysis we use iterative 2D radiation hydrostatic disk models which include irradiation by the star, and dust sublimation and deposition depending on the local temperature and vapor pressure. We perform a parameter study by varying the magnetized turbulence onset temperature, the accretion stress, the dust mass fraction, and the mass accretion rate. Our models feature a gas-only inner disk, a silicate sublimation front and dust rim starting at around 0.08 au, an ionization transition zone with a corresponding density jump, and a pressure maximum which acts as a pebble trap at around 0.12 au. Migration torque maps show Earth- and super-Earth-mass planets halt in our model disks at orbital periods ranging from 10 to 22 days. Such periods are in good agreement with both the inferred location of the innermost planets in multiplanetary systems, and the break in planet occurrence rates from the Kepler sample at 10 days. In particular, models with small grains depleted produce a trap located at a 10-day orbital period, while models with a higher abundance of small grains present a trap at around a 17-day orbital period. The snow line lies at 1.6 au, near where the occurrence rate of the giant planets peaks. We conclude that the dust sublimation zone is crucial for forming close-in planets, especially when considering tightly packed super-Earth systems.

Statistical analyses from exoplanet surveys around low-mass stars indicate that super-Earth and Neptune-mass planets are more frequent than gas giants around such stars, in agreement with core accretion theory of planet formation. Using precise radial velocities derived from visual and near-infrared spectra, we report the discovery of a giant planet with a minimum mass of 0.46 Jupiter masses in an eccentric 204-day orbit around the very low-mass star GJ 3512. Dynamical models show that the high eccentricity of the orbit is most likely explained from planet-planet interactions. The reported planetary system challenges current formation theories and puts stringent constraints on the accretion and migration rates of planet formation and evolution models, indicating that disc instability may be more efficient in forming planets than previously thought.

Convergent migration involving multiple planets embedded in a viscous protoplanetary disc is expected to produce a chain of planets in mean motion resonances, but the multiplanet systems observed by the Kepler spacecraft are generally not in resonance. We demonstrate that under equivalent conditions, where in a viscous disc convergent migration will form a long-term stable system of planets in a chain of mean motion resonances, migration in an inviscid disc often produces a system which is highly dynamically unstable. In particular, if planets are massive enough to significantly perturb the disc surface density and drive vortex formation, the smooth capture of planets into mean motion resonances is disrupted. As planets pile up in close orbits, not protected by resonances, close encounters increase the probability of planet-planet collisions, even while the gas disc is still present. While inviscid discs often produce unstable non-resonant systems, stable, closely packed, non-resonant systems can also be formed. Thus, when examining the expectation for planet migration to produce planetary systems in mean motion resonances, the effective turbulent viscosity of the protoplanetary disc is a key parameter.

Gas-giant planets, like Jupiter and Saturn, acquire massive gaseous envelopes during the approximately 3 Myr-long lifetimes of protoplanetary discs. In the core accretion scenario, the formation of a solid core of around 10 Earth masses triggers a phase of rapid gas accretion. Previous 3D grid-based hydrodynamical simulations found runaway gas accretion rates corresponding to approximately 10 to 100 Jupiter masses per Myr. Such high accretion rates would result in all planets with larger-than-10-Earth-mass cores forming Jupiter-like planets, in clear contrast to the ice giants in the Solar System and the observed exoplanet population. In this work, we use 3D hydrodynamical simulations, that include radiative transfer, to model the growth of the envelope on planets with different masses. We find that gas flows rapidly through the outer part of the envelope, but this flow does not drive accretion. Instead, gas accretion is the result of quasi-static contraction of the inner envelope, which can be orders of magnitude smaller than the mass flow through the outer atmosphere. For planets smaller than Saturn, we measure moderate gas accretion rates that are below 1 Jupiter mass per Myr. Higher mass planets, however, accrete up to 10 times faster and do not reveal a self-driven mechanism that can halt gas accretion. Therefore, the reason for the final masses of Saturn and Jupiter remains difficult to understand, unless their completion coincided with the dissipation of the Solar Nebula.

We present the highest resolution study to date of super-Earths migrating in inviscid and low-viscosity discs, motivated by the connection to laminar, wind-driven models of protoplanetary discs. Our models unveil the critical role of vortices in determining the migration behaviour for partial gap-opening planets. Vortices form in pressure maxima at gap edges, and prevent the disc-feedback stopping of migration for intermediate planets in low-viscosity and inviscid discs, contrary to the concept of the `inertial limit' or `disc feedback' halting predicted from analytical models. Vortices may also form in the corotation region, and can dramatically modify migration behaviour through direct gravitational interaction with the planet. These features become apparent at high resolution, and for all but the highest viscosities there exist significant difficulties in obtaining numerically converged results. The migration of partial gap-opening planets, however, clearly becomes chaotic for sufficiently low viscosities. At moderate viscosity, a smooth disc-feedback regime is found in which migration can slow substantially, and the migration time-scale observed corresponds to migration being driven by diffusive relaxation of the gap edges. At high viscosity classical Type I migration is recovered. For Jupiter-analogue planets in inviscid discs, a wide, deep gap is formed. Transient Type II migration occurs over radial length-scales corresponding to the gap width, beyond which migration can stall. Finally, we examine the particle trapping driven by structures left in inviscid discs by a migrating planet, and find that particle traps in the form of multiple rings and vortices can persist long after the planet has passed. In this case, the observation of particle traps by submillimetre interferometers such as ALMA cannot be used to infer the current presence of an adjacent planet.

At a distance of 1.8 parsecs, Barnard's star (Gl 699) is a red dwarf with the largest apparent motion of any known stellar object. It is the closest single star to the Sun, second only to the alpha Centauri triple stellar system. Barnard's star is also among the least magnetically active red dwarfs known and has an estimated age older than our Solar System. Its properties have made it a prime target for planet searches employing techniques such as radial velocity, astrometry, and direct imaging, all with different sensitivity limits but ultimately leading to disproved or null results. Here we report that the combination of numerous measurements from high-precision radial velocity instruments reveals the presence of a low-amplitude but significant periodic signal at 233 days. Independent photometric and spectroscopic monitoring, as well as the analysis of instrumental systematic effects, show that this signal is best explained as arising from a planetary companion. The candidate planet around Barnard's star is a cold super-Earth with a minimum mass of 3.2 Earth masses orbiting near its snow-line. The combination of all radial velocity datasets spanning 20 years additionally reveals a long-term modulation that could arise from a magnetic activity cycle or from a more distant planetary object. Because of its proximity to the Sun, the proposed planet has a maximum angular separation of 220 milli-arcseconds from Barnard's star, making it an excellent target for complementary direct imaging and astrometric observations.

We consider an impact on an asteroid that is energetic enough to cause resurfacing by seismic reverberation and just below the catastrophic disruption threshold, assuming that seismic waves are not rapidly attenuated. In asteroids with diameter less than 1 km we identify a regime where rare energetic impactors can excite seismic waves with frequencies near those of the asteroid's slowest normal modes. In this regime, the distribution of seismic reverberation is not evenly distributed across the body surface. With mass-spring model elastic simulations, we model impact excitation of seismic waves with a force pulse exerted on the surface and using three different asteroid shape models. The simulations exhibit antipodal focusing and normal mode excitation. If the impulse excited vibrational energy is long lasting, vibrations are highest at impact point, its antipode and at high surface elevations such as an equatorial ridge. A near equatorial impact launches a seismic impulse on a non-spherical body that can be focused on two additional points on an the equatorial ridge. We explore simple flow models for the morphology of vibration induced surface slumping. We find that the initial seismic pulse is unlikely to cause large shape changes. Long lasting seismic reverberation on Bennu caused by a near equatorial impact could have raised the height of its equatorial ridge by a few meters and raised two peaks on it, one near impact site and the other near its antipode.

Recent developments in non-ideal magnetohydrodynamic simulations of protoplanetary disks suggest that instead of being traditional turbulent (viscous) accretion disks, they have a largely laminar flow with accretion driven by large-scale wind torques. These disks are possibly threaded by Hall-effect generated large-scale horizontal magnetic fields. We have examined the dynamics of the corotation region of a low mass planet embedded in such a disk and the evolution of the associated migration torque. These disks lack strong turbulence and associated turbulent diffusion, and the presence of a magnetic field and radial gas flow presents a situation outside the applicability of previous corotation torque theory. We summarize the analytical analysis of the corotation torque, give details on the numerical methods used, and in particular the relative merits of different numerical schemes for the inviscid problem.

The known exoplanet population displays a great diversity of orbital architectures, and explaining the origin of this is a major challenge for planet formation theories. The gravitational interaction between young planets and their protoplanetary disks provides one way in which planetary orbits can be shaped during the formation epoch. Disk-planet interactions are strongly influenced by the structure and physical processes that drive the evolution of the protoplanetary disk. In this review we focus on how disk-planet interactions drive the migration of planets when different assumptions are made about the physics of angular momentum transport, and how it drives accretion flows in protoplanetary disk models. In particular, we consider migration in discs where: (i) accretion flows arise because turbulence diffusively transports angular momentum; (ii) laminar accretion flows are confined to thin, ionised layers near disk surfaces and are driven by the launching of magneto-centrifugal winds, with the midplane being completely inert; (iii) laminar accretion flows pervade the full column density of the disc, and are driven by a combination of large scale horizontal and vertical magnetic fields.

We examine the migration of low mass planets in laminar protoplanetary discs, threaded by large scale magnetic fields in the dead zone that drive radial gas flows. As shown in Paper I, a dynamical corotation torque arises due to the flow-induced asymmetric distortion of the corotation region and the evolving vortensity contrast between the librating horseshoe material and background disc flow. Using simulations of laminar torqued discs containing migrating planets, we demonstrate the existence of the four distinct migration regimes predicted in Paper I. In two regimes, the migration is approximately locked to the inward or outward radial gas flow, and in the other regimes the planet undergoes outward runaway migration that eventually settles to fast steady migration. In addition, we demonstrate torque and migration reversals induced by midplane magnetic stresses, with a bifurcation dependent on the disc surface density. We develop a model for fast migration, and show why the outward runaway saturates to a steady speed, and examine phenomenologically its termination due to changing local disc conditions. We also develop an analytical model for the corotation torque at late times that includes viscosity, for application to discs that sustain modest turbulence. Finally, we use the simulation results to develop torque prescriptions for inclusion in population synthesis models of planet formation.

Although most of the circumbinary planets detected by the Kepler spacecraft are on orbits that are closely aligned with the binary orbital plane, the systems Kepler-413 and Kepler-453 exhibit small misalignments of $\sim 2.5^\circ$. One possibility is that these planets formed in a circumbinary disc whose midplane was inclined relative to the binary orbital plane. Such a configuration is expected to lead to a warped and twisted disc, and our aim is to examine the inclination evolution of planets embedded in these discs. We employed 3D hydrodynamical simulations that examine the disc response to the presence of a modestly inclined binary with parameters that match the Kepler-413 system, as a function of disc parameters and binary inclinations. The discs all develop slowly varying warps, and generally display very small amounts of twist. Very slow solid body precession occurs because a large outer disc radius is adopted. Simulations of planets embedded in these discs resulted in the planet aligning with the binary orbit plane for disc masses close to the minimum mass solar nebular, such that nodal precession of the planet was controlled by the binary. For higher disc masses, the planet maintains near coplanarity with the local disc midplane. Our results suggest that circumbinary planets born in tilted circumbinary discs should align with the binary orbit plane as the disc ages and loses mass, even if the circumbinary disc remains misaligned from the binary orbit. This result has important implications for understanding the origins of the known circumbinary planets.

Planets are born in protostellar disks, which are now observed with enough resolution to address questions about internal gas flows. Candidates for driving the flows include magnetic forces, but ionization state estimates suggest much of the gas mass decouples from magnetic fields. Thus, hydrodynamical instabilities could play a major role. We investigate disk dynamics under conditions typical for a T Tauri system, using global 3D radiation hydrodynamics simulations with embedded particles and a resolution of 70 cells per scale height. Stellar irradiation heating is included with realistic dust opacities. The disk starts in joint radiative balance and hydrostatic equilibrium. The vertical shear instability (VSI) develops into turbulence that persists up to at least 1600 inner orbits (143 outer orbits). Turbulent speeds are a few percent of the local sound speed at the midplane, increasing to 20%, or 100 m/s, in the corona. These are consistent with recent upper limits on turbulent speeds from optically thin and thick molecular line observations of TW Hya and HD 163296. The predominantly vertical motions induced by the VSI efficiently lift particles upwards. Grains 0.1 and 1 mm in size achieve scale heights greater than expected in isotropic turbulence. We conclude that while kinematic constraints from molecular line emission do not directly discriminate between magnetic and nonmagnetic disk models, the small dust scale heights measured in HL Tau and HD 163296 favor turbulent magnetic models, which reach lower ratios of the vertical kinetic energy density to the accretion stress.

Motivated by models suggesting that the inner planet forming regions of protoplanetary discs are predominantly lacking in viscosity-inducing turbulence, and are possibly threaded by Hall-effect generated large-scale horizontal magnetic fields, we examine the dynamics of the corotation region of a low-mass planet in such an environment. The corotation torque in an inviscid, isothermal, dead zone ought to saturate, with the libration region becoming both symmetrical and of a uniform vortensity, leading to fast inward migration driven by the Lindblad torques alone. However, in such a low viscosity situation, the material on librating streamlines essentially preserves its vortensity. If there is relative radial motion between the disc gas and the planet, the librating streamlines will no longer be symmetrical. Hence, if the gas is torqued by a large scale magnetic field so that it undergoes a net inflow or outflow past the planet, driving evolution of the vortensity and inducing asymmetry of the corotation region, the corotation torque can grow, leading to a positive torque. In this paper we treat this effect by applying a symmetry argument to the previously studied case of a migrating planet in an inviscid disc. Our results show that the corotation torque due to a laminar Hall-induced magnetic field in a dead zone behaves quite differently from that studied previously for a viscous disc. Furthermore, the magnetic field induced corotation torque and the dynamical corotation torque in a low viscosity disc can be regarded as one unified effect.

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 present the results of hydrodynamic simulations examining migration and growth of planets embedded in self-gravitating circumbinary discs. The binary star parameters are chosen to mimic those of the Kepler-16, -34 and -35 systems; the aim of this study is to examine the role of disc mass in determining the stopping locations of migrating planets at the edge of the cavity created by the central binary. Disc self-gravity can cause significant shrinkage of the cavity for disc masses in excess of 5--10 $\times$ the minimum mass solar nebula model. Planets forming early in the disc lifetime can migrate through the disc and stall at locations closer to the central star than is normally the case for lower mass discs, resulting in closer agreement between simulated and observed orbital architecture. The presence of a planet orbiting in the cavity of a massive disc can prevent the cavity size from expanding to the size of a lower mass disc. As the disc mass reduces over long time scales, this indicates that circumbinary planet systems retain memory of their initial conditions. Our simulations produce planetary orbits in good agreement with Kepler-16b without the need for self-gravity; Kepler-34 analogue systems produce wide and highly eccentric cavities, and self-gravity improves the agreement between simulations and data. Kepler-35b is more difficult to explain in detail due to it's relatively low mass, which results in the simulated stopping location being at a larger radius than that observed.

The types of instability in the interacting binary stars are reviewed. The project "Inter-Longitude Astronomy" is a series of smaller projects on concrete stars or groups of stars. It has no special funds, and is supported from resources and grants of participating organizations, when informal working groups are created. Totally we studied 1900+ variable stars of different types. The characteristic timescale is from seconds to decades and (extrapolating) even more. The monitoring of the first star of our sample AM Her was initiated by Prof. V.P. Tsesevich (1907-1983). Since more than 358 ADS papers were published. Some highlights of our photometric and photo-polarimetric monitoring and mathematical modelling of interacting binary stars of different types are presented: classical, asynchronous, intermediate polars and magnetic dwarf novae (DO Dra) with 25 timescales corresponding to different physical mechanisms and their combinations (part "Polar"); negative and positive superhumpers in nova-like and many dwarf novae stars ("Superhumper"); eclipsing "non-magnetic" cataclysmic variables; symbiotic systems ("Symbiosis"); super-soft sources (SSS, QR And); spotted (and not spotted) eclipsing variables with (and without) evidence for a current mass transfer ("Eclipser") with a special emphasis on systems with a direct impact of the stream into the gainer star's atmosphere, or V361 Lyr-type stars. Other parts of the ILA project are "Stellar Bell" (interesting pulsating variables of different types and periods - M, SR, RV Tau, RR Lyr, Delta Sct) and "Novice"(="New Variable") discoveries and classification with a subsequent monitoring for searching and studying possible multiple components of variability. Special mathematical methods have been developed to create a set of complementary software for statistically optimal modelling of variable stars of different types.

We present a study of four different formation scenarios that may be able to explain the origin of the recently announced planet (`Proxima b') orbiting Proxima Centauri. The aim is to examine how the formation scenarios differ in their predictions for the multiplicity of the Proxima system, the water/volatile content of Proxima b and its orbital eccentricity, so that these can be tested by future observations. A scenario of in situ formation via giant impacts from a locally enhanced disc of planetary embryos/planetesimals predicts that Proxima b will be in a multiplanet system with a measurably finite eccentricity. Assuming that the local solid enhancement needed to form a Proxima b analogue arises because of the inwards drift of solids in the form of small planetesimals/boulders, this scenario also results in Proxima b analogues being only moderately endowed with water/volatiles. A scenario in which multiple embryos form, migrate and mutually collide within a gas disc results in Proxima b being a member of a multiple system, possibly displaying mean motion resonances, but where the constituent members are Ocean planets due to accretion occurring mainly outside of the snowline. A scenario in which a single accreting embryo forms outside the snowline, and migrates inwards while accreting planetesimals/pebbles results in Proxima b being an isolated Ocean planet on a circular orbit. A scenario in which Proxima b formed via pebble accretion interior to the snowline produces a dry planet on a circular orbit. Future observations that characterise the physical and orbital properties of Proxima b, and any additional planets in the system, will provide valuable insights into the formation history of this neighbouring planetary system.

We have recently shown that spiral density waves propagating in accretion disks can undergo a parametric instability by resonantly coupling with and transferring energy into pairs of inertial waves (or inertial-gravity waves when buoyancy is important). In this paper, we perform inviscid three-dimensional global hydrodynamic simulations to examine the growth and consequence of this instability operating on the spiral waves driven by a Jupiter-mass planet in a protoplanetary disk. We find that the spiral waves are destabilized via the spiral wave instability (SWI), generating hydrodynamic turbulence and sustained radially-alternating vertical flows that appear to be associated with long wavelength inertial modes. In the interval $0.3~R_p \leq R \leq 0.7~R_p$, where $R_p$ denotes the semi-major axis of the planetary orbit (assumed to be 5~au), the estimated vertical diffusion rate associated with the turbulence is characterized by $\alpha_{\rm diff} \sim (0.2-1.2) \times 10^{-2}$. For the disk model considered here, the diffusion rate is such that particles with sizes up to several centimeters are vertically mixed within the first pressure scale height. This suggests that the instability of spiral waves launched by a giant planet can significantly disperse solid particles and trace chemical species from the midplane. In planet formation models where the continuous local production of chondrules/pebbles occurs over Myr time scales to provide a feedstock for pebble accretion onto these bodies, this stirring of solid particles may add a time constraint: planetary embryos and large asteroids have to form before a gas giant forms in the outer disk, otherwise the SWI will significantly decrease the chondrule/pebble accretion efficiency.

We present the results of 2-dimensional hydrodynamic simulations of self-gravitating circumbinary discs around binaries whose parameters match those of the circumbinary planet-hosting systems Kepler-16, -34 and -35. Previous work has shown that non-self-gravitating discs in these systems form an eccentric precessing inner cavity due to tidal truncation by the binary, and planets which form at large radii migrate until stalling at this cavity. While this scenario appears to provide a natural explanation for the observed orbital locations of the circumbinary planets, previous simulations have failed to match the observed planet orbital parameters. The aim of this work is to examine the role of self-gravity in modifying circumbinary disc structure as a function of disc mass, prior to considering the evolution of embedded circumbinary planets. In agreement with previous work, we find that for disc masses between 1--5 times the minimum mass solar nebula (MMSN), disc self-gravity affects modest changes in the structure and evolution of circumbinary discs. Increasing the disc mass to 10 or 20 MMSN leads to two dramatic changes in disc structure. First, the scale of the inner cavity shrinks substantially, bringing its outer edge closer to the binary. Second, in addition to the eccentric inner cavity, additional precessing eccentric ring-like features develop in the outer regions of the discs. If planet formation starts early in the disc lifetime, these changes will have a significant impact on the formation and evolution of planets and precursor material.

Recent high angular resolution ($\sim$0.2") ALMA observations of the 1.1 mm continuum and of HCO+ J=3-2 and HCN J=3-2 gas towards the binary protostar IRS 43 reveal multiple Keplerian disks which are significantly misaligned ($\gt$ 60$^\circ$), both in inclination and position angle and also with respect to the binary orbital plane. Each stellar component has an associated circumstellar disk while the binary is surrounded by a circumbinary disk. Together with archival VLA measurements of the stellar positions over 25 years, and assuming a circular orbit, we use our continuum measurements to determine the binary separation, a = 74 $\pm$ 4 AU, and its inclination, i $\lt$ 30$^\circ$. The misalignment in this system suggests that turbulence has likely played a major role in the formation of IRS 43.

At a distance of 1.295 parsecs, the red-dwarf Proxima Centauri ($\alpha$ Centauri C, GL 551, HIP 70890, or simply Proxima) is the Sun's closest stellar neighbor and one of the best studied low-mass stars. It has an effective temperature of only $\sim$ 3050 K, a luminosity of $\sim$0.1 per cent solar, a measured radius of 0.14 R$_\odot$ and a mass of about 12 per cent the mass of the Sun. Although Proxima is considered a moderately active star, its rotation period is $\sim$ 83 days, and its quiescent activity levels and X-ray luminosity are comparable to the Sun's. New observations reveal the presence of a small planet orbiting Proxima with a minimum mass of 1.3~Earth masses and an orbital period of $\sim$11.2 days. Its orbital semi-major axis is $\sim0.05$ AU, with an equilibrium temperature in the range where water could be liquid on its surface.

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

In this paper we present a combined photometric, spectroscopic and orbital period study of three early-type eclipsing binary systems: XZ Aql, UX Her, and AT Peg. As a result, we have derived the absolute parameters of their components and, on that basis, we discuss their evolutionary states. Furthermore, we compare their parameters with those of other binary systems and with the theoretical models. An analysis of all available up-to-date times of minima indicated that all three systems studied here show cyclic orbital changes, their origin is discussed in detail. Finally, we performed a frequency analysis for possible pulsational behavior and as a result we suggest that XZ Aql hosts a {\delta} Scuti component.

We present results from a suite of three-dimensional global hydrodynamic simulations which show that spiral density waves propagating in circumstellar disks are unstable to the growth of a parametric instability that leads to break-down of the flow into turbulence. This spiral wave instability (SWI) arises from a resonant interaction between pairs of inertial waves, or inertial-gravity waves, and the background spiral wave. The development of the instability in the linear regime involves the growth of a broad spectrum of inertial modes, with growth rates on the order of the orbital time, and results in a nonlinear saturated state in which turbulent velocity perturbations are of a similar magnitude to those induced by the spiral wave. The turbulence induces angular momentum transport, and vertical mixing, at a rate that depends locally on the amplitude of the spiral wave (we obtain a stress parameter $\alpha \sim 5 \times 10^{-4}$ in our reference model). The instability is found to operate in a wide-range of disk models, including those with isothermal or adiabatic equations of state, and in viscous disks where the dimensionless kinematic viscosity $\nu\le 10^{-5}$. This robustness suggests that the instability will have applications to a broad range of astrophysical disk-related phenomena, including those in close binary systems, planets embedded in protoplanetary disks (including Jupiter in our own Solar System) and FU Orionis outburst models. Further work is required to determine the nature of the instability, and to evaluate its observational consequences, in physically more complete disk models than we have considered in this paper.

Our recent N-body simulations of planetary system formation, incorporating models for the main physical processes thought to be important during the building of planets (i.e. gas disc evolution, migration, planetesimal/boulder accretion, gas accretion onto cores, etc.), have been successful in reproducing some of the broad features of the observed exoplanet population (e.g. compact systems of low mass planets, hot Jupiters), but fail completely to form any surviving cold Jupiters. The primary reason for this failure is rapid inward migration of growing protoplanets during the gas accretion phase, resulting in the delivery of these bodies onto orbits close to the star. Here, we present the results of simulations that examine the formation of gas giant planets in protoplanetary discs that are radially structured due to spatial and temporal variations in the effective viscous stresses, and show that such a model results in the formation of a population of cold gas giants. Furthermore, when combined with models for disc photoevaporation and a central magnetospheric cavity, the simulations reproduce the well-known hot-Jupiter/cold-Jupiter dichotomy in the observed period distribution of giant exoplanets, with a period valley between 10-100 days.

OJ287 is a quasi-periodic quasar with roughly 12 year optical cycles. It displays prominent outbursts which are predictable in a binary black hole model. The model predicted a major optical outburst in December 2015. We found that the outburst did occur within the expected time range, peaking on 2015 December 5 at magnitude 12.9 in the optical R-band. Based on Swift/XRT satellite measurements and optical polarization data, we find that it included a major thermal component. Its timing provides an accurate estimate for the spin of the primary black hole, chi = 0.313 +- 0.01. The present outburst also confirms the established general relativistic properties of the system such as the loss of orbital energy to gravitational radiation at the 2 % accuracy level and it opens up the possibility of testing the black hole no-hair theorem with a 10 % accuracy during the present decade.

We present the results of planet formation N-body simulations based on a comprehensive physical model that includes planetary mass growth through mutual embryo collisions and planetesimal/boulder accretion, viscous disc evolution, planetary migration and gas accretion onto planetary cores. The main aim of this study is to determine which set of model parameters leads to the formation of planetary systems that are similar to the compact low mass multi-planet systems that have been discovered by radial velocity surveys and the Kepler mission. We vary the initial disc mass, solids-to-gas ratio and the sizes of the boulders/planetesimals, and for a restricted volume of the parameter space we find that compact systems containing terrestrial planets, super-Earths and Neptune-like bodies arise as natural outcomes of the simulations. Disc models with low values of the solids-to-gas ratio can only form short-period super-Earths and Neptunes when small planetesimals/boulders provide the main source of accretion, since the mobility of these bodies is required to overcome the local isolation masses for growing embryos. The existence of short-period super-Earths around low metallicity stars provides strong evidence that small, mobile bodies (planetesimals, boulders or pebbles) played a central role in the formation of the observed planets.

We present the results of 2D and 3D hydrodynamic simulations of idealized protoplanetary discs that examine the formation and evolution of vortices by the vertical shear instability (VSI). In agreement with recent work, we find that discs with radially decreasing temperature profiles and short thermal relaxation time-scales, are subject to the axisymmetric VSI. In three dimensions, the resulting velocity perturbations give rise to quasi-axisymmetric potential vorticity perturbations that break-up into discrete vortices, in a manner that is reminiscent of the Rossby wave instability. Discs with very short thermal evolution time-scales (i.e. {\tau}<0.1 local orbit periods) develop strong vorticity perturbations that roll up into vortices that have small aspect ratios ({\chi}<2) and short lifetimes (~ a few orbits). Longer thermal time-scales give rise to vortices with larger aspect ratios (6<{\chi}<10), and lifetimes that depend on the entropy gradient. A steeply decreasing entropy profile leads to vortex lifetimes that exceed the simulation run times of hundreds of orbital periods. Vortex lifetimes in discs with positive or weakly decreasing entropy profiles are much shorter, being 10s of orbits at most, suggesting that the subcritical baroclinic instability plays an important role in sustaining vortices against destruction through the elliptical instability. Applied to the outer regions of protoplanetary discs, where the VSI is most likely to occur, our results suggest that vortices formed by the VSI are likely to be short lived structures.

We measure the photometric rotation periods of the components of multiple systems in young stellar associations to investigate the causes of the observed rotation period dispersion. We present the case of the wide binary AG Tri in the 23-Myr young beta Pictoris Association consisting of K4 + M1 dwarfs. Our multi-band, multi-season photometric monitoring allowed us to measure the rotation periods of both components P_A = 12.4d and P_B = 4.66d, to detect a prominent magnetic activity in the photosphere, likely responsible for the measured radial velocity variations, and for the first time, a flare event on the M1 component AG Tri B. We investigate either the possibility that the faster rotating component may have suffered an enhanced primordial disc dispersal, starting its PMS spin-up earlier than the slower rotating component, or the possibility that the formation of a debris disc may have prevented AG Tri A from gaining part of the angular momentum from the accreting disc.

We present the results of the spectroscopic, photometric and orbital period variation analyses of the detached eclipsing binary \astrobj{V482~Per}. We derived the absolute parameters of the system (M$_{1}$ = 1.51 M$_{\odot}$, M$_{2}$ = 1.29 M$_{\odot}$, R$_{1}$ = 2.39 R$_{\odot}$, R$_{2}$ = 1.45 R$_{\odot}$, L$_{1}$ = 10.15 L$_{\odot}$, L$_{2}$ = 3.01 L$_{\odot}$) for the first time in literature, based on an analysis of our own photometric and spectroscopic observations. We confirm the nature of the variations observed in the system's orbital period, suggested to be periodic by earlier works. A light time effect due to a physically bound, star-sized companion (M$_{3}$ = 2.14 M$_{\odot}$) on a highly eccentric (e = 0.83) orbit, seems to be the most likely cause. We argue that the companion can not be a single star but another binary instead. We calculated the evolutionary states of the system's components, and we found that the primary is slightly evolving after the Main Sequence, while the less massive secondary lies well inside it.

The Vertical Shear Instability is one of two known mechanisms potentially active in the so-called dead zones of protoplanetary accretion disks. A recent analysis indicates that a subset of unstable modes shows unbounded growth - both as resolution is increased and when the nominal lid of the atmosphere is extended, possibly indicating ill-posedness in previous attempts of linear analysis. The reduced equations governing the instability are revisited and the generated solutions are examined using both the previously assumed separable forms and an improved non-separable solution form that is herewith introduced. Analyzing the reduced equations using the separable form shows that, while the low-order body modes have converged eigenvalues and eigenfunctions (as both the vertical boundaries of the atmosphere are extended and with increased radial resolution), it is also confirmed that the corresponding high-order body modes and the surface modes do indeed show unbounded growth rates. However, the energy contained in both the higher-order body modes and surface modes diminishes precipitously due to the disk's Gaussian density profile. Most of the energy of the instability is contained in the low-order modes. An inseparable solution form is introduced which filters out the inconsequential surface modes leaving only body modes (both low and high-order ones). The analysis predicts a fastest growing mode with a specific radial length scale. The growth rates associated with the fundamental corrugation and breathing modes matches the growth and length scales observed in previous nonlinear studies of the instability.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

While the Standard Model of particle physics does not include free particles with fractional charge, experimental searches have not ruled out their existence. We report results from the Cryogenic Dark Matter Search (CDMS II) experiment that give the first direct-detection limits for cosmogenically-produced relativistic particles with electric charge lower than $e$/6. A search for tracks in the six stacked detectors of each of two of the CDMS II towers found no candidates, thereby excluding new parameter space for particles with electric charges between $e$/6 and $e$/200.

Protoplanetary disks are believed to accrete onto their central T Tauri star because of magnetic stresses. Recently published shearing box simulations indicate that Ohmic resistivity, ambipolar diffusion and the Hall effect all play important roles in disk evolution. In the presence of a vertical magnetic field, the disk remains laminar between 1-5au, and a magnetocentrifugal disk wind forms that provides an important mechanism for removing angular momentum. Questions remain, however, about the establishment of a true physical wind solution in the shearing box simulations because of the symmetries inherent in the local approximation. We present global MHD simulations of protoplanetary disks that include Ohmic resistivity and ambipolar diffusion, where the time-dependent gas-phase electron and ion fractions are computed under FUV and X-ray ionization with a simplified recombination chemistry. Our results show that the disk remains laminar, and that a physical wind solution arises naturally in global disk models. The wind is sufficiently efficient to explain the observed accretion rates. Furthermore, the ionization fraction at intermediate disk heights is large enough for magneto-rotational channel modes to grow and subsequently develop into belts of horizontal field. Depending on the ionization fraction, these can remain quasi-global, or break-up into discrete islands of coherent field polarity. The disk models we present here show a dramatic departure from our earlier models including Ohmic resistivity only. It will be important to examine how the Hall effect modifies the evolution, and to explore the influence this has on the observational appearance of such systems, and on planet formation and migration.

We report on the results of a search for a Weakly Interacting Massive Particle (WIMP) signal in low-energy data of the Cryogenic Dark Matter Search (CDMS~II) experiment using a maximum likelihood analysis. A background model is constructed using GEANT4 to simulate the surface-event background from $^{210}$Pb decay-chain events, while using independent calibration data to model the gamma background. Fitting this background model to the data results in no statistically significant WIMP component. In addition, we perform fits using an analytic ad hoc background model proposed by Collar and Fields, who claimed to find a large excess of signal-like events in our data. We confirm the strong preference for a signal hypothesis in their analysis under these assumptions, but excesses are observed in both single- and multiple-scatter events, which implies the signal is not caused by WIMPs, but rather reflects the inadequacy of their background model.

We improve on our previous treatments of long-term evolution of protostellar disks by explicitly solving disk self-gravity in two dimensions. The current model is an extension of the one-dimensional layered accretion disk model of Bae et al. We find that gravitational instability (GI)-induced spiral density waves heat disks via compressional heating (i.e. $P\rm{d}V$ work), and can trigger accretion outbursts by activating the magnetorotational instability (MRI) in the magnetically inert disk dead-zone. The GI-induced spiral waves propagate well inside of gravitationally unstable region before they trigger outbursts at $R \lesssim 1$ AU where GI cannot be sustained. This long-range propagation of waves cannot be reproduced with the previously used local $\alpha$ treatments for GI. In our standard model where zero dead-zone residual viscosity ($\alpha_{\rm rd}$) is assumed, the GI-induced stress measured at the onset of outbursts is locally as large as $0.01$ in terms of the generic $\alpha$ parameter. However, as suggested in our previous one-dimensional calculations, we confirm that the presence of a small but finite $\alpha_{\rm rd}$ triggers thermally-driven bursts of accretion instead of the GI + MRI-driven outbursts that are observed when $\alpha_{\rm rd}=0$. The inclusion of non-zero residual viscosity in the dead-zone decreases the importance of GI soon after mass feeding from the envelope cloud ceases. During the infall phase while the central protostar is still embedded, our models stay in a quiescent accretion phase with $\dot{M}_{acc}\sim10^{-8}-10^{-7}~M_{\odot}~{\rm yr^{-1}}$ over $60~\%$ of the time and spend less than $15~\%$ of the infall phase in accretion outbursts.

We present N-body simulations of planetary system formation in thermally-evolving, viscous disc models. The simulations incorporate type I migration (including corotation torques and their saturation), gap formation, type II migration, gas accretion onto planetary cores, and gas disc dispersal through photoevaporation. The aim is to examine whether or not the oligarchic growth scenario, when combined with self-consistent disc models and up-to-date prescriptions for disc-driven migration, can produce planetary systems similar to those that have been observed. The results correlate with the initial disc mass. Low mass discs form close-packed systems of terrestrial-mass planets and super-Earths. Higher mass discs form multiple generations of planets, with masses in the range 10 < mp < 45M_Earth. These planets generally type I migrate into the inner disc, because of corotation torque saturation, where they open gaps and type II migrate into the central star. Occasionally, a final generation of low-to-intermediate mass planets forms and survives due to gas disc dispersal. No surviving gas giants were formed in our simulations. Analysis shows that these planets can only survive migration if a core forms and experiences runaway gas accretion at orbital radii r > 10 au prior to the onset of type II migration. We conclude that planet growth above masses mp > 10M_Earth during the gas disc life time leads to corotation torque saturation and rapid inward migration, preventing the formation and survival of gas giants. This result is in contrast to the success in forming gas giant planets displayed by some population synthesis models. This discrepancy arises, in part, because the type II migration prescription adopted in the population synthesis models causes too large a reduction in the migration speed when in the planet dominated regime.

Over the last twenty years, the search for extrasolar planets revealed us the rich diversity of the outcomes of the formation and evolution of planetary systems. In order to fully understand how these extrasolar planets came to be, however, the orbital and physical data we possess are not enough, and they need to be complemented with information on the composition of the exoplanets. Ground-based and space-based observations provided the first data on the atmospheric composition of a few extrasolar planets, but a larger and more detailed sample is required before we can fully take advantage of it. The primary goal of the Exoplanet Characterization Observatory (EChO) is to fill this gap, expanding the limited data we possess by performing a systematic survey of hundreds of extrasolar planets. The full exploitation of the data that EChO and other space-based and ground-based facilities will provide in the near future, however, requires the knowledge of what are the sources and sinks of the chemical species and molecules that will be observed. Luckily, the study of the past history of the Solar System provides several indications on the effects of processes like migration, late accretion and secular impacts, and on the time they occur in the life of planetary systems. In this work we will review what is already known about the factors influencing the composition of planetary atmospheres, focusing on the case of gaseous giant planets, and what instead still need to be investigated.

First photometric study for the newly discovered systems USNO-B1.0 1091-0130715 and GSC-03449-0680 were applied by means of recent windows interface version of the Wilson and Devinney code based on model atmospheres by Kurucz (1993). The accepted models reveal some absolute parameters for both systems, which used in adopting the spectral type of the systems components and follow their evolutionary status. Distances to each systems are and physical properties are estimated. Comparisons of the computed physical parameters with stellar models are discussed.

We report a first search for weakly interacting massive particles (WIMPs) using the background rejection capabilities of SuperCDMS. An exposure of 577 kg-days was analyzed for WIMPs with mass < 30 GeV/c2, with the signal region blinded. Eleven events were observed after unblinding. We set an upper limit on the spin-independent WIMP-nucleon cross section of 1.2e-42 cm2 at 8 GeV/c2. This result is in tension with WIMP interpretations of recent experiments and probes new parameter space for WIMP-nucleon scattering for WIMP masses < 6 GeV/c2.

PLATO 2.0 has recently been selected for ESA's M3 launch opportunity (2022/24). Providing accurate key planet parameters (radius, mass, density and age) in statistical numbers, it addresses fundamental questions such as: How do planetary systems form and evolve? Are there other systems with planets like ours, including potentially habitable planets? The PLATO 2.0 instrument consists of 34 small aperture telescopes (32 with 25 sec readout cadence and 2 with 2.5 sec candence) providing a wide field-of-view (2232 deg2) and a large photometric magnitude range (4-16 mag). It focusses on bright (4-11 mag) stars in wide fields to detect and characterize planets down to Earth-size by photometric transits, whose masses can then be determined by ground-based radial-velocity follow-up measurements. Asteroseismology will be performed for these bright stars to obtain highly accurate stellar parameters, including masses and ages. The combination of bright targets and asteroseismology results in high accuracy for the bulk planet parameters: 2%, 4-10% and 10% for planet radii, masses and ages, respectively. The planned baseline observing strategy includes two long pointings (2-3 years) to detect and bulk characterize planets reaching into the habitable zone (HZ) of solar-like stars and an additional step-and-stare phase to cover in total about 50% of the sky. PLATO 2.0 will observe up to 1,000,000 stars and detect and characterize hundreds of small planets, and thousands of planets in the Neptune to gas giant regime out to the HZ. It will therefore provide the first large-scale catalogue of bulk characterized planets with accurate radii, masses, mean densities and ages. This catalogue will include terrestrial planets at intermediate orbital distances, where surface temperatures are moderate. Coverage of this parameter range with statistical numbers of bulk characterized planets is unique to PLATO 2.0.

We present new radial velocities and photometry of the short-period Algol TT Herculis. Previous attempts to model the light curves of the system have met with limited success, primarily because of the lack of a reliable mass ratio. Our spectroscopic observations are the first to result in radial velocities for the secondary star, and thus provide a spectroscopic mass ratio. Simultaneous analysis of the radial velocities and new photometry shows that the system is a double contact binary, with a rapidly rotating primary that fills its limiting lobe.

SuperCDMS is an experiment designed to directly detect Weakly Interacting Massive Particles (WIMPs), a favored candidate for dark matter ubiquitous in the Universe. In this paper, we present WIMP-search results using a calorimetric technique we call CDMSlite, which relies on voltage- assisted Luke-Neganov amplification of the ionization energy deposited by particle interactions. The data were collected with a single 0.6 kg germanium detector running for 10 live days at the Soudan Underground Laboratory. A low energy threshold of 170 eVee (electron equivalent) was obtained, which allows us to constrain new WIMP-nucleon spin-independent parameter space for WIMP masses below 6 GeV/c2.

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.

We present global hydrodynamic and magnetohydrodynamic (MHD) simulations with mesh refinement of accreting planets embedded in protoplanetary disks (PPDs). The magnetized disk includes Ohmic resistivity that depends on the overlying mass column, leading to turbulent surface layers and a dead zone near the midplane. The main results are: (i) The accretion flow in the Hill sphere is intrinsically 3D for hydrodynamic and MHD models. Net inflow toward the planet is dominated by high latitude flows. A circumplanetary disk (CPD) forms. Its midplane flows outward in a pattern whose details differ between models. (ii) Gap opening magnetically couples and ignites the dead zone near the planet, leading to stochastic accretion, a quasi-turbulent flow in the Hill sphere and a CPD whose structure displays high levels of variability. (iii) Advection of magnetized gas onto the rotating CPD generates helical fields that launch magnetocentrifugally driven outflows. During one specific epoch a highly collimated, one-sided jet is observed. (iv) The CPD's surface density $\sim30{\rm\,g\,cm^{-2}}$, small enough for significant ionization and turbulence to develop. (v) The accretion rate onto the planet in the MHD simulation reaches a steady value $8 \times 10^{-3} {\rm M_\oplus yr^{-1}}$, and is similar in the viscous hydrodynamic runs. Our results suggest that gas accretion onto a forming giant planet within a magnetized PPD with dead zone allows rapid growth from Saturnian to Jovian masses. As well as being relevant for giant planet formation, these results have important implications for the formation of regular satellites around gas giant planets.

We report results of a search for Weakly Interacting Massive Particles (WIMPS) with the silicon detectors of the CDMS II experiment. This blind analysis of 140.2 kg-days of data taken between July 2007 and September 2008 revealed three WIMP-candidate events with a surface-event background estimate of 0.41^{+0.20}_{-0.08}(stat.)^{+0.28}_{-0.24}(syst.). Other known backgrounds from neutrons and 206Pb are limited to < 0.13 and <0.08 events at the 90% confidence level, respectively. The exposure of this analysis is equivalent to 23.4 kg-days for a recoil energy range of 7-100 keV for a WIMP of mass 10 GeV/c2. The probability that the known backgrounds would produce three or more events in the signal region is 5.4%. A profile likelihood ratio test of the three events that includes the measured recoil energies gives a 0.19% probability for the known-background-only hypothesis when tested against the alternative WIMP+background hypothesis. The highest likelihood occurs for a WIMP mass of 8.6 GeV/c2 and WIMP-nucleon cross section of 1.9e-41 cm2.

The SuperCDMS experiment in the Soudan Underground Laboratory searches for dark matter with a 9-kg array of cryogenic germanium detectors. Symmetric sensors on opposite sides measure both charge and phonons from each particle interaction, providing excellent discrimination between electron and nuclear recoils, and between surface and interior events. Surface event rejection capabilities were tested with two $^{210}$Pb sources producing $\sim$130 beta decays/hr. In $\sim$800 live hours, no events leaked into the 8--115 keV signal region, giving upper limit leakage fraction $1.7 \times 10^{-5}$ at 90% C.L., corresponding to $< 0.6$ surface event background in the future 200-kg SuperCDMS SNOLAB experiment.

We present the results of high resolution 2D simulations of low mass planets on fixed eccentric orbits embedded in protoplanetary discs. The aim of this study is to determine how the strength of the sustained, non-linear corotation torque experienced by embedded planets varies as a function of orbital eccentricity, disc parameters, and planetary mass. In agreement with previous work we find that the corotation torque diminishes as orbital eccentricity, $e$, increases. Analysis of the time-averaged streamlines in the disc demonstrates that the width of the horseshoe region narrows as the eccentricity increases, and we suggest that this narrowing largely explains the observed decrease in the corotation torque. We employ three distinct methods for estimating the strength of the unsaturated corotation torque from our simulations, and provide an empirical fit to these results. We find that a simple model where the corotation torque, $\Gamma_C$, decreases exponentially with increasing eccentricity (i.e. $\Gamma_C \propto \exp{(-e/e_f)}$) provides a good global fit to the data with an e-folding eccentricity, $e_f$, that scales linearly with the disc scale height at the planet location. We confirm that this model provides a good fit for planet masses of 5 and 10 $\Me$ in our simulations. The formation of planetary systems is likely to involve significant planet-planet interactions that will excite eccentric orbits, and this is likely to influence disc-driven planetary migration through modification of the corotation torque. Our results suggest that high fidelity models of planetary formation should account for these effects.

We report results of a search for Weakly Interacting Massive Particles (WIMPs) with the Si detectors of the CDMS II experiment. This report describes a blind analysis of the first data taken with CDMS II's full complement of detectors in 2006-2007; results from this exposure using the Ge detectors have already been presented. We observed no candidate WIMP-scattering events in an exposure of 55.9 kg-days before analysis cuts, with an expected background of ~1.1 events. The exposure of this analysis is equivalent to 10.3 kg-days over a recoil energy range of 7-100 keV for an ideal Si detector and a WIMP mass of 10 GeV/c2. These data set an upper limit of 1.7x10-41 cm2 on the WIMP-nucleon spin-independent cross section of a 10 GeV/c2 WIMP. These data exclude parameter space for spin-independent WIMP-nucleon elastic scattering that is relevant to recent searches for low-mass WIMPs.

Several circumbinary planets have been detected by the Kepler mission. Recent work has emphasized the difficulty of forming these planets at their observed locations. It has been suggested that these planets formed further out in their discs and migrated in to locations where they are observed. We examine the orbital evolution of planets embedded in circumbinary disc models for the three systems Kepler-16, Kepler-34 and Kepler-35. The aims are: to explore the plausibility of a formation scenario in which cores form at large distances from the binaries and undergo inward migration and gas accretion as the gas disc disperses; to determine which sets of disc parameters lead to planets whose final orbits provide reasonable fits to the observed systems. We performed simulations of a close binary system interacting with circumbinary discs with differing aspect ratios, and viscous stress parameters. Once the binary+disc system reaches quasi-equilibrium we embed a planet in the disc and examine its evolution under the action of binary and disc forces. We consider fully-formed planets with masses equal to those inferred from Kepler data, and low-mass cores that migrate and accrete gas while the gas disc is being dispersed. A typical outcome for all systems is stalling of inward migration as the planet enters the inner cavity formed by the binary system. The circumbinary disc becomes eccentric and the disc eccentricity forces the planet into a noncircular orbit. For each of the Kepler-16b, Kepler-34b and Kepler-35b systems we obtain planets whose parameters agree reasonably well with the observational data. The simulations presented here provide support for a formation scenario in which a core forms, migrates inward and accretes gas, but accurate fitting of the observed Kepler systems is likely to require disc models that are significantly more sophisticated in terms of their input physics.

We analyse the dynamics of the multiple planet system HD 181433. This consists of two gas giant planets (c and d) with msin i = 0.65 MJup and 0.53 MJup orbiting with periods 975 and 2468 days, respectively. The two planets appear to be in a 5:2 mean motion resonance, as this is required for the system to be dynamically stable. A third planet with mass m_b sin i = 0.023 MJup orbits close to the star with orbital period 9.37 days. Each planet orbit is significantly eccentric, with current values estimated to be e_b = 0.39, e_c = 0.27 and e_d = 0.47. In this paper we assess different scenarios that may explain the origin of these eccentric orbits, with particular focus on the innermost body, noting that the large eccentricity of planet b cannot be explained through secular interaction with the outer pair. We consider a scenario in which the system previously contained an additional giant planet that was ejected during a period of dynamical instability among the planets. N-body simulations are presented that demonstrate that during scattering and ejection among the outer planets a close encounter between a giant and the inner body can raise e_b to its observed value. We consider the possibility that an undetected planet in the system increases the secular forcing of planet b by the exterior giant planets, but we find that the resulting eccentricity is not large enough to agree with the observed one. We also consider a scenario in which the spin-down of the central star causes the system to pass through secular resonance. Spin-down rates below the critical value lead to longterm capture of planet b in secular resonance, driving the eccentricity toward unity. If additional short-period low mass planets are present in the system, however, we find that mutual scattering can release planet b from the secular resonance, leading to a system with orbital parameters similar to those observed today.

(Abridged) We analyse the stability and evolution of power-law accretion disc models. These have midplane densities that follow radial power-laws, and have either temperature or entropy distributions that are power-law functions of cylindrical radius. We employ two different hydrodynamic codes to perform 2D-axisymmetric and 3D simulations that examine the long-term evolution of the disc models as a function of the power-law indices of the temperature or entropy, the thermal relaxation time of the fluid, and the viscosity. We present a stability analysis of the problem that we use to interpret the simulation results. We find that disc models whose temperature or entropy profiles cause the equilibrium angular velocity to vary with height are unstable to the growth of modes with wavenumber ratios |k_R/k_Z| >> 1 when the thermodynamic response of the fluid is isothermal, or the thermal evolution time is comparable to or shorter than the local dynamical time scale. These discs are subject to the Goldreich-Schubert-Fricke (GSF) or `vertical shear' linear instability. Development of the instability involves excitation of vertical breathing and corrugation modes in the disc, with the corrugation modes in particular being a feature of the nonlinear saturated state. Instability operates when the dimensionless disc kinematic viscosity nu < 10^{-6} (Reynolds numbers Re>H c_s/nu > 2500). In 3D the instability generates a quasi-turbulent flow, and the Reynolds stress produces a fluctuating effective viscosity coefficient whose mean value reaches alpha ~ 6 x 10^{-4} by the end of the simulation. The vertical shear instability in disc models which include realistic thermal physics has yet to be examined. Should it occur, however, our results suggest that it will have significant consequences for their internal dynamics, transport properties, and observational appearance.

By carrying out two-dimensional two-fluid global simulations, we have studied the response of dust to gap formation by a single planet in the gaseous component of a protoplanetary disk - the so-called "dust filtration" mechanism. We have found that a gap opened by a giant planet at 20 AU in a \alpha=0.01, \dot{M}=10^{-8} Msun/yr disk can effectively stop dust particles larger than 0.1 mm drifting inwards, leaving a sub-millimeter dust cavity/hole. However, smaller particles are difficult to filter by a planet-induced gap due to 1) dust diffusion, and 2) a high gas accretion velocity at the gap edge. An analytic model is also derived to understand what size particles can be filtered by the gap edge. Finally, with our updated understanding of dust filtration, we have computed Monte-Carlo radiative transfer models with variable dust size distributions to generate the spectral energy distributions (SEDs) of disks with gaps. By comparing with transitional disk observations (e.g. GM Aur), we have found that dust filtration alone has difficulties to deplete small particles sufficiently to explain the near-IR deficit of transitional disks, except under some extreme circumstances. The scenario of gap opening by multiple planets studied previously suffers the same difficulty. One possible solution is by invoking both dust filtration and dust growth in the inner disk. In this scenario, a planet induced gap filters large dust particles in the disk, and the remaining small dust particles passing to the inner disk can grow efficiently without replenishment from fragmentation of large grains. Predictions for ALMA have also been made based on all these scenarios. We conclude that dust filtration with planet(s) in the disk is a promising mechanism to explain submm observations of transitional disks but it may need to be combined with other processes (e.g. dust growth) to explain the near-IR deficit.

We present the results of the study of the contact binary system BO CVn. We have obtained physical parameters of the components based on combined analysis of new, multi-color light curves and spectroscopic mass ratio. This is the first time the latter has been determined for this object. We derived the contact configuration for the system with a very high filling factor of about 88 percent. We were able to reproduce the observed light curve, namely the flat bottom of the secondary minimum, only if a third light has been added into the list of free parameters. The resulting third light contribution is significant, about 20-24 percent, while the absolute parameters of components are: M1=1.16, M2=0.39, R1=1.62 and R2=1.00 (in solar units). The O-C diagram shows an upward parabola which, under the conservative mass transfer assumption, would correspond to a mass transfer rate of dM/dt = 6.3 \times 10-8M\odot/yr, matter being transferred from the less massive component to the more massive one. No cyclic, short-period variations have been found in the O-C diagram (but longer-term variations remain a possibility)

As planets form and grow within gaseous protoplanetary disks, the mutual gravitational interaction between the disk and planet leads to the exchange of angular momentum, and migration of the planet. We review current understanding of disk-planet interactions, focussing in particular on physical processes that determine the speed and direction of migration. We describe the evolution of low mass planets embedded in protoplanetary disks, and examine the influence of Lindblad and corotation torques as a function of the disk properties. The role of the disk in causing the evolution of eccentricities and inclinations is also discussed. We describe the rapid migration of intermediate mass planets that may occur as a runaway process, and examine the transition to gap formation and slower migration driven by the viscous evolution of the disk for massive planets. The roles and influence of disk self-gravity and magnetohydrodynamic turbulence are discussed in detail, as a function of the planet mass, as is the evolution of multiple planet systems. Finally, we address the question of how well global models of planetary formation that include migration are able to match observations of extrasolar planets.

(Abridged) Planetesimals embedded in a protoplanetary disc are stirred by gravitational torques exerted by density fluctuations in the surrounding turbulence. In particular, planetesimals in a disc supporting fully developed magneto-rotational turbulence are readily excited to velocity dispersions above the threshold for catastrophic disruption, halting planet formation. We aim to examine the stirring of planetesimals lying instead in a magnetically-decoupled midplane dead zone, stirred only by spiral density waves propagating out of the disc's magnetically-coupled turbulent surface layers. We extend previous studies to include a wider range of disc models, and explore the effects of varying the disc column density and external magnetic field strength. [...] The strength of the stirring is found to be independent of the gas surface density, which is contrary to the increase with disc mass expected from a simple linear wave picture. The discrepancy arises from the shearing out of density waves as they propagate into the dead zone, resulting in density structures near the midplane that exert weaker stochastic torques on average. We provide a simple analytic fit to our numerically obtained torque amplitudes that accounts for this effect. The stirring on the other hand depends sensitively on the net vertical magnetic flux, up to a saturation level above which magnetic forces dominate in the turbulent layers. For the majority of our models, the equilibrium planetesimal velocity dispersions lie between the thresholds for disrupting strong and weak aggregates, suggesting that collision outcomes will depend on material properties. However, discs with relatively weak magnetic fields yield reduced stirring, and their dead zones provide safe-havens even for the weakest planetesimals against collisional destruction.

(Abridged) We present the results of N-body simulations of planetary systems formation in radiatively-inefficient disc models, where positive corotation torques may counter the rapid inward migration of low mass planets driven by Lindblad torques. The aim of this work is to examine the nature of planetary systems that arise from oligarchic growth in such discs. We adapt the commonly-used Mercury-6 symplectic integrator by including simple prescriptions for planetary migration (types I and II), planetary atmospheres that enhance the probability of planetesimal accretion by protoplanets, gas accretion onto forming planetary cores, and gas disc dispersal. We perform a suite of simulations for a variety of disc models with power-law surface density and tempera- ture profiles, with a focus on models in which unsaturated corotation torques can drive outward migration of protoplanets. In some models we account for the quenching of corotation torques that arises when planetary orbits become eccentric. Approximately half of our simulations lead to the successful formation of gas giant planets with a broad range of masses and semimajor axes. We conclude that convergent migration induced by corotation torques operating during planet formation can enhance the growth rate of planetary cores, but these often migrate into the central star because corotation torques saturate. Outward migration of planetary cores of modest mass can lead to the formation of gas giant planets at large distances from the central star, similar to those observed recently through direct imaging surveys. The excitation of planetary eccentricities through planet-planet scat- tering during oligarchic growth may quench the effects of corotation torques, however, such that inward migration is driven by Lindblad torques.

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

We present two-dimensional hydrodynamic simulations of self-gravitating protostellar disks subject to axisymmetric infall from envelopes and irradiation from the central star, to explore disk fragmentation due to gravitational instability (GI), and the fragmented clump evolution. We assume that the disk is built gradually and smoothly by the infall, resulting in good numerical convergence. We confirm that for disks around solar-mass stars, infall at high rates at radii beyond ~50 AU leads to disk fragmentation. At lower infall rates <1e-5 Msun/yr, however, irradiation suppresses fragmentation. We find that, once formed, the fragments or clumps migrate inward on typical type-I time scales of ~2e3 yr initially, but later migration deviates from the type-I time scale when the clump becomes more massive than the local disk mass, and/or when they starts to open gaps. As they migrate, the clumps accrete from the disk at a rate 1e-3 to 1e-1 MJupiter/yr, consistent with analytic estimates that assume a 1-2 Hill radii cross section. Most clumps can grow to their isolation masses >0.1 Msun quickly. The eventual fates of these clumps, however, diverges depending on the migration speed: 3 out of 13 clumps become massive enough (brown dwarf mass) to open gaps in the disk and essentially stop migrating; 4 are tidally destroyed during inward migration; 6 migrate across the inner simulated disk boundary. A simple analytic model for clump evolution is derived to explain these different fates. Overall, our results indicate that fast migration, accretion, and tidal destruction of the clumps pose challenges to the scenario of giant planet formation by GI in situ, but may provide a formation mechanism for close binary systems.

The migration of low-mass planets is driven by the differential Lindblad torque and the corotation torque in non-magnetic viscous models of protoplanetary discs. The corotation torque has recently received detailed attention as it may slow down, stall, or reverse migration. In laminar viscous disc models, the long-term evolution of the corotation torque is intimately related to viscous and thermal diffusion processes in the planet's horseshoe region. This paper examines the properties of the corotation torque in discs where MHD turbulence develops as a result of the magnetorotational instability, considering a weak initial toroidal magnetic field. We present results of 3D MHD simulations carried out with two different codes. Non-ideal MHD effects and the disc's vertical stratification are neglected, and locally isothermal disc models are considered. The running time-averaged torque exerted by the disc on a fixed planet is evaluated in three disc models. We first present results with an inner disc cavity (planet trap). As in viscous disc models, the planet is found to experience a positive running time-averaged torque over several hundred orbits, which highlights the existence of an unsaturated corotation torque maintained in the long term in MHD turbulent discs. Two disc models with initial power-law density and temperature profiles are also adopted, in which the time-averaged torque is found to be in decent agreement with its counterpart in laminar viscous disc models with similar viscosity at the planet location. Detailed analysis of the averaged torque density distributions indicates that the differential Lindblad torque takes very similar values in MHD turbulent and laminar viscous discs, and there exists an unsaturated corotation torque in MHD turbulent discs. This analysis also reveals the existence of an additional corotation torque in weakly magnetized discs.

(abridged) Accretion in protoplanetary discs is thought to be driven by [...] turbulence via the magnetorotational instability (MRI). Recent work has shown that a planetesimal swarm embedded in a fully turbulent disc is subject to strong excitation of the velocity dispersion, leading to collisional destruction of bodies with radii R_p < 100 km. Significant diffusion of planetesimal semimajor axes also arises, leading to large-scale spreading of the planetesimal population throughout the inner regions of the protoplanetary disc, in apparent contradiction of constraints provided by the distribution of asteroids within the asteroid belt. In this paper, we examine the dynamics of planetesimals embedded in vertically stratified turbulent discs, with and without dead zones. Our main aims are to examine the turbulent excitation of the velocity dispersion, and the radial diffusion, of planetesimals in these discs. We employ three dimensional MHD simulations [...], along with an equilibrium chemistry model [...] We find that planetesimals in fully turbulent discs develop large random velocities that will lead to collisional destruction/erosion for bodies with sizes below 100 km, and undergo radial diffusion on a scale \sim 2.5 au over a 5 Myr disc life time. But planetesimals in a dead zone experience a much reduced excitation of their random velocities, and equilibrium velocity dispersions lie between the disruption thresholds for weak and strong aggregates for sizes R_p < 100 km. We also find that radial diffusion occurs over a much reduced length scale \sim 0.25 au over the disc life time, this being consistent with solar system constraints. We conclude that planetesimal growth via mutual collisions between smaller bodies cannot occur in a fully turbulent disc. By contrast, a dead zone may provide a safe haven in which km-sized planetesimals can avoid mutual destruction through collisions.

Context. Abridged. Many stars are members of binary systems. During early phases when the stars are surrounded by discs, the binary orbit and disc midplane may be mutually inclined. The discs around T Tauri stars will become mildly warped and undergo solid body precession around the angular momentum vector of the binary system. It is unclear how planetesimals in such a disc will evolve and affect planet formation. Aims. We investigate the dynamics of planetesimals embedded in discs that are perturbed by a binary companion on a circular, inclined orbit. We examine collisional velocities of the planetesimals to determine when they can grow through accretion. We vary the binary inclination, binary separation, D, disc mass, and planetesimal radius. Our standard model has D=60 AU, inclination=45 deg, and a disc mass equivalent to the MMSN. Methods. We use a 3D hydrodynamics code to model the disc. Planetesimals are test particles which experience gas drag, the gravitational force of the disc, the companion star gravity. Planetesimal orbit crossing events are detected and used to estimate collisional velocities. Results. For binary systems with modest inclination (25 deg), disc gravity prevents planetesimal orbits from undergoing strong differential nodal precession (which occurs in absence of the disc), and forces planetesimals to precess with the disc on average. For bodies of different size the orbit planes become modestly mutually inclined, leading to collisional velocities that inhibit growth. For larger inclinations (45 degrees), the Kozai effect operates, leading to destructively large relative velocities. Conclusions. Planet formation via planetesimal accretion is difficult in an inclined binary system with parameters similar to those considered in this paper. For systems in which the Kozai mechanism operates, the prospects for forming planets are very remote.

We use 2D hydrodynamic simulations of viscous disks to examine whether dynamically-interacting multiple giant planets can explain the large gaps (spanning over one order of magnitude in radius) inferred for the transitional and pre-transitional disks around T Tauri stars. In the absence of inner disk dust depletion, we find that it requires three to four giant planets to open up large enough gaps to be consistent with inferences from spectral energy distributions, because the gap width is limited by the tendency of the planets to be driven together into 2:1 resonances. With very strong tidal torques and/or rapid planetary accretion, fewer planets can also generate a large cavity interior to the locally formed gap(s) by preventing outer disk material from moving in. In these cases, however, the reduction of surface density produces a corresponding reduction in the inner disk accretion rate onto the star; this makes it difficult to explain the observed accretion rates of the pre/transitional disks. We find that even with four planets in disks, additional substantial dust depletion is required to explain observed disk gaps/holes. Substantial dust settling and growth, with consequent significant reductions in optical depths, is inferred for typical T Tauri disks in any case, and an earlier history of dust growth is consistent with the hypothesis that pre/transitional disks are explained by the presence of giant planets. We conclude that the depths and widths of gaps, and disk accretion rates in pre/transitional disks cannot be reproduced by a planet-induced gap opening scenario alone. Significant dust depletion is also required within the gaps/holes. Order of magnitude estimates suggest the mass of small dust particles (<1 micron) relative to the gas must be depleted to 1e-5 -- 1e-2 of the interstellar medium value, implying a very efficient mechanism of small dust removal or dust growth.

(abridged) Angular momentum transport and accretion in protoplanetary discs are generally believed to be driven by MHD turbulence via the magneto-rotational instability (MRI). The dynamics of solid bodies embedded in such discs (dust grains, boulders, planetesimals and planets) may be strongly affected by the turbulence, such that the formation pathways for planetary systems are determined in part by the strength and spatial distribution of the turbulent flow. We examine the dynamics of planetesimals, with radii between 1m \^a 10 km, embedded in turbulent protoplanetary discs, using three dimensional MHD simulations. The planetesimals experience gas drag and stochastic gravitational forces due to the turbulent disc. We use, and compare the results from, local shearing box simulations and global models in this study. The main aims of this work are to examine: the growth, and possible saturation, of the velocity dispersion of embedded planetesimals as a function of their size and disc parameters; the rate of radial migration and diffusion of planetesimals; the conditions under which the results from shearing box and global simulations agree. We find good agreement between local and global simulations when shearing boxes of dimension 4H x 16H x 2H are used (H being the local scale height). The magnitude of the density fluctuations obtained is sensitive to the box size, due to the excitation and propagation of spiral density waves. This affects the stochastic forcing experienced by planetesimals. [...] Our models show that fully developed MHD turbulence in protoplanetary discs would have a destructive effect on embedded planetesimals. Relatively low levels of turbulence are required for traditional models of planetesimal accretion to operate, this being consistent with the existence of a dead zone in protoplanetary discs.

We present the results of hydrodynamic simulations of the growth and orbital evolution of giant planets embedded in a protoplanetary disk with a dead-zone. The aim is to examine to what extent the presence of a dead-zone affects the rates of mass accretion and migration for giant planets. We performed 3D numerical simulations using a grid-based hydrodynamics code. In these simulations of non-magnetised disks, the dead-zone is treated as a region where the vertical profile of the viscosity depends on the distance from the equatorial plane. We consider dead-zones with vertical sizes, H_dz, ranging from 0 to H_dz=2.3H, where H is the disk scale-height. For all models, the vertically integrated viscous stress, and the related mass flux through the disk, have the same value, such that the simulations test the dependence of planetary mass accretion and migration on the vertical distribution of the viscous stress. For each model, an embedded 30 earth-masses planet on a fixed circular orbit is allowed to accrete gas from the disk. Once the planet mass becomes equal to that of Saturn or Jupiter, we allow the planet orbit to evolve due to gravitational interaction with the disk. We find that the time scale over which a protoplanet grows to become a giant planet is essentially independent of the dead-zone size, and depends only on the total rate at which the disk viscously supplies material to the planet. For Saturn-mass planets, the migration rate depends only weakly on the size of the dead-zone for H_dz< 1.5H, but becomes slower when H_dz=2.3H. This effect is due to the desaturation of corotation torques which originate from residual material in the partial-gap region. For Jupiter-mass planets, there is a clear tendency for the migration to proceed more slowly as the size of the dead-zone increases.

With references to recent observational results on the nebula IC59, we applied our previously developed Smoothed Particle Hydrodynamics (SPH) code which is based on Radiative Driven Implosion (RDI) model, to investigate the possible formation mechanism for the observed morphological structures of differently shaped BRCs. The simulation results confirmed the existence of the 4th type morphology of BRCs -- type M BRC. We are able to find the necessary condition for the appearance of the type M BRC based on the fact that the simulated physical properties of the cloud are consistent with observations on IC59. More importantly, the prospect of RDI triggered star formation by RDI model in all of the observed type M BRCs is ruthlessly eliminated.

We aim to examine the detailed disc structure that arises in a misaligned binary system as a function of the disc aspect ratio h, viscosity parameter alpha, disc outer radius R, and binary inclination angle gamma_F. We also aim to examine the conditions that lead to an inclined disc being disrupted by strong differential precession. We use a grid-based hydrodynamic code to perform 3D simulations. This code has a relatively low numerical viscosity compared with the SPH schemes that have been used previously to study inclined discs. This allows the influence of viscosity on the disc evolution to be tightly controlled. We find that for thick discs (h=0.05) with low alpha, efficient warp communication in the discs allows them to precess as rigid bodies with very little warping or twisting. Such discs are observed to align with the binary orbit plane on the viscous evolution time. Thinner discs with higher viscosity, in which warp communication is less efficient, develop significant twists before achieving a state of rigid-body precession. Under the most extreme conditions we consider (h=0.01, alpha=0.005 and alpha=0.1), we find that discs can become broken or disrupted by strong differential precession. Discs that become highly twisted are observed to align with the binary orbit plane on timescales much shorter than the viscous timescale, possibly on the precession time. We find agreement with previous studies that show that thick discs with low viscosity experience mild warping and precess rigidly. We also find that as h is decreased substantially, discs may be disrupted by strong differential precession, but for disc thicknesses that are significantly less (h=0.01) than those found in previous studies (h=0.03).

The primary aim of this work is to examine the effect of parabolic stellar encounters on the evolution of a Jovian-mass giant planet forming within a protoplanetary disc. We consider the effect on both the mass accretion and the migration history as a function of encounter distance. We use a grid-based hydrodynamics code to perform 2D simulations of a system consisting of a giant planet embedded within a gaseous disc orbiting around a star, which is perturbed by a passing star on a prograde, parabolic orbit. The disc model extends out to 50 AU, and parabolic encounters are considered with impact parameters ranging from 100 - 250 AU. In agreement with previous work, we find that the disc is significantly tidally truncated for encounters < 150 AU, and the removal of angular momentum from the disc by the passing star causes a substantial inflow of gas through the disc. The gap formed by the embedded planet becomes flooded with gas, causing the gas accretion rate onto the planet to increase abruptly. Gas flow through the gap, and into the inner disc, causes the positive inner disc torques exerted on the planet to increase, resulting in a sustained period of outward migration. For weaker interactions, corresponding to an encounter distance of > 250 AU, we find that the planet-disc system experiences minimal perturbation. Our results indicate that stellar fly-bys in young clusters may significantly modify the masses and orbital parameters of giant planets forming within protostellar discs. Planets that undergo such encounters are expected to be more massive, and to orbit with larger semimajor axes, than planets in systems which have not experienced parabolic encounters.

The discovery of short-period Neptune-mass objects, now including the remarkable system HD69830 (Lovis et al. 2006) with three Neptune analogues, raises difficult questions about current formation models which may require a global treatment of the protoplanetary disc. Several formation scenarios have been proposed, where most combine the canonical oligarchic picture of core accretion with type I migration (e.g. Terquem & Papaloizou 2007) and planetary atmosphere physics (e.g. Alibert et al. 2006). To date, published studies have considered only a small number of progenitors at late times. This leaves unaddressed important questions about the global viability of the models. We seek to determine whether the most natural model -- namely, taking the canonical oligarchic picture of core accretion and introducing type I migration -- can succeed in forming objects of 10 Earth masses and more in the innermost parts of the disc. This problem is investigated using both traditional semianalytic methods for modelling oligarchic growth as well as a new parallel multi-zone N-body code designed specifically for treating planetary formation problems with large dynamic range (McNeil & Nelson 2009). We find that it is extremely difficult for oligarchic tidal migration models to reproduce the observed distribution. Even under many variations of the typical parameters, we form no objects of mass greater than 8 Earth masses. By comparison, it is relatively straightforward to form icy super-Earths. We conclude that either the initial conditions of the protoplanetary discs in short-period Neptune systems were substantially different from the standard disc models we used, or there is important physics yet to be understood.

We examine the effect of giant planet migration on the formation of inner terrestrial planet systems. We consider situations in which the giant planet halts migration at semi-major axes in the range 0.13 - 1.7 AU due to gas disk dispersal. An N-body code is employed that is linked to a viscous gas disk algorithm capable of simulating: gas loss via accretion onto the central star and photoevaporation; gap formation by the giant planet; type II migration of the giant; optional type I migration of protoplanets; gas drag on planetesimals. We find that most of the inner system planetary building blocks survive the passage of the giant planet, either by being shepherded inward or scattered into exterior orbits. Systems of one or more hot-Earths are predicted to form and remain interior to the giant planet, especially if type II migration has been limited, or where type I migration has affected protoplanetary dynamics. Habitable planets in low eccentricity warm-Jupiter systems appear possible if the giant planet makes a limited incursion into the outer regions of the habitable zone (HZ), or traverses its entire width and ceases migrating at a radial distance of less than half that of the HZ's inner edge. We conclude that Type II migration does not prevent terrestrial planet formation.

The aim of this paper is to study the vertical profile of small dust particles in protoplanetary discs in which angular momentum transport is due to MHD turbulence driven by the magnetorotational instability. We consider particle sizes that range from approximately 1 micron up to a few millimeters.We use a grid--based MHD code to perform global two-fluid simulations of turbulent protoplanetary discs which contain dust grains of various sizes. In quasi--steady state, the gravitational settling of dust particles is balanced by turbulent diffusion. Simple and standard models of this process fail to describe accurately the vertical profile of the dust density. The disagreement is larger for small dust particles (of a few microns in size), especially in the disc upper layers ($Z>3H$, where $H$ is the scale-height). Here there can be orders of magnitude in the disagreement between the simple model predictions and the simulation results. This is because MHD turbulence is not homogeneous in accretion discs, since velocity fluctuations increase significantly in the disc upper layer where a strongly magnetized corona develops. We provide an alternative model that gives a better fit to the simulations. In this model, dust particles are diffused away from the midplane by MHD turbulence, but the diffusion coefficient varies vertically and is everywhere proportional to the square of the local turbulent vertical velocity fluctuations. The spatial distribution of dust particles can be used to trace the properties of MHD turbulence in protoplanetary discs, such as the amplitude of the velocity fluctuations. In the future, detailed and direct comparison between numerical simulations and observations should prove a useful tool for constraining the properties of turbulence in protoplanetary discs.

We investigate the stability of those low-mass Trojan planets that form in a protoplanetary disc and subsequently accrete gas to become gas giants. We calculate their evolution before, during, and after gas disc dispersal. A two-dimensional hydrodynamics code combined with an N-body solver is used to evolve the system of disc and planets. Gas disc dispersal is simulated in a simple manner by assuming global exponential decay of the disc mass, leading to the stalling of migration after semi-major axes have approximately halved from their initial values. We consider Trojan pairs with different initial masses and gas accretion rates and gas disc models with different masses and viscosities. An N-body code, adapted to model disc forces, is used to examine large-scale migration and the formation of very short period Trojan planets. For each combination of planetary pair and disc model that we consider in our hydrodynamic simulations, each Trojan system remains stable before, during, and after disc dispersal. The long-term stability of these systems in the absence of gas is tested using N-body simulations, and all systems remain stable for those evolution times equal to 10^9 years. Eccentricities remain low (e<0.02) in all cases. Increases in the amplitude of libration about the L4/L5 Lagrange points accompany the inward migration, and during very large-scale migration Trojan systems may be disrupted prior to the onset of disc dispersal. The stability of Trojan pairs during rapid type I migration, during the transition to type II migration with the accompanying gap formation in the gas disc, and during gas loss when the disc disperses, indicates that isolated Trojan planet systems are very stable. If a common mechanism exists for their formation, we suggest they may be readily observed in nature.

Protoplanet eccentricities of e >~ H/r can slow or reverse migration, but previous 2D studies have shown that gravitational scattering cannot maintain significant planet eccentricities against disc-induced damping. We simulate the evolution of low-mass protoplanetary swarms in three dimensions. The aim is to examine both protoplanet survival rates and the dynamical structure of the resulting planetary systems, and to compare them with 2D simulations. We present results from a 3D hydrodynamic simulation of eight protoplanets embedded in a protoplanetary disc. We also present a suite of simulations performed using an N-body code, modified to include prescriptions for planetary migration and for eccentricity and inclination damping. These prescriptions were obtained by fitting analytic formulae to hydrodynamic simulations of planets embedded in discs with initially eccentric and/or inclined orbits. As was found in two dimensions, differential migration produces groups of protoplanets in stable, multiple mean-motion resonances that migrate in lockstep, preventing prolonged periods of gravitational scattering. In almost all simulations, this leads to large-scale migration of the protoplanet swarm into the central star in the absence of a viable stopping mechanism. The evolution involves mutual collisions, occasional instances of large-scale scattering, and the frequent formation of the long-lived, co-orbital planet systems that arise in > 30% of all runs. Disc-induced damping overwhelms eccentricity and inclination growth due to planet-planet interactions. Co-orbital planets are a natural outcome of dynamical relaxation in a strongly dissipative environment, and if observed in nature would imply that such a period of evolution commonly arises during planetary formation.

Modern N-body techniques for planetary dynamics are generally based on symplectic algorithms specially adapted to the Kepler problem. These methods have proven very useful in studying planet formation, but typically require the timestep for all objects to be set to a small fraction of the orbital period of the innermost body. This computational expense can be prohibitive for even moderate particle number for many physically interesting scenarios, such as recent models of the formation of hot exoplanets, in which the semimajor axis of possible progenitors can vary by orders of magnitude. We present new methods which retain most of the benefits of the standard symplectic integrators but allow for radial zones with distinct timesteps. These approaches should make simulations of planetary accretion with large dynamical range tractable. As proof of concept we present preliminary science results from an implementation of the algorithm as applied to an oligarchic migration scenario for forming hot Neptunes.

A new Radiative Driven Implosion (RDI) model based on Smoothed Particle Hydrodynamics (SPH) technique is developed and applied to investigate the morphological evolutions of molecular clouds under the effect of ionising radiation. This model self-consistently includes the self-gravity of the cloud in the hydrodynamical evolution, the UV radiation component in the radiation transfer equations, the relevant heating and cooling mechanisms in the energy evolution and a comprehensive chemical network. The simulation results reveal that under the effect of ionising radiation, a molecular cloud may evolve through different evolutionary sequences. Dependent on its initial gravitational state, the evolution of a molecular cloud does not necessarily follow a complete morphological evolution sequence from type A to B to C, as described by previous RDI models. When confronted with observations, the simulation results provide satisfactory physical explanations for a series of puzzles derived from Bright-Rimmed Clouds (BRCs) observations. The consistency of the modelling results with observations shows that the self-gravity of a molecular cloud should not be neglected in any investigation on the dynamical evolution of molecular clouds when they are exposed to ionising radiation.

Over 30 planetary systems have been discovered to reside in binary stars. For small separations gravitational perturbation of the secondary star has a strong influence on the planet formation process. It truncates the protoplanetary disk, may shortens its lifetime, and stirs up the embedded planetesimals. Due to its small semi-major axis (18.5 AU) and large eccentricity (e=0.35) the binary $\gamma$ Cephei represents a particularly challenging example. In the present study we model the orbital evolution and growth of embedded protoplanetary cores of about 30 earth masses in the putative protoplanetary disk surrounding the primary star in the $\gamma$ Cep system. We assume coplanarity of the disk, binary and planet and perform two-dimensional hydrodynamic simulations of embedded cores in a protoplanetary disk. The presence of the eccentric secondary star perturbs the disk periodically and generates strong spiral arms at periapse which propagate toward the disk centre. The disk also becomes slightly eccentric (with e_d = 0.1-0.15), and displays a slow retrograde precession in the inertial frame. For all initial separations (2.5 to 3.5 AU) we find inward migration of the cores. For initial semi-major axes (a_p \gsim 2.7), we find a strong increase in the planetary eccentricity despite the presence of inward migration. Only cores which are initially far from the disk outer edge have a bounded orbital eccentricity which converges, roughly to the value of the planet observed in the $\gamma$ Cep system. We have shown that under the condition protoplanetary cores can form at around 2.5 AU, it is possible to evolve and grow such a core to form a planet with final outcome similar to that observed.

We present the results of hydrodynamic simulations of the formation and subsequent orbital evolution of giant planets embedded in a circumbinary disc. We assume that a 20 earth masses core has migrated to the edge of the inner cavity formed by the binary where it remains trapped by corotation torques. This core is then allowed to accrete gas from the disc, and we study its orbital evolution as it grows in mass. For each of the two accretion time scales we considered, we performed three simulations. In two of the three simulations, we stop the accretion onto the planet once its mass becomes characteristic of that of Saturn or Jupiter. In the remaining case, the planet can accrete disc material freely in such a way that its mass becomes higher than Jupiter's. The simulations show different outcomes depending on the final mass m_p of the giant. For m_p=1 M_S (where M_S is Saturn's mass), we find that the planet migrates inward through its interaction with the disc until its eccentricity becomes high enough to induce a torque reversal. The planet then migrates outward, and the system remains stable on long time scales. For m_p > 1 M_J (where M_J is Jupiter's mass) we observed two different outcomes. In each case the planet enters the 4:1 resonance with the binary, and resonant interaction drives up the eccentricity of the planet until it undergoes a close encounter with the secondary star. The result can either be ejection from the system or scattering out into the disc followed by a prolonged period of outward migration. These results suggest that circumbinary planets are more likely to be quite common in the Saturn-mass range. Jupiter-mass circumbinary planets are likely to be less common because of their less stable evolution.

Protostellar accretion discs have cool, dense midplanes where externally originating ionisation sources such as X-rays or cosmic rays are unable to penetrate. This suggests that for a wide range of radii, MHD turbulence can only be sustained in the surface layers where the ionisation fraction is sufficiently high. A dead zone is expected to exist near the midplane, such that active accretion only occurs near the upper and lower disc surfaces. Recent work, however, suggests that under suitable conditions the dead zone may be enlivened by turbulent transport of ions from the surface layers into the dense interior. In this paper we present a suite of simulations that examine where, and under which conditions, a dead zone can be enlivened by turbulent mixing. We use three-dimensional, multifluid shearing box MHD simulations, which include vertical stratification, ionisation chemistry, ohmic resistivity, and ionisation due to X-rays from the central protostar. We compare the results of the MHD simulations with a simple reaction-diffusion model. The simulations show that in the absence of gas-phase heavy metals, such as magnesium, turbulent mixing has essentially no effect on the dead zone. The addition of a relatively low abundance of magnesium, however, increases the recombination time and allows turbulent mixing of ions to enliven the dead zone completely beyond a distance of 5 AU from the central star, for our particular disc model. During the late stages of protoplanetary disc evolution, when small grains have been depleted and the disc surface density has decreased below its high initial value, the structure of the dead zone may be significantly altered by the action of turbulent transport.

We investigate the evolution of two-planet systems embedded in a protoplanetary disc, which are composed of a Jupiter-mass planet plus another body located further out in the disc. We consider outermost planets with masses ranging from 10 earth masses to 1 M_J. We also examine the case of outermost bodies with masses < 10 earth masses (M_E). Differential migration of the planets due to disc torques leads to different evolution outcomes depending on the mass of the outer protoplanet. For planets with mass < 3.5 M_E the type II migration rate of the giant exceeds the type I migration rate of the outer body, resulting in divergent migration. Outer bodies with masses in the range 3.5 < m_o < 20 M_E become trapped at the edge of the gap formed by the giant planet, because of corotation torques. Higher mass planets are captured into resonance with the inner planet. If 30 < m_o < 40 M_E or m_o=1 M_J, then the 2:1 resonance is established. If 80 < m_o < 100 M_E, the 3:2 resonance is favoured. Simulations of gas-accreting protoplanets of mass m_o > 20 M_E, trapped initially at the edge of the gap, or in the 2:1 resonance, also result in eventual capture in the 3:2 resonance as the planet mass grows to become close to the mass of Saturn. Our results suggest that there is a theoretical lower limit to the mass of an outer planet that can be captured into resonance with an inner Jovian planet, which is relevant to observations of extrasolar multiplanet systems. Furthermore, capture of a Saturn-like planet into the 3:2 resonance with a Jupiter-like planet is a very robust outcome of simulations. This result is relevant to recent scenarios of early Solar System evolution which require Saturn to have existed interior to the 2:1 resonance with Jupiter prior to the onset of the Late Heavy Bombardment.

Previous work has shown that the tidal interaction between a binary system and a circumbinary disc leads to the formation of a large inner cavity in the disc. Subsequent formation and inward migration of a low mass planet causes it to become trapped at the cavity edge, where it orbits until further mass growth or disc dispersal. The question of how systems of multiple planets in circumbinary discs evolve has not yet been addressed. Here, we present a suite of simulations which study the evolution of pairs of planets migrating in the disc. We also present the results of hydrodynamic simulations of five-planet systems, and study their long term evolution after disc dispersal using a N-body code. For the two-planet simulations we assume that the innermost planet has migrated to the edge of the inner cavity and remains trapped there, and study the subsequent evolution of the system as the outermost planet migrates inward. We find that the outcomes largely depend on the mass ratio q=m_i/m_o, where m_i (m_o) is the mass of the innermost (outermost) planet. For q<1, planets usually undergo dynamical scattering or orbital exchange. For values of q>1 the systems reach equilibrium configurations in which the planets are locked into mean motion resonances, and remain trapped at the edge of the inner cavity without further migration. Most simulations of five-planet systems we performed resulted in collisions and scattering events, such that only a single planet remained in orbit about the binary. In one case however, a multiplanet resonant system was found to be dynamically stable over long time scales, suggesting that such systems may be observed in planet searches focussed on close binary systems.

A well-known bottleneck for the core-accretion model of giant-planet formation is the loss of the cores into the star by Type-I migration, due to the tidal interactions with the gas disk. It has been shown that a steep surface-density gradient in the disk, such as the one expected at the boundary between an active and a dead zone, acts as a planet trap and prevents isolated cores from migrating down to the central star. We study the relevance of the planet trap concept for the accretion and evolution of systems of multiple planetary embryos/cores. We performed hydrodynamical simulations of the evolution of systems of multiple massive objects in the vicinity of a planet trap. The planetary embryos evolve in 3 dimensions, whereas the disk is modeled with a 2D grid. Synthetic forces are applied onto the embryos to mimic the damping effect that the disk has on their inclinations. Systems with two embryos tend to acquire stable, separated and non-migrating orbits, with the more massive embryo placed at the planet trap and the lighter one farther out in the disk. Systems of multiple embryos are intrinsically unstable. Consequently, a long phase of mutual scattering can lead to accreting collisions among embryos; some embryos are injected into the inner part of the disk, where they can be evacuated into the star by Type I migration. The system can resume a stable, non-migrating configuration only when the number of surviving embryos decreases to a small value (~2-4). This can explain the limited number of giant planets in our solar system. These results should apply in general to any case in which the Type-I migration of the inner embryo is prevented by some mechanism, and not solely to the planet trap scenario.

Models of terrestrial planet formation in the presence of a migrating giant planet have challenged the notion that hot-Jupiter systems lack terrestrial planets. We briefly review this issue and suggest that hot-Jupiter systems should be prime targets for future observational missions designed to detect Earth-sized and potentially habitable worlds.

Context: Our previous models of a giant planet migrating through an inner protoplanet/planetesimal disk find that the giant shepherds a portion of the material it encounters into interior orbits, whilst scattering the rest into external orbits. Scattering tends to dominate, leaving behind abundant material that can accrete into terrestrial planets. Aims: We add to the possible realism of our model by simulating type I migration forces which cause an inward drift, and strong eccentricity and inclination damping of protoplanetary bodies. This extra dissipation might be expected to enhance shepherding at the expense of scattering, possibly modifying our previous conclusions. Methods: We employ an N-body code that is linked to a viscous gas disk algorithm capable of simulating: gas accretion onto the central star; gap formation in the vicinity of the giant planet; type II migration of the giant planet; type I migration of protoplanets; and the effect of gas drag on planetesimals. We use the code to re-run three scenarios from a previous work where type I migration was not included. Results: The additional dissipation introduced by type I migration enhances the inward shepherding of material but does not severely reduce scattering. We find that > 50% of the solids disk material still survives the migration in scattered exterior orbits: most of it well placed to complete terrestrial planet formation at < 3 AU. The shepherded portion of the disk accretes into hot-Earths, which survive in interior orbits for the duration of our simulations. Conclusions: Water-rich terrestrial planets can form in the habitable zones of hot-Jupiter systems and hot-Earths and hot-Neptunes may also be present. These systems should be targets of future planet search missions.

We present the results of hydrodynamical simulations of low mass protoplanets embedded in circumbinary accretion disks. The aim is to examine the migration and long term orbital evolution of the protoplanets, in order to establish the stability properties of planets that form in circumbinary disks. Simulations were performed using a grid--based hydrodynamics code. First we present a set of calculations that study how a binary interacts with a circumbinary disk. We evolve the system for 10^5 binary orbits, which is the time needed for the system to reach a quasi-equilibrium state. From this time onward the apsidal lines of the disk and the binary are aligned, and the binary eccentricity remains essentially unchanged with a value of e_b ~ 0.08. Once this stationary state is obtained, we embed a low mass protoplanet in the disk and let it evolve under the action of the binary and disk forces. We consider protoplanets with masses of 5, 10 and 20 Earth masses. In each case, we find that inward migration of the protoplanet is stopped at the edge of the tidally truncated cavity formed by the binary. This effect is due to positve corotation torques, which can counterbalance the net negative Lindblad torques in disk regions where the surface density profile has a sufficiently large positive gradient. Halting of migration occurs in a region of long-term stability, suggesting that low mass circumbinary planets may be common, and that gas giant circumbinary planets should be able to form in circumbinary disks.

Young planets embedded in their protoplanetary disk interact gravitationally with it leading to energy and angular momentum exchange. This interaction determines the evolution of the planet through changes to the orbital parameters. We investigate changes in the orbital elements of a 20 Earth--mass planet due to the torques from the disk. We focus on the non-linear evolution of initially non-vanishing eccentricity $e$ and/or inclination $i$. We treat the disk as a two- or three-dimensional viscous fluid and perform hydrodynamical simulations with an embedded planet. We find rapid exponential decay of the planet orbital eccentricity and inclination for small initial values of $e$ and $i$, in agreement with linear theory. For larger values of $e > 0.1$ the decay time increases and the decay rate scales as $\dot{e} \propto e^{-2}$, consistent with existing theoretical models. For large inclinations ($i$ > 6 deg) the inclination decay rate shows an identical scaling $di/dt \propto i^{-2}$. We find an interesting dependence of the migration on the eccentricity. In a disk with aspect ratio $H/r=0.05$ the migration rate is enhanced for small non-zero eccentricities ($e < 0.1$), while for larger values we see a significant reduction by a factor of $\sim 4$. We find no indication for a reversal of the migration for large $e$, although the torque experienced by the planet becomes positive when $e \simeq 0.3$. This inward migration is caused by the persisting energy loss of the planet. For non gap forming planets, eccentricity and inclination damping occurs on a time scale that is very much shorter than the migration time scale. The results of non linear hydrodynamic simulations are in very good agreement with linear theory for small $e$ and $i$.

Among the extrasolar planetary systems about 30 are located in a stellar binary orbiting one of the stars, preferably the more massive primary. The dynamical influence of the second companion alters firstly the orbital elements of the forming protoplanet directly and secondly the structure of the disk from which the planet formed which in turn will modify the planet's evolution. We present detailed analysis of these effects and present new hydrodynamical simulations of the evolution of protoplanets embedded in circumstellar disks in the presence of a companion star, and compare our results to the system $\gamma$ Cep. To analyse the early formation of planetary embryos, we follow the evolution of a swarm of planetesimals embedded in a circumstellar disk. Finally, we study the evolution of planets embedded in circumbinary disks.

We present a series of calculations aimed at examining how an inner system of planetesimals/protoplanets, undergoing terrestrial planet formation, evolves under the influence of a giant planet undergoing inward type II migration through the region bounded between 5 - 0.1 AU. We find that > 60% of the solids disk survives by being scattered by the giant planet into external orbits. Planetesimals are scattered outward almost as efficiently as protoplanets, resulting in the regeneration of a solids disk where dynamical friction is strong and terrestrial planet formation is able to resume. A simulation extended for a few Myr after the migration of the giant planet halted at 0.1 AU, resulted in an apparently stable planet of ~ 2 Earth masses forming in the habitable zone. Migration-induced mixing of volatile-rich material from beyond the `snowline' into the inner disk regions means that terrestrial planets that form there are likely to be water-rich. We predict that hot--Jupiter systems are likely to harbor water-rich terrestrial planets in their habitable zones. These planets may be detected by future planet search missions.

We present the results of global 3-D MHD simulations of stratified and turbulent protoplanetary disc models. The aim of this work is to develop thin disc models capable of sustaining turbulence for long run times, which can be used for on-going studies of planet formation in turbulent discs. The results are obtained using two codes written in spherical coordinates: GLOBAL and NIRVANA. Both are time--explicit and use finite differences along with the Constrained Transport algorithm to evolve the equations of MHD. In the presence of a weak toroidal magnetic field, a thin protoplanetary disc in hydrostatic equilibrium is destabilised by the magnetorotational instability (MRI). When the resolution is large enough (25 vertical grid cells per scale height), the entire disc settles into a turbulent quasi steady-state after about 300 orbits. Angular momentum is transported outward such that the standard alpha parameter is roughly 4-6*10^{-3}. We find that the initial toroidal flux is expelled from the disc midplane and that the disc behaves essentially as a quasi-zero net flux disc for the remainder of the simulation. As in previous studies, the disc develops a dual structure composed of an MRI--driven turbulent core around its midplane, and a magnetised corona stable to the MRI near its surface. By varying disc parameters and boundary conditions, we show that these basic properties of the models are robust. The high resolution disc models we present in this paper achieve a quasi--steady state and sustain turbulence for hundreds of orbits. As such, they are ideally suited to the study of outstanding problems in planet formation such as disc--planet interactions and dust dynamics.

We perform numerical simulations of a disc-planet system using various grid-based and smoothed particle hydrodynamics (SPH) codes. The tests are run for a simple setup where Jupiter and Neptune mass planets on a circular orbit open a gap in a protoplanetary disc during a few hundred orbital periods. We compare the surface density contours, potential vorticity and smoothed radial profiles at several times. The disc mass and gravitational torque time evolution are analyzed with high temporal resolution. There is overall consistency between the codes. The density profiles agree within about 5% for the Eulerian simulations while the SPH results predict the correct shape of the gap although have less resolution in the low density regions and weaker planetary wakes. The disc masses after 200 orbital periods agree within 10%. The spread is larger in the tidal torques acting on the planet which agree within a factor 2 at the end of the simulation. In the Neptune case the dispersion in the torques is greater than for Jupiter, possibly owing to the contribution from the not completely cleared region close to the planet.

Context. Recent observations of the X-ray emission from T Tauri stars in the Orion nebula have shown that they undergo frequent outbursts in their X-ray luminosity. These X-ray flares are characterised by increases in luminosity by two orders of magnitude, a typical duration of less than one day, and a significant hardening of the X-ray spectrum. Aims. It is unknown what effect these X-ray flares will have on the ionisation fraction and dead-zone structure in protoplanetary disks. We present the results of calculations designed to address this question. Methods. We have performed calculations of the ionisation fraction in a standard $\alpha$-disk model using two different chemical reaction networks. We include in our models ionisation due to X-rays from the central star, and calculate the time-dependent ionisation fraction and dead--zone structure for the inner 10 AU of a protoplanetary disk model. Results. We find that the disk response to X-ray flares depends on whether the plasma temperature increases during flares and/or whether heavy metals (such as magnesium) are present in the gas phase. Under favourable conditions the outer disk dead--zone can disappear altogether,and the dead-zone located between 0.5 < R < 2 AU can disappear and reappear in phase with the X-ray luminosity. Conclusions. X-ray flares can have a significant effect on the dead-zone structure in protoplanetary disks. Caution is required in interpreting this result as the duration of X-ray bursts is considerably shorter than the growth time of MHD turbulence due to the magnetorotational instability.

A three dimensional Smoothed Particle Hydrodynamical (SPH) model has been extended to study the radiative driven implosion effect of massive stars on the dynamical evolutions of surrounding molecular clouds. The new elements in the upgraded code are the inclusion of Lyman continuum in the incident radiation flux and the treatment of hydrogen ionisation process; introducing ionisation heating & recombination cooling effects; and adding a proper description of the magnetic and turbulent pressures to the internal pressure of the molecular cloud. The application of this newly developed model to the structure of the middle Eagle Nebula finger suggests that the shock induced by the ionising radiation at the front side of the head precedes an ionisation front moving towards the center of the core, and that the core at the fingertip is at transition stage evolving toward a state of induced star formation. The dynamical evolution of the velocity field of the simulated cloud structure is discussed to illustrate the role of the self-gravity and the different cloud morphologies which appear at different stages in the evolutionary process of the cloud. The modelled gas evaporation rate is consistent with that of current other models and the density, temperature and chemical profiles are agreement with the observed values. The relative lifetimes of different simulated cloud morphologies suggests a possible answer to the question of why more bright-rimmed clouds are observed to possess a flat-core than an elongated-core morphology. [Abridged]

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.

We calculate the ionisation fraction in protostellar disk models using a number of different chemical reaction networks, including gas-phase and gas-grain reaction schemes. The disk models we consider are conventional alpha-disks, which include viscous heating and radiative cooling. The primary source of ionisation is assumed to be X-ray irradiation from the central star. We consider a number of gas-phase chemical networks. In general we find that the simple models predict higher fractional ionisation levels and more extensive active zones than the more complex models. When heavy metal atoms are included the simple models predict that the disk is magnetically active throughout. The complex models predict that extensive regions of the disk remain magnetically uncoupled even with a fractional abundance of magnesium of 10(-8). The addition of submicron sized grains with a concentration of 10(-12) causes the size of the dead zone to increase dramatically for all kinetic models considered. We find that the simple and complex gas-grain reaction schemes agree on the size and structure of the resulting dead zone. We examine the effects of depleting the concentration of small grains as a crude means of modeling the growth of grains during planet formation. We find that a depletion factor of 10(-4) causes the gas-grain chemistry to converge to the gas-phase chemistry when heavy metals are absent. 10(-8) is required when magnesium is included. This suggests that efficient grain growth and settling will be required in protoplanetary disks, before a substantial fraction of the disk mass in the planet forming zone between 1 - 10 AU becomes magnetically active and turbulent.

We study the migration of solid bodies in turbulent protoplanetary accretion discs by means of global MHD simulations. The bodies range in size from 5 centimetres up to 1 metre, and so include objects whose migration is expected to be the most rapid due to gas drag interaction with the disc. As they drift inward through the disc, some of them are trapped in regions where gas pressure maxima are created by long lived anticyclonic vortices. This accumulation is very efficient, locally increasing the dust--to--gas ratio by a factor > 100 in some cases. We discuss the possible implications of this result for theories of planet formation.

We calculate the ionisation fraction in protostellar disk models using two different gas-phase chemical networks, and examine the effect of turbulent mixing by modelling the diffusion of chemical species vertically through the disk. The aim is to determine in which regions of the disk gas can couple to a magnetic field and sustain MHD turbulence. We find that the effect of diffusion depends crucially on the elemental abundance of heavy metals (magnesium) included in the chemical model. In the absence of heavy metals, diffusion has essentially no effect on the ionisation structure of the disks, as the recombination time scale is much shorter than the turbulent diffusion time scale. When metals are included with an elemental abundance above a threshold value, the diffusion can dramatically reduce the size of the magnetically decoupled region, or even remove it altogther. For a complex chemistry the elemental abundance of magnesium required to remove the dead zone is 10(-10) - 10(-8). We also find that diffusion can modify the reaction pathways, giving rise to dominant species when diffusion is switched on that are minor species when diffusion is absent. This suggests that there may be chemical signatures of diffusive mixing that could be used to indirectly detect turbulent activity in protoplanetary disks. We find examples of models in which the dead zone in the outer disk region is rendered deeper when diffusion is switched on. Overall these results suggest that global MHD turbulence in protoplanetary disks may be self-sustaining under favourable circumstances, as turbulent mixing can help maintain the ionisation fraction above that necessary to ensure good coupling between the gas and magnetic field.

Using 2D MHD numerical simulations performed with two different finite difference Eulerian codes, we analyze the effect that a toroidal magnetic field has on low mass planet migration in nonturbulent protoplanetary discs. The presence of the magnetic field modifies the waves that can propagate in the disc. In agreement with a recent linear analysis (Terquem 2003), we find that two magnetic resonances develop on both sides of the planet orbit, which contribute to a significant global torque. In order to measure the torque exerted by the disc on the planet, we perform simulations in which the latter is either fixed on a circular orbit or allowed to migrate. For a 5 earth mass planet, when the ratio \beta between the square of the sound speed and that of the Alfven speed at the location of the planet is equal to 2, we find inward migration when the magnetic field B_{\phi} is uniform in the disc, reduced migration when B_{\phi} decreases as r^{-1} and outward migration when B_{\phi} decreases as r^{-2}. These results are in agreement with predictions from the linear analysis. Taken as a whole, our results confirm that even a subthermal stable field can stop inward migration of an earth--like planet.

(Abridged).We present the results of MHD simulations of low mass protoplanets interacting with turbulent disks. We calculate the orbital evolution of `planetesimals' and protoplanets with masses in the range 0 < m_p < 30 M_Earth. Planetesimals and protoplanets undergo stochastic migration due to interaction with turbulent density fluctuations. Over run times of ~ 150 planet orbits, stochastic migration dominates over type I migration for many models. Fourier analysis of the torques experienced by planets indicates that the torque fluctuations contain components with significant power whose time scales of variation are similar to the simulation run times. These low frequency fluctuations partly explain the dominance of stochastic torques, and may provide a powerful means of counteracting the type I migration of some planets in turbulent disks. Turbulence is a source of eccentricity driving. Planetesimals attained eccentricities in the range 0.02 < e < 0.14, m_p=1 M_Earth planets attained eccentricities 0.02 < e < 0.08, and m_p=10 M_Earth protoplanets reached 0.02 < e < 0.03. This is in basic agreement with a model in which turbulence drives e-growth, and interaction with disk material at coorbital Lindblad resonances causes e-damping. These results are significant for planet formation. Stochastic migration may prevent some planet cores migrating into their star via type I before becoming gas giants. The growth of planetary cores may be enhanced by preventing isolation. Eccentricity excitation by turbulence, however, may reduce growth rates of planetary cores during the runaway and oligarchic growth stages, and cause collisions between planetesimals to become destructive.

We present the results of N--body simulations which examine the effect that gas giant planet migration has on the formation of terrestrial planets. The models incorporate a 0.5 Jupiter mass planet undergoing type II migration through an inner protoplanet--planetesimal disk, with gas drag included. Each model is initiated with the inner disk being at successively increased levels of maturity, so that it is undergoing either oligarchic or giant impact style growth as the gas giant migrates. In all cases, a large fraction of the disk mass survives the passage of the giant, either by accreting into massive terrestrial planets shepherded inward of the giant, or by being scattered into external orbits. Shepherding is favored in younger disks where there is strong dynamical friction from planetesimals and gas drag is more influential, whereas scattering dominates in more mature disks where dissipation is weaker. In each scenario, sufficient mass is scattered outward to provide for the eventual accretion of a set of terrestrial planets in external orbits, including within the system's habitable zone. An interesting result is the generation of massive, short period, terrestrial planets from compacted material pushed ahead of the giant. These planets are reminiscent of the short period Neptune mass planets discovered recently, suggesting that such `hot Neptunes' could form locally as a by-product of giant planet migration.

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

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.

We present new spectroscopic and photometric data on the early-type overcontact binary TU Muscae. The analysis of the spectroscopic data shows that the line of sight to the system crosses three kinematically sharp and well-separated interstellar reddening sources and that the stars rotate synchronously. We present new radial velocities that are in good agreement with earlier optical velocities and, thus, do not confirm the systematically smaller velocities obtained from IUE spectra. The optical velocities are analyzed simultaneously with the photometric data to derive accurate absolute dimesions for the binary components.The results show that TU Mus consists of an O7.5 primary with M_1=23.5 +/- 0.8 M_sun, R_1=7.48 +/- 0.08 R_sun and an O9.5 secondary with M_2=15.3 +/- 0.4 M_sun, R_2=6.15 +/- 0.07 R_sun in an overcontact configuration and that the orbital period has remained constant over the three decades covered by the observations. These results might imply that the mass transfer seen in late-type overcontact binaries does not occur in their early-type counterparts.

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

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

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

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

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

We present near-infrared medium-resolution ($R \sim 875$) spectra of the microquasar GRS 1915+105 on 1997 August 13-15 UTC from the Hale 200-inch telescope. The spectra showed broad emission lines of He I (2.058 $\mu$m) and H I (2.166 $\mu$m - Br$\gamma$), consistent with previous work. On August 14 UTC, we took spectra with $\sim 6$-minute time resolution during infrared flaring events similar to those reported in Eikenberry et al. (1998a), which appear to reveal plasma ejection from the system. During the flares, the emission line fluxes varied in approximately linear proportionality to the IR continuum flux, implying that the lines are radiatively pumped by the flares. We also detected a weak He II (2.189 $\mu$m) emission line on August 14 UTC. The nature of the line variability and the presence of the He II feature indicate that the emission lines in GRS 1915+105 arise in an accretion disk around the compact object, rather than in the circumstellar disk of a proposed Oe/Be companion. The radiative line pumping also implies that the flare emission originates from ejecta which have moved out of the accretion disk plane.

We study the vertical heating and thickening of galaxy disks due to accretion of small satellites. Our simulations are restricted to axial symmetry, which largely eliminates numerical evolution of the target galaxy but requires the trajectory of the satellite to be along the symmetry axis of the target. We find that direct heating of disk stars by the satellite is not important because the satellite's gravitational perturbation has little power at frequencies resonant with the vertical stellar orbits. The satellite does little damage to the disk until its decaying orbit resonantly excites large-scale disk bending waves. Bending waves can damp through dynamical friction from the halo or internal wave-particle resonances; we find that wave-particle resonances dominate the damping. The principal vertical heating mechanism is therefore dissipation of bending waves at resonances with stellar orbits in the disk. Energy can thus be deposited some distance from the point of impact of the satellite. The net heating from a tightly bound satellite can be substantial, but satellites that are tidally disrupted before they are able to excite bending waves do not thicken the disk.

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.

We report simultaneous X-ray and infrared (IR) observations of the Galactic microquasar GRS1915+105 using XTE and the Palomar 200-inch telescope on August 13-15, 1997 UTC. During the last two nights, the microquasar GRS 1915+105 exhibited quasi-regular X-ray/infrared (IR) flares with a spacing of $\sim 30$ minutes. While the physical mechanism triggering the flares is currently unknown, the one-to-one correspondence and consistent time offset between the X-ray and IR flares establish a close link between the two. At late times in the flares the X-ray and IR bands appear to ``decouple'', with the X-ray band showing large-amplitude fast oscillations while the IR shows a much smoother, more symmetrical decline. In at least three cases, the IR flare has returned to near its minimum while the X-rays continue in the elevated oscillatory state, ruling out thermal reprocessing of the X-ray flux as the source of IR flare. Furthermore, observations of similar IR and radio flares by Fender et al. (1997) imply that the source of the IR flux in such flares is synchrotron emission. The common rise and subsequent decoupling of the X-ray and IR flux and probable synchrotron origin of the IR emission is consistent with a scenario wherein the IR flux originates in a relativistic plasma which has been ejected from the inner accretion disk. In that case, these simultaneous X-ray/IR flares from a black-hole/relativistic-jet system are the first clear observational evidence linking of the time-dependent interaction of the jet and the inner disk in decades of quasar and microquasar studies.

Dramatic torque reversals between spin up and spin down have been observed in half of the persistent X-ray pulsars monitored by the BATSE all-sky monitor on CGRO. Theoretical models developed to explain early pulsar timing data can explain spin down torques via a disk-magnetosphere interaction if the star nearly corotates with the inner accretion disk. To produce the observed BATSE torque reversals, however, these equilibrium models require the disk to alternate between two mass accretion rates, with $\dot M_{\pm}$ producing accretion torques of similar magnitude, but always of opposite sign. Moreover, in at least one pulsar (GX 1+4) undergoing secular spin down the neutron star spins down faster during brief ($\sim 20$ day) hard X-ray flares -- this is opposite the correlation expected from standard theory, assuming BATSE pulsed flux increases with mass accretion rate. The $10$ day to 10 yr intervals between torque reversals in these systems are much longer than any characteristic magnetic or viscous time scale near the inner disk boundary and are more suggestive of a global disk phenomenon. We discuss possible explanations of the observed torque behavior. Despite the preferred sense of rotation defined by the binary orbit, the BATSE observations are surprisingly consistent with an earlier suggestion by Makishima \etal (1988) for GX~1+4: the disks in these systems somehow alternate between episodes of prograde and retrograde rotation. We are unaware of any mechanism that could produce a stable retrograde disk in a binary undergoing Roche-lobe overflow, but such flip-flop behavior does occur in numerical simulations of wind-fed systems. One possibility is that the disks in some of these binaries are fed by an X-ray excited wind.

We summarize five years of continuous monitoring of accretion-powered pulsars with the Burst and Transient Source Experiment (BATSE) on the Compton Gamma Ray Observatory. Our 20-70 keV observations have determined or refined the orbital parameters of 13 binaries, discovered 5 new transient accreting pulsars, measured the pulsed flux history during outbursts of 12 transients (GRO J1744-28, 4U 0115+634, GRO J1750-27, GS 0834-430, 2S 1417-624, GRO J1948+32, EXO 2030+375, GRO J1008-57, A 0535+26, GRO J2058+42, 4U 1145-619 and A 1118-616), and also measured the accretion torque history of during outbursts of 6 of those transients whose orbital parameters were also known. We have also continuously measured the pulsed flux and spin frequency for eight persistently accreting pulsars (Her X-1, Cen X-3, Vela X-1, OAO 1657-415, GX 301-2, 4U 1626-67, 4U 1538-52, and GX 1+4). Because of their continuity and uniformity over a long baseline, BATSE observations have provided new insights into the long-term behavior of accreting magnetic stars. We have found that all accreting pulsars show stochastic variations in their spin frequencies and luminosities, including those displaying secular spin-up or spin-down on long time scales, blurring the conventional distinction between disk-fed and wind-fed binaries. Pulsed flux and accretion torque are strongly correlated in outbursts of transient accreting pulsars, but uncorrelated, or even anticorrelated, in persistent sources.

We apply linear response theory to a general, inhomogeneous, stationary stellar system, with particular emphasis on dissipative processes analogous to Landau damping. Assuming only that the response is causal, we show that the irreversible work done by an external perturber is described by the anti-Hermitian part of a linear response operator, and damping of collective modes is described by the anti-Hermitian part of a related polarization operator. We derive an exact formal expression for the response operator, which is the classical analog of a well-known result in quantum statistical physics. When the self-gravity of the response can be ignored, and the ensemble-averaged gravitational potential is integrable, the expressions for the mode energy, damping rate, and polarization operator reduce to well-known formulae derived from perturbation theory in action-angle variables. In this approximation, dissipation occurs only via resonant interaction with stellar orbits or collective modes. For stellar systems in thermal equilibrium, the anti-Hermitian part of the response operator is directly related to the correlation function of the fluctuations. Thus dissipative properties of the system are completely determined by the spectrum of density fluctuations---the fluctuation-dissipation theorem. In particular, we express the coefficient of dynamical friction for an orbiting test particle in terms of the fluctuation spectrum; this reduces to the known Chandrasekhar formula in the restrictive case of an infinite homogeneous system with a Maxwellian velocity distribution.

Over five years of daily hard X-ray (>20 keV) monitoring of the 2-min accretion-powered pulsar GX 1+4 with the Compton Gamma Ray Observatory/BATSE large-area detectors has found nearly continuous rapid spin-down, interrupted by a bright 200-d spin-up episode. During spin-down, the torque becomes more negative as the luminosity increases (assuming that the 20-60 keV pulsed flux traces bolometric luminosity), the opposite of what is predicted by standard accretion torque theory. No changes in the shape of the 20-100 keV pulsed energy spectrum were detected, so that a very drastic change in the spectrum below 20 keV or the pulsed fraction would be required to make the 20-60 keV pulsed flux a poor luminosity tracer. These are the first observations which flatly contradict standard magnetic disk accretion theory, and they may have important implications for understanding the spin evolution of X-ray binaries, cataclysmic variables, and protostars. We briefly discuss the possibility that GX 1+4 may be accreting from a retrograde disk during spin-down, as previously suggested.

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.

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.

Large-scale warps in the outer parts of spiral galaxy discs have been observed for almost forty years, but their origin remains obscure. We review the dynamics of warped galaxy discs. We identify several mechanisms that could excite warps, all involving the gravitational interaction between the disc and the dark-matter halo.

Estimates indicate there may be $\simgreat 10^3$ low luminosity X-ray pulsars ($L \simless 10^{34} \rom{erg\;s^{-1}}$) in the Galaxy undergoing ``low-state'' wind accretion in Be/X-ray binary systems, and $\sim 10^8-10^9$ isolated neutron stars which may be accreting directly from the interstellar medium. Despite their low effective temperatures ($kT_e \simless 300 \rom{eV}$), low luminosity accreting neutron stars with magnetic fields $B \sim (0.7-7) \times 10^{12}$ G could emit a substantial fraction ($0.5 - 5\%$) of their total luminosity in a moderately broadened ($E/\Delta E \sim 2-4$) cyclotron {\it emission} line which peaks in the energy range $ \sim 5-20\rom{keV}$. The bulk of the thermal emission from these stars will be in the extreme ultraviolet/soft X-ray regime, which is subject to strong interstellar absorption and would be difficult to distinguish from spectra of other types of objects. In sharp contrast, the {\it nonthermal} cyclotron component predicted here will not be strongly absorbed, and consequently it may be the only distinguishing signature for the bulk of these low luminosity sources. We propose a search for this cyclotron emission feature in long pointed observations of the newly discovered candidate isolated neutron star MS0317.7-6477, and the Be/X-ray transient

We investigate the dynamical interaction of galactic warps with the surrounding dark matter halo, using analytic perturbation theory. A precessing warp induces a density wake in the collisionless dark matter, which acts back on the original warp, transferring energy and angular momentum between the warp and halo (dynamical friction). In most cases dynamical friction damps the warp, but in unusual circumstances (such as a halo that rotates in the same direction as the precession of the warp, or a warp in the equatorial plane of an axisymmetric prolate halo) friction can excite a warp. The damping/excitation time is usually short compared to the Hubble time for realistic systems. Thus most warps cannot be primordial; they must be maintained by some ongoing excitation mechanism.

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