Showing 1–43 of 43 results for author: Jones, C E

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

Angular Momentum Loss Rates in Be Stars Determined by the Viscous Decretion Disc Model

Authors: M. R. Ghoreyshi, C. E. Jones, A. Granada

Abstract: Circumstellar discs around Be stars are formed by the material ejected by the central star. This process removes excess angular momentum from the star as viscosity facilitates the mass and angular momentum transfer within the disc and its growth. The angular momentum loss rates (AMLR) of Be stars is a subject of debate in the literature. Through the modelling of the disc formation and dissipation… ▽ More Circumstellar discs around Be stars are formed by the material ejected by the central star. This process removes excess angular momentum from the star as viscosity facilitates the mass and angular momentum transfer within the disc and its growth. The angular momentum loss rates (AMLR) of Be stars is a subject of debate in the literature. Through the modelling of the disc formation and dissipation phases observed from Be stars, their average AMLR can be determined and this is the goal of this work. We use the viscous decretion disc (VDD) model to provide a range of the average AMLR for Be stars and compare these rates with predicted values from the literature. We explore the reasons for discrepancies between the predicted values of average AMLR using the VDD and Geneva stellar evolution (GSE) models that were previously reported in literature and find that the largest differences occur when Be stars are rotating below their critical speeds. We show that the time over which the mass reservoir builds up is inversely proportional to the average AMLR. Also, we determine a revised value of the average AMLR for the Galactic Be star omega CMa of 4.7x10^36 g cm^2/s^2, which is in better agreement with the values expected for a typical B2 type star. Finally, the effect of disc truncation due to the presence of a companion star is investigated and we find that this has a minimal effect on the average AMLR. △ Less

Submitted 26 October, 2022; originally announced October 2022.

UV Spectropolarimetry with Polstar: Massive Star Binary Colliding Winds

Authors: Nicole St-Louis, Kenneth Gayley, D. John Hillier, Richard Ignace, Carol E. Jones, Alexandre David-Uraz, Noel D. Richardson, Jorick S. Vink, Geraldine J. Peters, Jennifer L. Hoffman, Yael, Nazé, Heloise Stevance, Tomer Shenar, Andrew G. Fullard, Jaimie R. Lomax, Paul A. Scowen

Abstract: The winds of massive stars are important for their direct impact on the interstellar medium, and for their influence on the final state of a star prior to it exploding as a supernova. However, the dynamics of these winds is understood primarily via their illumination from a single central source. The Doppler shift seen in resonance lines is a useful tool for inferring these dynamics, but the mappi… ▽ More The winds of massive stars are important for their direct impact on the interstellar medium, and for their influence on the final state of a star prior to it exploding as a supernova. However, the dynamics of these winds is understood primarily via their illumination from a single central source. The Doppler shift seen in resonance lines is a useful tool for inferring these dynamics, but the mapping from that Doppler shift to the radial distance from the source is ambiguous. Binary systems can reduce this ambiguity by providing a second light source at a known radius in the wind, seen from orbitally modulated directions. From the nature of the collision between the winds, a massive companion also provides unique additional information about wind momentum fluxes. Since massive stars are strong ultraviolet (UV) sources, and UV resonance line opacity in the wind is strong, UV instruments with a high resolution spectroscopic capability are essential for extracting this dynamical information. Polarimetric capability also helps to further resolve ambiguities in aspects of the wind geometry that are not axisymmetric about the line of sight, because of its unique access to scattering direction information. We review how the proposed MIDEX-scale mission Polstar can use UV spectropolarimetric observations to critically constrain the physics of colliding winds, and hence radiatively-driven winds in general. We propose a sample of 20 binary targets, capitalizing on this unique combination of illumination by companion starlight, and collision with a companion wind, to probe wind attributes over a range in wind strengths. Of particular interest is the hypothesis that the radial distribution of the wind acceleration is altered significantly, when the radiative transfer within the winds becomes optically thick to resonance scattering in multiple overlapping UV lines. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: 26 pages, 12 figures, Review in a topical collection series of Astrophysics and Space Sciences on the proposed Polstar satellite. arXiv admin note: substantial text overlap with arXiv:2111.11552

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

Ultraviolet Spectropolarimetry:Conservative and Nonconservative Mass Transfer in OB Interacting Binaries

Authors: Geraldine J. Peters, Ken Gayley, Richard Ignace, Carol E. Jones, Yael Naze, Nicole St-Louis, Heloise Stevance, Jorick S. Vink, Noel D. Richardson, Jennifer L. Hoffman, Jamie R. Lomax, Tomer Shenar, Andrew G. Fullard, Paul A. Scowen

Abstract: One objective of the Polstar spectropolarimetry mission is to characterize the degree of nonconservative mass transfer that occurs at various stages of binary evolution, from the initial mass reversal to the late Algol phase. The proposed instrument combines spectroscopic and polarimetric capabilities, where the spectroscopy can resolve Doppler shifts in UV resonance lines with 10 km/s precision,… ▽ More One objective of the Polstar spectropolarimetry mission is to characterize the degree of nonconservative mass transfer that occurs at various stages of binary evolution, from the initial mass reversal to the late Algol phase. The proposed instrument combines spectroscopic and polarimetric capabilities, where the spectroscopy can resolve Doppler shifts in UV resonance lines with 10 km/s precision, and polarimetry can resolve linear polarization with 1e-3 precision or better. The spectroscopy will identify absorption by mass streams seen in projection against the stellar disk as a function of orbital phase, hot accretion spots, as well as scattering from extended splash structures, circumbinary disks, and other flows in and above/below the orbital plane (e.g. jets) that fail to be transferred conservatively. The polarimetry affects more the light coming from material not seen against the stellar disk, allowing the geometry of the scattering to be tracked, resolving ambiguities left by the spectroscopy and light-curve information. For example, nonconservative mass streams ejected in the polar direction will produce polarization of the opposite sign from conservative transfer accreting in the orbital plane. Also, time domain coverage over a range of phases of the binary orbit are well supported by the Polstar observing strategy. Combining these elements will significantly improve our understanding of the mass transfer process and the amount of mass that can escape from the system, an important channel for changing the final mass, and ultimate supernova, of the large number of massive stars found in binaries at close enough separation to undergo interaction. △ Less

