Showing 1–50 of 99 results for author: Carciofi, A

On the Origin of Fast Rotating Stars. I. Photometric calibration and results of AO-assisted BVRI+Halpha imaging of NGC330 with SAMI/SOAR

Authors: Felipe Navarete, Pedro Ticiani dos Santos, Alex Cavaliéri Carciofi, André Luiz Figueiredo

Abstract: H$α$ emission is a clear indicator of circumstellar activity in Be stars, historically employed to assess the classical Be star (CBe) population in young open clusters (YOCs). The YOC NGC330 in the Small Magellanic Cloud exhibits a large known fraction of {CBe} stars and was selected for a pilot study to establish a comprehensive methodology for identifying H$α$ emitters in the Magellanic Clouds,… ▽ More H$α$ emission is a clear indicator of circumstellar activity in Be stars, historically employed to assess the classical Be star (CBe) population in young open clusters (YOCs). The YOC NGC330 in the Small Magellanic Cloud exhibits a large known fraction of {CBe} stars and was selected for a pilot study to establish a comprehensive methodology for identifying H$α$ emitters in the Magellanic Clouds, encompassing the entire B-type spectral range. Using the SOAR Adaptative Module Imager (SAMI), we investigated the stellar population of NGC330 using multi-band BVRI+H$α$ imaging. We identified H$α$ emitters within the entire V-band range covered by SAMI/SOAR observations ($V\lesssim22$), comprising the complete B-type stellar population and offering a unique opportunity to explore the Be phenomenon across all spectral sub-classes. The stellar radial distribution shows a clear bimodal pattern between the most massive (B5 or earlier) and the lower-mass main-sequence objects (later than B6) within the cluster. The former is concentrated towards the cluster center (showing a dispersion of $σ=4.26\pm0.20$ pc), whereas the latter extends across larger radii ($σ=10.83\pm0.65$ pc), indicating mass stratification within NGC330. The total fraction of emitters is $4.4\pm0.5\%$, notably smaller than previous estimates from flux- or seeing-limited observations. However, a higher fraction of H$α$ emitters is observed among higher-mass stars ($32.8\pm3.4\%$) than within lower-mass ($4.4\pm0.9\%$). Consequently, the putative CBe population exhibits distinct dynamical characteristics compared to the bulk of the stellar population in NGC330. These findings highlight the significance of the current observations in providing a complete picture of the CBe population in NGC330. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Comments: Accepted for publication in the Astrophysical Journal. 24 pages, 17 figures, 3 tables

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.

Comments: 20 pages, 13 figures, accepted for publication in ApJ

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.

Comments: 8 pages, 12 figures

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.

Comments: 20 pages, 23 figures

Dynamical masses of the primary Be star and the secondary sdB star in the single-lined binary kappa Dra (B6 IIIe)

Authors: R. Klement, D. Baade, Th. Rivinius, D. R. Gies, L. Wang, J. Labadie-Bartz, P. Ticiani Dos Santos, J. D. Monnier, A. C. Carciofi, A. Mérand, N. Anugu, G. H. Schaefer, J. -B. Le Bouquin, C. L. Davies, J. Ennis, T. Gardner, S. Kraus, B. R. Setterholm, A. Labdon

