-
Evolution of the disk in the Be binary $δ$ Scorpii probed during three periastron passages
Authors:
R. G. Rast,
C. E. Jones,
A. C. Carciofi,
M. W. Suffak,
A. C. F. Silva,
G. W. Henry,
C. Tycner
Abstract:
We examine the evolution of the disk surrounding the Be star in the highly eccentric binary system $δ$ Scorpii over its three most recent periastron passages. $V$-band and $B-V$ photometry, along with H$α$ spectroscopy are combined with a new set of extensive multi-band polarimetry data to produce a detailed comparison of the disk's physical conditions during the time periods surrounding each clos…
▽ More
We examine the evolution of the disk surrounding the Be star in the highly eccentric binary system $δ$ Scorpii over its three most recent periastron passages. $V$-band and $B-V$ photometry, along with H$α$ spectroscopy are combined with a new set of extensive multi-band polarimetry data to produce a detailed comparison of the disk's physical conditions during the time periods surrounding each closest approach of the secondary star. We use the three-dimensional Monte Carlo radiative transfer code \textsc{HDUST} and smoothed particle hydrodynamics (\textsc{SPH}) code to support our observations with models of disk evolution, discussing the behaviour of the H$α$ and He\,\textsc{i} 6678 lines, $V$-band magnitude, and polarization degree. We compare the characteristics of the disk immediately before each periastron passage to create a baseline for the unperturbed disk. We find that the extent of the H$α$ emitting region increased between each periastron passage, and that transient asymmetries in the disk become more pronounced with each successive encounter. Asymmetries of the H$α$ and He\,\textsc{i} 6678 lines in 2011 indicate that perturbations propagate inward through the disk near periastron. When the disk's direction of orbit is opposite to that of the secondary, the parameters used in our models do not produce spiral density enhancements in the H$α$ emitting region because the tidal interaction time is short due to the relative velocities of the disk particles with the secondary. The effects of the secondary star on the disk are short-lived and the disk shows independent evolution between each periastron event.
△ Less
Submitted 22 April, 2024;
originally announced April 2024.
-
Disc Tearing in a Be Star: Predicted 3D Observations
Authors:
M. W. Suffak,
C. E. Jones,
A. C. Carciofi
Abstract:
We build on our previous work involving smoothed particle hydrodynamic simulations of Be stars, by using the model that exhibited disc tearing as input into the three-dimensional Monte Carlo radiative transfer code HDUST to predict observables from a variety of viewing angles throughout the disc tearing process. We run one simulation at the start of each orbital period from 20 to 72 orbital period…
▽ More
We build on our previous work involving smoothed particle hydrodynamic simulations of Be stars, by using the model that exhibited disc tearing as input into the three-dimensional Monte Carlo radiative transfer code HDUST to predict observables from a variety of viewing angles throughout the disc tearing process. We run one simulation at the start of each orbital period from 20 to 72 orbital periods, which covers two complete disc tearing events. The resulting trends in observables are found to be dependent on the relative position of the observer and the tearing disc. The $\rm Hα$ equivalent width, $V$ magnitude, and polarization can all increase or decrease in any combination depending on the viewpoint of the observer. The $\rm Hα$ line profile also displays changes in strength and peak separation throughout the tearing process. We show how the outer disc of the torn system can have a large effect on the $\rm Hα$ line profile, and also contributes to a wavelength-dependent polarization position angle, resulting in a similar sawtooth shape to the polarization percentage. Finally, we compare our predictions to Pleione (28 Tau) where evidence has suggested that a disc tearing event has occurred in the past. We find that our tearing disc model can broadly match the trends seen in Pleione's observables, as well as produce the two-component $\rm Hα$ lines observed in Pleione. This is the strongest evidence, thus far, of Pleione's disc having indeed experienced a tearing event.
△ Less
Submitted 23 November, 2023;
originally announced November 2023.
-
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.
-
Non-LTE Monte Carlo Radiative Transfer. III. The thermal properties of Tilted and Warped Be Star Discs
Authors:
M. W. Suffak,
C. E. Jones,
A. C. Carciofi,
T. H. de Amorim
Abstract:
We use the three-dimensional Monte Carlo radiative transfer code HDUST to model Be stars where the disc is tilted from the equatorial plane of the star. We compute 128 models across 4 spectral types, B0, B2, B5 and B8, tilting the disc by $0^o$, $10^o$, $20^o$, and $40^o$, while varying disc density according to spectral type. We also compute every model for an average and high stellar rotation ra…
▽ More
We use the three-dimensional Monte Carlo radiative transfer code HDUST to model Be stars where the disc is tilted from the equatorial plane of the star. We compute 128 models across 4 spectral types, B0, B2, B5 and B8, tilting the disc by $0^o$, $10^o$, $20^o$, and $40^o$, while varying disc density according to spectral type. We also compute every model for an average and high stellar rotation rate. We first discuss non-tilted disc temperatures and show its non-linear dependence on stellar and disc parameters. We find that tilting the disc minimally affects the density-weighted average disc temperature, but tilting does create a temperature asymmetry in disc cross sections, which is more pronounced for a faster rotation rate. We also investigate the effect tilting has on $V$-band magnitude, polarization, and the H$α$ line. Tilting the disc does affect these observables, but the changes are entirely dependent on the position of the observer relative to the direction of tilt. We find the observables that distinguish tilting from a change in density or geometry are the H$α$ line shape, where it can transition between single-peaked and double-peaked, and the polarization position angle, whose value is dependent on the projected major elongation axis of the disc on the sky. We also present one early and one late-type model with warped discs. We find their temperature structure varies a small amount from the uniformly tilted models, and the different observables correspond to different tilt angles, consistent with their expected volume of origin within the disc.
△ Less
Submitted 9 September, 2023;
originally announced September 2023.
-
Strong Variability in AzV 493, an Extreme Oe-Type Star in the SMC
Authors:
M. S. Oey,
N. Castro,
M. Renzo,
I. Vargas-Salazar,
M. W. Suffak,
M. Ratajczak,
J. D. Monnier,
M. K. Szymanski,
G. D. Phillips,
N. Calvet,
A. Chiti,
G. Micheva,
K. C. Rasmussen,
R. H. D. Townsend
Abstract:
We present 18 years of OGLE photometry together with spectra obtained over 12 years, revealing that the early Oe star AzV 493 shows strong photometric (Delta I < 1.2 mag) and spectroscopic variability with a dominant, 14.6-year pattern and ~40-day oscillations. We estimate stellar parameters T_eff = 42000 K, log L/L_sun = 5.83 +/- 0.15, M/M_sun = 50 +/- 9, and vsini = 370 +/- 40 km/s. Direct spect…
▽ More
We present 18 years of OGLE photometry together with spectra obtained over 12 years, revealing that the early Oe star AzV 493 shows strong photometric (Delta I < 1.2 mag) and spectroscopic variability with a dominant, 14.6-year pattern and ~40-day oscillations. We estimate stellar parameters T_eff = 42000 K, log L/L_sun = 5.83 +/- 0.15, M/M_sun = 50 +/- 9, and vsini = 370 +/- 40 km/s. Direct spectroscopic evidence shows episodes of both gas ejection and infall. There is no X-ray detection, and it is likely a runaway star. AzV 493 may have an unseen companion on a highly eccentric (e > 0.93) orbit. We propose that close interaction at periastron excites ejection of the decretion disk, whose variable emission-line spectrum suggests separate inner and outer components, with an optically thick outer component obscuring both the stellar photosphere and the emission-line spectrum of the inner disk at early phases in the photometric cycle. It is plausible that AzV 493's mass and rotation have been enhanced by binary interaction followed by the core-collapse supernova explosion of the companion, which now could be either a black hole or neutron star. This system in the Small Magellanic Cloud can potentially shed light on OBe decretion disk formation and evolution, massive binary evolution, and compact binary progenitors.
△ Less
Submitted 26 January, 2023;
originally announced January 2023.
-
Ultraviolet Spectropolarimetry: on the origin of rapidly rotating B stars
Authors:
C. E. Jones,
J. Labadie-Bartz,
D. V. Cotton,
Y. Nazé,
G. J. Peters,
D. J. Hillier,
C. Neiner,
N. D. Richardson,
J. L. Hoffman,
A. C. Carciofi,
J. P. Wisniewski,
K. G. Gayley,
M. W. Suffak,
R. Ignace,
P. A. Scowen
Abstract:
UV spectroscopy and spectropolarimetry hold the key to understanding certain aspects of massive stars that are largely inaccessible with optical or longer wavelength observations. This is especially true for the rapidly-rotating Be and Bn stars, owing to their high temperatures, geometric asymmetries, binary properties, and evolutionary history. UV spectropolarimetric observations are extremely se…
▽ More
UV spectroscopy and spectropolarimetry hold the key to understanding certain aspects of massive stars that are largely inaccessible with optical or longer wavelength observations. This is especially true for the rapidly-rotating Be and Bn stars, owing to their high temperatures, geometric asymmetries, binary properties, and evolutionary history. UV spectropolarimetric observations are extremely sensitive to the photospheric consequences of rapid rotation (i.e. oblateness, temperature, and surface gravity gradients). Our polarized radiative-transfer modelling predicts that with low-resolution UV spectropolarimetry covering 120 -- 300 nm the inclination angle of a rapid rotator can be determined to within 5 degrees, and the rotation rate to within 1%. The origin of rapid rotation in Be/n stars can be explained by either single-star or binary evolution, but their relative importance is largely unknown. Some Be stars have hot sub-luminous (sdO) companions, which at an earlier phase transferred their envelope (and with it mass and angular momentum) to the present-day rapid rotator. Through spectral modelling of a wide range of simulated Be/n+sdO configurations, we demonstrate that high-resolution high-SNR UV spectroscopy can detect an sdO star even when $\sim$1,000 times fainter in the UV than its Be/n star companion. This degree of sensitivity is needed to more fully explore the parameter space of Be/n+sdO binaries, which so far has been limited to about a dozen systems with relatively luminous sdO stars. We suggest that a UV spectropolarimetric survey of Be/n stars is the next step forward in understanding this population. Such a dataset would, when combined with population synthesis models, allow for the determination of the relative importance of the possible evolutionary pathways traversed by these stars, which is also crucial for understanding their future evolution and fate.
△ Less
Submitted 27 December, 2022; v1 submitted 15 November, 2021;
originally announced November 2021.
-
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.