Submitted 12 October, 2022; v1 submitted 27 November, 2021; originally announced November 2021.

Comments: Refereed paper in Ap&SS 367:Topical Collection, 23 pages, 10 figures

Ultraviolet Spectropolarimetry with Polstar: Massive Star Binary Colliding Winds

Authors: Nicole St-Louis, Kenneth Gayley, Desmond John Hillier, Richard Ignace, Carol E. Jones, Alexandre David-Uraz, Noel D. Richardson, Jorick S. Vink, Geraldine J. Peters, Jennifer L. Hoffman, Yael Naze, Heloise Stevance, Tomer Shenar, Andrew G. Fullard, Jamie R. Lomax, Paul A. Scowen

Abstract: As sources of chemical enrichment, ionizing radiation and energetic feedback, massive stars drive the ecology of their host galaxies despite their relative rarity, additionally to yielding compact remnants, which can generate gravitational waves. The evolution of massive stars is crucially informed by their detailed mass-loss history; however, wind structures on a variety of scales cause important… ▽ More As sources of chemical enrichment, ionizing radiation and energetic feedback, massive stars drive the ecology of their host galaxies despite their relative rarity, additionally to yielding compact remnants, which can generate gravitational waves. The evolution of massive stars is crucially informed by their detailed mass-loss history; however, wind structures on a variety of scales cause important uncertainties on their mass-loss rates. Binary systems can place further constraints on the mass-loss properties of massive stars, especially colliding-wind binaries. In this paper, we review how the proposed MIDEX-scale mission Polstar can critically constrain the physics of colliding winds (and hence radiatively-driven winds in general) with ultraviolet spectropolarimetric observations, providing an unprecedented improvement on the accuracy of the determination of both mass-loss rates and the velocity structure of the winds of massive stars. We propose a sample of 17 targets that will allow us to study a variety of wind-colliding systems spanning a large parameter space using the spatial information yielded by both spectroscopic and polarimetric data obtained with Polstar. △ Less

Submitted 16 December, 2021; v1 submitted 22 November, 2021; originally announced November 2021.

Comments: 21 pages, 12 figures 1 table, Whitepaper

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

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

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

Quiescent and active phases in Be stars: A WISE snapshot of young galactic open clusters

Authors: A. Granada, C. E. Jones, T. A. A. Sigut, T. Semaan, C. Georgy, G. Meynet, S. Ekström

Abstract: Through the modelling of near-infrared photometry of star plus disk systems with the codes BEDISK/BERAY, we successfully describe the WISE photometric characteristics of Be stars in five young open clusters, NGC 663, NGC 869, NGC 884, NGC 3766 and NGC 4755, broadly studied in the literature. WISE photometry allows previously known Be stars to be detected and to find new Be candidates which could b… ▽ More Through the modelling of near-infrared photometry of star plus disk systems with the codes BEDISK/BERAY, we successfully describe the WISE photometric characteristics of Be stars in five young open clusters, NGC 663, NGC 869, NGC 884, NGC 3766 and NGC 4755, broadly studied in the literature. WISE photometry allows previously known Be stars to be detected and to find new Be candidates which could be confirmed spectroscopically. The location of Be stars in the WISE colour-magnitude diagram, separates them in two groups; active (Be stars hosting a developed circumstellar disk) and quiescent objects (Be stars in a diskless phase), and this way, we can explore how often stars are observed in these different stages. The variability observed in most active variable Be stars is compatible with a disk dissipation phase. We find that 50 percent of Be stars in the studied open clusters are in an active phase. We can interpret this as Be stars having a developed circumstellar disk half of the time. The location of Be stars with a developed disk in the CMD require mass loss rates in agreement with values recently reported in the literature. For these objects, we expect to have a tight relation between the equivalent width of the Hα line and the mass of the disk, if the inclination is known. Also, near-infrared photometry of Be stars in stellar clusters has the potential of being useful to test whether there is a preferential viewing angle. △ Less

Submitted 14 December, 2017; originally announced December 2017.

Comments: 25 pages, 14 figures, Accepted for publication in AJ

Computing the Photometric and Polarimetric Variability of Be Stars

Authors: Keegan C. Marr, C. E. Jones, R. J. Halonen

Abstract: We investigate variations in the linear polarization as well as in the V-band and B-band colour-magnitudes for classical Be star disks. We present two models: disks with enhanced disk density and disks that are tilted or warped from the stellar equatorial plane. In both cases, we predict variation in observable properties of the system as the disk rotates. We use a non-LTE radiative transfer code… ▽ More We investigate variations in the linear polarization as well as in the V-band and B-band colour-magnitudes for classical Be star disks. We present two models: disks with enhanced disk density and disks that are tilted or warped from the stellar equatorial plane. In both cases, we predict variation in observable properties of the system as the disk rotates. We use a non-LTE radiative transfer code BEDISK (Sigut & Jones) in combination with a Monte Carlo routine that includes multiple scattering (Halonen et al.) to model classical Be star systems. We find that a disk with an enhanced density region that is one order of magnitude denser than the disk's base density shows as much as $\sim$0.2% variability in the polarization while the polarization position angle varies by $\sim$8°. The $Δ$V magnitude for the same system shows variations of up to $\sim$0.4 magnitude while the $Δ$(B-V) colour varies by at most $\sim$0.01 magnitude. We find that disks tilted from the equatorial plane at small angles of $\sim$30°more strongly reflect the values of polarization and colour-magnitudes reported in the literature than disks tilted at larger angles. For this model, the linear polarization varies by $\sim$0.3%, the polarization position angle varies by $\sim$60°, the $Δ$V magnitude varies up to 0.35 magnitude, and the $Δ$(B-V) colour varies up to 0.1 magnitude. We find that the enhanced disk density models show ranges of polarization and colour-magnitudes that are commensurate with what is reported in the literature for all sizes of the density enhanced regions. From this, we cannot determine any preference for small or large density enhanced regions. △ Less

