The Be Star 66 Ophiuchi: 60 Years of Disk Evolution

We use a time-dependent hydrodynamic code and a non-LTE Monte Carlo code to model disk dissipation for the Be star 66 Ophiuchi. We compiled 63 years of observations from 1957 to 2020 to encompass the complete history of the growth and subsequent dissipation of the star's disk. Our models are constrained by new and archival photometry, spectroscopy, and polarization observations, allowing us to model the disk dissipation event. Using Markov Chain Monte Carlo methods, we find that the properties of 66 Oph are consistent with those of a standard B2Ve star. We computed a grid of 61,568 Be star disk models to constrain the density profile of the disk before dissipation using observations of the H alpha line profile and spectral energy distribution. We find at the onset of dissipation the disk has a base density of 2.5 x 10(-11) g cm(-3) with a radial power-law index of n = 2.6. Our models indicate that after 21 yr of disk dissipation 66 Oph's outer disk remained present and bright in the radio. We find an isothermal disk with constant viscosity with an alpha = 0.4 and an outer disk radius of similar to 115 stellar radii best reproduces the rate of 66 Oph's disk dissipation. We determined the interstellar polarization in the direction of the star in the V band is p = 0.63 +/- 0.02% with a polarization position angle of theta (IS) approximate to 85.degrees 7 +/- 0.degrees 7. Using the Stokes QU diagram, we find the intrinsic polarization position angle of 66 Oph's disk is theta (int) approximate to 98 degrees +/- 3 degrees.