We first present a summary of our numerical work on accretion discs in close binary systems. Our recent studies on numerical simulations of the surface flow on the mass-losing star in a close binary star is then reviewed.
Discs in binaries have a complex behavior because of the perturbations of the companion star. Planet formation in binary-star systems both depend on the companion star parameters and on the properties of the circumstellar disc. An eccentric disc may increase the impact velocity of planetesimals and therefore jeopardize the accumulation process. We model the evolution of discs in close binaries including the effects of self-gravity and adopting different prescriptions to model the discs radiative properties. We focus on the dynamical properties and evolutionary tracks of the discs. We use the hydrodynamical code FARGO and we include in the energy equation heating and cooling effects. Radiative discs have a lower disc eccentricity compared to locally isothermal discs with same temperature profile. As a consequence, we do not observe the formation of an internal elliptical low density region as in locally isothermal disc models. However, the disc eccentricity depends on the disc mass through the opacities. Akin to locally isothermal disc models, self-gravity forces the discs longitude of pericenter to librate about a fixed orientation with respect to the binary apsidal line ($pi$). The discs radiative properties play an important role in the evolution of discs in binaries. A radiative disc has an overall shape and internal structure that are significantly different compared to a locally isothermal disc with same temperature profile. This is an important finding both for describing the evolutionary track of the disc during its progressive mass loss, and for planet formation since the internal structure of the disc is relevant for planetesimals growth in binary systems. The non-symmetrical distribution of mass in these discs causes large eccentricities for planetesimals that may affect their growth.
We study the excitation and damping of tides in close binary systems, accounting for the leading order nonlinear corrections to linear tidal theory. These nonlinear corrections include two distinct effects: three-mode nonlinear interactions and nonlinear excitation of modes by the time-varying gravitational potential of the companion. This paper presents the formalism for studying nonlinear tides and studies the nonlinear stability of the linear tidal flow. Although the formalism is applicable to binaries containing stars, planets, or compact objects, we focus on solar type stars with stellar or planetary companions. Our primary results include: (1) The linear tidal solution often used in studies of binary evolution is unstable over much of the parameter space in which it is employed. More specifically, resonantly excited gravity waves are unstable to parametric resonance for companion masses M > 10-100 M_Earth at orbital periods P = 1-10 days. The nearly static equilibrium tide is, however, parametrically stable except for solar binaries with P < 2-5 days. (2) For companion masses larger than a few Jupiter masses, the dynamical tide causes waves to grow so rapidly that they must be treated as traveling waves rather than standing waves. (3) We find a novel form of parametric instability in which a single parent wave excites a very large number of daughter waves (N = 10^3[P / 10 days]) and drives them as a single coherent unit with growth rates that are ~N times faster than the standard three wave parametric instability. (4) Independent of the parametric instability, tides excite a wide range of stellar p-modes and g-modes by nonlinear inhomogeneous forcing; this coupling appears particularly efficient at draining energy out of the dynamical tide and may be more important than either wave breaking or parametric resonance at determining the nonlinear dissipation of the dynamical tide.
Short period binary systems containing magnetic Ap stars are anomalously rare. This apparent anomaly may provide insight into the origin of the magnetic fields in theses stars. As an early investigation of this, we observed three close binary systems that have been proposed to host Ap stars. Two of these systems (HD 22128 and HD 56495) we find contain Am stars, but not Ap stars. However, for one system (HD 98088) we find the primary is indeed an Ap star, while the secondary is an Am star. Additionally, the Ap star is tidally locked to the secondary, and the predominately dipolar magnetic field of the Ap star is roughly aligned with the secondary. Further investigations of HD 98088 are planned by the BinaMIcS collaboration.
We simulate numerically the surface flow of a gas-supplying companion star in a semi-detached binary system. Calculations are carried out for a region including only the mass-losing star, thus not the mass accreting star. The equation of state is that of an ideal gas characterized by a specific heat ratio gamma, and the case with gamma=5/3 is mainly studied. A system of eddies appears on the surface of the companion star: an eddy in the low pressure region near the L1 point, one around the high pressure at the north pole, and one or two eddies around the low pressure at the opposite side of the L1 point. Gas elements starting near the pole region rotate clockwise around the north pole (here the binary system rotates counter-clockwise as seen from the north pole). Because of viscosity, the gas drifts to the equatorial region, switches to the counter-clockwise eddy near the L1 point and flows through the L1 point to finally form the L1 stream. The flow field in the L1 region and the structure of the L1 stream are also considered.
We study the development of finite eccentricity in accretion disks in close binary systems using a two-dimensional grid-based numerical scheme. We perform detailed parameter studies to explore the dependence on viscosity, disk aspect ratio, the inclusion of a mass-transfer stream and the role of the boundary conditions. We consider mass ratios 0.05<q<0.3 appropriate to superoutbursting cataclysmic binary systems. Instability to the formation of a precessing eccentric disk that attains a quasi-steady state with mean eccentricity in the range 0.3-0.5 occurs readily. The shortest growth times are ~15 binary orbits for the largest viscosities and the instability mechanism is for the most part consistent with the mode-coupling mechanism associated with the 3:1 resonance proposed by Lubow. However, the results are sensitive to the treatment of the inner boundary and to the incorporation of the mass-transfer stream. In the presence of a stream we found a critical viscosity below which the disk remains circular. Incorporation of a mass-transfer stream tends to impart stability for small enough viscosity (or, equivalently, mass-transfer rate through the disk) and does assist in obtaining a prograde precession rate that is in agreement with observations. For the larger q the location of the 3:1 resonance is pushed outwards towards the Roche lobe where higher-order mode couplings and nonlinearity occur. It is likely that three-dimensional simulations that properly resolve the disks vertical structure are required to make significant progress in this case.
Log in to be able to interact and post comments
comments
Fetching comments
Sorry, something went wrong while fetching comments!