Submitted 7 December, 2017; originally announced December 2017.

Comments: Accepted in ApJ. 11 pages, 9 figures. The University of Western Ontario

$Ω$-Slow Solutions and Be Star Disks

Authors: I. Araya, C. E. Jones, M. Curé, J. Silaj, L. Cidale, A. Granada, A. Jiménez

Abstract: As the disk formation mechanism(s) in Be stars is(are) as yet unknown, we investigate the role of rapidly rotating radiation-driven winds in this process. We implemented the effects of high stellar rotation on m-CAK models accounting for: the shape of the star, the oblate finite disk correction factor, and gravity darkening. For a fast rotating star, we obtain a two-component wind model, i.e., a f… ▽ More As the disk formation mechanism(s) in Be stars is(are) as yet unknown, we investigate the role of rapidly rotating radiation-driven winds in this process. We implemented the effects of high stellar rotation on m-CAK models accounting for: the shape of the star, the oblate finite disk correction factor, and gravity darkening. For a fast rotating star, we obtain a two-component wind model, i.e., a fast, thin wind in the polar latitudes and an $Ω$-slow, dense wind in the equatorial regions. We use the equatorial mass densities to explore H$α$ emission profiles for the following scenarios: 1) a spherically symmetric star, 2) an oblate shaped star with constant temperature, and 3) an oblate star with gravity darkening. One result of this work is that we have developed a novel method for solving the gravity darkened, oblated m-CAK equation of motion. Furthermore, from our modeling we find a) the oblate finite disk correction factor, for the scenario considering the gravity darkening, can vary by at least a factor of two between the equatorial and polar directions, influencing the velocity profile and mass-loss rate accordingly, b) the H$α$ profiles predicted by our model are in agreement with those predicted by a standard power-law model for following values of the line-force parameters: $1.5 \lesssim k \lesssim 3$, $ \, α\sim 0.6$ and $\, δ\gtrsim 0.1$, and c) the contribution of the fast wind component to the H$α$ emission line profile is negligible; therefore, the line profiles arise mainly from the equatorial disks of Be stars. △ Less

Submitted 1 August, 2017; originally announced August 2017.

Comments: 20 pages and 14 figures. Accepted for publication in The Astrophysical Journal

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

The disk physical conditions of 48 Persei constrained by contemporaneous Hα spectroscopy and interferometry

Authors: C. E. Jones, T. A. A. Sigut, B. J. Grzenia, C. Tycner, R. T. Zavala

Abstract: We utilize a multi-step modelling process to produce synthetic interferometric and spectroscopic observables, which are then compared to their observed counterparts. Our extensive set of interferometric observations of the Be star 48 Per, totaling 291 data points, were obtained at the Navy Precision Optical Interferometer from 2006 Nov 07 to 2006 Nov 23. Our models were further constrained by comp… ▽ More We utilize a multi-step modelling process to produce synthetic interferometric and spectroscopic observables, which are then compared to their observed counterparts. Our extensive set of interferometric observations of the Be star 48 Per, totaling 291 data points, were obtained at the Navy Precision Optical Interferometer from 2006 Nov 07 to 2006 Nov 23. Our models were further constrained by comparison with contemporaneous H$α$ line spectroscopy obtained at the John S. Hall Telescope at Lowell Observatory on 2006 Nov 1. Theoretical spectral energy distributions, SEDs, for 48~Per were confirmed by comparison with observations over a wavelength regime of 0.4 to 60 microns from Touhami et al. (2010) and Vieira et al. (2017). Our best-fitting combined model from H$α$ spectroscopy, H$α$ interferometry and SED fitting has a power law density fall off, $n$, of 2.3 and an initial density at the stellar surface of $ρ_0 = 1.0 \times 10^{-11} \rm{g\,cm}^{-3}$ with a inclination constrained by H$α$ spectroscopy and interferometry of 45$^{\circ} \pm 5^{\circ}$. The position angle for the system, measured east from north, is 121$^{\circ}$ $\pm$ 1$^{\circ}$. Our best-fit model shows that the disk emission originates in a moderately large disk with a radius of 25 R$_{*}$ which is consistent with a disk mass of approximately 5 $\times$ 10$^{24}$ g or 3 $\times$ 10$^{-10}$ M$_{*}$. Finally, we compare our results with previous studies of 48~Per by Quirrenbach et al. (1997) with agreement but find, with our much larger data set, that our disk size contradicts the findings of Delaa et al. (2011). △ Less

Submitted 27 April, 2017; originally announced April 2017.

Comments: 26 pages, 9 figures

Evidence for Different Disk Mass Distributions Between Early and Late-Type Be Stars in the BeSOS Survey

Authors: C. Arcos, C. E. Jones, T. A. A. Sigut, S. Kanaan, M. Curé

Abstract: The circumstellar disk density distributions for a sample of 63 Be southern stars from the BeSOS survey were found by modelling their H$α$ emission line profiles. These disk densities were used to compute disk masses and disk angular momenta for the sample. Average values for the disk mass are 3.4$\times$10$^{-9}$ and 9.5$\times$10$^{-10}$ $M_{\star}$ for early (B0-B3) and late (B4-B9) spectral ty… ▽ More The circumstellar disk density distributions for a sample of 63 Be southern stars from the BeSOS survey were found by modelling their H$α$ emission line profiles. These disk densities were used to compute disk masses and disk angular momenta for the sample. Average values for the disk mass are 3.4$\times$10$^{-9}$ and 9.5$\times$10$^{-10}$ $M_{\star}$ for early (B0-B3) and late (B4-B9) spectral types, respectively. We also find that the range of disk angular momentum relative to the star are between 150-200 and 100-150 $J_{\star}/M_{\star}$, again for early and late-type Be stars respectively. The distributions of the disk mass and disk angular momentum are different between early and late-type Be stars at a 1\% level of significance. Finally, we construct the disk mass distribution for the BeSOS sample as a function of spectral type and compare it to the predictions of stellar evolutionary models with rapid rotation. The observed disk masses are typically larger than the theoretical predictions, although the observed spread in disk masses is typically large. △ Less

Submitted 26 April, 2017; originally announced April 2017.

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

