We use data from the Hipparcos Catalog and the Barbier-Brossat & Figon (2000) catalog of stellar radial velocities to test the hypothesis that the Beta Pic planetesimal disk was disrupted by a close stellar encounter. We trace the space motions of 21,497 stars and discover 18 that have passed within 5 pc of Beta Pic in the past 1 Myr. Beta Pic's closest encounter is with the K2III star HIP 27628 (0.6 pc), but dynamically the most important encounter is with the F7V star HIP 23693 (0.9 pc). We calculate the velocity and eccentricity changes induced by the 18 perturbations and conclude that they are dynamically significant if planetesimals exist in a Beta Pic Oort cloud. We provide a first-order estimate for the evolutionary state of a Beta Pic Oort cloud and conclude that the primary role of these stellar perturbations would be to help build a comet cloud rather than destroy a pre-existing structure. The stellar sample is 20% complete and motivates future work to identify less common close interactions that would significantly modify the observed circumstellar disk. For future radial velocity study we identify six stars in the Hipparcos Catalog that may have approached Beta Pic to within 0.1 pc and therefore remain as candidate disk perturbers.

We numerically investigate the dynamics of how a close stellar fly-by encounter of a symmetrical circumstellar planetesimal disc can give rise to the many kinds of asymmetries and substructures attributed to the edge-on dusty disc of Beta Pic. In addition we present new optical coronagraphic observations of the outer parts of Beta Pic's disc, and report that the radial extent is significantly greater than was found in previous measurements. The northeasterly extension of the disc's midplane is now measured out to 1835au from the star; the southwesterly component is measured out to 1450au. Hence we use the length asymmetry induced in a distribution of simulation test particles as the principal diagnostic feature when modelling the disc response, in order to constrain fly-by parameters. In particular we favour a low inclination prograde and near-parabolic orbit perturber of mass approximately 0.5 Solar masses. These initial conditions suggest that the perturber could have been physically associated with Beta Pic prior to the encounter. Thus we also consider the possibility that the perturber could be bound to Beta Pic: a consideration also of general interest where dust discs are known to exist in binary star systems. In some of our models, we can relate groupings of perturbed particles to the large-scale structure of the Beta Pic disc. The groupings correspond to: high eccentricity and inclination particles that reach apocentre and maximum height in the southwest, moderately eccentric and low inclination particles that reach apocentre in the northeast, and relatively unperturbed particles inside approximately 200au radius.

We have numerically investigated the dynamics of how a close stellar fly-by encounter of a symmetrical circumstellar planetesimal disc can give rise to the many kinds of asymmetries and substructures attributed to the dusty disc of Beta Pic. We find three recognizable groupings of test particles that can be related to the morphology of the Beta Pic disc. These are: highly eccentric and inclined orbit particles that reach apocentre in the southwest, moderately eccentric and inclined orbit particles that reach apocentre in the northeast, and a relatively unperturbed region inside ~200au radius.

The nearby main sequence star Beta Pictoris is surrounded by an edge-on disk of dust produced by the collisional erosion of larger planetesimals. Here we report the discovery of substructure within the northeast extension of the disk midplane that may represent an asymmetric ring system around Beta Pic. We present a dynamical model showing that a close stellar flyby with a quiescient disk of planetesimals can create such rings, along with previously unexplained disk asymmetries. Thus we infer that Beta Pic's planetesimal disk was highly disrupted by a stellar encounter in the last hundred thousand years.

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

We have investigated effects of early stellar encounters on a protoplanetary disk (planetesimal disk) and found that they can explain the high eccentricities and inclinations observed in the outer part ($>42$AU) of the Edgeworth-Kuiper Belt (EKB). The proto-sun is considered as a member of a stellar aggregation that undergoes dissolution on a timescale $\sim 10^8$yrs, such that a planetesimal disk experiences a flyby encounter at pericenter distance ($q$) on the order of 100AU. We simulated dynamical evolution of a planetesimal (test particle) disk perturbed by a passing star. We show that the stellar encounter pumps the velocity dispersion in the disk in the outer parts ($> 0.25q$). Planet formation is forestalled in that region. We also find that a stellar encounter with pericenter distance $q \sim 100-200$AU could have pumped up the velocity dispersion of EKB objects outside 42AU to the observed magnitude while preserving that inside Neptune's 3:2 mean-motion resonance (located at 39.5AU), which allows for the efficient capture of objects by the sweeping of the 3:2 resonance during orbital migration by proto-Neptune.

In this letter recent results on the nodal precession of accretion discs in close binaries are applied to the discs in some X-ray binary systems. The ratio between the tidally forced precession period and the binary orbital period is given, as well as the condition required for the rigid precession of gaseous Keplerian discs. Hence the minimum precessional period that may be supported by a fluid Keplerian disc is determined. It is concluded that near rigid body precession of tilted accretion discs can occur and generally reproduce observationally inferred precession periods, for reasonable system parameters. In particular long periods in SS433, Her X-1, LMC X-4 and SMC X-1 can be fit by the tidal model. It is also found that the precession period that has been tentatively put forward for Cyg X-2 cannot be accomodated by a tidally precessing disc model for any realistic choice of system parameters.

In this paper we revisit the problem of the tidal interaction occuring between a protostellar accretion disc and a secondary point mass following a parabolic trajectory. We model the disc response analytically and we compare our results with three-dimensional SPH simulations. Inviscid as well as viscous hydrodynamics is considered. We show that in a viscous system the response derived from inviscid considerations is predominant even for the highest estimates of an anomalous disc shear viscosity. The angular momentum lost from the disc during the encounter is derived from linear theory, for distant fly-bys, as well as the changes to the disc orientation expected in non-coplanar encounters. It is shown that the target discs can become warped and precess by a small amount during non-coplanar encounters. This small precession is shown to give rise to a relative tilt of the disc which is always more important for determining its final orientation than is the change to the orbital inclination. We discuss the implications of our results for protostellar accretion discs and planetary systems.

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

We describe some recent nonlinear three dimensional hydrodynamic simulations of accretion discs in binary systems where the orbit is circular and not necessarily coplanar with the disc midplane. The calculations are relevant to a number of observed astrophysical phenomena, including the precession of jets associated with young stars, the high spectral index of some T Tauri stars, and the light curves of X-ray binaries such as Hercules X-1 which suggest the presence of precessing accretion discs.

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

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