Abstract: Because many classical Be stars may owe their nature to mass and angular-momentum transfer in a close binary, the present masses, temperatures, and radii of their components are of high interest for comparison to stellar evolution models. Kappa Dra is a 61.5-day single-lined binary with a B6 IIIe primary. With the CHARA Array instruments MIRC/MIRC-X and MYSTIC, we detected the secondary at (approx… ▽ More Because many classical Be stars may owe their nature to mass and angular-momentum transfer in a close binary, the present masses, temperatures, and radii of their components are of high interest for comparison to stellar evolution models. Kappa Dra is a 61.5-day single-lined binary with a B6 IIIe primary. With the CHARA Array instruments MIRC/MIRC-X and MYSTIC, we detected the secondary at (approximately photospheric) flux ratios of 1.49 +- 0.10% and 1.63 +- 0.09% in the H and K band, respectively. From a large and diverse optical spectroscopic database only the radial velocity curve of the Be star could be extracted. However, employing the parallaxes from Hipparcos and Gaia, which agree within their nominal 1-sigma errors, we could derive the total mass and found component masses of 3.65 +- 0.48 Msun and 0.426 +- 0.043 Msun for the Be star and the companion, respectively. Previous cross-correlation of the observed far-UV spectrum with sdO spectral model templates had not detected a companion belonging to the hot O-type subdwarf (sdO) population known from ~20 earlier-type Be stars. Guided by our full 3D orbital solution, we found a strong cross-correlation signal for a stripped subdwarf B-type companion (far-UV flux ratio of 2.3 +- 0.5%), enabling the first firm characterization of such a star, and making kappa Dra the first mid- to late-type Be star with a directly-observed subdwarf companion. △ Less

Submitted 6 October, 2022; originally announced October 2022.

Comments: Accepted to ApJ. Figure sets available on request

The binary system of the spinning-top Be star Achernar

Authors: P. Kervella, S. Borgniet, A. Domiciano de Souza, A. Mérand, A. Gallenne, Th. Rivinius, S. Lacour, A. Carciofi, D. Moser Faes, J. -B. Le Bouquin, M. Taormina, B. Pilecki, J. -Ph. Berger, Ph. Bendjoya, R. Klement, F. Millour, E. Janot-Pacheco, A. Spang, F. Vakili

Abstract: Achernar, the closest and brightest classical Be star, presents rotational flattening, gravity darkening, occasional emission lines due to a gaseous disk, and an extended polar wind. It is also a member of a close binary system with an early A-type dwarf companion. We aim to determine the orbital parameters of the Achernar system and to estimate the physical properties of the components. We monito… ▽ More Achernar, the closest and brightest classical Be star, presents rotational flattening, gravity darkening, occasional emission lines due to a gaseous disk, and an extended polar wind. It is also a member of a close binary system with an early A-type dwarf companion. We aim to determine the orbital parameters of the Achernar system and to estimate the physical properties of the components. We monitored the relative position of Achernar B using a broad range of high angular resolution instruments of the VLT/VLTI (VISIR, NACO, SPHERE, AMBER, PIONIER, GRAVITY, and MATISSE) over a period of 13 years (2006-2019). These astrometric observations are complemented with a series of more than 700 optical spectra for the period from 2003 to 2016. We determine that Achernar B orbits the Be star on a seven-year period, eccentric orbit (e = 0.7255 +/- 0.0014) which brings the two stars within 2 au at periastron. The mass of the Be star is found to be mA = 6.0 +/- 0.6 Msun for a secondary mass of mB = 2.0 +/- 0.1 Msun. We find a good agreement of the parameters of Achernar A with the evolutionary model of a critically rotating star of 6.4 Msun at an age of 63 million years. We also identify a resolved comoving low-mass star, which leads us to propose that Achernar is a member of the Tucana-Horologium moving group. Achernar A is presently in a short-lived phase of its evolution following the turn-off, during which its geometrical flattening ratio is the most extreme. Considering the orbital parameters, no significant interaction occurred between the two components, demonstrating that Be stars may form through a direct, single-star evolution path without mass transfer. Since component A will enter the instability strip in a few hundred thousand years, Achernar appears to be a promising progenitor of the Cepheid binary systems. △ Less

Submitted 15 September, 2022; originally announced September 2022.

Comments: 27 pages, Astronomy & Astrophysics, in press. The full set of continuum normalized high resolution spectra of Achernar is available at https://doi.org/10.5281/zenodo.6977303

Journal ref: A&A 667, A111 (2022)

The Role of Disk Tearing and Precession in the Observed Variability of Pleione

Authors: Keegan C. Marr, Carol E. Jones, Chris Tycner, Alex C. Carciofi, Ariane C. Fonseca Silva

Abstract: We acquired H$α$ spectroscopic observations from 2005 to 2019 showing Pleione has transitioned from a Be phase to a Be-shell phase during this period. Using the radiative transfer code \hdust\ we created a grid of $\sim100,000$ disk models for Pleione. We successfully reproduced the observed transition with a disk model that varies in inclination while maintaining an equatorial density of… ▽ More We acquired H$α$ spectroscopic observations from 2005 to 2019 showing Pleione has transitioned from a Be phase to a Be-shell phase during this period. Using the radiative transfer code \hdust\ we created a grid of $\sim100,000$ disk models for Pleione. We successfully reproduced the observed transition with a disk model that varies in inclination while maintaining an equatorial density of $ρ_0(r) = 3\times 10^{-11} (r/R_{eq})^{-2.7}~\rm{g~cm^{-3}}$, and an H$α$ emitting region extending to $15~\rm{R_{eq}}$. We use a precessing disk model to extrapolate the changing disk inclination over $120$ years and follow the variability in archival observations. The best-fit disk model precesses over a line of sight inclination between $\sim25\rm{^{\circ}}$ and $\sim144\rm{^{\circ}}$ with a precessional period of $\sim80.5$ years. Our precessing models match some of the observed variability but fail to reproduce all of the historical data available. Therefore, we propose an ad-hoc model based on our precessing disk model inspired by recent SPH simulations of similar systems, where the disk tears due to the tidal influence of a companion star. In this model, a single disk is slowly tilted to an angle of $30^{\circ}$ from the stellar equator over $34$ years. Then, the disk is torn by the companion's tidal torque, with the outer region separating from the innermost disk. The small inner disk returns to the stellar equator as mass injection remains constant. The outer disk precesses for $\sim15$ years before gradually dissipating. The process repeats every $34$ years and reproduces all trends in Pleione's variability. △ Less

Submitted 11 February, 2022; originally announced February 2022.

Comments: 19 pages, 15 figures

Untangling magnetic massive star properties with linear polarization variability and the Analytic Dynamical Magnetosphere model

Authors: M. S. Munoz, G. A. Wade, D. M. Faes, A. C. Carciofi, J. Labadie-Bartz

Abstract: Light scattered off particles can become linearly polarized. Stars surrounded by oblique, co-rotating envelopes are therefore expected to manifest periodic linear polarimetric variations. The electron scattering magnetospheres of magnetic massive stars are expected to be suitable candidates to observe this effect. In this paper, we present the first semi-analytical model capable of synthesizing th… ▽ More Light scattered off particles can become linearly polarized. Stars surrounded by oblique, co-rotating envelopes are therefore expected to manifest periodic linear polarimetric variations. The electron scattering magnetospheres of magnetic massive stars are expected to be suitable candidates to observe this effect. In this paper, we present the first semi-analytical model capable of synthesizing the continuum polarimetric signatures of magnetic O-type stars in an optically thin, single electron scattering limit. The purpose of this investigation is to improve our general understanding of magnetic hot stars by characterizing their polarimetric behaviour. Our linear polarization model is constructed by combining the analytical expressions for the polarimetric variations of an obliquely rotating envelope with the Analytic Dynamical Magnetosphere model to represent a physical model for the envelope density structure. We compute grids of model Stokes $Q$ and $U$ curves and show that their shapes are unique to the choice of inclination and obliquity angles. We apply our model to HD 191612, a prototypical Of?p-type star, having both polarimetric and photometric observations. We find that the polarimetric modulations are best reproduced with $i=19^{+12}_{-3}$$^\circ$, $β=71^{+3}_{-9}$$^\circ$, and $\log \dot{M}_{B=0}=-6.11^{+0.12}_{-0.06}$ [M$_{\odot}$ yr$^{-1}$]. These results agree with previous investigations of this star. By combining both polarimetric and photometric synthesis tools, we simultaneously model the observations thus adding further refinement of the wind and magnetic properties of HD 191612. △ Less

Submitted 22 December, 2021; originally announced December 2021.

Comments: Accepted by MNRAS

Interferometric detections of sdO companions orbiting three classical Be stars

Authors: R. Klement, G. H. Schaefer, D. R. Gies, L. Wang, D. Baade, Th. Rivinius, A. Gallenne, A. C. Carciofi, J. D. Monnier, A. Mérand, N. Anugu, S. Kraus, C. L. Davies, C. Lanthermann, T. Gardner, P. Wysocki, J. Ennis, A. Labdon, B. R. Setterholm, J. Le Bouquin

Abstract: Classical Be stars are possible products of close binary evolution, in which the mass donor becomes a hot, stripped O or B-type subdwarf (sdO/sdB), and the mass gainer spins up and grows a disk to become a Be star. While several Be+sdO binaries have been identified, dynamical masses and other fundamental parameters are available only for a single Be+sdO system, limiting the confrontation with bina… ▽ More Classical Be stars are possible products of close binary evolution, in which the mass donor becomes a hot, stripped O or B-type subdwarf (sdO/sdB), and the mass gainer spins up and grows a disk to become a Be star. While several Be+sdO binaries have been identified, dynamical masses and other fundamental parameters are available only for a single Be+sdO system, limiting the confrontation with binary evolution models. In this work, we present direct interferometric detections of the sdO companions of three Be stars 28 Cyg, V2119 Cyg, and 60 Cyg, all of which were previously found in UV spectra. For two of the three Be+sdO systems, we present first orbits and preliminary dynamical masses of the components, revealing that one of them could be the first identified progenitor of a Be/X-ray binary with a neutron star companion. These results provide new sets of fundamental parameters that are crucially needed to establish the evolutionary status and origin of Be stars. △ Less

Submitted 9 December, 2021; originally announced December 2021.

Comments: accepted to ApJ

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.

Comments: 38 pages, 16 figures, 1 table, topical collection in Astrophysics and Space Science: UV Spectropolarimetry for Stellar, Interstellar, and Exoplanetary Astrophysics with Polstar

Journal ref: Astrophysics and Space Science, Volume 367, Issue 12, article id.124, December 2022

Growth and Dissipation of Be Star Discs in Misaligned Binary Systems

Authors: M. Suffak, C. E. Jones, A. C. Carciofi

Abstract: We use a three-dimensional smoothed particle hydrodynamics code to simulate growth and dissipation of Be star discs in systems where the binary orbit is misaligned with respect to the spin axis of the primary star. We investigate six different scenarios of varying orbital period and misalignment angle, feeding the disc at a constant rate for 100 orbital periods, and then letting the disc dissipate… ▽ More We use a three-dimensional smoothed particle hydrodynamics code to simulate growth and dissipation of Be star discs in systems where the binary orbit is misaligned with respect to the spin axis of the primary star. We investigate six different scenarios of varying orbital period and misalignment angle, feeding the disc at a constant rate for 100 orbital periods, and then letting the disc dissipate for 100 orbital periods. During the disc growth phase, we find that the binary companion tilts the disc away from its initial plane at the equator of the primary star before settling to a constant orientation after 40 to 50 orbital periods. While the mass-injection into the disc is ongoing, the tilting of the disc can cause material to reaccrete onto the primary star prematurely. Once disc dissipation begins, usually the disc precesses about the binary companion's orbital axis with precession periods ranging from 20 to 50 orbital periods. In special cases we detect phenomena of disc tearing, as well as Kozai-Lidov oscillations of the disc. These oscillations reach a maximum eccentricity of about 0.6, and a minimum inclination of about \ang{20} with respect to the binary's orbit. We also find the disc material to have highly eccentric orbits beyond the transition radius, where the disc changes from being dominated by viscous forces, to heavily controlled by the companion star, in contrast to its nearly circular motion inward of the transition radius. Finally, we offer predictions to how these changes will affect Be star observables. △ Less

Submitted 15 October, 2021; originally announced October 2021.

Comments: 14 pages, 14 figures

Outbursts and stellar properties of the classical Be star HD 6226

Authors: Noel D. Richardson, Olivier Thizy, Jon E. Bjorkman, Alex Carciofi, Amanda C. Rubio, Joshua D. Thomas, Karen S. Bjorkman, Jonathan Labadie-Bartz, Matheus Genaro, John P. Wisniewski, Luqian Wang, Douglas R. Gies, S. Drew Chojnowski, Andrea Daly, Thompson Edwards, Carlie Fowler, Allison D. Gullingsrud, Nolan Habel, David J. James, Emily Kehoe, Heidi Kuchta, Alexis Lane, Anatoly Miroshnichenko, Ashish Mishra, Herbert Pablo , et al. (36 additional authors not shown)

Abstract: The bright and understudied classical Be star HD 6226 has exhibited multiple outbursts in the last several years during which the star grew a viscous decretion disk. We analyze 659 optical spectra of the system collected from 2017-2020, along with a UV spectrum from the Hubble Space Telescope and high cadence photometry from both TESS and the KELT survey. We find that the star has a spectral type… ▽ More The bright and understudied classical Be star HD 6226 has exhibited multiple outbursts in the last several years during which the star grew a viscous decretion disk. We analyze 659 optical spectra of the system collected from 2017-2020, along with a UV spectrum from the Hubble Space Telescope and high cadence photometry from both TESS and the KELT survey. We find that the star has a spectral type of B2.5IIIe, with a rotation rate of 74% of critical. The star is nearly pole-on with an inclination of $13.4$ degree. We confirm the spectroscopic pulsational properties previously reported, and report on three photometric oscillations from KELT photometry. The outbursting behavior is studied with equivalent width measurements of H$α$ and H$β$, and the variations in both of these can be quantitatively explained with two frequencies through a Fourier analysis. One of the frequencies for the emission outbursts is equal to the difference between two photometric oscillations, linking these pulsation modes to the mass ejection mechanism for some outbursts. During the TESS observation time period of 2019 October 7 to 2019 November 2, the star was building a disk. With a large dataset of H$α$ and H$β$ spectroscopy, we are able to determine the timescales of dissipation in both of these lines, similar to past work on Be stars that has been done with optical photometry. HD 6226 is an ideal target with which to study the Be disk-evolution given its apparent periodic nature, allowing for targeted observations with other facilities in the future. △ Less

Submitted 22 September, 2021; originally announced September 2021.

Comments: 48 pages including appendices, accepted to MNRAS. MNRAS online version has a 3D printed dynamical spectrum in the appendix. Contact for this section (due to size limitations)

Spiral Density Enhancements in Be Binary Systems

Authors: Isabelle H. Cyr, Carol E. Jones, Alex C. Carciofi, Camryn Steckel, Christopher Tycner, Atsuo T. Okazaki

Abstract: We use a smoothed particle hydrodynamics (SPH) code to examine the effects of a binary companion on a Be star disk for a range of disk viscosities and misalignment angles, i.e. the angle between the orbital plane and the primary's spin axis. The density structures in the disk due to the tidal interaction with the binary companion are investigated. Expanding on our previous work, the shape and dens… ▽ More We use a smoothed particle hydrodynamics (SPH) code to examine the effects of a binary companion on a Be star disk for a range of disk viscosities and misalignment angles, i.e. the angle between the orbital plane and the primary's spin axis. The density structures in the disk due to the tidal interaction with the binary companion are investigated. Expanding on our previous work, the shape and density structure of density enhancements due to the binary companion are analyzed and the changes in observed interferometric features due to these orbiting enhancements are also predicted. We find that larger misalignment angles and viscosity values result in more tightly wound spiral arms with densities that fall-off more slowly with radial distance from the central star. We show that the orbital phase has very little effect on the structure of the spiral density enhancements. We demonstrate that these spiral features can be detected with an interferometer in H$α$ and K-band emission. We also show that the spiral features affect the axis ratios determined by interferometry depending on the orientation of these features and the observer. For example, our simulations show that the axis ratios can vary by 20% for our co-planar binary disk system depending on the location of the disk density enhancements. △ Less

Submitted 13 July, 2021; originally announced July 2021.

Comments: 13 pages, 14 figures

Journal ref: MNRAS, 497, 3525 (2020)

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.

Comments: 20 pages, 16 figures

A Multi-observing-technique Study of the Dynamical Evolution of the Viscous Disk around the Be Star $ω$ CMa

Authors: Mohammad R. Ghoreyshi, Alex C. Carciofi, Carol E. Jones, Daniel M. Faes, Dietrich Baade, Thomas Rivinius

Abstract: The observed emission lines of Be stars originate from a circumstellar Keplerian disk that are generally well explained by the Viscous Decretion Disk model. In an earlier work we performed the modeling of the full light curve of the bright Be star $ω$ CMa (Ghoreyshi et al. 2018) with the 1-D time-dependent hydrodynamics code SINGLEBE and the Monte Carlo radiative-transfer code HDUST. We used the… ▽ More The observed emission lines of Be stars originate from a circumstellar Keplerian disk that are generally well explained by the Viscous Decretion Disk model. In an earlier work we performed the modeling of the full light curve of the bright Be star $ω$ CMa (Ghoreyshi et al. 2018) with the 1-D time-dependent hydrodynamics code SINGLEBE and the Monte Carlo radiative-transfer code HDUST. We used the V -band light curve that probes the inner disk through four disk formation and dissipation cycles. This new study compares predictions of the same set of model parameters with time-resolved photometry from the near UV through the mid-infrared, comprehensive series of optical spectra, and optical broad-band polarimetry, that overall represent a larger volume of the disk. Qualitatively, the models reproduce the trends in the observed data due to the growth and decay of the disk. However, quantitative differences exist, e.g., an overprediction of the flux increasing with wavelength, too slow decreases in Balmer emission-line strength that are too slow during disk dissipation, and the discrepancy between the range of polarimetric data and the model. We find that a larger value of the viscosity parameter alone, or a truncated disk by a companion star, reduces these discrepancies by increasing the dissipation rate in the outer regions of the disk. △ Less

Submitted 8 February, 2021; originally announced February 2021.

Comments: 25 pages, 19 figures

Short-term variability and mass loss in Be stars VI. Frequency groups in $γ$ Cas detected by TESS

Authors: Jonathan Labadie-Bartz, Dietrich Baade, Alex C. Carciofi, Amanda Rubio, Thomas Rivinius, Camilla C. Borre, Christophe Martayan, Robert J. Siverd

Abstract: In photometry of $γ$ Cas (B0.5 IVe) from the SMEI and BRITE-Constellation satellites, indications of low-order non-radial pulsation have recently been found, which would establish an important commonality with the class of classical Be stars at large. New photometry with the TESS satellite has detected three frequency groups near 1.0 ($g1$), 2.4 ($g2$), and 5.1 ($g3$) d$^{-1}$, respectively. Some… ▽ More In photometry of $γ$ Cas (B0.5 IVe) from the SMEI and BRITE-Constellation satellites, indications of low-order non-radial pulsation have recently been found, which would establish an important commonality with the class of classical Be stars at large. New photometry with the TESS satellite has detected three frequency groups near 1.0 ($g1$), 2.4 ($g2$), and 5.1 ($g3$) d$^{-1}$, respectively. Some individual frequencies are nearly harmonics or combination frequencies but not exactly so. Frequency groups are known from roughly three quarters of all classical Be stars and also from pulsations of $β$ Cep, SPB, and $γ$ Dor stars and, therefore, firmly establish $γ$ Cas as a non-radial pulsator. The total power in each frequency group is variable. An isolated feature exists at 7.57 d$^{-1}$ and, together with the strongest peaks in the second and third groups ordered by increasing frequency ($g2$ and $g3$), is the only one detected in all three TESS sectors. The former long-term 0.82 d$^{-1}$ variability would fall into $g1$ and has not returned at a significant level, questioning its attribution to rotational modulation. Low-frequency stochastic variability is a dominant feature of the TESS light curve, possibly caused by internal gravity waves excited at the core-envelope interface. These are known to be efficient at transporting angular momentum outward, and may also drive the oscillations that constitute $g1$ and $g2$. The hard X-ray flux of $γ$ Cas is the only remaining major property that distinguishes this star from the class of classical Be stars. △ Less

Submitted 11 December, 2020; originally announced December 2020.

Comments: 19 pages, 11 figures. To be published in MNRAS

Classifying Be star variability with TESS I: the southern ecliptic

Authors: Jonathan Labadie-Bartz, Alex C. Carciofi, Tajan Henrique de Amorim, Amanda Rubio, André Luiz, Pedro Ticiani dos Santos, Keegan Thomson-Paressant

Abstract: TESS photometry is analyzed for 432 classical Be stars observed in the first year of the mission. The often complex and diverse variability of each object in this sample is classified to obtain an understanding of the behavior of this class as a population. 98\% of the systems are variable above the noise level, with timescales spanning nearly the entire range of what is accessible with TESS, from… ▽ More TESS photometry is analyzed for 432 classical Be stars observed in the first year of the mission. The often complex and diverse variability of each object in this sample is classified to obtain an understanding of the behavior of this class as a population. 98\% of the systems are variable above the noise level, with timescales spanning nearly the entire range of what is accessible with TESS, from tens of minutes to tens of days. The variability seen with TESS is summarized as follows. Nearly every system contains multiple periodic signals in the frequency regime between about 0.5 -- 4 d$^{-1}$. One or more groups of closely-spaced frequencies is the most common feature, present in 85\% of the sample. Among the Be stars with brightening events that are characteristic of mass ejection episodes (17\% of the full sample, or 30\% of early-type stars), all have at least one frequency group, and the majority of these (83\%) show a concurrent temporary amplitude enhancement in one or more frequency groups. About one third of the sample is dominated by low frequency ($f < 0.5$ d$^{-1}$, and often much lower) variability. Stochastic signals are prominent in about 26\% of the sample, with varying degrees of intensity. Higher frequency signals ($6 < f < 15$ d$^{-1}$) are sometimes seen (in 14\% of the sample) and in most cases likely reflect p mode pulsation. In rare cases ($\sim$3\%), even higher frequencies beyond the traditional p mode regime ($f > 15$ d$^{-1}$) are observed. △ Less

Submitted 26 October, 2020; originally announced October 2020.

Comments: Submitted to AAS Journals. 52 pgs, 12 figures, 2 tables

The photometric and polarimetric variability of magnetic O-type stars

Authors: Melissa Munoz, Gregg Wade, Daniel Faes, Alex Carciofi

Abstract: Massive star winds are important contributors to the energy, momentum and chemical enrichment of the interstellar medium. Strong, organized and predominantly dipolar magnetic fields have been firmly detected in a small subset of massive O-type stars. Magnetic massive stars are known to exhibit phase-locked variability of numerous observable quantities that is hypothesized to arise due to the prese… ▽ More Massive star winds are important contributors to the energy, momentum and chemical enrichment of the interstellar medium. Strong, organized and predominantly dipolar magnetic fields have been firmly detected in a small subset of massive O-type stars. Magnetic massive stars are known to exhibit phase-locked variability of numerous observable quantities that is hypothesized to arise due to the presence of an obliquely rotating magnetosphere formed via the magnetic confinement of their strong outflowing winds. Analyzing the observed modulations of magnetic O-type stars is thus a key step towards the better understanding of the physical processes that occur within their magnetospheres. The dynamical processes that lead to the formation of a magnetosphere are formally solved utilizing complex MHD simulations. Recently, an Analytic Dynamical Magnetosphere (ADM) model has been developed that can quickly be employed to compute the time-averaged density, temperature and velocity gradients within a dynamical magnetosphere. Here, we exploit the ADM model to compute photometric and polarimetric observables of magnetic Of?p stars, to test geometric models inferred from magnetometry. We showcase important results on the prototypical Of?p-type star HD 191612, that lead to a better characterization of massive star wind and magnetic properties. △ Less

Submitted 28 April, 2020; originally announced April 2020.

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.

Comments: 24 pages, 17 figures, accepted for publication in Astronomy & Astrophysics on 06/02/2020. Abstract shortened to fit within arXiv character limit

Journal ref: A&A 636, A110 (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.

Comments: 16 pages, 12 figures, 4 tables, accepted for publishing in ApJ

Prevalence of SED turndown among classical Be stars: Are all Be stars close binaries?

Authors: Robert Klement, A. C. Carciofi, T. Rivinius, R. Ignace, L. D. Matthews, K. Torstensson, D. Gies, R. G. Vieira, N. D. Richardson, A. Domiciano de Souza, J. E. Bjorkman, G. Hallinan, D. M. Faes, B. Mota, A. D. Gullingsrud, C. de Breuck, P. Kervella, M. Curé, D. Gunawan

Abstract: Rapid rotation is a fundamental characteristic of classical Be stars and a crucial property allowing for the formation of their circumstellar disks. Past evolution in a mass and angular momentum transferring binary system offers a plausible solution to how Be stars attained their fast rotation. Although the subdwarf remnants of mass donors in such systems should exist in abundance, only a few have… ▽ More Rapid rotation is a fundamental characteristic of classical Be stars and a crucial property allowing for the formation of their circumstellar disks. Past evolution in a mass and angular momentum transferring binary system offers a plausible solution to how Be stars attained their fast rotation. Although the subdwarf remnants of mass donors in such systems should exist in abundance, only a few have been confirmed due to tight observational constraints. An indirect method of detecting otherwise hidden companions is offered by their effect on the outer parts of Be star disks, which are expected to be disrupted or truncated. In the context of the IR and radio continuum excess radiation originating in the disk, the disk truncation can be revealed by a turndown in the spectral energy distribution due to reduced radio flux levels. In this work we search for signs of spectral turndown in a sample of 57 classical Be stars with radio data, which include new data for 23 stars and the longest wavelength detections so far ($λ\approx$ 10\,cm) for 2 stars. We confidently detect the turndown for all 26 stars with sufficient data coverage (20 of which are not known to have close binary companions). For the remaining 31 stars, data are inconclusive as to whether the turndown is present or not. The analysis suggests that many if not all Be stars have close companions influencing their outer disks. If confirmed to be subdwarf companions, the mass transfer spin-up scenario might explain the existence of the vast majority of classical Be stars. △ Less

Submitted 11 December, 2019; v1 submitted 26 September, 2019; originally announced September 2019.

Journal ref: 2019ApJ...885..147K

Modelling the periodical variations in multiband polarisation and photometry for discs of binary Be stars

Authors: Despina Panoglou, Marcelo Borges Fernandes, Dietrich Baade, Daniel M. Faes, Thomas Rivinius, Alex C. Carciofi, Atsuo T. Okazaki

Abstract: The tidal interaction of a Be star with a binary companion forms two spiral arms that cause orbital modulation of the Be disc structure. The aim of this work is to identify observables in which this modulation is apparent. The structure of a Be disc in a coplanar circular binary system is computed with a smoothed-particle hydrodynamics code, and a radiation transfer code calculates the spectral en… ▽ More The tidal interaction of a Be star with a binary companion forms two spiral arms that cause orbital modulation of the Be disc structure. The aim of this work is to identify observables in which this modulation is apparent. The structure of a Be disc in a coplanar circular binary system is computed with a smoothed-particle hydrodynamics code, and a radiation transfer code calculates the spectral energy distribution. Line depolarisation was confirmed, with polarisation profiles nearly reverse to emission-line profiles. The continuum flux maximizes for pole-on discs, but photometric variability maximizes for edge-on discs. The linear polarisation exhibits one or two maxima per orbital cycle. While polarisation variability in visible passbands is important only at low inclinations, infrared bands may demonstrate high orbital variability even at large inclinations. More evident is the modulation in the polarisation angle (PA) for low inclinations. The latter can be used to track azimuthal asymmetries for pole-on discs, where the spectroscopic variability in the violet-to-red (V/R) emission-component ratio disappears. PA reversals coincide with phases where V/R=1, tracking lines of sight directed towards regions where the approaching and receding arms overlap. Continuum flux and polarisation are mostly in phase for neighbouring wavelength regions. It is suggested that studies of non-symmetric discs distorted by tidal forces from a secondary star may be used to study disc variabilities of other origins. △ Less

Submitted 18 April, 2019; originally announced April 2019.

Comments: 20 pages, 19 figures; accepted for publication in MNRAS

Short-term variability and mass loss in Be stars IV. Two groups of closely spaced, approximately equidistant frequencies in three decades of space photometry of $ν$ Puppis (B7-8 IIIe)

Authors: D. Baade, A. Pigulski, Th. Rivinius, L. Wang, Ch. Martayan, G. Handler, D. Panoglou, A. C. Carciofi, R. Kuschnig, A. Mehner, A. F. J. Moffat, H. Pablo, S. M. Rucinski, G. A. Wade, W. W. Weiss, K. Zwintz

Abstract: In early-type Be stars, groups of nonradial pulsation (NRP) modes with numerically related frequencies may be instrumental for the release of excess angular momentum through mass-ejection events. Difference and sum/harmonic frequencies often form additional groups. The goal of this study is to find out whether a similar frequency pattern occurs in the cooler third-magnitude B7-8\,IIIe shell star… ▽ More In early-type Be stars, groups of nonradial pulsation (NRP) modes with numerically related frequencies may be instrumental for the release of excess angular momentum through mass-ejection events. Difference and sum/harmonic frequencies often form additional groups. The goal of this study is to find out whether a similar frequency pattern occurs in the cooler third-magnitude B7-8\,IIIe shell star $ν$ Pup. Time-series analyses are performed of space photometry with BRITE-Constellation (2015, 2016/17, and 2017/18), SMEI (2003--011), and Hipparcos (1989-1993). Two IUE SWP and 27 optical echelle spectra spanning 20 years were retrieved from various archives. The optical spectra exhibit no anomalies or well-defined variabilities. A magnetic field was not detected. All three photometry satellites recorded variability near 0.656 c/d which is resolved into three features separated by ~0.0021 c/d. First harmonics form a second frequency group, also spaced by ~0.0021 c/d. The frequency spacing is very nearly but not exactly equidistant. Variability near 0.0021 c/d was not detected. The long-term frequency stability could be used to derive meaningful constraints on the properties of a putative companion star. The IUE spectra do not reveal the presence of a hot subluminous secondary. $ν$\,Pup is another Be star exhibiting an NRP variability pattern with long-term constancy and underlining the importance of combination frequencies and frequency groups. The star is a good target for efforts to identify an effectively single Be star. △ Less

Submitted 30 August, 2018; originally announced August 2018.

Comments: 16 pages, 12 figures, 5 tables, submitted to Astronomy & Astrophysics

Journal ref: A&A 620, A145 (2018)

The life cycles of Be viscous decretion discs: The case of ω CMa

Authors: M. R. Ghoreyshi, A. C. Carciofi, L. R. Rimulo, R. G. Vieira, D. M. Faes, D. Baade, J. E. Bjorkman, S. Otero, Th. Rivinius

Abstract: We analyzed V-band photometry of the Be star ω CMa, obtained during the last four decades, during which the star went through four complete cycles of disc formation and dissipation. The data were simulated by hydrodynamic models based on a time-dependent implementation of the viscous decretion disc (VDD) paradigm, in which a disc around a fast-spinning Be star is formed by material ejected by the… ▽ More We analyzed V-band photometry of the Be star ω CMa, obtained during the last four decades, during which the star went through four complete cycles of disc formation and dissipation. The data were simulated by hydrodynamic models based on a time-dependent implementation of the viscous decretion disc (VDD) paradigm, in which a disc around a fast-spinning Be star is formed by material ejected by the star and driven to progressively larger orbits by means of viscous torques. Our simulations offer a good description of the photometric variability during phases of disc formation and dissipation, which suggests that the VDD model adequately describes the structural evolution of the disc. Furthermore, our analysis allowed us to determine the viscosity parameter α, as well as the net mass and angular momentum (AM) loss rates. We find that α is variable, ranging from 0.1 to 1.0, not only from cycle to cycle but also within a given cycle. Additionally, build-up phases usually have larger values of α than the dissipation phases. Furthermore, during dissipation the outward AM flux is not necessarily zero, meaning that ω CMa does not experience a true quiescence but, instead, switches between a high to a low AM loss rate during which the disc quickly assumes an overall lower density but never zero. We confront the average AM loss rate with predictions from stellar evolution models for fast-rotating stars, and find that our measurements are smaller by more than one order of magnitude. △ Less

Submitted 11 June, 2018; originally announced June 2018.

Comments: 16 pages, 16 figures

Orbital solution leading to an acceptable interpretation for the enigmatic gamma-ray binary HESS J0632+057

Authors: Yuki Moritani, Takafumi Kawano, Sho Chimasu, Akiko Kawachi, Hiromitsu Takahashi, Jumpei Takata, Alex C. Carciofi

Abstract: High-dispersion spectroscopic monitoring of HESS J0632+057 has been carried out over four orbital cycles in order to search for orbital modulation, covering the entire orbital phase. We have measured radial velocity of H$α$ emission line with the method introduced by Shafter et al. (1986), which has been successfully applied to some Be stars. The velocity is seen to increase much earlier than expe… ▽ More High-dispersion spectroscopic monitoring of HESS J0632+057 has been carried out over four orbital cycles in order to search for orbital modulation, covering the entire orbital phase. We have measured radial velocity of H$α$ emission line with the method introduced by Shafter et al. (1986), which has been successfully applied to some Be stars. The velocity is seen to increase much earlier than expected for the orbital period of 315 days, and much more steeply than expected at around "apastron". The period of the H$α$ modulation is found to be as $308^{+26}_{-23}$ days. We have also analyzed Swift/XRT data from 2009 to 2015 to study the orbital modulation, selecting the data with good statistics ($\geq$ 30 counts). With additional two-year data to the previous works, the orbital period has been updated to $313^{+11}_{-8}$ days, which is consistent with the previous X-ray periods and the spectroscopic one. The past XMM-Newton and Chandra observations prefer the period of 313 days. With the new period, assuming that H$α$ velocities accurately trace the motion of the Be star, we have derived a new set of the orbital parameters. In the new orbit, which is less eccentric ($e \sim 0.6$), two outbursts occur after apastron, and just after periastron. Besides, the column density in bright phase ($4.7^{+0.9}_{-08}\times10^{21}\;\mathrm{cm^{-2}}$) is higher than in faint phase ($2.2\pm0.5\times10^{21}\;\mathrm{cm^{-2}}$). These facts suggest that outbursts occur when the compact object passes nearby/through the Be disk. The mass function implies that mass of the compact object is less than 2.5 $\mathrm{M_{sun}}$ assuming that the mass of the Be star is 13.2--18.2 $\mathrm{M_{sun}}$ (Aragona at al. 2010) unless the inclination is extremely small. The photon index indicates that the spectra becomes softer when the system is bright. These suggest that the compact object is a pulsar. △ Less

Submitted 11 April, 2018; originally announced April 2018.

Comments: 15 pages, 7 figures, accepted to PASJ

The life cycles of Be viscous decretion discs: fundamental disc parameters of 54 SMC Be stars

Authors: Leandro Rocha Rímulo, Alex Cavaliéri Carciofi, Rodrigo Georgetti Vieira, Thomas Rivinius, Daniel Moser Faes, André Luiz Figueiredo, Jon Eric Bjorkman, Cyril Georgy, Mohammad Reza Ghoreyshi, Igor Soszyński

Abstract: Be stars are main-sequence massive stars with emission features in their spectrum, which originates in circumstellar gaseous discs. Even though the viscous decretion disc (VDD) model can satisfactorily explain most observations, two important physical ingredients, namely the magnitude of the viscosity ($α$) and the disk mass injection rate, remain poorly constrained. The light curves of Be stars t… ▽ More Be stars are main-sequence massive stars with emission features in their spectrum, which originates in circumstellar gaseous discs. Even though the viscous decretion disc (VDD) model can satisfactorily explain most observations, two important physical ingredients, namely the magnitude of the viscosity ($α$) and the disk mass injection rate, remain poorly constrained. The light curves of Be stars that undergo events of disc formation and dissipation offer an opportunity to constrain these quantities. A pipeline was developed to model these events that uses a grid of synthetic light curves, computed from coupled hydrodynamic and radiative transfer calculations. A sample of 54 Be stars from the OGLE survey of the Small Magellanic Cloud (SMC) was selected for this study. Because of the way our sample was selected (bright stars with clear disc events), it likely represents the densest discs in the SMC. Like their siblings in the Galaxy, the mass of the disc in the SMC increases with the stellar mass. The typical mass and angular momentum loss rates associated with the disk events are of the order of $\sim$$10^{-10}\, M_\odot\,\mathrm{yr^{-1}}$ and $\sim$$5\times 10^{36}\, \mathrm{g\, cm^{2}\, s^{-2}}$, respectively. The values of $α$ found in this work are typically of a few tenths, consistent with recent results in the literature and with the ones found in dwarf novae, but larger than current theory predicts. Considering the sample as a whole, the viscosity parameter is roughly two times larger at build-up ($\left\langleα_\mathrm{bu}\right\rangle = 0.63$) than at dissipation ($\left\langleα_\mathrm{d}\right\rangle = 0.26$). Further work is necessary to verify whether this trend is real or a result of some of the model assumptions. △ Less

Submitted 21 February, 2018; originally announced February 2018.

BRITEning up the Be Phenomenon

Authors: D. Baade, Th. Rivinius, A. Pigulski, D. Panoglou, A. Carciofi, G. Handler, R. Kuschnig, Ch. Martayan, A. Mehner, A. F. J. Moffat, H. Pablo, A. Popowicz, S. M. Rucinski, G. A. Wade, W. W. Weiss, K. Zwintz

Abstract: Observations of 25 Ori much expand the picture derived of other early-type Be stars with BRITE and SMEI. Two instead of one difference frequencies rule the variability: (a) The lower one, 0.0129 c/d, is the frequency of events with full amplitudes of 100-200 mmag which may signal mass loss possibly driven by the higher one, 0.1777 c/d. (b) Much of the entire power spectrum is a tightly woven netwo… ▽ More Observations of 25 Ori much expand the picture derived of other early-type Be stars with BRITE and SMEI. Two instead of one difference frequencies rule the variability: (a) The lower one, 0.0129 c/d, is the frequency of events with full amplitudes of 100-200 mmag which may signal mass loss possibly driven by the higher one, 0.1777 c/d. (b) Much of the entire power spectrum is a tightly woven network of combination frequencies: (i) Below 0.25\,c/d, numerous frequencies are difference frequencies. (ii) Many frequencies above 2.5 c/d can be represented as sum frequencies and in a few cases as harmonics. (iii) Many frequencies between 1.1 and 1.75\,c/d can be portrayed as parents of combination frequencies. The number and fraction of combination frequencies increases steeply with decreasing amplitude and and accuracy of the frequency matching. △ Less

Submitted 28 August, 2017; originally announced August 2017.

Comments: 8 pages, 0 figures, to appear in Proc. 3rd BRITE Science Conference, Lac Taureau, Quebec, Aug. 7-10, 2017 (eds. G.A. Wade et al.); the page footer is wrong (bug in template)

Short-term variability and mass loss in Be stars III. BRITE and SMEI satellite photometry of 28 Cygni

Authors: D. Baade, A. Pigulski, Th. Rivinius, A. C. Carciofi, D. Panoglou, M. Ghoreyshi, G. Handler, R. Kuschnig, A. F. J. Moffat, H. Pablo, A. Popowicz, G. A. Wade, W. W. Weiss, K. Zwintz

Abstract: The BRITE Constellation of nanosatellites obtained mmag photometry of 28 Cygni for 11 months in 2014-2016. Observations with the Solar Mass Ejection Imager in 2003-2010 and 118 H$α$ line profiles were added. For decades, 28 Cyg has exhibited four large-amplitude frequencies: two closely spaced frequencies of spectroscopically confirmed $g$ modes near 1.5 c/d, one slightly lower exophotospheric (… ▽ More The BRITE Constellation of nanosatellites obtained mmag photometry of 28 Cygni for 11 months in 2014-2016. Observations with the Solar Mass Ejection Imager in 2003-2010 and 118 H$α$ line profiles were added. For decades, 28 Cyg has exhibited four large-amplitude frequencies: two closely spaced frequencies of spectroscopically confirmed $g$ modes near 1.5 c/d, one slightly lower exophotospheric (Stefl) frequency, and at 0.05 c/d the difference frequency between the two g modes. This top-level framework is indistinguishable from eta Cen (Paper I), which is also very similar in spectral type, rotation rate, and viewing angle. The Stefl frequency is the only one that does not seem to be affected by the difference frequency. The amplitude of the latter undergoes large variations; around maximum the amount of near-circumstellar matter is increased, and the amplitude of the Stefl frequency grows by some factor. During such brightenings dozens of transient spikes appear in the frequency spectrum, concentrated in three groups. Only eleven frequencies were common to all years of BRITE observations. Be stars seem to be controlled by several coupled clocks, most of which are not very regular on timescales of weeks to months but function for decades. The combination of g modes to the low difference frequency and/or the atmospheric response to it appears significantly nonlinear. Like in eta Cen, the difference-frequency variability seems the main responsible for the modulation of the star-to-disc mass transfer in 28 Cyg. A hierarchical set of difference frequencies may reach the longest timescales known of the Be phenomenon. △ Less

Submitted 24 August, 2017; originally announced August 2017.

Comments: 17 pages, 21 figures, submitted to Astronomy & Astrophysics

Journal ref: A&A 610, A70 (2018)

Outbursts and Disk Variability in Be Stars

Authors: Jonathan Labadie-Bartz, S. Drew Chojnowski, David G. Whelan, Joshua Pepper, M. Virginia McSwain, Marcelo Borges Fernandes, John P. Wisniewski, Guy S. Stringfellow, Alex C. Carciofi, Robert J. Siverd, Amy L. Glazier, Sophie G. Anderson, Anthoni J. Caravello, Keivan G. Stassun, Michael B. Lund, Daniel J. Stevens, Joseph E. Rodriguez, David J. James, Rudolf B. Kuhn

Abstract: In order to study the growth and evolution of circumstellar disks around classical Be stars, we analyze optical time-series photometry from the KELT survey with simultaneous infrared and visible spectroscopy from the APOGEE survey and BeSS database for a sample of 160 Galactic classical Be stars. The systems studied here show variability including transitions from a diskless to a disk-possessing s… ▽ More In order to study the growth and evolution of circumstellar disks around classical Be stars, we analyze optical time-series photometry from the KELT survey with simultaneous infrared and visible spectroscopy from the APOGEE survey and BeSS database for a sample of 160 Galactic classical Be stars. The systems studied here show variability including transitions from a diskless to a disk-possessing state (and vice versa), and persistent disks that vary in strength, being replenished at either regularly or irregularly occurring intervals. We detect disk-building events (outbursts) in the light curves of 28\% of our sample. Outbursts are more commonly observed in early- (57\%), compared to mid- (27\%) and late-type (8\%) systems. A given system may show anywhere between 0 -- 40 individual outbursts in its light curve, with amplitudes ranging up to $\sim$0.5 mag and event durations between $\sim$2 -- 1000 days. We study how both the photometry and spectroscopy change together during active episodes of disk growth or dissipation, revealing details about the evolution of the circumstellar environment. We demonstrate that photometric activity is linked to changes in the inner disk, and show that, at least in some cases, the disk growth process is asymmetrical. Observational evidence of Be star disks both growing and clearing from the inside out is presented. The duration of disk buildup and dissipation phases are measured for 70 outbursts, and we find that the average outburst takes about twice as long to dissipate as it does to build up in optical photometry. Our analysis hints that dissipation of the inner disk occurs relatively slowly for late-type Be stars. △ Less

Submitted 23 January, 2018; v1 submitted 8 August, 2017; originally announced August 2017.

Comments: Published in The Astronomical Journal. 23 pages, 25 figures, 2 tables

Journal ref: Labadie-Bartz, Jonathan et al. 2018 AJ 155 53

The viscous disk properties of 80 Be stars

Authors: Rodrigo G. Vieira, Alex C. Carciofi

Abstract: We investigated the properties of $80$ Be star disks, using the time dependent viscous decretion disk (VDD) model to interpret their mid-infrared continuum fluxes. Besides estimating the typical parameter ranges for such disks, this approach allowed us to relate the disk structure to its evolutionary status. In this contribution, we present a brief summary of the recently published results, and di… ▽ More We investigated the properties of $80$ Be star disks, using the time dependent viscous decretion disk (VDD) model to interpret their mid-infrared continuum fluxes. Besides estimating the typical parameter ranges for such disks, this approach allowed us to relate the disk structure to its evolutionary status. In this contribution, we present a brief summary of the recently published results, and discuss some properties and time-scales of the VDD model solution. △ Less

Submitted 10 July, 2017; originally announced July 2017.

Comments: 6 pages, 3 figures

Journal ref: ASP Conference Series, Vol. 508. San Francisco: Astronomical Society of the Pacific, 2017, p. 33

The life cycles of Be viscous decretion discs: time-dependent modelling of infrared continuum observations

Authors: Rodrigo G. Vieira, Alex C. Carciofi, Jon E. Bjorkman, Thomas Rivinius, Dietrich Baade, Leandro R. Rímulo

Abstract: We apply the viscous decretion disc (VDD) model to interpret the infrared disc continuum emission of 80 Be stars observed in different epochs. In this way, we determined 169 specific disc structures, namely their density scale, $ρ_0$, and exponent, $n$. We found that the $n$ values range mainly between $1.5$ and $3.5$, and $ρ_0$ varies between $10^{-12}$ and $10^{-10}\,\mathrm{g\,cm^{-3}}$, with a… ▽ More We apply the viscous decretion disc (VDD) model to interpret the infrared disc continuum emission of 80 Be stars observed in different epochs. In this way, we determined 169 specific disc structures, namely their density scale, $ρ_0$, and exponent, $n$. We found that the $n$ values range mainly between $1.5$ and $3.5$, and $ρ_0$ varies between $10^{-12}$ and $10^{-10}\,\mathrm{g\,cm^{-3}}$, with a peak close to the lower value. Our large sample also allowed us to firmly establish that the discs around early-type stars are denser than in late-type stars. Additionally, we estimated the disc mass decretion rates and found that they range between $10^{-12}$ and $10^{-9}\,\mathrm{M_{\odot}\,yr^{-1}}$. These values are compatible with recent stellar evolution models of fast-rotating stars. One of the main findings of this work is a correlation between the $ρ_0$ and $n$ values. In order to find out whether these relations can be traced back to the evolution of discs or have some other origin, we used the VDD model to calculate temporal sequences under different assumptions for the time profile of the disc mass injection. The results support the hypothesis that the observed distribution of disc properties is due to a common evolutionary path. In particular, our results suggest that the timescale for disc growth, during which the disc is being actively fed by mass injection episodes, is shorter than the timescale for disc dissipation, when the disc is no longer fed by the star and dissipates as a result of the viscous diffusion of the disc material. △ Less

Submitted 10 July, 2017; originally announced July 2017.

Comments: 21 pages, 15 figures, 4 tables, 2016 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society

Be discs in binary systems: II. Misaligned orbits

Authors: Isabelle H. Cyr, Carol E. Jones, Despina Panoglou, Alex C. Carciofi, Atsuo T. Okazaki

Abstract: We use a smoothed particle hydrodynamics (SPH) code to examine the effects of misaligned binary companions on Be star discs. We systematically vary the degree of misalignment between the disc and the binary orbit, as well as the disc viscosity and orbital period to study their effects on the density in the inner and outer parts of the disc. We find that varying the degree of misalignment, the visc… ▽ More We use a smoothed particle hydrodynamics (SPH) code to examine the effects of misaligned binary companions on Be star discs. We systematically vary the degree of misalignment between the disc and the binary orbit, as well as the disc viscosity and orbital period to study their effects on the density in the inner and outer parts of the disc. We find that varying the degree of misalignment, the viscosity, and the orbital period affects both the truncation radius and the density structure of the outer disc, while the inner disc remains mostly unaffected. We also investigate the tilting of the disc in the innermost part of the disc and find the tilt increases with radius until reaching a maximum around 5 stellar radii. The direction of the line of nodes, with respect to the equator of the central star, is found to be offset compared to the orbital line of nodes, and to vary periodically in time, with a period of half a orbital phase. We also compare the scale height of our discs with the analytical scale height of an isothermal disc, which increases with radius as $r^{1.5}$. We find that this formula reproduces the scale height well for both aligned and misaligned systems but underestimates the scale height in regions of the disc where density enhancements develop. △ Less

Submitted 21 June, 2017; originally announced June 2017.

Comments: 11 pages, 10 figures, to be published in MNRAS

Be discs in coplanar circular binaries: Phase-locked variations of emission lines

Authors: Despina Panoglou, Daniel M. Faes, Alex C. Carciofi, Atsuo T. Okazaki, Dietrich Baade, Thomas Rivinius, Marcelo Borges Fernandes

Abstract: The first results of radiative transfer calculations on decretion discs of binary Be stars are presented. A smoothed particle hydrodynamics code computes the structure of Be discs in coplanar circular binary systems for a range of orbital and disc parameters. The resulting disc configuration consists of two spiral arms, and can be given as input into a Monte Carlo code, which calculates the radiat… ▽ More The first results of radiative transfer calculations on decretion discs of binary Be stars are presented. A smoothed particle hydrodynamics code computes the structure of Be discs in coplanar circular binary systems for a range of orbital and disc parameters. The resulting disc configuration consists of two spiral arms, and can be given as input into a Monte Carlo code, which calculates the radiative transfer along the line of sight for various observational coordinates. Making use of the property of steady disc structure in coplanar circular binaries, observables are computed as functions of the orbital phase. Orbital-phase series of line profiles are given for selected parameter sets under various viewing angles, to allow comparison with observations. Flat-topped profiles with and without superimposed multiple structures are reproduced, showing, for example, that triple-peaked profiles do not have to be necessarily associated with warped discs and misaligned binaries. It is demonstrated that binary tidal effects give rise to phase-locked variability of the violet-to-red (V/R) ratio of hydrogen emission lines. The V/R ratio exhibits two maxima per cycle; in certain cases those maxima are equal, leading to a clear new V/R cycle every half orbital period. This study opens a way in identifying binaries and in constraining the parameters of binary systems that exhibit phase-locked variations induced by tidal interaction with a companion star. △ Less

Submitted 22 September, 2017; v1 submitted 22 April, 2017; originally announced April 2017.

Comments: 12 pages, 16 figures; accepted for publication by MNRAS

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.

Comments: 20 pages, 8 figures

Journal ref: A&A 601, A74 (2017)

The viscous decretion disk model of the classical Be star $β$ CMi revisited

Authors: Robert Klement, Alex Cavalieri Carciofi, Thomas Rivinius

Abstract: We revisit the viscous decretion disk (VDD) model of the classical Be star $β$ CMi as presented by \citet{klement} using an updated version of the radiative transfer code {\ttfamily HDUST}. A software bug was causing the mean intensities to be slightly underestimated in the equatorial region of the disk, with small but detectable effects on the disk temperature and emergent spectrum. The new model… ▽ More We revisit the viscous decretion disk (VDD) model of the classical Be star $β$ CMi as presented by \citet{klement} using an updated version of the radiative transfer code {\ttfamily HDUST}. A software bug was causing the mean intensities to be slightly underestimated in the equatorial region of the disk, with small but detectable effects on the disk temperature and emergent spectrum. The new model fixes an unexplained feature of the original model, which was able to reproduce the observations only when considering a dual density behavior: a steep density fall-off in the very inner parts of the disk followed by a shallower density profile. The new model is able to reproduce all the observables reasonably well using a single power-law for the density profile throughout the whole disk, as predicted by the VDD model. All the other original conclusions, most importantly the reported truncation of the disk at the distance of 35 stellar equatorial radii from the central star, remain unchanged. △ Less

Submitted 21 March, 2017; originally announced March 2017.

Comments: 6 pages, 3 figures, 1 table, conference proceedings ASPC 508

Journal ref: The B[e] Phenomenon: Forty Years of Studies, Edited by Anatoly Miroshnichenko, Sergey Zharikov, Daniela Korčáková and Marek Wolf. ASP Conference Series, Vol. 508. San Francisco: Astronomical Society of the Pacific, 2017, p.39

VLTI/PIONIER images the Achernar disk swell

Authors: G. Dalla Vedova, F. Millour, A. Domiciano de Souza, R. G. Petrov, D. Moser Faes, A. C. Carciofi, P. Kervella, T. Rivinius

Abstract: Context. The mechanism of disk formation around fast-rotating Be stars is not well understood. In particular, it is not clear which mechanisms operate, in addition to fast rotation, to produce the observed variable ejection of matter. The star Achernar is a privileged laboratory to probe these additional mechanisms because it is close, presents B-Be phase variations on timescales ranging from 6 yr… ▽ More Context. The mechanism of disk formation around fast-rotating Be stars is not well understood. In particular, it is not clear which mechanisms operate, in addition to fast rotation, to produce the observed variable ejection of matter. The star Achernar is a privileged laboratory to probe these additional mechanisms because it is close, presents B-Be phase variations on timescales ranging from 6 yr to 15 yr, a companion star was discovered around it, and probably presents a polar wind or jet. Aims. Despite all these previous studies, the disk around Achernar was never directly imaged. Therefore we seek to produce an image of the photosphere and close environment of the star. Methods. We used infrared long-baseline interferometry with the PIONIER/VLTI instrument to produce reconstructed images of the photosphere and close environment of the star over four years of observations. To study the disk formation, we compared the observations and reconstructed images to previously computed models of both the stellar photosphere alone (normal B phase) and the star presenting a circumstellar disk (Be phase). Results. The observations taken in 2011 and 2012, during the quiescent phase of Achernar, do not exhibit a disk at the detection limit of the instrument. In 2014, on the other hand, a disk was already formed and our reconstructed image reveals an extended H-band continuum excess flux. Our results from interferometric imaging are also supported by several H-alpha line profiles showing that Achernar started an emission-line phase sometime in the beginning of 2013. The analysis of our reconstructed images shows that the 2014 near-IR flux extends to 1.7 - 2.3 equatorial radii. Our model-independent size estimation of the H-band continuum contribution is compatible with the presence of a circumstellar disk, which is in good agreement with predictions from Be-disk models. △ Less

Submitted 8 March, 2017; originally announced March 2017.

Science cases for a visible interferometer

Authors: Philippe Stee, France Allard, Myriam Benisty, Lionel Bigot, Nicolas Blind, Henri Boffin, Marcelo Borges Fernandes, Alex Carciofi, Andrea Chiavassa, Orlagh Creevey, Pierre Cruzalebes, Willem-Jan de Wit, Armando Domiciano de Souza, Martin Elvis, Nicolas Fabas, Daniel Faes, Alexandre Gallenne, Carlos Guerrero Pena, Michel Hillen, Sebastian Hoenig, Michael Ireland, Pierre Kervella, Makoto Kishimoto, Nadia Kostogryz, Stefan Kraus , et al. (32 additional authors not shown)

Abstract: High spatial resolution is the key for the understanding various astrophysical phenomena. But even with the future E-ELT, single dish instruments are limited to a spatial resolution of about 4 mas in the visible. For the closest objects within our Galaxy most of the stellar photosphere remains smaller than 1 mas. With the success of long baseline interferometry these limitations were soom overcome… ▽ More High spatial resolution is the key for the understanding various astrophysical phenomena. But even with the future E-ELT, single dish instruments are limited to a spatial resolution of about 4 mas in the visible. For the closest objects within our Galaxy most of the stellar photosphere remains smaller than 1 mas. With the success of long baseline interferometry these limitations were soom overcome. Today low and high resolution interferometric instruments on the VLTI and CHARA offer an immense range of astrophysical studies. Combining more telescopes and moving to visible wavelengths broadens the science cases even more. With the idea of developing strong science cases for a future visible interferometer, we organized a science group around the following topics: pre-main sequence and main sequence stars, fundamental parameters, asteroseismology and classical pulsating stars, evolved stars, massive stars, active galactic nuclei (AGNs) and imaging techniques. A meeting was organized on the 15th and 16th of January, 2015 in Nice with the support of the Action Specific in Haute Resolution Angulaire (ASHRA), the Programme National en Physique Stellaire (PNPS), the Lagrange Laboratory and the Observatoire de la Cote d'Azur, in order to present these cases and to discuss them further for future visible interferometers. This White Paper presents the outcome of the exchanges. This book is dedicated to the memory of our colleague Olivier Chesneau who passed away at the age of 41. △ Less

Submitted 21 March, 2017; v1 submitted 7 March, 2017; originally announced March 2017.

Comments: White Paper prospective. This book is dedicated to the memory of our colleague Olivier Chesneau who passed away at the age of 41. v2 includes some corrections to text

An analysis of the V-band light curve of the Be star omega CMa with the viscous decretion disk model

Authors: Mohammad Reza Ghoreyshi, Alex Cavalieri Carciofi

Abstract: We analyze the V-band photometry data of the Be star omega CMa, observed during the last four decades. The data is fitted by hydrodynamic models based on the viscous decretion disk (VDD) theory, in which a disk around a fast-spinning Be star is formed by material ejected by the central star and driven to progressively wider orbits by means of viscous torques. For the first time, we apply the model… ▽ More We analyze the V-band photometry data of the Be star omega CMa, observed during the last four decades. The data is fitted by hydrodynamic models based on the viscous decretion disk (VDD) theory, in which a disk around a fast-spinning Be star is formed by material ejected by the central star and driven to progressively wider orbits by means of viscous torques. For the first time, we apply the model for both the disk build up and the dissipation phases. Our simulations offer a good description of the photometric variability in both phases, which suggests that the VDD model adequately describes the disk structural evolution. Furthermore, our analysis allowed us to determine the viscosity parameter (alpha) of the gas, as well as the net mass loss rate. We find that alpha is variable, ranging from 0.1 to 1.0, and that build up phases have larger values of alpha than the dissipation phases. Additionally, we find that, contrary to what is generally assumed, even during quiescence the outward mass flux is never zero, suggesting that the star alternates between a high mass-loss phase (outburst) and a low mass-loss phase (quiescence). △ Less

Submitted 22 February, 2017; originally announced February 2017.

A Spectroscopic Orbit for the late-type Be star $β$ CMi

Authors: Nick Dulaney, Noel D. Richardson, Cody J. Gerhartz, J. E. Bjorkman, K. S. Bjorkman, Alex C. Carciofi, Robert Klement, Luqian Wang, Nancy D. Morrison, Allison D. Bratcher, Jennifer J. Greco, Kevin K. Hardegree-Ullman, Ludwik Lembryk, Wayne L. Oswald, Jesica L. Trucks

Abstract: The late-type Be star $β$ CMi is remarkably stable compared to other Be stars that have been studied. This has led to a realistic model of the outflowing Be disk by Klement et al. These results showed that the disk is likely truncated at a finite radius from the star, which Klement et al.~suggest is evidence for an unseen binary companion in orbit. Here we report on an analysis of the Ritter Obser… ▽ More The late-type Be star $β$ CMi is remarkably stable compared to other Be stars that have been studied. This has led to a realistic model of the outflowing Be disk by Klement et al. These results showed that the disk is likely truncated at a finite radius from the star, which Klement et al.~suggest is evidence for an unseen binary companion in orbit. Here we report on an analysis of the Ritter Observatory spectroscopic archive of $β$ CMi to search for evidence of the elusive companion. We detect periodic Doppler shifts in the wings of the H$α$ line with a period of 170 d and an amplitude of 2.25 km s$^{-1}$, consistent with a low-mass binary companion ($M\approx 0.42 M_\odot$). We then compared the small changes in the violet-to-red peak height changes ($V/R$) with the orbital motion. We find weak evidence that it does follow the orbital motion, as suggested by recent Be binary models by Panoglou et al. Our results, which are similar to those for several other Be stars, suggest that $β$ CMi may be a product of binary evolution where Roche lobe overflow has spun up the current Be star, likely leaving a hot subdwarf or white dwarf in orbit around the star. Unfortunately, no direct sign of this companion star is found in the very limited archive of {\it International Ultraviolet Explorer} spectra. △ Less

Submitted 18 January, 2017; originally announced January 2017.

Comments: Accepted to AAS Journals, 2 tables, 3 figures, 11 pages

Pulsations and outbursts in Be stars: Small differences - big impacts

Authors: D. Baade, Th. Rivinius, A. Pigulski, A. Carciofi, G. Handler, R. Kuschnig, Ch. Martayan, A. Mehner, A. F. J. Moffat, H. Pablo, A. Popowicz, S. M. Rucinski, G. A. Wade, W. W. Weiss, K. Zwintz

Abstract: New high-cadence observations with BRITE covering many months confirm that coupled pairs of nonradial pulsation modes are widespread among early-type Be stars. With the difference frequency between the parental variations they may form a roughly sinusoidal variability or the amplitude may cyclicly vary. A first - amplified - beat pattern is also found. In all three cases the amplitudes of differen… ▽ More New high-cadence observations with BRITE covering many months confirm that coupled pairs of nonradial pulsation modes are widespread among early-type Be stars. With the difference frequency between the parental variations they may form a roughly sinusoidal variability or the amplitude may cyclicly vary. A first - amplified - beat pattern is also found. In all three cases the amplitudes of difference frequencies can exceed the amplitude sum of the base frequencies, and modulations of the star-to-circumstellar-disk mass-transfer rate may be associated with these slow variations. This suggests more strongly than any earlier observations that significant dissipation of pulsational energy in the atmosphere may be a cause of mass ejections from Be stars. A unifying interpretative concept is presented. △ Less

Submitted 3 November, 2016; originally announced November 2016.

Comments: 10 pages, 3 figures, preliminary formatting. To appear in Proc. 2nd BRITE-Constellation Science Conference: 'small satellites - big science', Innsbruck, Aug. 2016, eds.: G. Stachowski, E. Poretti, and J. Matthews

Radiative transfer on decretion discs of Be binaries

Authors: Despina Panoglou, Daniel M. Faes, Alex C. Carciofi, Atsuo T. Okazaki, Thomas Rivinius

Abstract: In this work we explore the effect of binarity in the decretion disc of Be stars, in order to explain their variability. To this aim, we performed smoothed particle hydrodynamics (SPH) simulations on Be binary systems, following the matter ejected isotropically from the equator of the Be star towards the base of an isothermal decretion disc. We let the system evolve for time long enough to be cons… ▽ More In this work we explore the effect of binarity in the decretion disc of Be stars, in order to explain their variability. To this aim, we performed smoothed particle hydrodynamics (SPH) simulations on Be binary systems, following the matter ejected isotropically from the equator of the Be star towards the base of an isothermal decretion disc. We let the system evolve for time long enough to be considered at steady state, and focus on the effect of viscosity for coplanar prograde binary orbits. The disc structure is found to be locked to the orbital phase, exhibiting also a dependence on the azimuthal angle. Additionally, we present the first results from detailed non-local thermodynamic equilibrium (non-LTE) radiative transfer calculations of the disc structure computed with the SPH code. This is achieved by the use of the three-dimensional (3D) Monte Carlo code HDUST, which can produce predictions with respect to a series of observables. △ Less

Submitted 28 October, 2016; originally announced October 2016.

Comments: 6 pages, 3 figures; to appear in the proceedings of the conference "The B[e] phenomenon: Forty years of studies" held in Prague, 2016

Be Stars Seen by Space Photometry

Authors: Thomas Rivinius, Dietrich Baade, Alex C. Carciofi

Abstract: Classical Be stars are introduced as object class and their particular potential for space based photometry is highlighted. A brief summary of the various types of variability observed in Be stars makes clear that an interpretation of every single frequency as a pulsation mode falls short, instead there are as well purely circumstellar variations and those that originate in the immediate stellar t… ▽ More Classical Be stars are introduced as object class and their particular potential for space based photometry is highlighted. A brief summary of the various types of variability observed in Be stars makes clear that an interpretation of every single frequency as a pulsation mode falls short, instead there are as well purely circumstellar variations and those that originate in the immediate stellar to circumstellar interaction region. In particular the latter offer great potential, as they are linked to one of the few remaining great riddles of Be stars, namely how they feed their disks. △ Less

Submitted 23 October, 2016; originally announced October 2016.

Comments: 8 pages, 4 figures, invited talk to 2nd BRITE-Constellation Science Conference held in Innsbruck, to be published in the Proceedings of the Polish Astronomical Society, eds: G. Stachowski, E. Poretti, and J. Matthews

Difference frequencies of nonradial pulsation modes and repetitive outbursts in classical Be stars

Authors: D. Baade, Th. Rivinius, A. Pigulski, A. Carciofi, BEST

Abstract: Space photometry has revealed rich pulsation spectra in classical Be stars. Often frequency pairs can be found that differ by nearly the frequency of a third genuine stellar variability. The lowest currently known of these so-called difference (or $Δ$) frequencies reach 0.03 c/d, which is the upper limit of outburst frequencies. The potential of such slow variations for pulsation-assisted mass los… ▽ More Space photometry has revealed rich pulsation spectra in classical Be stars. Often frequency pairs can be found that differ by nearly the frequency of a third genuine stellar variability. The lowest currently known of these so-called difference (or $Δ$) frequencies reach 0.03 c/d, which is the upper limit of outburst frequencies. The potential of such slow variations for pulsation-assisted mass loss is discussed. △ Less

Submitted 7 October, 2016; originally announced October 2016.

Comments: 6 pages, 0 figures, Proc. 'THE B[e] PHENOMENOM: FORTY YEARS OF STUDIES', Prgaue

The Complex Circumstellar and Circumbinary Environment of V356 Sgr

Authors: Jamie R. Lomax, Andrew G. Fullard, Michael A. Malatesta, Brian Babler, Daniel Bednarski, Jodi R. Berdis, Karen S. Bjorkman, Jon E. Bjorkman, Alex C. Carciofi, James W. Davidson Jr., Marcus Keil, Marilyn R. Meade, Kenneth Nordsieck, Matt Scheffler, Jennifer L. Hoffman, John P. Wisniewski

Abstract: We analyze 45 spectropolarimetric observations of the eclipsing, interacting binary star V356 Sgr, obtained over a period of 21 years, to characterize the geometry of the system's circumstellar material. After removing interstellar polarization from these data, we find the system exhibits a large intrinsic polarization signature arising from electron scattering. In addition, the lack of repeatable… ▽ More We analyze 45 spectropolarimetric observations of the eclipsing, interacting binary star V356 Sgr, obtained over a period of 21 years, to characterize the geometry of the system's circumstellar material. After removing interstellar polarization from these data, we find the system exhibits a large intrinsic polarization signature arising from electron scattering. In addition, the lack of repeatable eclipses in the polarization phase curves indicates the presence of a substantial pool of scatterers not occulted by either star. We suggest that these scatterers form either a circumbinary disk coplanar with the gainer's accretion disk or an elongated structure perpendicular to the orbital plane of V356 Sgr, possibly formed by bipolar outflows. We also observe small-scale, cycle-to-cycle variations in the magnitude of intrinsic polarization at individual phases, which we interpret as evidence of variability in the amount of scattering material present within and around the system. This may indicate a mass transfer or mass loss rate that varies on the time-scale of the system's orbital period. Finally, we compare the basic polarimetric properties of V356 Sgr with those of the well studied beta Lyr system; the significant differences observed between the two systems suggests diversity in the basic circumstellar geometry of Roche-lobe overflow systems △ Less

Submitted 23 September, 2016; originally announced September 2016.

Comments: 14 pages, 6 figures, 5 tables, accepted to MNRAS, Tables 1 and 3 will be available in full at CDS after publication

Short-term variability and mass loss in Be stars II. Physical taxonomy of photometric variability observed by the Kepler spacecraft

Authors: Thomas Rivinius, Dietrich Baade, Alex Carciofi

Abstract: Context: [abbreviated] Aims: Kepler data of three known Be stars are re-visited to establish their pulsational nature and assess the properties of additional, non-pulsational variations. The three program stars turned out to be one inactive Be star, one active, continuously outbursting Be star, and one Be star transiting from a non-outbursting into an outbursting phase, thus forming an excellent… ▽ More Context: [abbreviated] Aims: Kepler data of three known Be stars are re-visited to establish their pulsational nature and assess the properties of additional, non-pulsational variations. The three program stars turned out to be one inactive Be star, one active, continuously outbursting Be star, and one Be star transiting from a non-outbursting into an outbursting phase, thus forming an excellent sample to distill properties of Be stars in the various phases of their life-cycle. Methods: [abbreviated] Results: The short-term photometric variability of Be stars must be disentangled into a stellar and a circumstellar part. The stellar part is on the whole not different from what is seen in non-Be stars. However, some of the observed phenomena might be to be due to resonant mode coupling, a mechanism not typically considered for B-type stars. Short-term circumstellar variability comes in the form of either a group of relatively well-defined, short-lived frequencies during outbursts, which are called Štefl frequencies, and broad bumps in the power spectra, indicating aperiodic variability on a time scale similar to typical low-order $g$-mode pulsation frequencies, rather than true periodicity. Conclusions: From a stellar pulsation perspective, Be stars are rapidly rotating SPB stars, that is they pulsate in low order $g$-modes, even if the rapid rotation can project the observed frequencies into the traditional high-order $p$-mode regime above about 4 c/d. However, when a circumstellar disk is present, Be star power spectra are complicated by both cyclic, or periodic, and aperiodic circumstellar phenomena, possibly even dominating the power spectrum. △ Less

Submitted 9 August, 2016; originally announced August 2016.

Comments: Accepted by Astronomy and Astrophysics

Journal ref: A&A 593, A106 (2016)

Investigating the Circumstellar Disk of the Be Shell Star 48 Librae

Authors: J. Silaj, C. E. Jones, A. C. Carciofi, C. Escolano, A. T. Okazaki, C. Tycner, T. Rivinius, R. Klement, D. Bednarski

Abstract: A global disk oscillation implemented in the viscous decretion disk (VDD) model has been used to reproduce most of the observed properties of the well known Be star $ζ$ Tau. 48 Librae shares several similarities with $ζ$ Tau -- they are both early-type Be stars, they display shell characteristics in their spectra, and they exhibit cyclic $V/R$ variations -- but has some marked differences as well,… ▽ More A global disk oscillation implemented in the viscous decretion disk (VDD) model has been used to reproduce most of the observed properties of the well known Be star $ζ$ Tau. 48 Librae shares several similarities with $ζ$ Tau -- they are both early-type Be stars, they display shell characteristics in their spectra, and they exhibit cyclic $V/R$ variations -- but has some marked differences as well, such as a much denser and more extended disk, a much longer $V/R$ cycle, and the absence of the so-called triple-peak features. We aim to reproduce the photometric, polarimetric, and spectroscopic observables of 48 Librae with a self-consistent model, and to test the global oscillation scenario for this target. Our calculations are carried out with the three-dimensional NLTE radiative transfer code HDUST. We employ a rotationally deformed, gravity-darkened central star, surrounded by a disk whose unperturbed state is given by the VDD model. A two-dimensional global oscillation code is then used to calculate the disk perturbation, and superimpose it on the unperturbed disk. A very good, self-consistent fit to the time-averaged properties of the disk is obtained with the VDD. The calculated perturbation has a period $P = 12$ yr, which agrees with the observed period, and the behaviour of the $V/R$ cycle is well reproduced by the perturbed model. The perturbed model improves the fit to the photometric data and reproduces some features of the observed spectroscopic data. Some suggestions to improve the synthesized spectroscopy in a future work are given. △ Less

Submitted 2 August, 2016; originally announced August 2016.

Journal ref: 2016ApJ...826...81S

Be discs in binary systems I. Coplanar orbits

Authors: Despina Panoglou, Alex C. Carciofi, Rodrigo G. Vieira, Isabelle H. Cyr, Carol E. Jones, Atsuo T. Okazaki, Thomas Rivinius

Abstract: Be stars are surrounded by outflowing circumstellar matter structured in the form of decretion discs. They are often members of binary systems, where it is expected that the decretion disc interacts both radiatively and gravitationally with the companion. In this work we study how various orbital (period, mass ratio and eccentricity) and disc (viscosity) parameters affect the disc structure in cop… ▽ More Be stars are surrounded by outflowing circumstellar matter structured in the form of decretion discs. They are often members of binary systems, where it is expected that the decretion disc interacts both radiatively and gravitationally with the companion. In this work we study how various orbital (period, mass ratio and eccentricity) and disc (viscosity) parameters affect the disc structure in coplanar systems. We simulate such binaries with the use of a smoothed particle hydrodynamics code. The main effects of the secondary on the disc are its truncation and the accumulation of material inwards of truncation. We find two cases with respect to the effects of eccentricity: (i) In circular or nearly circular prograde orbits, the disc maintains a rotating, constant in shape, configuration, which is locked to the orbital phase. The disc is smaller in size, more elongated and more massive for low viscosity parameter, small orbital separation and/or high mass ratio. (ii) Highly eccentric orbits are more complex, with the disc structure and total mass strongly dependent on the orbital phase and the distance to the secondary. We also study the effects of binarity in the disc continuum emission. Since the infrared and radio SED are sensitive to the disc size and density slope, the truncation and matter accumulation result in considerable modifications in the emergent spectrum. We conclude that binarity can serve as an explanation for the variability exhibited in observations of Be stars, and that our model can be used to detect invisible companions. △ Less

Submitted 21 June, 2016; v1 submitted 21 May, 2016; originally announced May 2016.

Comments: 15 pages, 30 figures

Non-isothermal effects on Be disks

Authors: Rodrigo G. Vieira, Alex C. Carciofi, Jon E. Bjorkman

Abstract: In the last decade, the viscous decretion disk model has emerged as the new paradigm for Be star disks. In this contribution, we propose a simple analytical model to estimate the continuum infrared excess arising from these circumstellar disks, in the light of the currently accepted scenario. We demonstrate that the disk can be satisfactorily described by a two component system: an inner optically… ▽ More In the last decade, the viscous decretion disk model has emerged as the new paradigm for Be star disks. In this contribution, we propose a simple analytical model to estimate the continuum infrared excess arising from these circumstellar disks, in the light of the currently accepted scenario. We demonstrate that the disk can be satisfactorily described by a two component system: an inner optically thick region, which we call the pseudo-photosphere, and a diffuse outer part. In particular, a direct connexion between the disk brightness profile and the thermal structure is derived, and then confronted to realistic numerical simulations. This result quantifies how the non-isothermality of the disk ultimately affects both infrared measured fluxes and visibilities. △ Less

Submitted 13 April, 2016; originally announced April 2016.

Comments: To appear in conference proceedings "Bright Emissaries: Be Stars as Messengers of Star-Disk Physics", eds. A. Sigut and C. Jones (11 Aug 2014 - 13 Aug 2014, London, Ontario, Canada)

The Variable Line Width of Achernar

Authors: Th. Rivinius, R. H. D. Townsend, D. Baade, A. C. Carciofi, N. Leister, S. Štefl

Abstract: Spectroscopic observations of Achernar over the past decades, have shown the photospheric line width, as measured by the rotational parameter $v \sin i$, to vary in correlation with the emission activity. Here we present new observations, covering the most recent activity phase, and further archival data collected from the archives. The $v \sin i$ variation is confirmed. On the basis of the availa… ▽ More Spectroscopic observations of Achernar over the past decades, have shown the photospheric line width, as measured by the rotational parameter $v \sin i$, to vary in correlation with the emission activity. Here we present new observations, covering the most recent activity phase, and further archival data collected from the archives. The $v \sin i$ variation is confirmed. On the basis of the available data it cannot be decided with certainty whether the increased line width precedes the emission activity, i.e. is a signature of the ejection mechanism, or postdates is, which would make it a signature of re-accretion of some of the disk-material. However, the observed evidence leans towards the re-accretion hypothesis. Two further stars showing the effect of variable line width in correlation with emission activity, namely 66 Oph and $π$ Aqr, are presented as well. △ Less

Submitted 10 February, 2016; originally announced February 2016.

Comments: To appear in conference proceedings "Bright Emissaries: Be Stars as Messengers of Star-Disk Physics", eds. A. Sigut and C. Jones (11 Aug 2014 - 13 Aug 2014, London, Ontario, Canada)