Statistical Analysis of Interferometric Measurements of Axis Ratios for Classical Be Stars

Authors: Richard P. Cyr, Carol E. Jones, Christopher Tycner

Abstract: This work presents a novel method to estimate the effective opening angle of CBe star disks from projected axis ratio measurements, obtained by interferometry using Bayesian statistics. A Monte Carlo scheme was used to generate a large set of theoretical axis ratios from disk models using different distributions of disk densities and opening angles. These theoretical samples were then compared to… ▽ More This work presents a novel method to estimate the effective opening angle of CBe star disks from projected axis ratio measurements, obtained by interferometry using Bayesian statistics. A Monte Carlo scheme was used to generate a large set of theoretical axis ratios from disk models using different distributions of disk densities and opening angles. These theoretical samples were then compared to observational samples, using a two-sample Kolmogorov-Smirnov test, to determine which theoretical distribution best reproduces the observations. The results suggest that the observed ratio distributions in the K-, H-, and N-band can best be explained by the presence of thin disks, with opening half-angles of the order of 0.15$^{\circ}$ to 4.0$^{\circ}$. Results for measurements over the H$α$ line point toward slightly thicker disks, 3.7$^{\circ}$ to 14$^{\circ}$, which is consistent with a flaring disk predicted by the viscous disk model. △ Less

Submitted 27 November, 2014; originally announced November 2014.

Comments: Accepted for publication in ApJ, 10 pages, 13 figures

Line-Driven Winds Revisited in the Context of Be Stars: $Ω$-slow Solutions with High $k$ Values

Authors: J. Silaj, M. Cure, C. E. Jones

Abstract: The standard, or fast, solutions of m-CAK line-driven wind theory cannot account for slowly outflowing disks like the ones that surround Be stars. It has been previously shown that there exists another family of solutions --- the $Ω$-slow solutions --- that is characterized by much slower terminal velocities and higher mass-loss rates. We have solved the one-dimensional m-CAK hydrodynamical equati… ▽ More The standard, or fast, solutions of m-CAK line-driven wind theory cannot account for slowly outflowing disks like the ones that surround Be stars. It has been previously shown that there exists another family of solutions --- the $Ω$-slow solutions --- that is characterized by much slower terminal velocities and higher mass-loss rates. We have solved the one-dimensional m-CAK hydrodynamical equation of rotating radiation-driven winds for this latter solution, starting from standard values of the line force parameters ($α$, $k$, and $δ$), and then systematically varying the values of $α$ and $k$. Terminal velocities and mass-loss rates that are in good agreement with those found in Be stars are obtained from the solutions with lower $α$ and higher $k$ values. Furthermore, the equatorial densities of such solutions are comparable to those that are typically assumed in ad hoc models. For very high values of $k$, we find that the wind solutions exhibit a new kind of behavior. △ Less

Submitted 5 November, 2014; originally announced November 2014.

Journal ref: ApJ, 795, 78, 2014

Using Photometry to Probe the Circumstellar Environment of delta Scorpii

Authors: C. E. Jones, Paul Wiegert, C. Tycner, G. W. Henry, R. P. Cyr, R. J. Halonen, M. W. Muterspaugh

Abstract: We acquired Johnson BV photometry of the binary Be disk system delta~Scorpii during its 2009, 2010, 2011, and 2012 observing seasons and used it to probe the innermost regions of the disk. We found that several disk building events have occurred during this time, resulting in an overall brightening in the V-band and reddening of the system. In addition to these long-term trends, we found cyclical… ▽ More We acquired Johnson BV photometry of the binary Be disk system delta~Scorpii during its 2009, 2010, 2011, and 2012 observing seasons and used it to probe the innermost regions of the disk. We found that several disk building events have occurred during this time, resulting in an overall brightening in the V-band and reddening of the system. In addition to these long-term trends, we found cyclical variability in each observing season on timescales between 60 and 100 days. We were able to reproduce the changes in the magnitude and colour of delta Sco using our theoretical models and found that variable mass-loss rates in the range 2.5-7.0x10^{-9} M_{sun}/yr over ~35 days can reproduce the observed increase in brightness. △ Less

Submitted 4 June, 2014; originally announced June 2014.

Comments: 15 pages, emulateapj

Journal ref: AJ, 145, 142, 2013

On the Intrinsic Continuum Linear Polarization of Classical Be Stars: The Effects of Metallicity and One-Armed Density Perturbations

Authors: Robbie J. Halonen, Carol E. Jones

Abstract: We report on the effects of two disk properties on the continuum linear polarization signature of classical Be stars. First, we investigate the effect of including metallicity in computing the thermal structure of the circumstellar gas on the resulting polarimetric Balmer jump. The Balmer jump is a distinguishing feature of the polarization signature in these objects and, as such, can be used as a… ▽ More We report on the effects of two disk properties on the continuum linear polarization signature of classical Be stars. First, we investigate the effect of including metallicity in computing the thermal structure of the circumstellar gas on the resulting polarimetric Balmer jump. The Balmer jump is a distinguishing feature of the polarization signature in these objects and, as such, can be used as a tool for differentiating classical Be stars from similar H alpha-emitters identified through conventional photometric techniques. We find that although low-metallicity environments will have hotter disk temperatures on average, the temperature change alone cannot account for the discrepancy in the frequency of Balmer jumps between low-metallicity and solar-metallicity stellar populations. Second, we investigate the effect of including a global one-armed oscillation in the gas density distribution of the modeled disk. We find that a non-axisymmetric perturbation pattern yields discernible variations in the predicted polarization level. If these density oscillations are present in the inner region of classical Be star disks, the polarimetric variations should produce a periodic signature which can help characterize the dynamical nature of the gas near the star. △ Less

Submitted 23 July, 2013; originally announced July 2013.

Comments: Accepted to ApJS on July 22, 2013

On the Intrinsic Continuum Linear Polarization of Classical Be Stars during Disk Growth and Dissipation

Authors: Robbie J. Halonen, Carol E. Jones

