-
Bayesian sampling with BeAtlas, a grid of synthetic Be star spectra I. Recovering the fundamental parameters of αEri and βCMi
Authors:
A. C. Rubio,
A. C. Carciofi,
P. Ticiani,
B. C. Mota,
R. G. Vieira,
D. M. Faes,
M. Genaro,
T. H. de Amorim,
R. Klement,
I. Araya,
C. Arcos,
M. Curé,
A. Domiciano de Souza,
C. Georgy,
C. E. Jones,
M. W. Suffak,
A. C. F. Silva
Abstract:
Classical Be stars are fast rotating, near main sequence B-type stars. The rotation and the presence of circumstellar discs profoundly modify the observables of active Be stars. Our goal is to infer stellar and disc parameters, as well as distance and interstellar extinction, using the currently most favoured physical models for these objects. We present BeAtlas, a grid of 61.600 NLTE radiative tr…
▽ More
Classical Be stars are fast rotating, near main sequence B-type stars. The rotation and the presence of circumstellar discs profoundly modify the observables of active Be stars. Our goal is to infer stellar and disc parameters, as well as distance and interstellar extinction, using the currently most favoured physical models for these objects. We present BeAtlas, a grid of 61.600 NLTE radiative transfer models for Be stars, calculated with the HDUST code. The grid was coupled with a Monte Carlo Markov chain code to sample the posterior distribution. We test our method on two well-studied Be stars, $α$ Eri and $β$ CMi, using photometric, polarimetric and spectroscopic data as input to the code. We recover literature determinations for most of the parameters of the targets, in particular the mass and age of $α$ Eri, the disc parameters of $β$ CMi, and their distances and inclinations. The main discrepancy is that we estimate lower rotational rates than previous works. We confirm previously detected signs of disc truncation in $β$ CMi and note that its inner disc seems to have a flatter density slope than its outer disc. The correlations between the parameters are complex, further indicating that exploring the entire parameter space simultaneously is a more robust approach, statistically. The combination of BeAtlas and Bayesian-MCMC techniques proves successful, and a powerful new tool for the field: the fundamental parameters of any Be star can now be estimated in a matter of hours or days.
△ Less
Submitted 12 September, 2023;
originally announced September 2023.
-
The Be Star 66 Ophiuchi: 60 Years of Disk Evolution
Authors:
Keegan C. Marr,
Carol E. Jones,
Alex C. Carciofi,
Amanda C. Rubio,
Bruno C. Mota,
Mohammad R. Ghoreyshi,
Daniel W. Hatfield,
Leandro R. Rímulo
Abstract:
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 mode…
▽ More
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 66 Oph is consistent with standard B2Ve stellar properties. We computed a grid of 61568 Be star disk models to constrain the density profile of the disk before dissipation using observations of the H$α$ line profile and SED. We find at the onset of dissipation the disk has a base density of $2.5\times10^{-11}\ \rm{g\ cm^{-3}}$ with a radial power-law index of $n=2.6$. Our models indicate that after 21 years 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 $α= 0.4$ and an outer disk radius of $\sim$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 \pm 0.02\%$ with a polarization position angle of $θ_{IS}\approx85.7 \pm 0.7^\circ$. Using the Stokes QU diagram, we find the intrinsic polarization position angle of 66 Oph's disk is $θ_{int}\approx98 \pm 3^\circ$.
△ Less
Submitted 22 April, 2021; v1 submitted 11 March, 2021;
originally announced March 2021.
-
Visible and near-infrared spectro-interferometric analysis of the edge-on Be star $ο$ Aquarii
Authors:
E. S. G. de Almeida,
A. Meilland,
A. Domiciano de Souza,
P. Stee,
D. Mourard,
N. Nardetto,
R. Ligi,
I. Tallon-Bosc,
D. M. Faes,
A. C. Carciofi,
D. Bednarski,
B. C. Mota,
N. Turner,
T. A. ten Brummelaar
Abstract:
We present a detailed visible and near-IR spectro-interferometric analysis of the Be-shell star $ο$ Aquarii from quasi-contemporaneous CHARA/VEGA and VLTI/AMBER observations. We measured the stellar radius of $ο$ Aquarii as 4.0 $\pm$ 0.3 $\mathrm{R_{\odot}}$. We constrained the disk geometry and kinematics using a kinematic model and a MCMC fitting procedure. The disk sizes in H$α$ and Br$γ$ were…
▽ More
We present a detailed visible and near-IR spectro-interferometric analysis of the Be-shell star $ο$ Aquarii from quasi-contemporaneous CHARA/VEGA and VLTI/AMBER observations. We measured the stellar radius of $ο$ Aquarii as 4.0 $\pm$ 0.3 $\mathrm{R_{\odot}}$. We constrained the disk geometry and kinematics using a kinematic model and a MCMC fitting procedure. The disk sizes in H$α$ and Br$γ$ were found to be similar, at $\sim$10-12 $\mathrm{D_{\star}}$, which is uncommon since most results for Be stars show a larger extension in H$α$ than in Br$γ$. We found that the inclination angle $i$ derived from H$α$ is significantly lower ($\sim$15 deg) than the one derived from Br$γ$. The disk kinematics were found to be near to the Keplerian rotation in Br$γ$, but not in H$α$. After analyzing all our data using a grid of HDUST models (BeAtlas), we found a common physical description for the disk in both lines: $Σ_{0}$ = 0.12 g cm\textsuperscript{-2} and $m$ = 3.0. The stellar rotational rate was found to be very close ($\sim$96\%) to the critical value. Our analysis of multi-epoch H$α$ profiles and imaging polarimetry indicates that the disk has been stable for at least 20 years. Compared to Br$γ$, the data in H$α$ shows a substantially different picture that cannot fully be understood using the current physical models of Be star disks. $ο$ Aquarii presents a stable disk, but the measured $m$ is lower than the standard value in the VDD model for steady-state. Such long-term stability can be understood in terms of the high rotational rate for this star, the rate being a main source for the mass injection in the disk. Our results on the stellar rotation and disk stability are consistent with results in the literature showing that late-type Be stars are more likely to be fast rotators and have stable disks.
△ Less
Submitted 21 February, 2020;
originally announced February 2020.
-
The Large Scale Behaviour in the Disk of $δ$ Scorpii from 2000-2018
Authors:
M. W. Suffak,
C. E. Jones,
C. Tycner,
G. W. Henry,
A. C. Carciofi,
B. C. Mota,
A. C. Rubio
Abstract:
We model the circumstellar disk of $δ$ Sco using the 3-dimensional Monte Carlo radiative transfer code HDUST in order to quantify the large scale changes in the disk through the years 2000 to 2018, and to see if these changes can be attributed to the secondary star affecting the disk throughout its orbit. We determine our best-fitting models through matching simulated observations to actual H…
▽ More
We model the circumstellar disk of $δ$ Sco using the 3-dimensional Monte Carlo radiative transfer code HDUST in order to quantify the large scale changes in the disk through the years 2000 to 2018, and to see if these changes can be attributed to the secondary star affecting the disk throughout its orbit. We determine our best-fitting models through matching simulated observations to actual H$\rm α$ spectroscopy and V-band photometric observations. Our modelling results confirm previous findings that the disk of $δ$ Sco was forming early in the century. We also find a period of disk dissipation when the companion is at apastron, as well as a significant growth of the disk between 2009 and 2011, prior to the periastron of 2011. Due to the steady-state nature of the disk after 2011, it is difficult to say whether the variations seen are due to the effect of the close passage of the binary companion.
△ Less
Submitted 14 January, 2020;
originally announced January 2020.
-
Revealing the structure of the outer disks of Be stars
Authors:
R. Klement,
A. C. Carciofi,
T. Rivinius,
L. D. Matthews,
R. G. Vieira,
R. Ignace,
J. E. Bjorkman,
B. C. Mota,
D. M. Faes,
A. D. Bratcher,
M. Curé,
S. Štefl
Abstract:
Context. The structure of the inner parts of Be star disks (20 stellar radii) is well explained by the viscous decretion disk (VDD) model, which is able to reproduce the observable properties of most of the objects studied so far. The outer parts, on the ther hand, are not observationally well-explored, as they are observable only at radio wavelengths. A steepening of the spectral slope somewhere…
▽ More
Context. The structure of the inner parts of Be star disks (20 stellar radii) is well explained by the viscous decretion disk (VDD) model, which is able to reproduce the observable properties of most of the objects studied so far. The outer parts, on the ther hand, are not observationally well-explored, as they are observable only at radio wavelengths. A steepening of the spectral slope somewhere between infrared and radio wavelengths was reported for several Be stars that were previously detected in the radio, but a convincing physical explanation for this trend has not yet been provided. Aims. We test the VDD model predictions for the extended parts of a sample of six Be disks that have been observed in the radio to address the question of whether the observed turndown in the spectral energy distribution (SED) can be explained in the framework of the VDD model, including recent theoretical development for truncated Be disks in binary systems. Methods. We combine new multi-wavelength radio observations from the Karl. G. Jansky Very Large Array (JVLA) and Atacama Pathfinder Experiment (APEX) with previously published radio data and archival SED measurements at ultraviolet, visual, and infrared wavelengths. The density structure of the disks, including their outer parts, is constrained by radiative transfer modeling of the observed spectrum using VDD model predictions. In the VDD model we include the presumed effects of possible tidal influence from faint binary companions. Results. For 5 out of 6 studied stars, the observed SED shows strong signs of SED turndown between far-IR and radio wavelengths. A VDD model that extends to large distances closely reproduces the observed SEDs up to far IR wavelengths, but fails to reproduce the radio SED. ... (abstract continues but did not fit here)
△ Less
Submitted 21 March, 2017;
originally announced March 2017.
-
Disk-Loss and Disk-Renewal Phases in Classical Be Stars II. Contrasting with Stable and Variable Disks
Authors:
Zachary H. Draper,
John P. Wisniewski,
Karen S. Bjorkman,
Marilyn R. Meade,
Xavier Haubois,
Bruno C. Mota,
Alex C. Carciofi,
Jon E. Bjorkman
Abstract:
Recent observational and theoretical studies of classical Be stars have established the utility of polarization color diagrams (PCD) in helping to constrain the time-dependent mass decretion rates of these systems. We expand on our pilot observational study of this phenomenon, and report the detailed analysis of a long-term (1989-2004) spectropolarimetric survey of 9 additional classical Be stars,…
▽ More
Recent observational and theoretical studies of classical Be stars have established the utility of polarization color diagrams (PCD) in helping to constrain the time-dependent mass decretion rates of these systems. We expand on our pilot observational study of this phenomenon, and report the detailed analysis of a long-term (1989-2004) spectropolarimetric survey of 9 additional classical Be stars, including systems exhibiting evidence of partial disk-loss/disk-growth episodes as well as systems exhibiting long-term stable disks. After carefully characterizing and removing the interstellar polarization along the line of sight to each of these targets, we analyze their intrinsic polarization behavior. We find that many steady-state Be disks pause at the top of the PCD, as predicted by theory. We also observe sharp declines in the Balmer jump polarization for later spectral type, near edge-on steady-state disks, again as recently predicted by theory, likely caused when the base density of the disk is very high, and the outer region of the edge-on disk starts to self absorb a significant number of Balmer jump photons. The intrinsic $V$-band polarization and polarization position angle of $γ$ Cas exhibits variations that seem to phase with the orbital period of a known one-armed density structure in this disk, similar to the theoretical predictions of Halonen & Jones. We also observe stochastic jumps in the intrinsic polarization across the Balmer jump of several known Be+sdO systems, and speculate that the thermal inflation of part of the outer region of these disks could be responsible for producing this observational phenomenon. Finally, we estimate the base densities of this sample of stars to be between $\approx 8\times 10^{-11}$ to $\approx 4 \times 10^{-12}\,\rm g cm^{-3}$ during quasi steady state periods given there maximum observed polarization.
△ Less
Submitted 21 February, 2014;
originally announced February 2014.
-
Dynamical Evolution of Viscous Disks around Be Stars. II : polarimetry
Authors:
Xavier Haubois,
Bruno C. Mota,
Alex C. Carciofi,
Zachary H. Draper,
John P. Wisniewski,
Daniel Bednarski,
Thomas Rivinius
Abstract:
Be stars exhibit variability for a great number of observables. Putting the pieces together of the disk dynamics is not an easy task and requires arduous modeling before achieving a good fit of the observational data. In order to guide the modeling process and make it more efficient, it is very instructive to investigate reference dynamical cases. This paper focuses on continuum polarimetric quant…
▽ More
Be stars exhibit variability for a great number of observables. Putting the pieces together of the disk dynamics is not an easy task and requires arduous modeling before achieving a good fit of the observational data. In order to guide the modeling process and make it more efficient, it is very instructive to investigate reference dynamical cases. This paper focuses on continuum polarimetric quantities and is the second of a series that aims to demonstrate the capacity of deriving the dynamical history and fundamental parameters of a classical Be star through the follow-up of various observables. After a detailed study of the different opacities at play in the formation of polarized spectra, we investigate predictions of polarimetric observables in the continuum for different dynamical scenarios. Our models are based on a coupling of a hydrodynamic viscous decretion simulations in a disk and a 3-D non-LTE radiative transfer code. Introducing the polarization color diagram (PCD), we show that certain combinations of polarimetric observables exhibit features that are characteristic of a mass loss history. This diagram also enables to estimate fundamental parameters such as the inclination angle, the disk density scale and the $α$ viscous diffusion parameter. We present the PCD as a powerful diagnosis tool to track the dynamical phases of a Be star such as disk building-up, dissipation, periodic and episodic outbursts. Finally we confront our models with observations of 4 Be stars that exhibited long-term polarimetric activity.
△ Less
Submitted 27 February, 2014; v1 submitted 9 February, 2014;
originally announced February 2014.
