The Formation and Role of Vortices in Protoplanetary Disks

We carry out a two-dimensional, compressible, simulation of a disk, including dust particles, to study the formation and role of vortices in protoplanetary disks. We find that anticyclonic vortices can form out of an initial random perturbation of the vorticity field. Vortices have a typical decay time of the order of 50 orbital periods (for a viscosity parameter alpha=0.0001 and a disk aspect ratio of H/r = 0.15). If vorticity is continuously generated at a constant rate in the flow (e.g. by convection), then a large vortex can form and be sustained (due to the merger of vortices). We find that dust concentrates in the cores of vortices within a few orbital periods, when the drag parameter is of the order of the orbital frequency. Also, the radial drift of the dust induces a significant increase in the surface density of dust particles in the inner region of the disk. Thus, vortices may represent the preferred location for planetesimal formation in protoplanetary disks. We show that it is very difficult for vortex mergers to sustain a relatively coherent outward flux of angular momentum.Comment: Sumitted to the Astrophysical Journal, October 20, 199

V3885 Sagittarius: a Comparison with a Range of Standard Model Accretion Disks

A $\widetilde{\chi}^2$ analysis of standard model accretion disk synthetic spectrum fits to combined $FUSE$ and STIS spectra of V3885 Sagittarius, on an absolute flux basis, selects a model that accurately represents the observed SED. Calculation of the synthetic spectrum requires the following system parameters. The cataclysmic variable secondary star period-mass relation calibrated by Knigge in 2007 sets the secondary component mass. A mean white dwarf (WD) mass from the same study, that is consistent with an observationally-determined mass ratio, sets the adopted WD mass of $0.7M_{\odot}$, and the WD radius follows from standard theoretical models. The adopted inclination, i=65{\arcdeg}, is a literature consensus, and is subsequently supported by $\widetilde{\chi}^2$ analysis. The mass transfer rate is the remaining parameter to set the accretion disk $T_{\rm eff}$ profile, and the $Hipparcos$ parallax constrains that parameter to $\dot{M}=5.0{\pm}2.0{\times}10^{-9} M_{\odot} {\rm yr}^{-1}$ by a comparison with observed spectra. The fit to the observed spectra adopts the contribution of a $57,000{\pm}5000$K WD. The model thus provides realistic constraints on $\dot{M}$ and $T_{\rm eff}$ for a large $\dot{M}$ system above the period gap.Comment: 41 pages, 7 figures, 9 tables. Astrophysical Journal (accepted

The anomalous accretion disk of the Cataclysmic Variable RW Sextantis

Synthetic spectra covering the wavelength range 900\AA~to 3000\AA~provide an accurate fit, established by a ${\chi}_{\nu}^2$ analysis, to a combined observed spectrum of RW Sextantis. Two separately calibrated distances to the system establish the synthetic spectrum comparison on an absolute flux basis but with two alternative scaling factors, requiring alternative values of $\dot{M}$ for final models. Based on comparisons for a range of $\dot{M}$ values, the observed spectrum does not follow the standard model. Rather than the exponent 0.25 in the expression for the radial temperature profile, a value close to 0.125 produces a synthetic spectrum with an accurate fit to the combined spectrum. A study of time-series $FUSE$ spectra shows that a proposed warped or tilted disk is not supported by the data; an alternative proposal is that an observed non-axisymmetric wind results from an interaction with the mass transfer stream debris.Comment: 56 pages, 15 figures, 11 tables. Accepted for The Astrophysical Journa

Spiral Structure in IP Peg: Confronting Theory and Observations

The first convincing piece of evidence of spiral structure in the accretion disc in IP Pegasi was found by Steeghs et al. (1997). We performed two kinds of 2D hydrodynamic simulations, a SFS finite volume scheme and a SPH scheme, with a mass ratio of 0.5. Both results agreed well with each other. We constructed Doppler maps and line flux-binary phase relations based on density distributions, the results agreeing well with those obtained by observation.Comment: 4 pages, LaTeX with 2 ps figures using crckapb.sty. To appear in the Proceedings of Numerical Astrophysics 1998, Tokyo, Japan, 10-13 March, 1998, eds. S. M. Miyama, K. Tomisaka and T. Hanawa (Kluwer Academic Publishers