Abstract: We investigate the intrinsic continuum linear polarization from axisymmetric density distributions of gas surrounding classical Be stars during the formation and dissipation of their circumstellar disks. We implement a Monte Carlo calculation of the Stokes parameters with the use of the non-LTE radiative transfer code of Sigut & Jones (2007) to reproduce the continuous polarimetric spectra of clas… ▽ More We investigate the intrinsic continuum linear polarization from axisymmetric density distributions of gas surrounding classical Be stars during the formation and dissipation of their circumstellar disks. We implement a Monte Carlo calculation of the Stokes parameters with the use of the non-LTE radiative transfer code of Sigut & Jones (2007) to reproduce the continuous polarimetric spectra of classical Be stars. The scattering of light in the nonspherical circumstellar envelopes of classical Be stars produces a distinct polarization signature that can be used to study the physical nature of the scattering environment. In this paper, we highlight the utility of polarimetric measurements as important diagnostics in the modeling of these systems. We illustrate the effects of using self-consistent calculation of the thermal structure of the circumstellar gas on the characteristic wavelength-dependence of the polarization spectrum. In showing that the principal features of the polarization spectrum originate from different parts of the disk, we emphasize the capability of polarimetric observations to trace the evolution of the disk on critical scales. We produce models that approximate the disk formation and dissipation periods and illustrate how the polarimetric properties of these systems can have a pivotal role in determining the mechanism for mass decretion from the central star. △ Less

Submitted 23 July, 2013; originally announced July 2013.

Comments: Accepted to ApJ on January 15, 2013

Journal ref: The Astrophysical Journal, 765:17, 2013 March 1

The Effect of Density on the Thermal Structure of Gravitationally-Darkened Be Star Disks

Authors: M. A. McGill, T. A. A. Sigut, C. E. Jones

Abstract: The effects of gravitational darkening on the thermal structure of Be star disks of differing densities are systematically examined. Gravitational darkening is the decrease of the effective temperature near the equator and the corresponding increase near the poles of a star caused by rapid rotation. We also include the rotational distortion of the star using the Roche Model. Increasing the disk de… ▽ More The effects of gravitational darkening on the thermal structure of Be star disks of differing densities are systematically examined. Gravitational darkening is the decrease of the effective temperature near the equator and the corresponding increase near the poles of a star caused by rapid rotation. We also include the rotational distortion of the star using the Roche Model. Increasing the disk density increases the optical depths in the equatorial plane, resulting in the formation of an inner cool region near the equatorial plane of the disk. High rotation rates result in disks that have temperatures similar to those of a denser disk, namely cooler overall. However the effect of increasing rotation produces additional heating in the upper disk due to the hotter stellar pole. Cool regions in the equatorial plane normally associated with high density are seen in low density models at high rotation rates. Gravitational darkening increases the amount of very cool and very hot material in the disk and decreases the amount of disk material at moderate temperatures. We also present models which study the effect of gravitational darkening on hydrostatically-converged disks, in which the temperature structure is consistent with vertical hydrostatic equilibrium. Because the equatorial regions become cooler, hydrostatically converged models that include gravity darkening have smaller vertical scale heights, and $H/R$ is smaller by as much as 56% near $v_{\rm crit}$. Finally we explore differences in disk temperatures when alternate formulations of gravitational darkening, which lower the temperature difference between the pole and the equator, are used. △ Less

Submitted 29 October, 2012; originally announced October 2012.

The Thermal Structure of Gravitationally-Darkened Classical Be Star Disks

Authors: M. A. McGill, T. A. A. Sigut, C. E. Jones

Abstract: The effect of gravitational darkening on models of the thermal structure of Be star disks is systematically studied for a wide range of Be star spectral types and rotation rates. Gravitational darkening causes a reduction of the stellar effective temperature towards the equator and a redirection of energy towards the poles. It is an important physical effect in these star-disk systems because the… ▽ More The effect of gravitational darkening on models of the thermal structure of Be star disks is systematically studied for a wide range of Be star spectral types and rotation rates. Gravitational darkening causes a reduction of the stellar effective temperature towards the equator and a redirection of energy towards the poles. It is an important physical effect in these star-disk systems because the photoionizing radiation from the central B star is the main energy source for the disk. We have added gravitational darkening to the bedisk code to produce circumstellar disk models that include both the variation in the effective temperature with latitude and the non-spherical shape of the star in the calculation of the stellar photoionizing radiation field. The effect of gravitational darkening on global measures of disk temperature is generally significant for rotation rates above 80% of critical rotation. For example, a B0V model rotating at 95% of critical has a density-averaged disk temperature \approx 2500 K cooler than its nonrotating counterpart. However, detailed differences in the temperature structure of disks surrounding rotating and non-rotating stars reveal a complex pattern of heating and cooling. Spherical gravitational darkening, an approximation that ignores the changing shape of the star, gives good results for disk temperatures for rotation rates less than \approx 80% of critical. However for the highest rotation rates, the distortion of the stellar surface caused by rotation becomes important. △ Less

Submitted 30 August, 2011; v1 submitted 29 August, 2011; originally announced August 2011.

Comments: 35 pages, 19 figures

The Variability of Halpha Equivalent Widths in Be Stars

Authors: C. E. Jones, C. Tycner, A. D. Smith

Abstract: Focusing on B-emission stars, we investigated a set of H$α$ equivalent widths calculated from observed spectra acquired over a period of about 4 years from 2003 to 2007. During this time, changes in equivalent width for our program stars were monitored. We have found a simple statistical method to quantify these changes in our observations. This statistical test, commonly called the F ratio, invol… ▽ More Focusing on B-emission stars, we investigated a set of H$α$ equivalent widths calculated from observed spectra acquired over a period of about 4 years from 2003 to 2007. During this time, changes in equivalent width for our program stars were monitored. We have found a simple statistical method to quantify these changes in our observations. This statistical test, commonly called the F ratio, involves calculating the ratio of the external and internal error. We show that the application of this technique can be used to place bounds on the degree of variability of Be stars. This observational tool provides a quantitative way to find Be stars at particular stages of variability requiring relatively little observational data. △ Less

Submitted 8 March, 2011; originally announced March 2011.

Comments: 8 figures accepted by The Astronomical Journal

Be Star Disk Models in Consistent Vertical Hydrostatic Equilibrium

Authors: T. A. A. Sigut, M. A. McGill, C. E. Jones

Abstract: A popular model for the circumstellar disks of Be stars is that of a geometrically thin disk with a density in the equatorial plane that drops as a power law of distance from the star. It is usually assumed that the vertical structure of such a disk (in the direction parallel to the stellar rotation axis) is governed by the hydrostatic equilibrium set by the vertical component of the star's grav… ▽ More A popular model for the circumstellar disks of Be stars is that of a geometrically thin disk with a density in the equatorial plane that drops as a power law of distance from the star. It is usually assumed that the vertical structure of such a disk (in the direction parallel to the stellar rotation axis) is governed by the hydrostatic equilibrium set by the vertical component of the star's gravitational acceleration. Previous radiative equilibrium models for such disks have usually been computed assuming a fixed density structure. This introduces an inconsistency as the gas density is not allowed to respond to temperature changes and the resultant disk model is not in vertical, hydrostatic equilibrium. In this work, we modify the {\sc bedisk} code of \citet{sig07} so that it enforces a hydrostatic equilibrium consistent with the temperature solution. We compare the disk densities, temperatures, H$α$ line profiles, and near-IR excesses predicted by such models with those computed from models with a fixed density structure. We find that the fixed models can differ substantially from the consistent hydrostatic models when the disk density is high enough that the circumstellar disk develops a cool ($T\lesssim10,000 $K) equatorial region close to the parent star. Based on these new hydrostatic disks, we also predict an approximate relation between the (global) density-averaged disk temperature and the $T_{\rm eff}$ of the central star, covering the full range of central Be star spectral types. △ Less

Submitted 8 May, 2009; originally announced May 2009.

Comments: 25 pages; 11 figures

Journal ref: Astrophys.J.699:1973-1981,2009

The case for optical interferometric polarimetry

Authors: Nicholas M. Elias II, Carol E. Jones, Henrique R. Schmitt, Anders M. Jorgensen, Michael J. Ireland, Karine Perraut

Abstract: Within the last 10 years, long-baseline optical interferometry (LBOI) has benefited significantly from increased sensitivity, spatial resolution, and spectral resolution, e.g., measuring the diameters and asymmetries of single stars, imaging/fitting the orbits of multiple stars, modeling Be star disks, and modeling AGN nuclei. Similarly, polarimetry has also yielded excellent astrophysical resul… ▽ More Within the last 10 years, long-baseline optical interferometry (LBOI) has benefited significantly from increased sensitivity, spatial resolution, and spectral resolution, e.g., measuring the diameters and asymmetries of single stars, imaging/fitting the orbits of multiple stars, modeling Be star disks, and modeling AGN nuclei. Similarly, polarimetry has also yielded excellent astrophysical results, e.g., characterizing the atmospheres and shells of red giants/supergiants, modeling the envelopes of AGB stars, studying the morphology of Be stars, and monitoring the short- and long- term behavior of AGNs. The next logical evolutionary step in instrumentation is to combine LBOI with polarimetry, which is called optical interferometric polarimetry (OIP). In other words, measurements of spatial coherence are performed simultaneously with measurements of coherence between orthogonal polarization states. △ Less

Submitted 19 November, 2008; originally announced November 2008.

Comments: This is a white paper

Constraining the Physical Parameters of the Circumstellar Disk of chi Ophiuchi

Authors: C. Tycner, C. E. Jones, T. A. A. Sigut, H. R. Schmitt, J. A. Benson, D. J. Hutter, R. T. Zavala

Abstract: We present a numerical model describing a circularly symmetric gaseous disk around the Be star chi Ophiuchi. The model is constrained by long-baseline interferometric observations that are sensitive to the H-alpha Balmer line emission from the disk. For the first time our interferometric observations spatially resolve the inner region of the circumstellar disk around chi Ophiuchi and we use thes… ▽ More We present a numerical model describing a circularly symmetric gaseous disk around the Be star chi Ophiuchi. The model is constrained by long-baseline interferometric observations that are sensitive to the H-alpha Balmer line emission from the disk. For the first time our interferometric observations spatially resolve the inner region of the circumstellar disk around chi Ophiuchi and we use these results to place a constraint on the physical extent of the H-alpha-emitting region. We demonstrate how this in turn results in very specific constraints on the parameters that describe the variation of the gas density as a function of radial distance from the central star. △ Less

Submitted 24 July, 2008; originally announced July 2008.

Comments: Accepted for publication in ApJ

A Parameter Study of Classical Be Star Disk Models Constrained by Optical Interferometry

Authors: C. E. Jones, C. Tycner, T. A. A. Sigut, J. A. Benson, D. J. Hutter

Abstract: We have computed theoretical models of circumstellar disks for the classical Be stars $κ$ Dra, $β$ Psc, and $\upsilon$ Cyg. Models were constructed using a non-LTE radiative transfer code developed by \citet{sig07} which incorporates a number of improvements over previous treatments of the disk thermal structure, including a realistic chemical composition. Our models are constrained by direct co… ▽ More We have computed theoretical models of circumstellar disks for the classical Be stars $κ$ Dra, $β$ Psc, and $\upsilon$ Cyg. Models were constructed using a non-LTE radiative transfer code developed by \citet{sig07} which incorporates a number of improvements over previous treatments of the disk thermal structure, including a realistic chemical composition. Our models are constrained by direct comparison with long baseline optical interferometric observations of the H$α$ emitting regions and by contemporaneous H$α$ line profiles. Detailed comparisons of our predictions with H$α$ interferometry and spectroscopy place very tight constraints on the density distributions for these circumstellar disks. △ Less

Submitted 9 July, 2008; originally announced July 2008.

Comments: 10 figures,28 pages, accepted by ApJ

Infrared Excess in the Be Star Delta Scorpii

Authors: R. J. Halonen, C. E. Jones, T. A. A. Sigut, R. T. Zavala, C. Tycner, S. E. Levine, C. B. Luginbuhl, F. J. Vrba, N. Vlieg

Abstract: We present infrared photometric observations of the Be binary system delta Scorpii obtained in 2006. The J,H and K magnitudes are the same within the errors compared to observations taken 10 months earlier. We derive the infrared excess from the observation and compare this to the color excess predicted by a radiative equilibrium model of the primary star and its circumstellar disk. We use a non… ▽ More We present infrared photometric observations of the Be binary system delta Scorpii obtained in 2006. The J,H and K magnitudes are the same within the errors compared to observations taken 10 months earlier. We derive the infrared excess from the observation and compare this to the color excess predicted by a radiative equilibrium model of the primary star and its circumstellar disk. We use a non-LTE computational code to model the gaseous envelope concentrated in the star's equatorial plane and calculate the expected spectral energy distribution and Halpha emission profile of the star with its circumstellar disk. Using the observed infrared excess of delta Sco, as well as Halpha spectroscopy bracketing the IR observations in time, we place constraints on the radial density distribution in the circumstellar disk. Because the disk exhibits variability in its density distribution, this work will be helpful in understanding its dynamics. △ Less

Submitted 7 May, 2008; originally announced May 2008.

Comments: 12 pages, 14 figures, to be published in PASP May 2008

The Thermal Structure of the Circumstellar Disk Surrounding the Classical Be Star gamma Cassiopeia

Authors: T. A. A. Sigut, C. E. Jones

Abstract: We have computed radiative equilibrium models for the gas in the circumstellar envelope surrounding the hot, classical Be star $γ$Cassiopeia. This calculation is performed using a code that incorporates a number of improvements over previous treatments of the disk's thermal structure by \citet{mil98} and \citet{jon04}; most importantly, heating and cooling rates are computed with atomic models f… ▽ More We have computed radiative equilibrium models for the gas in the circumstellar envelope surrounding the hot, classical Be star $γ$Cassiopeia. This calculation is performed using a code that incorporates a number of improvements over previous treatments of the disk's thermal structure by \citet{mil98} and \citet{jon04}; most importantly, heating and cooling rates are computed with atomic models for H, He, CNO, Mg, Si, Ca, & Fe and their relevant ions. Thus, for the first time, the thermal structure of a Be disk is computed for a gas with a solar chemical composition as opposed to assuming a pure hydrogen envelope. We compare the predicted average disk temperature, the total energy loss in H$α$, and the near-IR excess with observations and find that all can be accounted for by a disk that is in vertical hydrostatic equilibrium with a density in the equatorial plane of $ρ(R)\approx 3$ to $5\cdot 10^{-11} (R/R_*)^{-2.5} \rm g cm^{-3}$. We also discuss the changes in the disk's thermal structure that result from the additional heating and cooling processes available to a gas with a solar chemical composition over those available to a pure hydrogen plasma. △ Less

Submitted 27 June, 2007; originally announced June 2007.

Comments: 11 pages, 8 figures high resolution figures available at http://inverse.astro.uwo.ca/sig_jon07.html

Constraining Disk Parameters of Be Stars using Narrowband H-alpha Interferometry with the NPOI

Authors: C. Tycner, G. C. Gilbreath, R. T. Zavala, J. T. Armstrong, J. A. Benson, A. R. Hajian, D. J. Hutter, C. E. Jones, T. A. Pauls, N. M. White

Abstract: Interferometric observations of two well-known Be stars, gamma Cas and phi Per, were collected and analyzed to determine the spatial characteristics of their circumstellar regions. The observations were obtained using the Navy Prototype Optical Interferometer equipped with custom-made narrowband filters. The filters isolate the H-alpha emission line from the nearby continuum radiation, which res… ▽ More Interferometric observations of two well-known Be stars, gamma Cas and phi Per, were collected and analyzed to determine the spatial characteristics of their circumstellar regions. The observations were obtained using the Navy Prototype Optical Interferometer equipped with custom-made narrowband filters. The filters isolate the H-alpha emission line from the nearby continuum radiation, which results in an increased contrast between the interferometric signature due to the H-alpha-emitting circumstellar region and the central star. Because the narrowband filters do not significantly attenuate the continuum radiation at wavelengths 50 nm or more away from the line, the interferometric signal in the H-alpha channel is calibrated with respect to the continuum channels. The observations used in this study represent the highest spatial resolution measurements of the H-alpha-emitting regions of Be stars obtained to date. These observations allow us to demonstrate for the first time that the intensity distribution in the circumstellar region of a Be star cannot be represented by uniform disk or ring-like structures, whereas a Gaussian intensity distribution appears to be fully consistent with our observations. △ Less

Submitted 3 February, 2006; originally announced February 2006.

Comments: 23 pages, 14 figures, accepted for publication in AJ

A minimum hypothesis explanation for an IMF with a lognormal body and power law tail

Authors: Shantanu Basu, C. E. Jones

Abstract: We present a minimum hypothesis model for an IMF that resembles a lognormal distribution at low masses but has a distinct power-law tail. Even if the central limit theorem ensures a lognormal distribution of condensation masses at birth, a power-law tail in the distribution arises due to accretion from the ambient cloud, coupled with a non-uniform (exponential) distribution of accretion times. We present a minimum hypothesis model for an IMF that resembles a lognormal distribution at low masses but has a distinct power-law tail. Even if the central limit theorem ensures a lognormal distribution of condensation masses at birth, a power-law tail in the distribution arises due to accretion from the ambient cloud, coupled with a non-uniform (exponential) distribution of accretion times. △ Less

Submitted 15 November, 2004; originally announced November 2004.

Comments: 2 pages, 1 figure, to appear in IMF@50, eds. E. Corbelli, F. Palla, and H. Zinnecker, Kluwer, Astrophysics and Space Science Library

On the Power-Law Tail in the Mass Function of Protostellar Condensations and Stars

Authors: Shantanu Basu, C. E. Jones

Abstract: We explore the idea that the power-law tail in the mass function of protostellar condensations and stars arises from the accretion of ambient cloud material on to a condensation, coupled with a nonuniform (exponential) distribution of accretion lifetimes. This model allows for the generation of power-law distributions in all star-forming regions, even if condensations start with a lognormal mass… ▽ More We explore the idea that the power-law tail in the mass function of protostellar condensations and stars arises from the accretion of ambient cloud material on to a condensation, coupled with a nonuniform (exponential) distribution of accretion lifetimes. This model allows for the generation of power-law distributions in all star-forming regions, even if condensations start with a lognormal mass distribution, as may be expected from the central limit theorem, and supported by some recent numerical simulations of turbulent molecular clouds. For a condensation mass $m$ with growth rate $dm/dt \propto m$, an analytic three-parameter probability density function is derived; it resembles a lognormal at low mass and has a pure power-law high-mass tail. An approximate power-law tail is also expected for other growth laws, and we calculate the distribution for the plausible case $dm/dt \propto m^{2/3}$. Furthermore, any single time snapshot of the masses of condensations that are still accreting (and are of varying ages) also yields a distribution with a power-law tail similar to that of the IMF. △ Less

Submitted 16 November, 2003; originally announced November 2003.

Comments: 5 pages, 4 figures, to appear in MNRAS (Letters section)

Journal ref: Mon.Not.Roy.Astron.Soc. 347 (2004) L47

The peculiar circumstellar environment of NGC2024-IRS2

Authors: A. Lenorzer, A. Bik, A. de Koter, S. E. Kurtz, L. B. F. M. Waters, L. Kaper, C. E. Jones, T. R. Geballe

Abstract: We re-examine the nature of NGC2024-IRS2 in light of the recent discovery of the late O-type star, IRS2b, located 5 arcsec from IRS2. Using L-band spectroscopy, we set a lower limit of Av = 27.0 mag on the visual extinction towards IRS2. Arguments based on the nature of the circumstellar material, favor an Av of 31.5 mag. IRS2 is associated with the UCHII region G206.543-16.347 and the infrared… ▽ More We re-examine the nature of NGC2024-IRS2 in light of the recent discovery of the late O-type star, IRS2b, located 5 arcsec from IRS2. Using L-band spectroscopy, we set a lower limit of Av = 27.0 mag on the visual extinction towards IRS2. Arguments based on the nature of the circumstellar material, favor an Av of 31.5 mag. IRS2 is associated with the UCHII region G206.543-16.347 and the infrared source IRAS 05393-0156. We show that much of the mid-infrared emission towards IRS2, as well as the far infrared emission peaking at ~ 100 micron, do not originate in the direct surroundings of IRS2, but instead from an extended molecular cloud. Using new K-, L- and L'-band spectroscopy and a comprehensive set of infrared and radio continuum measurements from the literature, we apply diagnostics based on the radio slope, the strength of the infrared hydrogen recombination lines, and the presence of CO band-heads to constrain the nature and spatial distribution of the circumstellar material of IRS2. Using simple gaseous and/or dust models of prescribed geometry, we find strong indications that the infrared flux originating in the circumstellar material of IRS2 is dominated by emission from a dense gaseous disk with a radius of about 0.6 AU. At radio wavelengths the flux density distribution is best described by a stellar wind recombining at a radius of about 100 AU. Although NGC2024/IRS2 shares many similarities with BN-like objects, we do not find evidence for the presence of a dust shell surrounding this object. Therefore, IRS2 is likely more evolved. △ Less

Submitted 9 October, 2003; originally announced October 2003.

Comments: 15 pages, 10 figures. Accepted for publication in Astronomy & Astrophysics

Journal ref: Astron.Astrophys. 414 (2004) 245-259

The Intrinsic Shapes of Molecular Cloud Fragments over a Range of Length Scales

Authors: C. E. Jones, Shantanu Basu

Abstract: We decipher intrinsic three-dimensional shape distributions of molecular clouds, cloud cores, Bok globules, and condensations using recently compiled catalogues of observed axis ratios for these objects mapped in carbon monoxide, ammonia, through optical selection, or in continuum dust emission. We apply statistical techniques to compare assumed intrinsic axis ratio distributions with observed p… ▽ More We decipher intrinsic three-dimensional shape distributions of molecular clouds, cloud cores, Bok globules, and condensations using recently compiled catalogues of observed axis ratios for these objects mapped in carbon monoxide, ammonia, through optical selection, or in continuum dust emission. We apply statistical techniques to compare assumed intrinsic axis ratio distributions with observed projected axis ratio distributions. Intrinsically triaxial shapes produce projected distributions which agree with observations. Molecular clouds mapped in $^{12}$CO are intrinsically triaxial but more nearly prolate than oblate, while the smaller cloud cores, Bok globules, and condensations are also intrinsically triaxial but more nearly oblate than prolate. △ Less

Submitted 19 December, 2001; originally announced December 2001.

Comments: 12 pages, 11 figures. Version with color figures can be found at http://www.astro.uwo.ca/~cjones/ or http://www.astro.uwo.ca/~basu/. To appear in ApJ, 10 April 2002, v. 569, no. 1

Intrinsic Shapes of Molecular Cloud Cores

Authors: C. E. Jones, Shantanu Basu, John Dubinski

Abstract: We conduct an analysis of the shapes of molecular cloud cores using recently compiled catalogs of observed axis ratios of individual cores mapped in ammonia or through optical selection. We apply both analytical and statistical techniques to deproject the observed axis ratios in order to determine the true distribution of cloud core shapes. We find that neither pure oblate nor pure prolate cores… ▽ More We conduct an analysis of the shapes of molecular cloud cores using recently compiled catalogs of observed axis ratios of individual cores mapped in ammonia or through optical selection. We apply both analytical and statistical techniques to deproject the observed axis ratios in order to determine the true distribution of cloud core shapes. We find that neither pure oblate nor pure prolate cores can account for the observed distribution of core shapes. Intrinsically triaxial cores produce distributions which agree with observations. The best-fit triaxial distribution contains cores which are more nearly oblate than prolate. △ Less

Submitted 8 January, 2001; originally announced January 2001.

Comments: 20 pages, 10 figures. To appear in ApJ (2001 April 1). Color figures available at http://www.astro.uwo.ca/~cjones/ or http://www.astro.uwo.ca/~basu/pub.html