Showing 1–50 of 72 results for author: Curé, M

Evolution of rotating massive stars adopting a newer, self-consistent wind prescription at SMC metallicity

Authors: Alex Camilo Gormaz-Matamala, Jorge Cuadra, Sylvia Ekström, Georges Meynet, Michel Curé, Krzysztof Belczynski

Abstract: We use Geneva-evolution-code to run evolutionary tracks for stellar masses ranging from $20$ to $85$ $M_\odot$ at SMC metallicity ($Z=0.002$). We upgrade the recipe for stellar winds by adopting our self-consistent m-CAK prescription, which reduces the value of mass-loss rate by a factor between 2 and 6 depending on the mass range. The impact of our new winds is wide, and it can be divided betwe… ▽ More We use Geneva-evolution-code to run evolutionary tracks for stellar masses ranging from $20$ to $85$ $M_\odot$ at SMC metallicity ($Z=0.002$). We upgrade the recipe for stellar winds by adopting our self-consistent m-CAK prescription, which reduces the value of mass-loss rate by a factor between 2 and 6 depending on the mass range. The impact of our new winds is wide, and it can be divided between direct and indirect impact. For the most massive models ($60$ and $85$ $M_\odot$) with $\dot M\gtrsim2\times10^{-7}$ $M_\odot$ yr$^{-1}$, the impact is direct because lower mass loss make stars remove less envelope and therefore remain more massive and less chemically enriched at their surface at the end of their MS phase. For the less massive models ($20$ and $25$ $M_\odot$) with $\dot M\lesssim2\times10^{-8}$ $M_\odot$ yr$^{-1}$, the impact is indirect because lower mass loss make the stars keep high rotational velocities for a longer period of time, then extending the H-core burning lifetime and reaching the end of the MS with higher surface enrichment. Given that the conditions at the H-depletion change, the stars will lose more mass during their He-core burning stages anyways. For $M_\text{zams}=20$ to $40$ $M_\odot$, our models predict stars will evolve through the Hertzsprung gap, from O-type supergiants to BSG and finally RSG, with larger mass fractions of helium compared to old evolution models. New models also set down to $M_\text{zams}=85\,M_\odot$ the minimal initial mass required for a single star to become WR at metallicity $Z=0.002$. New values for $\dot M$ need to be complemented with upgrades in additional features such as convective core overshooting and distribution of rotational velocities, besides more detailed observations from projects such as XShootU, in order to provide a robust framework for the study of massive stars at low metallicity environments. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 687, A290 (2024)

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.

Radiation-Driven Wind Hydrodynamics of Massive Stars: A Review

Authors: Michel Cure, Ignacio Araya

Abstract: Mass loss from massive stars plays a determining role in their evolution through the upper Hertzsprung-Russell diagram. The hydrodynamic theory that describes their steady-state winds is the line-driven wind theory (m-CAK). From this theory, the mass loss rate and the velocity profile of the wind can be derived, and estimating these properly will have a profound impact on quantitative spectroscopy… ▽ More Mass loss from massive stars plays a determining role in their evolution through the upper Hertzsprung-Russell diagram. The hydrodynamic theory that describes their steady-state winds is the line-driven wind theory (m-CAK). From this theory, the mass loss rate and the velocity profile of the wind can be derived, and estimating these properly will have a profound impact on quantitative spectroscopy analyses from the spectra of these objects. Currently, the so-called beta-law, which is an approximation for the fast solution, is widely used instead of m-CAK hydrodynamics, and when the derived value is beta greater than 1.2, there is no hydrodynamic justification for these values. This review focuses on (1) a detailed topological analysis of the equation of motion (EoM), (2) solving the EoM numerically for all three different (fast and two slow) wind solutions, (3) deriving analytical approximations for the velocity profile via the LambertW function and (4) presenting a discussion of the applicability of the slow solutions. △ Less

Submitted 19 May, 2023; originally announced May 2023.

Comments: Published in Galaxies Journal (33 pages, 18 figures)

Journal ref: Galaxies 2023, 11(3), 68

Evolution of rotating massive stars with new hydrodynamic wind models

Authors: Alex Camilo Gormaz-Matamala, Jorge Cuadra, Georges Meynet, Michel Curé

Abstract: Mass loss due to line-driven winds is central to our understanding of the evolution of massive stars. We extend the evolution models introduced in Paper I, where the mass loss recipe is based on the simultaneous calculation of the wind hydrodynamics and the line-acceleration, by incorporating the effects of stellar rotation. We introduce a grid of self-consistent line-force parameters for a set… ▽ More Mass loss due to line-driven winds is central to our understanding of the evolution of massive stars. We extend the evolution models introduced in Paper I, where the mass loss recipe is based on the simultaneous calculation of the wind hydrodynamics and the line-acceleration, by incorporating the effects of stellar rotation. We introduce a grid of self-consistent line-force parameters for a set of standard evolutionary tracks. With that, we generate a new set of evolutionary tracks with rotation for $M_\text{ZAMS}=25,40,70,$ and $120\,M_\odot$, and metallicities $Z=0.014$ and $0.006$. The self-consistent approach gives lower mass loss rates than the standard values adopted in previous evolution models. This decrease impacts strongly on the tracks of the most massive models. Weaker winds allow the star to retain more mass, but also more angular momentum. As a consequence, weaker wind models rotate faster and show a less efficient mixing in their inner stellar structure. The new tracks predict an evolution of the rotational velocities through the MS in close agreement with the range of $\varv\sin i$ values found by recent surveys of Galactic O-type stars. As subsequent implications, the weaker winds from self-consistent models suggest a reduction of the contribution of the isotope $^{26}$Al to the ISM due to stellar winds of massive stars during the MS phase. Moreover, the higher luminosities found for the self-consistent evolutionary models suggest that some populations of massive stars might be less massive than previously thought, as in the case of Ofpe stars at the Galactic Centre. Therefore, this study opens a wide range of consequences for further research based on the evolution of massive stars. △ Less

Submitted 23 March, 2023; originally announced March 2023.

Comments: Accepted for publication in Astronomy & Astrophysics journal

Journal ref: A&A 673, A109 (2023)

V4142 Sgr: a Double Periodic Variable with an accretor surrounded by the accretion-disk's atmosphere

Authors: J. A. Rosales, R. E. Mennickent, G. Djurašević, I. Araya, M. Curé, D. R. G. Schleicher, J. Petrović

Abstract: Context: A detailed study of the close interacting binary V4142\,Sgr based on photometric and spectroscopic analysis is presented.This system belongs to the enigmatic class of Algol-like variables showing a long photometric cycle of unknown nature. Aims: Performing photometric data-mining and spectroscopic observations covering the orbital cycle, we obtain the orbital parameters and the stellar pr… ▽ More Context: A detailed study of the close interacting binary V4142\,Sgr based on photometric and spectroscopic analysis is presented.This system belongs to the enigmatic class of Algol-like variables showing a long photometric cycle of unknown nature. Aims: Performing photometric data-mining and spectroscopic observations covering the orbital cycle, we obtain the orbital parameters and the stellar properties of the binary system, along with the physical properties of the accretion disk located around the hot star. Insights on the evolutive path of the system are obtained. Methods: The light curve was modeled through an inverse modeling method using a theoretical light curve of the binary system, considering the light curve contribution of both stars and the accretion disk of the hot star to obtain the fundamental parameters. To constrain the main stellar parameters the mass ratio was fixed, as well as the donor temperature using the obtained values from our spectroscopic analysis including deblending methods to isolate the spectral lines of the stellar components. The system parameters were compared with a grid of binary star evolutive models in order to get insights on the evolutionary history of the system. Results: The orbital period and the long cycle were re-calculated and found to be of $30.633 \pm 0.002 ~\mathrm{days}$ and $1201 \pm 14 ~\mathrm{days}$. The spectral analysis reveals H$α$ double emission with a persistent $V \leq R$ asymmetry which is considered evidence of a possible wind emergin from the hotspot region... △ Less

Submitted 21 February, 2023; originally announced February 2023.

Journal ref: Astronomy & Astrophysics,670,2023,A94

Evolution of massive stars with new hydrodynamic wind models

Authors: A. C. Gormaz-Matamala, M. Curé, G. Meynet, J. Cuadra, J. H. Groh, L. J. Murphy

Abstract: Here we present evolutionary models for a set of massive stars, introducing a new prescription for the mass-loss rate obtained from hydrodynamical calculations in which the wind velocity profile, $v(r)$, and the line-acceleration, $g_\text{line}$, are obtained in a self consistently way. Replacing mass-loss rates at the Main Sequence stage from the standard Vink's formula by our new recipe, we gen… ▽ More Here we present evolutionary models for a set of massive stars, introducing a new prescription for the mass-loss rate obtained from hydrodynamical calculations in which the wind velocity profile, $v(r)$, and the line-acceleration, $g_\text{line}$, are obtained in a self consistently way. Replacing mass-loss rates at the Main Sequence stage from the standard Vink's formula by our new recipe, we generate a new set of evolutionary tracks for $M_\text{ZAMS}=25,40,70$ and $120\,M_\odot$ and metallicities $Z=0.014$ (Galactic), $Z=0.006$ (LMC), and $Z=0.002$ (SMC). Our new derived formula for mass-loss rate predicts a dependence $\dot M\propto Z^a$, where $a$ is not longer constant but dependent on the stellar mass: ranging from $a\sim0.53$ when $M_*\sim120\;M_\odot$, to $a\sim1.02$ when $M_*\sim25\;M_\odot$. We found that models adopting the new recipe for $\dot M$ retain more mass during their evolution, which is expressed in larger radii and consequently more luminous tracks over the Hertzsprung-Russell diagram. These differences are more prominent for the cases of $M_\text{ZAMS}=70$ and 120 $M_\odot$ at solar metallicity, where we found self-consistent tracks are $\sim0.1$ dex brighter and keep extra mass up to 20 $M_\odot$, compared with the classical models using the previous formulation for mass-loss rate. Moreover, we observed remarkable differences for the evolution of the radionuclide isotope $^{26}$Al in the core and the surface of the star. Since $\dot M_\text{sc}$ are weaker than the commonly adopted values for evolutionary tracks, self-consistent tracks predict a later modification in the abundance number of $^{26}$Al in the stellar winds. This new behaviour could provide useful information about the real contribution of this isotope from massive stars to the Galactic interstellar medium. △ Less

Submitted 11 July, 2022; originally announced July 2022.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 665, A133 (2022)

Revisiting Viscous Transonic Decretion Disks of Be Stars

Authors: Michel Curé, Rodrigo Meneses, Ignacio Araya, Catalina Arcos, Greco Peña, Natalia Machuca, Abigali Rodriguez

Abstract: In the context of Be stars, we restudied the viscous transonic decretion disk model of these stars. This model is driven by a radiative force due to an ensemble of optically-thin lines and viscosity considering the Shakura Sunyaev prescription. The non-linear equation of motion presents a singularity (sonic point) and an eigenvalue, which is also the initial condition at the stellar surface. Then,… ▽ More In the context of Be stars, we restudied the viscous transonic decretion disk model of these stars. This model is driven by a radiative force due to an ensemble of optically-thin lines and viscosity considering the Shakura Sunyaev prescription. The non-linear equation of motion presents a singularity (sonic point) and an eigenvalue, which is also the initial condition at the stellar surface. Then, to obtain this eigenvalue, we set it as a radial quantity and perform a detailed topological analysis. Thereafter, we describe a numerical method for solving either Nodal and Saddle transonic solutions. The value of the viscosity,"alpha", barely determine the location of the sonic point, but it determines the topology of the solution. We found two Nodal solutions, which are almost indistinguishable between them. Saddle solutions are founded for lower values of "alpha" than the required of the Nodal solutions. In addition, rotational velocity do not play a determine role in the velocity (and density) profile, because viscosity effects collapse all the solutions to almost a unique one in a small region above the stellar surface. A suitable combination of line force parameters and/or disk temperature, give location of the sonic point lower than 50 stellar radii, describing a truncated disk. This could explain the SED turndown observed in Be stars without needing a binary companion. △ Less

Submitted 17 June, 2022; originally announced June 2022.

Comments: Accepted in A&A

Journal ref: A&A 664, A185 (2022)

Application limit of the photocentre displacement to fundamental stellar parameters of fast rotators -- Illustration on the edge-on fast rotator Regulus

Authors: M. Hadjara, R. G. Petrov, S. Jankov, P. Cruzalèbes, A. Boskri, A. Spang, S. Lagarde, J. He, X. Chen, C. Nitschelm, E. S. G. de Almeida, G. Pereira, E. A. Michael, Q. Gao, W. Wang, I. Reyes, C. Arcos, I. Araya, M. Curé

Abstract: Differential Interferometry allows to obtain the differential visibility and phase, in addition to the spectrum. The differential phase contains important information about the structure and motion of stellar photosphere such as stellar spots and non-radial pulsations, and particularly the rotation. Thus, this interferometric observable strongly helps to constrain the stellar fundamental parameter… ▽ More Differential Interferometry allows to obtain the differential visibility and phase, in addition to the spectrum. The differential phase contains important information about the structure and motion of stellar photosphere such as stellar spots and non-radial pulsations, and particularly the rotation. Thus, this interferometric observable strongly helps to constrain the stellar fundamental parameters of fast rotators. The spectro-astrometry mainly uses the photocentre displacements, which is a first approximation of the differential phase, and is applicable only for unresolved or marginally objects. We study here the sensitivity of relevant stellar parameters to the simulated photocentres using the SCIROCCO code: a semi-analytical algorithm dedicated to fast rotators, applied to two theoretical modeling stars based on Achernar and Regulus, in order to classify the importance of these parameters and their impact on the modeling. We compare our simulations with published VLTI/AMBER data. This current work sets the limits of application of photocentre displacements to fast rotators, and under which conditions we can use the photocentres and/or the differential phase, through a pre-established physical criterion. To validate our theoretical study, we apply our method of analysis on observed data of the edge-on fast rotator Regulus. For unresolved targets, with a visibility $V\sim 1$, the photocentre can constrain the main stellar fundamental parameters of fast rotators, whereas from marginally resolved objects ($0.8 \leq V < 1$), mainly the rotation axis position angle ($\rm PA_{\rm rot}$) can be directly deduced from the vectorial photocentre displacement, which is very important for young cluster studies. △ Less

Submitted 19 March, 2022; originally announced March 2022.

Journal ref: MNRAS 511, 4724-4740 (2022)

New self-consistent wind parameters fitting optical spectra of O-type stars observed with HERMES spectrograph

Authors: Alex C. Gormaz-Matamala, Michel Curé, Alex Lobel, Jorge A. Panei, Jorge Cuadra, Ignacio Araya, Catalina Arcos, Felipe Figueroa-Tapia

Abstract: We perform spectral fitting for a set of O-type stars based on self-consistent wind solutions, which provide mass-loss rate and velocity profiles directly derived from the initial stellar parameters. The great advantage of this self-consistent spectral fitting is therefore the reduction of the number of free parameters to be tuned. Self-consistent values for the line-force parameters (k,alpha,de… ▽ More We perform spectral fitting for a set of O-type stars based on self-consistent wind solutions, which provide mass-loss rate and velocity profiles directly derived from the initial stellar parameters. The great advantage of this self-consistent spectral fitting is therefore the reduction of the number of free parameters to be tuned. Self-consistent values for the line-force parameters (k,alpha,delta) and subsequently for the mass-loss rate and terminal velocity are provided by the m-CAK prescription introduced in Paper I, updated in this work with improvements such as a temperature structure for the wind, self-consistently evaluated from the line-acceleration. Synthetic spectra are calculated using the radiative transfer code FASTWIND, replacing the classical beta-law for our new calculated velocity profiles. We found that self-consistent m-CAK solutions provide values for theoretical mass-loss rates on the order of the most recent predictions of other studies. From here, we generate synthetic spectra with self-consistent hydrodynamics to fit and obtain a new set of stellar and wind parameters for our sample of O-type stars whose spectra was taken with the high resolution echelle spectrograph HERMES (R=85000). We find a satisfactory global fit for our observations, with good accuracy for photospheric He I and He II lines and a quite acceptable fit for H lines. Although this self-consistent spectral analysis is currently constrained in the optical wavelength range only, this is an important step towards the determination of stellar and wind parameters without using a beta-law. Given these results, we expect that the values introduced here should be helpful for future studies about the stars constituting this sample, together with the prospective that the m-CAK self-consistent prescription be extended to numerous studies about massive stars in future. △ Less

Submitted 3 March, 2022; v1 submitted 17 February, 2022; originally announced February 2022.

Comments: Accepted for publication in Astronomy & Astrophysics, in 17 February 2022 In this updated version, we corrected information about funding data in Acknowledgements section

Journal ref: A&A 661, A51 (2022)

Automatic algorithm to obtain $v \sin i$ values via Fourier Transform in BeSOS database

Authors: Martín Solar, Catalina Arcos, Michel Curé, Ronaldo S. Levenhagen, Ignacio Araya

Abstract: Be stars are found to rotate close to their critical rotation and therefore they are considered as an important laboratory of study for stellar rotation. In this context, we obtain the projected rotational velocity of a sample of classical Be Southern stars in the BeSOS database via Fourier Transform in an automated way for several absorption lines at different epochs. A Gaussian profile is fitted… ▽ More Be stars are found to rotate close to their critical rotation and therefore they are considered as an important laboratory of study for stellar rotation. In this context, we obtain the projected rotational velocity of a sample of classical Be Southern stars in the BeSOS database via Fourier Transform in an automated way for several absorption lines at different epochs. A Gaussian profile is fitted to eight observed photospheric HeI lines in order to select automatically the spectral signal given by areas under the curve of 95.45%, 98.75% and 99.83% from the profile to obtain $v \sin i$ via Fourier Transform technique. The values obtained are in global agreement with the literature. Analysing only one line is not enough to set the $v \sin i$ value, depending on the line the value in most cases are underestimated with respect to $λ$4471. When gravity darkening effects are including, apparent values increases by $\sim10$%. The resolution of the instrument PUCHEROS used for BeSOS spectra ($R \sim 17\,000$) constrain a theoretical lower bound possible at $v \sin i \sim 100$ km s$^{-1}$. The procedure has limitations using a linear limb-darkening function with $\varepsilon = 0.6$ for classical Be stars rotating close to the break-up velocity without gravity-darkening corrections, which can't be negligible. Previous works measure $v \ sin i$ values using just one spectral line and here we demonstrate that with more lines the results can varies. This could be due to the photospheric distribution of atomic transitions on classical Be stars. △ Less

Submitted 21 January, 2022; originally announced January 2022.

Comments: 14 pages, 10 figures, MNRAS, accepted

Fundamental parameters of the eclipsing binary DD CMa and evidence for mass exchange

Authors: J. A. Rosales, R. E. Mennickent, G. Djurašević, J. F. González, I. Araya, M. Cabezas, D. R. G. Schleicher, M. Curé

Abstract: We present a detailed photometric and spectroscopic analysis of DD CMa, based on published survey photometry and new spectroscopic data. We find an improved orbital period of $P_\mathrm{o}= 2.0084530 \pm 0.0000006 ~\mathrm{d}$. Our spectra reveal H$β$ and H$α$ absorptions with weak emission shoulders and we also find color excess in the WISE multiband photometry, interpreted as signatures of circu… ▽ More We present a detailed photometric and spectroscopic analysis of DD CMa, based on published survey photometry and new spectroscopic data. We find an improved orbital period of $P_\mathrm{o}= 2.0084530 \pm 0.0000006 ~\mathrm{d}$. Our spectra reveal H$β$ and H$α$ absorptions with weak emission shoulders and we also find color excess in the WISE multiband photometry, interpreted as signatures of circumstellar matter. We model the $V$-band orbital light curve derived from the ASAS and ASAS-SN surveys, assuming a semidetached configuration and using the mass ratio and temperature of the hotter star derived from our spectroscopic analysis. Our model indicates that the system consists of a B 2.5 dwarf and a B 9 giant of radii 3.2 and 3.7 $\mathrm{R_{\odot}}$, respectively, orbiting in a circular orbit of radius 6.75 $\mathrm{R_{\odot}}$. We also found $M_{\mathrm{c}} = 1.7 \pm 0.1 ~\mathrm{M_{\odot}}$, $T_{\mathrm{c}} = 11350 \pm 100 ~\mathrm{K}$ and $M_{\mathrm{h}} = 6.4 \pm 0.1 ~\mathrm{M_{\odot}}$, $T_{\mathrm{h}} = 20000 \pm 500 ~\mathrm{K}$, for the cooler and hotter star, respectively. We find broad single emission peaks in H$α$ and H$β$ after subtracting the synthetic stellar spectra. Our results are consistent with mass exchange between the stars, and suggest the existence of a stream of gas being accreted onto the early B-type star. △ Less

Submitted 1 September, 2021; originally announced September 2021.

Comments: 14 pages with 13 figures

Journal ref: The Astronomical Journal, 161, 2012, 165

Stellar and accretion disk parameters of the close binary HD 50526

Authors: J. A. Rosales, R. E. Mennickent, G. Djurašević, D. R. G. Schleicher, S. Zharikov, I. Araya, L. Celedón, M. Curé

Abstract: We present a photometric and spectroscopic study of HD 50526, an ellipsoidal binary member of the group Double Periodic Variable stars. Performing data-mining in photometric surveys and conducting new spectroscopic observations with several spectrographs during 2008 to 2015, we obtained orbital and stellar parameters of the system. The radial velocities were analyzed with the genetic PIKAIA algori… ▽ More We present a photometric and spectroscopic study of HD 50526, an ellipsoidal binary member of the group Double Periodic Variable stars. Performing data-mining in photometric surveys and conducting new spectroscopic observations with several spectrographs during 2008 to 2015, we obtained orbital and stellar parameters of the system. The radial velocities were analyzed with the genetic PIKAIA algorithm, whereas Doppler tomography maps for the H$α$ and H$β$ lines were constructed with the Total Variation Minimization code. An optimized simplex-algorithm was used to solve the inverse-problem adjusting the light curve with the best stellar parameters for the system. We find an orbital period of $6.701 \pm 0.001 ~\mathrm{d}$ and a long photometric cycle of $191 \pm 2 ~\mathrm{d}$. We detected the spectral features of the coldest star, and modeled it with a $\log{g} = 2.79 \pm 0.02 ~\mathrm{dex}$ giant of mass $1.13 \pm 0.02 ~\mathrm{M_{\odot}}$ and effective temperature $10500 \pm 125 ~\mathrm{K}$. In addition, we determine a mass ratio $q= 0.206 \pm 0.033$ and that the hot star is a B-type dwarf of mass $5.48 \pm 0.02 ~\mathrm{M_{\odot}}$. The $V$-band orbital light curve can be modeled including the presence of an accretion disk around the hotter star. This fills the Roche lobe of the hotter star, and has a radius $14.74 \pm 0.02 ~\mathrm{R_{\odot}}$ and temperature at the outer edge $9400 ~\mathrm{K}$. Two bright spots located in the disk account for the global morphology of the light curve. The Doppler tomography maps of H$α$ and H$β$, reveal complex structures of mass fluxes in the system. △ Less

Submitted 1 September, 2021; originally announced September 2021.

Comments: 16 pages with 14 figures

Journal ref: The Astronomical Journal, 162, 2021, 66

New hydrodynamic solutions for line-driven winds of hot massive stars using Lambert $W$-function

Authors: Alex C. Gormaz-Matamala, Michel Curé, D. John Hillier, Francisco Najarro, Brankica Kubátová, Jiří Kubát

Abstract: Hot massive stars present strong stellar winds that are driven by absorption, scattering and re\-emission of photons by the ions of the atmosphere (\textit{line-driven winds}). A better comprehension of this phenomenon, and a more accurate calculation of hydrodynamics and radiative acceleration is required to reduce the number of free parameters in spectral fitting, to determine accurate wind para… ▽ More Hot massive stars present strong stellar winds that are driven by absorption, scattering and re\-emission of photons by the ions of the atmosphere (\textit{line-driven winds}). A better comprehension of this phenomenon, and a more accurate calculation of hydrodynamics and radiative acceleration is required to reduce the number of free parameters in spectral fitting, to determine accurate wind parameters such as mass-loss rates and velocity profiles. We use the non-LTE model-atmosphere code CMFGEN to numerically solve the radiative transfer equation in the stellar atmosphere and to calculate the radiative acceleration $g_\text{rad}(r)$. Under the assumption that the radiative acceleration depends only on the radial coordinate, we solve analytically the equation of motion by means of the Lambert $W$-function. An iterative procedure between the solution of the radiative transfer and the equation of motion is executed in order to obtain a final self-consistent velocity field that is no longer based on any $β$-law. We apply the Lambert-procedure to three O supergiant stars ($ζ$-Puppis, HD~165763 and $α$-Cam) and discuss the Lambert-solutions for the velocity profiles. It is found that, even without recalculation of the mass-loss rate, the Lambert-procedure allows the calculation of consistent velocity profiles that reduce the number of free parameters when a spectral fitting using CMFGEN is performed. Synthetic spectra calculated from our Lambert-solutions show significant differences compared to the initial $β$-law CMFGEN models. The results indicate the importance of consistent velocity profile calculation in the CMFGEN code and its usage in a fitting procedure and interpretation of observed spectra. △ Less

Submitted 3 March, 2022; v1 submitted 28 June, 2021; originally announced June 2021.

Comments: Published in the Astrophysical Journal (ApJ), in 13 October 2021

Journal ref: 2021 The Astrophysical Journal, Volume 920, Issue 1, id.64, 30 pp

Analytical Solutions for Radiation-Driven Winds in Massive Stars II: The $δ$-slow Regime

Authors: I. Araya, A. Christen, M. Curé, L. S. Cidale, R. O. J. Venero, C. Arcos, A. C. Gormaz-Matamala, M. Haucke, P. Escárate, H. Clavería

Abstract: Accurate mass-loss rates and terminal velocities from massive stars winds are essential to obtain synthetic spectra from radiative transfer calculations and to determine the evolutionary path of massive stars. From a theoretical point of view, analytical expressions for the wind parameters and velocity profile would have many advantages over numerical calculations that solve the complex non-linear… ▽ More Accurate mass-loss rates and terminal velocities from massive stars winds are essential to obtain synthetic spectra from radiative transfer calculations and to determine the evolutionary path of massive stars. From a theoretical point of view, analytical expressions for the wind parameters and velocity profile would have many advantages over numerical calculations that solve the complex non-linear set of hydrodynamic equations. In a previous work, we obtained an analytical description for the fast wind regime. Now, we propose an approximate expression for the line-force in terms of new parameters and obtain a velocity profile closed-form solution (in terms of the Lambert $W$ function) for the $δ$-slow regime. Using this analytical velocity profile, we were able to obtain the mass-loss rates based on the m-CAK theory. Moreover, we established a relation between this new set of line-force parameters with the known stellar and m-CAK line-force parameters. To this purpose, we calculated a grid of numerical hydrodynamical models and performed a multivariate multiple regression. The numerical and our descriptions lead to good agreement between their values. △ Less

Submitted 7 April, 2021; originally announced April 2021.

Comments: 8 pages, 1 figure, accepted for publication in MNRAS

A triple star system with a misaligned and warped circumstellar disk shaped by disk tearing

Authors: Stefan Kraus, Alexander Kreplin, Alison K. Young, Matthew R. Bate, John D. Monnier, Tim J. Harries, Henning Avenhaus, Jacques Kluska, Anna S. E. Laws, Evan A. Rich, Matthew Willson, Alicia N. Aarnio, Fred C. Adams, Sean M. Andrews, Narsireddy Anugu, Jaehan Bae, Theo ten Brummelaar, Nuria Calvet, Michel Curé, Claire L. Davies, Jacob Ennis, Catherine Espaillat, Tyler Gardner, Lee Hartmann, Sasha Hinkley , et al. (7 additional authors not shown)

Abstract: Young stars are surrounded by a circumstellar disk of gas and dust, within which planet formation can occur. Gravitational forces in multiple star systems can disrupt the disk. Theoretical models predict that if the disk is misaligned with the orbital plane of the stars, the disk should warp and break into precessing rings, a phenomenon known as disk tearing. We present observations of the triple… ▽ More Young stars are surrounded by a circumstellar disk of gas and dust, within which planet formation can occur. Gravitational forces in multiple star systems can disrupt the disk. Theoretical models predict that if the disk is misaligned with the orbital plane of the stars, the disk should warp and break into precessing rings, a phenomenon known as disk tearing. We present observations of the triple star system GWOrionis, finding evidence for disk tearing. Our images show an eccentric ring that is misaligned with the orbital planes and the outer disk. The ring casts shadows on a strongly warped intermediate region of the disk. If planets can form within the warped disk, disk tearing could provide a mechanism for forming wide-separation planets on oblique orbits. △ Less

Submitted 4 September, 2020; v1 submitted 2 April, 2020; originally announced April 2020.

Comments: 63 pages, 4+13 Figures, 6 Tables, published at Science

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

AutoRegressive Planet Search: Application to the Kepler Mission

Authors: Gabriel A. Caceres, Eric D. Feigelson, G. Jogesh Babu, Natalia Bahamonde, Alejandra Christen, Karine Bertin, Cristian Meza, Michel Curé

Abstract: The 4-year light curves of 156,717 stars observed with NASA's Kepler mission are analyzed using the AutoRegressive Planet Search (ARPS) methodology described by Caceres et al. (2019). The three stages of processing are: maximum likelihood ARIMA modeling of the light curves to reduce stellar brightness variations; constructing the Transit Comb Filter periodogram to identify transit-like periodic di… ▽ More The 4-year light curves of 156,717 stars observed with NASA's Kepler mission are analyzed using the AutoRegressive Planet Search (ARPS) methodology described by Caceres et al. (2019). The three stages of processing are: maximum likelihood ARIMA modeling of the light curves to reduce stellar brightness variations; constructing the Transit Comb Filter periodogram to identify transit-like periodic dips in the ARIMA residuals; Random Forest classification trained on Kepler Team confirmed planets using several dozen features from the analysis. Orbital periods between 0.2 and 100 days are examined. The result is a recovery of 76% of confirmed planets, 97% when period and transit depth constraints are added. The classifier is then applied to the full Kepler dataset; 1,004 previously noticed and 97 new stars have light curve criteria consistent with the confirmed planets, after subjective vetting removes clear False Alarms and False Positive cases. The 97 Kepler ARPS Candidate Transits mostly have periods $P<10$ days; many are UltraShort Period hot planets with radii $<1$% of the host star. Extensive tabular and graphical output from the ARPS time series analysis is provided to assist in other research relating to the Kepler sample. △ Less

Submitted 23 May, 2019; originally announced May 2019.

Comments: 66 pages with 26 figures and 5 tables to appear in the Astronomical Journal. A version with high-resolution graphics, machine readable tables and FigureSet for the 97KACTs is available at https://drive.google.com/drive/folders/107xZIAj3C0HHqsW66Xc-8FA6O30EncVt?usp=sharing

A method to deconvolve stellar rotational velocities III. The probability distribution function via Maximum Likelihood utilizing Finite Distribution Mixtures

Authors: Rafael Orellana, Pedro Escarate, Michel Cure, Alejandra Christen, Rodrigo Carvajal, Juan Carlos Agüero

Abstract: The study of accurate methods to estimate the distribution of stellar rotational velocities is important for understanding many aspects of stellar evolution. From such observations we obtain the projected rotational speed v sin(i) in order to recover the true distribution of the rotational velocity. To that end, we need to solve a difficult inverse problem that can be posed as a Fredholm integral… ▽ More The study of accurate methods to estimate the distribution of stellar rotational velocities is important for understanding many aspects of stellar evolution. From such observations we obtain the projected rotational speed v sin(i) in order to recover the true distribution of the rotational velocity. To that end, we need to solve a difficult inverse problem that can be posed as a Fredholm integral of the first kind. n this work we have used a novel approach based on Maximum likelihood (ML) estimation to obtain an approximation of the true rotational velocity probability density function expressed as a sum of known distribution families. In our proposal, the measurements have been treated as random variables drawn from the projected rotational velocity probability density function. We analyzed the case of Maxwellian sum approximation, where we estimated the parameters that define the sum of distributions. The performance of the proposed method is analyzed using Monte Carlo simulations considering two theoretical cases for the probability density function of the true rotational stellar velocities: i) an unimodal Maxwellian probability density distribution and ii) a bimodal Maxwellian probability density distribution. The results show that the proposed method yielded more accurate estimates in comparison with the Tikhonov regularization method, especially for small sample length N=50. Our proposal was evaluated using real data from three sets of measurements, and our findings were validated using three statistical tests. The ML approach with Maxwellian sum approximation is a accurate method to deconvolve the rotational velocity probability density function, even when the sample length is small (N= 50) △ Less

Submitted 8 March, 2019; originally announced March 2019.

Journal ref: A&A 623, A138 (2019)

Self-consistent solutions for line-driven winds of hot massive stars. The m-CAK procedure

Authors: Alex C. Gormaz-Matamala, Michel Curé, Lydia S. Cidale, Roberto O. J. Venero

Abstract: Massive stars present strong stellar that which are described by the radiation driven wind theory. Accurate mass-loss rates are necessary to properly describe the stellar evolution across the Hertzsprung--Russel Diagram. We present a self-consistent procedure that coupled the hydrodynamics with calculations of the line-force, giving as results the line-force parameters, the velocity field, and t… ▽ More Massive stars present strong stellar that which are described by the radiation driven wind theory. Accurate mass-loss rates are necessary to properly describe the stellar evolution across the Hertzsprung--Russel Diagram. We present a self-consistent procedure that coupled the hydrodynamics with calculations of the line-force, giving as results the line-force parameters, the velocity field, and the mass-loss rate. Our calculations contemplate the contribution to the line-force multiplier from more than $\sim 900,000$ atomic transitions, an NLTE radiation flux from the photosphere and a quasi-LTE approximation for the occupational numbers. A full set of line-force parameters for $T_\text{eff}\ge 32,000$ K, surface gravities higher than 3.4 dex for two different metallicities are presented, with their corresponding wind parameters (terminal velocities and mass-loss rates). The already known dependence of line-force parameters on effective temperature is enhanced by the dependence on $\log g$. The terminal velocities present a stepper scaling relation with respect to the escape velocity, this might explain the scatter values observed in the hot side of the bistability jump. Moreover, a comparison of self-consistent mass-loss rates with empirical values shows a good agreement. Self-consistent wind solutions are used as input in FASTWIND to calculate synthetic spectra. We show, comparing with the observed spectra for three stars, that varying the clumping factor, the synthetic spectra rapidly converge into the neighbourhood region of the solution. It is important to stress that our self-consistent procedure significantly reduces the number of free parameters needed to obtain a synthetic spectrum. △ Less

Submitted 1 March, 2019; originally announced March 2019.

Wind properties of variable B supergiants. Evidence of pulsations connected with mass-loss episodes

Authors: M. Haucke, L. S. Cidale, R. O. J. Venero, M. Curé, M. Kraus, S. Kanaan, C. Arcos

Abstract: Variable B supergiants (BSGs) constitute a heterogeneous group of stars with complex photometric and spectroscopic behaviours. They exhibit mass-loss variations and experience different types of oscillation modes, and there is growing evidence that variable stellar winds and photospheric pulsations are closely related. To discuss the wind properties and variability of evolved B-type stars, we deri… ▽ More Variable B supergiants (BSGs) constitute a heterogeneous group of stars with complex photometric and spectroscopic behaviours. They exhibit mass-loss variations and experience different types of oscillation modes, and there is growing evidence that variable stellar winds and photospheric pulsations are closely related. To discuss the wind properties and variability of evolved B-type stars, we derive new stellar and wind parameters for a sample of 19 Galactic BSGs by fitting theoretical line profiles of H, He, and Si to the observed ones and compare them with previous determinations. The mass-loss rate of three stars has been obtained for the first time. All stars follow the known empirical wind momentum-luminosity relationships, and the late BSGs show the trend of the mid BSGs. We find an empirical relationship that associates the amplitude of mass-loss variations with photometric/spectroscopic variability on timescales of tens of days. We find that stars located on the cool side of the bi-stability jump show a decrease in the terminal and escape velocities ratio, while their corresponding mass-loss rates are similar to or lower than the values found for stars on the hot side. Particularly, for those variable stars a decrease in the terminal and escape velocities ratio is accompanied by a decrease in the mass-loss rate. Our results also suggest that radial pulsation modes with periods longer than 6 days might be responsible for the wind variability in the mid/late-type BSGs. These radial modes might be identified with strange modes, which are known to facilitate (enhanced) mass loss. On the other hand, we propose that the wind behaviour of stars on the cool side of the bi-stability jump could fit with predictions of the delta-slow hydrodynamics solution for radiation-driven winds with highly variable ionization. △ Less

Submitted 4 February, 2019; originally announced February 2019.

Journal ref: Astronomy & Astrophysics, 06/2018, Volume 614, id.A91, 28 pp

AutoRegressive Planet Search: Methodology

Authors: Gabriel A. Caceres, Eric D. Feigelson, G. Jogesh Babu, Natalia Bahamonde, Alejandra Christen, Karine Bertin, Cristian Meza, Michel Curé

Abstract: The detection of periodic signals from transiting exoplanets is often impeded by extraneous aperiodic photometric variability, either intrinsic to the star or arising from the measurement process. Frequently, these variations are autocorrelated wherein later flux values are correlated with previous ones. In this work, we present the methodology of the Autoregessive Planet Search (ARPS) project whi… ▽ More The detection of periodic signals from transiting exoplanets is often impeded by extraneous aperiodic photometric variability, either intrinsic to the star or arising from the measurement process. Frequently, these variations are autocorrelated wherein later flux values are correlated with previous ones. In this work, we present the methodology of the Autoregessive Planet Search (ARPS) project which uses Autoregressive Integrated Moving Average (ARIMA) and related statistical models that treat a wide variety of stochastic processes, as well as nonstationarity, to improve detection of new planetary transits. Providing a time series is evenly spaced or can be placed on an evenly spaced grid with missing values, these low-dimensional parametric models can prove very effective. We introduce a planet-search algorithm to detect periodic transits in the residuals after the application of ARIMA models. Our matched-filter algorithm, the Transit Comb Filter (TCF), is closely related to the traditional Box-fitting Least Squares and provides an analogous periodogram. Finally, if a previously identified or simulated sample of planets is available, selected scalar features from different stages of the analysis -- the original light curves, ARIMA fits, TCF periodograms, and folded light curves -- can be collectively used with a multivariate classifier to identify promising candidates while efficiently rejecting false alarms. We use Random Forests for this task, in conjunction with Receiver Operating Characteristic (ROC) curves, to define discovery criteria for new, high fidelity planetary candidates. The ARPS methodology can be applied to both evenly spaced satellite light curves and densely cadenced ground-based photometric surveys. △ Less

Submitted 14 May, 2019; v1 submitted 15 January, 2019; originally announced January 2019.

Comments: 40 pages, 12 figures, to appear in the Astronomical Journal

Spectro-interferometric observations of a sample of Be stars. Setting limits to the geometry and kinematics of stable Be disks

Authors: Y. R. Cochetti, C. Arcos, S. Kanaan, A. Meilland, L. S. Cidale, M. Curé

Abstract: Be stars are rapid rotators surrounded by a gaseous disk envelope whose origin is still under debate. This envelope is responsible for observed emission lines and large infrared excess. To progress in the understanding of the physical processes involved in the disk formation, we estimate the disk parameters for a sample of Be stars and search for correlations between these parameters and stellar p… ▽ More Be stars are rapid rotators surrounded by a gaseous disk envelope whose origin is still under debate. This envelope is responsible for observed emission lines and large infrared excess. To progress in the understanding of the physical processes involved in the disk formation, we estimate the disk parameters for a sample of Be stars and search for correlations between these parameters and stellar properties. We performed spectro-interferometric observations of 26 Be stars in the region of the Br$γ$ line to study the kinematical properties of their disks through the Doppler effect. Observations were performed at the Paranal observatory with the VLTI/AMBER interferometer. This instrument provides high spectral and high spatial resolutions. We modeled 18 Be stars with emission in the Br$γ$ line. The disk kinematic is described by a quasi-Keplerian rotation law, with the exception of HD28497 that presents a one-arm density-wave structure. Using a combined sample, we derived a mean value for the velocity ratio V/Vc=0.75, and found that rotation axes are probably randomly distributed in the sky. Disk sizes in the line component model are in the range of 2-13 stellar radii and do not correlate with the effective temperature or spectral type. However, we found that the maximum size of a stable disk correlates with the rotation velocity at the inner part of the disk and the stellar mass. We found that, on average, the Be stars of our combined sample do not rotate at their critical velocity. However, the centrifugal force and mass of the star defines an upper limit size for a stable disk configuration. For a given rotation, high-mass Be stars tend to have more compact disks than their low-mass counterparts. It would be interesting to follow up the evolution of the disk size in variable stars to better understand the formation and dissipation processes of their circumstellar disks. △ Less

Submitted 27 November, 2018; originally announced November 2018.

Comments: 30 pages. This paper includes 4 tables and 36 Figures. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 621, A123 (2019)

Imaging the disc rim and a moving close-in companion candidate in the pre-transitional disc of V1247 Orionis

Authors: Matthew Willson, Stefan Kraus, Jacques Kluska, John D. Monnier, Michel Cure, Mike Sitko, Alicia Aarnio, Michael J. Ireland, Aaron Rizzuto, Edward Hone, Alexander Kreplin, Sean Andrews, Nuria Calvet, Catherine Espaillat, Misato Fukagawa, Tim J. Harries, Sasha Hinkley, Samer Kanaan, Takayuki Muto, David J. Wilner

Abstract: V1247 Orionis harbours a pre-transitional disc with a partially cleared gap. Earlier interferometric and polarimetric observations revealed strong asymmetries both in the gap region and in the outer disc. The presence of a companion was inferred to explain these asymmetric structures and the ongoing disc clearing. Using an extensive set of multi-wavelength and multi-epoch observations we aimed to… ▽ More V1247 Orionis harbours a pre-transitional disc with a partially cleared gap. Earlier interferometric and polarimetric observations revealed strong asymmetries both in the gap region and in the outer disc. The presence of a companion was inferred to explain these asymmetric structures and the ongoing disc clearing. Using an extensive set of multi-wavelength and multi-epoch observations we aimed to identify the origin of the previously detected asymmetries. We have observed V1247 Ori at three epochs spanning $\sim678$ days using sparse aperture masking interferometry with Keck/NIRC2 and VLT/NACO. In addition, we search for signs of accretion through VLT/SPHERE-ZIMPOL spectral differential imaging in H$α$ and R-band continuum. Our SMA sub-millimetre interferometry in 880 $μ$m continuum and in the CO(3-2) line allows us to constrain the orientation and direction of rotation of the outer disc. We find the L'-band emission to be dominated by static features which trace forward-scattered dust emission from the inner edge of the outer disc located to the north-east. In H- and K-band, we see evidence for a companion candidate that moved systematically by 45$^{\circ}$ within the first $\sim$345 days. The separation of the companion candidate is not well constrained, but the observed position angle change is consistent with Keplerian motion of a body located on a 6 au orbit. From the SMA CO moment map, the location of the disc rim, and the detected orbital motion, we deduced the three-dimensional orientation of the disc. We see no indication of accretion in H$α$ and set upper limits for an accreting companion. The measured contrast of the companion candidate in H and K is consistent with an actively accreting protoplanet. Hence, we identify V1247 Ori as a unique laboratory for studying companion-disc interactions and disc clearing. △ Less

Submitted 16 November, 2018; originally announced November 2018.

Comments: 17 pages, 14 figures

Journal ref: A&A 621, A7 (2019)

Probing the Inner Disk Emission of the Herbig Ae Stars HD 163296 and HD 190073

Authors: Benjamin R. Setterholm, John D. Monnier, Claire L. Davies, Alexander Kreplin, Stefan Kraus, Fabien Baron, Alicia Aarnio, Jean-Philippe Berger, Nuria Calvet, Michel Curé, Samer Kanaan, Brian Kloppenborg, Jean-Baptiste Le Bouquin, Rafael Millan-Gabet, Adam E. Rubinstein, Michael L. Sitko, Judit Sturmann, Theo A. ten Brummelaar, Yamina Touhami

Abstract: The physical processes occurring within the inner few astronomical units of proto-planetary disks surrounding Herbig Ae stars are crucial to setting the environment in which the outer planet-forming disk evolves and put critical constraints on the processes of accretion and planet migration. We present the most complete published sample of high angular resolution H- and K-band observations of the… ▽ More The physical processes occurring within the inner few astronomical units of proto-planetary disks surrounding Herbig Ae stars are crucial to setting the environment in which the outer planet-forming disk evolves and put critical constraints on the processes of accretion and planet migration. We present the most complete published sample of high angular resolution H- and K-band observations of the stars HD 163296 and HD 190073, including 30 previously unpublished nights of observations of the former and 45 nights of the latter with the CHARA long-baseline interferometer, in addition to archival VLTI data. We confirm previous observations suggesting significant near-infrared emission originates within the putative dust evaporation front of HD 163296 and show this is the case for HD 190073 as well. The H- and K-band sizes are the same within $(3 \pm 3)\%$ for HD 163296 and within $(6 \pm 10)\%$ for HD 190073. The radial surface brightness profiles for both disks are remarkably Gaussian-like with little or no sign of the sharp edge expected for a dust evaporation front. Coupled with spectral energy distribution analysis, our direct measurements of the stellar flux component at H and K bands suggest that HD 190073 is much younger (<400 kyr) and more massive (~5.6 M$_\odot$) than previously thought, mainly as a consequence of the new Gaia distance (891 pc). △ Less

Submitted 9 November, 2018; originally announced November 2018.

Comments: 19 pages, 6 figures

Resolving the kinematics of the disks around Galactic B[e] supergiants

Authors: Grigoris Maravelias, Michaela Kraus, Lydia S. Cidale, Marcelo Borges Fernandes, Maria L. Arias, Michel Curé, Georgios Vasilopoulos

Abstract: B[e] Supergiants are luminous evolved massive stars. The mass-loss during this phase creates a complex circumstellar environment with atomic, molecular, and dusty regions usually found in rings or disk-like structures. For a better comprehension of the mechanisms behind the formation of these rings, detailed knowledge about their structure and dynamics is essential. To address that, we obtained hi… ▽ More B[e] Supergiants are luminous evolved massive stars. The mass-loss during this phase creates a complex circumstellar environment with atomic, molecular, and dusty regions usually found in rings or disk-like structures. For a better comprehension of the mechanisms behind the formation of these rings, detailed knowledge about their structure and dynamics is essential. To address that, we obtained high-resolution optical and near-infrared spectra for 8 selected Galactic B[e] Supergiants, for which CO emission has been detected. Assuming Keplerian rotation for the disk, we combine the kinematics obtained from the CO bands in the near-IR with those obtained by fitting the forbidden emission [OI] $λ$5577, [OI] $λλ$6300,6363, and [CaII] $λλ$7291,7323 lines in the optical to probe the disk structure. We find that the emission originates from multiple ring structures around all B[e] Supergiants, with each one of them displaying a unique combination of rings regardless of whether the object is part of a binary system. The confirmed binaries display spectroscopic variations of their line intensities and profiles as well as photometric variability, whereas the ring structures around the single stars are stable. △ Less

Submitted 4 September, 2018; v1 submitted 2 July, 2018; originally announced July 2018.

Comments: 27 pages, 13 figures, 4 tables. Published in MNRAS. Minor corrections in this version: (a) adding references in the text to Fig. 6, tables A6 and A8, and footnote 5, (b) adding/modifying the acknowledgments section

Differences in the gas and dust distribution in the transitional disk of a sun-like young star, PDS 70

Authors: Zachary C. Long, Eiji Akiyama, Michael Sitko, Rachel B. Fernandes, Korash Assani, Carol A. Grady, Michel Cure, Ruobing Dong, Misato Fukagawa, Yasuhiro Hasegawa, Jun Hashimoto, Thomas Henning, Shu-Ichiro Inutsuka, Stefan Kraus, Jungmi Kwon, Carey M. Lisse, Hauyu Baobabu Liu, Satoshi Mayama, Takayuki Muto, Takao Nakagawa, Michihiro Takami, Motohide Tamura, Thayne Currie, John P. Wisniewski, Yi Yang

Abstract: We present ALMA 0.87 mm continuum, HCO+ J=4--3 emission line, and CO J=3--2 emission line data of the disk of material around the young, Sun-like star PDS 70. These data reveal the existence of a possible two component transitional disk system with a radial dust gap of 0."2 +/- 0."05, an azimuthal gap in the HCO+ J=4--3 moment zero map, as well as two bridge-like features in the gas data. Interest… ▽ More We present ALMA 0.87 mm continuum, HCO+ J=4--3 emission line, and CO J=3--2 emission line data of the disk of material around the young, Sun-like star PDS 70. These data reveal the existence of a possible two component transitional disk system with a radial dust gap of 0."2 +/- 0."05, an azimuthal gap in the HCO+ J=4--3 moment zero map, as well as two bridge-like features in the gas data. Interestingly these features in the gas disk have no analogue in the dust disk making them of particular interest. We modeled the dust disk using the Monte Carlo radiative transfer code HOCHUNK3D (Whitney et al. 2013) using a two disk components. We find that there is a radial gap that extends from 15-60 au in all grain sizes which differs from previous work. △ Less

Submitted 2 April, 2018; originally announced April 2018.

Co-existence and switching between fast and $Ω$-slow wind solutions in rapidly rotating massive stars

Authors: Ignacio Araya, Michel Curé, Asif ud-Doula, Alfredo Santillán, Lydia Cidale

Abstract: Most radiatively-driven massive star winds can be modelled with m-CAK theory resulting in so called fast solution. However, those most rapidly rotating among them, especially when the stellar rotational speed is higher than $\sim 75\%$ of the critical rotational speed, can adopt a different solution called $Ω$-slow solution characterized by a dense and slow wind. Here, in this work we study the tr… ▽ More Most radiatively-driven massive star winds can be modelled with m-CAK theory resulting in so called fast solution. However, those most rapidly rotating among them, especially when the stellar rotational speed is higher than $\sim 75\%$ of the critical rotational speed, can adopt a different solution called $Ω$-slow solution characterized by a dense and slow wind. Here, in this work we study the transition region of the solutions where the fast solution changes to the $Ω$-slow. Using both time-steady and time-dependent numerical codes, we study this transition region for different equatorial models of B-type stars. In all the cases, at certain range of rotational speeds, we found a region where the fast and $Ω$-slow solution can co-exist. We find that the type of solution obtained in this co-existence region depends heavily on the initial conditions of our models. We also test the stability of the solutions within the co-existence region by performing base density perturbations in the wind. We find that under certain conditions, the fast solution can switch to a $Ω$-slow solution, or vice versa. Such solution switching may be a possible contributor of material injected into the circumstellar environment of Be stars, without requiring rotational speeds near critical values. △ Less

Submitted 20 March, 2018; originally announced March 2018.

Comments: 12 pages, 14 figures, accepted for publication in MNRAS

Spectroscopic and photometric study of the eclipsing interacting binary V495 Centauri

Authors: J. A. Rosales, R. E. Mennickent, G. Djurašević, I. Araya, M. Curé

Abstract: Double Periodic Variables (DPV) are among the new enigmas of semi-detached eclipsing binaries. These are intermediate-mass binaries characterized by a long photometric period lasting on average 33 times the orbital period. We present a spectroscopic and photometric study of the DPV V495 Cen based on new high-resolution spectra and the ASAS V-band light curve. We have determined an improved orbital… ▽ More Double Periodic Variables (DPV) are among the new enigmas of semi-detached eclipsing binaries. These are intermediate-mass binaries characterized by a long photometric period lasting on average 33 times the orbital period. We present a spectroscopic and photometric study of the DPV V495 Cen based on new high-resolution spectra and the ASAS V-band light curve. We have determined an improved orbital period of $33.492 \pm 0.002$ d and a long period of 1283 d. We find a cool evolved star of $M_{2}=0.91\pm 0.2 M_{\odot}$, $T_{2}= 6000\pm 250 K$ and $R_{2}=19.3 \pm 0.5 R_{\odot}$ and a hot companion of $M_{1}= 5.76\pm 0.3 M_{\odot}$, $T_{1}=16960\pm 400 K$ and $R=4.5\pm0.2 R_{\odot}$. The mid-type B dwarf is surrounded by a concave and geometrically thick disc, of radial extension $R_{d}= 40.2\pm 1.3 R_{\odot}$ contributing $\sim$ 11 percent to the total luminosity of the system at the V band. The system is seen under inclination $84.\!\!^{\circ}8$ $\pm$ $0.\!\!^{\circ}6$ and it is at a distance $d= 2092 \pm 104.6$ pc. The light curve analysis suggests that the mass transfer stream impacts the external edge of the disc forming a hot region 11 % hotter than the surrounding disc material. The persistent $V<R$ asymmetry of the H$α$ emission suggests the presence of a wind and the detection of a secondary absorption component in He I lines indicates a possible wind origin in the hotspot region. △ Less

Submitted 29 January, 2018; v1 submitted 25 January, 2018; originally announced January 2018.

Journal ref: Monthly Notices of the Royal Astronomical Society, 476, 3039 (2018)

Stellar parameters and H$α$ line profile variability of Be stars in the BeSOS survey

Authors: C. Arcos, S. Kanaan, J. Chávez, L. Vanzi, I. Araya, M. Curé

Abstract: The Be phenomenon is present in about 20$\%$ of the B-type stars. Be stars show variability on a broad range of timescales, which in most cases is related to the presence of a circumstellar disk of variable size and structure. For this reason a time resolved survey is highly desirable in order to understand the mechanisms of disk formation which are still poorly understood. In addition, a complete… ▽ More The Be phenomenon is present in about 20$\%$ of the B-type stars. Be stars show variability on a broad range of timescales, which in most cases is related to the presence of a circumstellar disk of variable size and structure. For this reason a time resolved survey is highly desirable in order to understand the mechanisms of disk formation which are still poorly understood. In addition, a complete observational sample would improve the statistical significance of the study of the stellar and disk parameters. The "Be Stars Observation Survey" (BeSOS) is a survey containing reduced spectra obtained using the echelle spectrograph PUCHEROS with a spectral resolution of 17000 in the range of 4260-7300 $\textÅ$. BeSOS's main objective is to offer consistent spectroscopic and time resolved data obtained with one instrument. The user can download or plot the data and get the stellar parameters directly from the website. We also provide a star-by-star analysis based on photometric, spectroscopic and interferometric data as well as general information about the whole BeSOS sample. Recently, BeSOS led to the discovery of a new Be star HD 42167 and facilitated study of the V/R variation of HD 35165 and HD 120324, the steady disk of HD 110335 and the Be shell status of HD 127972. Optical spectra used in this work, as well as the derived stellar parameters are available online in \url{http://besos.ifa.uv.cl}. △ Less

Submitted 23 November, 2017; originally announced November 2017.

Comments: 16 pages, 13 Figures and 4 tables

Molecular gas in debris disks around young A-type stars

Authors: A. Moór, M. Curé, Á. Kóspál, P. Ábrahám, T. Csengeri, C. Eiroa, D. Gunawan, Th. Henning, A. M. Hughes, A. Juhász, N. Pawellek, M. Wyatt

Abstract: According to the current paradigm of circumstellar disk evolution, gas-rich primordial disks evolve into gas-poor debris disks compose of second-generation dust. To explore the transition between these phases, we searched for $^{12}$CO, $^{13}$CO, and C$^{18}$O emission in seven dust-rich debris disks around young A-type stars, using ALMA in Band 6. We discovered molecular gas in three debris disk… ▽ More According to the current paradigm of circumstellar disk evolution, gas-rich primordial disks evolve into gas-poor debris disks compose of second-generation dust. To explore the transition between these phases, we searched for $^{12}$CO, $^{13}$CO, and C$^{18}$O emission in seven dust-rich debris disks around young A-type stars, using ALMA in Band 6. We discovered molecular gas in three debris disks. In all these disks, the $^{12}$CO line was optically thick, highlighting the importance of less abundant molecules in reliable mass estimates. Supplementing our target list by literature data, we compiled a volume-limited sample of dust-rich debris disks around young A-type stars within 150 pc. We obtained a CO detection rate of 11/16 above a $^{12}$CO J=2$-$1 line luminosity threshold of $\sim 1.4 \times 10 ^4$ Jykms$^{-1}$pc$^2$ in the sample. This high incidence implies that the presence of CO gas in bright debris disks around young A-type stars is likely more the rule than the exception. Interestingly, dust-rich debris disks around young FG-type stars exhibit, with the same detectability threshold as for A-type stars, significantly lower gas incidence. While the transition from protoplanetary to debris phase is associated with a drop of dust content, our results exhibit a large spread in the CO mass in our debris sample, with peak values comparable to those in protoplanetary Herbig Ae disks. In the particularly CO-rich debris systems the gas may have primordial origin, characteristic of a hybrid disk. △ Less

Submitted 25 September, 2017; originally announced September 2017.

Comments: Accepted for publication in ApJ, 11 pages, 3 figures, 2 tables

Resolving the circumstellar environment of the Galactic B[e] supergiant star MWC 137 from large to small scales

Authors: Michaela Kraus, Tiina Liimets, Cristina E. Cappa, Lydia S. Cidale, Dieter H. Nickeler, Nicolas U. Duronea, Maria L. Arias, Diah S. Gunawan, Mary E. Oksala, Marcelo Borges Fernandes, Grigoris Maravelias, Michel Cure, Miguel Santander-Garcia

Abstract: The Galactic object MWC 137 was suggested to belong to the group of B[e] supergiants. However, with its large-scale optical bipolar ring nebula and the high velocity jet and knots, it is a rather atypical representative of this class. We performed multi-wavelength observations spreading from the optical to the radio regime. Based on optical imaging and long-slit spectroscopic data we found that th… ▽ More The Galactic object MWC 137 was suggested to belong to the group of B[e] supergiants. However, with its large-scale optical bipolar ring nebula and the high velocity jet and knots, it is a rather atypical representative of this class. We performed multi-wavelength observations spreading from the optical to the radio regime. Based on optical imaging and long-slit spectroscopic data we found that the northern parts of the large-scale nebula are predominantly blue-shifted, while the southern regions appear mostly red-shifted. We developed a geometrical model consisting of two double-cones. While various observational features can be approximated with such a scenario, the observed velocity pattern is more complex. Using near-infrared integral-field unit spectroscopy we studied the hot molecular gas in the close vicinity of the star. The emission from the hot CO gas arises in a small-scale disk revolving around the star on Keplerian orbits. While the disk itself cannot be spatially resolved, its emission is reflected by dust arranged in arc-like structures and clumps surrounding MWC 137 on small scales. In the radio regime we mapped the cold molecular gas in the outskirts of the optical nebula. We found that large amounts of cool molecular gas and warm dust embrace the optical nebula in the east, south and west. No cold gas or dust were detected in the north and north-western regions. Despite the new insights on the nebula kinematics gained from our studies, the real formation scenario of the large-scale nebula remains an open issue. △ Less

Submitted 19 September, 2017; originally announced September 2017.

Comments: 19 pages, 14 figures, Table 1 will be available in CDS only. Accepted for publication in AJ

$Ω$-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

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.

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)

TYC 8241 2652 1 and the case of the disappearing disk: no smoking gun yet

Authors: Hans Moritz Günther, Stefan Kraus, Carl Melis, Michel Curé, Tim Harries, Michael Ireland, Samer Kanaan, Katja Poppenhaeger, Aaron Rizzuto, David Rodriguez, Christian P. Schneider, Michael Sitko, Gerd Weigelt, Matthew Willson, Scott Wolk

Abstract: TYC 8241 2652 1 is a young star that showed a strong mid-infrared (mid-IR, 8-25 mu) excess in all observations before 2008 consistent with a dusty disk. Between 2008 and 2010 the mid-IR luminosity of this system dropped dramatically by at least a factor of 30 suggesting a loss of dust mass of an order of magnitude or more. We aim to constrain possible models including removal of disk material by s… ▽ More TYC 8241 2652 1 is a young star that showed a strong mid-infrared (mid-IR, 8-25 mu) excess in all observations before 2008 consistent with a dusty disk. Between 2008 and 2010 the mid-IR luminosity of this system dropped dramatically by at least a factor of 30 suggesting a loss of dust mass of an order of magnitude or more. We aim to constrain possible models including removal of disk material by stellar activity processes, the presence of a binary companion, or other explanations suggested in the literature. We present new X-ray observations, optical spectroscopy, near-IR interferometry, and mid-IR photometry of this system to constrain its parameters and further explore the cause of the dust mass loss. In X-rays TYC 8241 2652 1 has all properties expected from a young star: Its luminosity is in the saturation regime and the abundance pattern shows enhancement of O/Fe. The photospheric Ha line is filled with a weak emission feature, indicating chromospheric activity consistent with the observed level of coronal emission. Interferometry does not detect a companion and sets upper limits on the companion mass of 0.2, 0.35, 0.1 and 0.05 M_sun at projected physical separations of 0.1-4 AU,4-5 AU, 5-10 AU, and 10-30 AU, respectively (assuming a distance of 120.9 pc). Our mid-IR measurements, the first of the system since 2012, are consistent with the depleted dust level seen after 2009. The new data confirms that stellar activity is unlikely to destroy the dust in the disk and shows that scenarios where either TYC 8241 2652 1 heats the disk of a binary companion or a potential companion heats the disk of TYC 8241 2652 1 are unlikely. △ Less

Submitted 4 November, 2016; originally announced November 2016.

Comments: accepted by A&A

Journal ref: A&A 598, A82 (2017)

Clumpy molecular structures revolving the B[e] supergiant MWC 137

Authors: M. Kraus, L. S. Cidale, T. Liimets, C. E. Cappa, N. Duronea, D. S. Gunawan, M. E. Oksala, M. Santander-Garcia, M. L. Arias, D. H. Nickeler, G. Maravelias, M. Borges Fernandes, M. Cure

Abstract: The peculiar emission-line star MWC 137 with its extended optical nebula was recently classified as B[e] supergiant. To study the spatial distribution of its circumstellar molecular gas on small and large scales, we obtained near-infrared and radio observations using SINFONI and APEX, respectively. We find that the hot CO gas is arranged in moving clumpy ring and shell structures close to the star… ▽ More The peculiar emission-line star MWC 137 with its extended optical nebula was recently classified as B[e] supergiant. To study the spatial distribution of its circumstellar molecular gas on small and large scales, we obtained near-infrared and radio observations using SINFONI and APEX, respectively. We find that the hot CO gas is arranged in moving clumpy ring and shell structures close to the star, while a cold CO envelope is encircling the borders of the optical nebula from the south to the west. △ Less

Submitted 17 October, 2016; originally announced October 2016.

Comments: 2 pages, 1 figure, to appear in the proceedings of the workshop 'The B[e] Phenomenom: Forty Years of Studies'

Disk formation in oblate B[e] stars

Authors: I. Araya, C. Arcos, M. Curé

Abstract: We investigate the possible role of line-driven winds in the circumstellar envelope in B[e] stars, mainly the role of the $Ω$-slow wind solution, which is characterized by a slower terminal velocity and higher mass-loss rate, in comparison with the standard (m-CAK) wind solution. In this work, we assume two scenarios: 1) a spherically symmetric star and 2) a scenario that considers the oblate shap… ▽ More We investigate the possible role of line-driven winds in the circumstellar envelope in B[e] stars, mainly the role of the $Ω$-slow wind solution, which is characterized by a slower terminal velocity and higher mass-loss rate, in comparison with the standard (m-CAK) wind solution. In this work, we assume two scenarios: 1) a spherically symmetric star and 2) a scenario that considers the oblate shape, considering only the oblate correction factor. For certain values of the line force parameters (according to previous works), we obtain in both scenarios a density contrast $\gtrsim10^{2}$ between equatorial and polar densities, characterized for a fast polar wind and a slow and denser wind when the $Ω$-slow wind solution is obtained. All this properties are enhanced when the oblate correction factor is included in our calculations. △ Less

Submitted 13 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" (Prague, Czech Republic, 27 June - 1 July 2016)

A method to deconvolve stellar rotational velocities II

Authors: A. Christen, P. Escarate, M. Cure, D. F. Rial, J. Cassetti

Abstract: Knowing the distribution of stellar rotational velocities is essential for the understanding stellar evolution. Because we measure the projected rotational speed vsini, we need to solve an ill-posed problem given by a Fredholm integral of the first kind to recover the true rotational velocity distribution. After discretization of the Fredholm integral, we apply the Tikhonov regularization method t… ▽ More Knowing the distribution of stellar rotational velocities is essential for the understanding stellar evolution. Because we measure the projected rotational speed vsini, we need to solve an ill-posed problem given by a Fredholm integral of the first kind to recover the true rotational velocity distribution. After discretization of the Fredholm integral, we apply the Tikhonov regularization method to obtain directly the probability distribution function for stellar rotational velocities. We propose a simple and straightforward procedure to determine the Tikhonov parameter. We applied Monte Carlo simulations to prove that Tikhonov method is a consistent estimator and asymptotically unbiased. This method is applied to a sample of cluster stars. We obtain confidences intervals using bootstrap method. Our results are in good agreement with the one obtained using the Lucy method, in recovering the probability density distribution of rotational velocities. Furthermore, Lucy estimation lies inside our confidence interval. Tikhonov regularization is a very robust method that deconvolve the rotational velocity probability density function from a sample of vsini data straightforward without needing any convergence criteria. △ Less

Submitted 15 September, 2016; originally announced September 2016.

Comments: 15 pages, 7 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 595, A50 (2016)

Inhomogeneous molecular ring around the B[e] supergiant LHA 120-S 73

Authors: M. Kraus, L. S. Cidale, M. L. Arias, G. Maravelias, D. H. Nickeler, A. F. Torres, M. Borges Fernandes, A. Aret, M. Cure, R. Vallverdu, R. H. Barba

Abstract: We aim to improve our knowledge on the structure and dynamics of the circumstellar disk of the LMC B[e] supergiant LHA 120-S 73. High-resolution optical and near-IR spectroscopic data were obtained over a period of 16 and 7 years, respectively. The spectra cover the diagnostic emission lines from [CaII] and [OI], as well as the CO bands. These features trace the disk at different distances from th… ▽ More We aim to improve our knowledge on the structure and dynamics of the circumstellar disk of the LMC B[e] supergiant LHA 120-S 73. High-resolution optical and near-IR spectroscopic data were obtained over a period of 16 and 7 years, respectively. The spectra cover the diagnostic emission lines from [CaII] and [OI], as well as the CO bands. These features trace the disk at different distances from the star. We analyzed the kinematics of the individual emission regions by modeling their emission profiles. A low-resolution mid-infrared spectrum was obtained as well, which provides information on the composition of the dusty disk. All diagnostic emission features display double-peaked line profiles, which we interpret as due to Keplerian rotation. We find that LHA 120-S 73 is surrounded by at least four individual rings of material with alternating densities (or by a disk with strongly non-monotonic radial density distribution). Moreover, we find that the molecular ring must have gaps or at least strong density inhomogeneities, or in other words, a clumpy structure. The mid-infrared spectrum displays features of oxygen- and carbon-rich grain species, which indicates a long-lived, stable dusty disk. We cannot confirm the previously reported high value for the stellar rotation velocity. The line profile of HeI 5876 A is strongly variable in both width and shape and resembles of those seen in non-radially pulsating stars. A proper determination of the real underlying stellar rotation velocity is hence not possible. The existence of multiple stable and clumpy rings of alternating density recalls ring structures around planets. Although there is currently insufficient observational evidence, it is tempting to propose a scenario with one (or more) minor bodies or planets revolving around LHA 120-S 73 and stabilizing the ring system, in analogy to the shepherd moons in planetary systems. △ Less

Submitted 1 July, 2016; originally announced July 2016.

Comments: 14 pages, 13 figure, accepted for pulication in A&A

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

Improving prediction performance of stellar parameters using functional models

Authors: Sylvain Robbiano, Matthieu Saumard, Michel Curé

Abstract: This paper investigates the problem of prediction of stellar parameters, based on the star's electromagnetic spectrum. The knowledge of these parameters permits to infer on the evolutionary state of the star. From a statistical point of view, the spectra of different stars can be represented as functional data. Therefore, a two-step procedure decomposing the spectra in a functional basis combined… ▽ More This paper investigates the problem of prediction of stellar parameters, based on the star's electromagnetic spectrum. The knowledge of these parameters permits to infer on the evolutionary state of the star. From a statistical point of view, the spectra of different stars can be represented as functional data. Therefore, a two-step procedure decomposing the spectra in a functional basis combined with a regression method of prediction is proposed. We also use a bootstrap methodology to build prediction intervals for the stellar parameters. A practical application is also provided to illustrate the numerical performance of our approach. △ Less

Submitted 20 October, 2015; originally announced October 2015.

Discovery of molecular gas around HD 131835 in an APEX molecular line survey of bright debris disks

Authors: A. Moór, Th. Henning, A. Juhász, P. Ábrahám, Z. Balog, Á. Kóspál, I. Pascucci, Gy. M. Szabó, R. Vavrek, M. Curé, T. Csengeri, C. Grady, R. Güsten, Cs. Kiss

Abstract: Debris disks are considered to be gas-poor, but recent observations revealed molecular or atomic gas in several 10-40 Myr old systems. We used the APEX and IRAM 30m radiotelescopes to search for CO gas in 20 bright debris disks. In one case, around the 16 Myr old A-type star HD 131835, we discovered a new gas-bearing debris disk, where the CO 3-2 transition was successfully detected. No other indi… ▽ More Debris disks are considered to be gas-poor, but recent observations revealed molecular or atomic gas in several 10-40 Myr old systems. We used the APEX and IRAM 30m radiotelescopes to search for CO gas in 20 bright debris disks. In one case, around the 16 Myr old A-type star HD 131835, we discovered a new gas-bearing debris disk, where the CO 3-2 transition was successfully detected. No other individual system exhibited a measurable CO signal. Our Herschel Space Observatory far-infrared images of HD 131835 marginally resolved the disk both at 70 and 100$μ$m, with a characteristic radius of ~170 au. While in stellar properties HD 131835 resembles $β$ Pic, its dust disk properties are similar to those of the most massive young debris disks. With the detection of gas in HD 131835 the number of known debris disks with CO content has increased to four, all of them encircling young ($\leq$40 Myr) A-type stars. Based on statistics within 125 pc, we suggest that the presence of detectable amount of gas in the most massive debris disks around young A-type stars is a common phenomenon. Our current data cannot conclude on the origin of gas in HD 131835. If the gas is secondary, arising from the disruption of planetesimals, then HD 131835 is a comparably young and in terms of its disk more massive analogue of the $β$ Pic system. However, it is also possible that this system similarly to HD 21997 possesses a hybrid disk, where the gas material is predominantly primordial, while the dust grains are mostly derived from planetesimals. △ Less

Submitted 13 October, 2015; originally announced October 2015.

Comments: Accepted for publication in ApJ, 18 pages, 9 figures, 5 tables

Multitechnique testing of the viscous decretion disk model I. The stable and tenuous disk of the late-type Be star $β$ CMi

Authors: R. Klement, A. C. Carciofi, T. Rivinius, D. Panoglou, R. G. Vieira, J. E. Bjorkman, S. Štefl, C. Tycner, D. M. Faes, D. Korčáková, A. Müller, R. T. Zavala, M. Curé

Abstract: The viscous decretion disk (VDD) model is able to explain most of the currently observable properties of the circumstellar disks of Be stars. However, more stringent tests, focusing on reproducing multitechnique observations of individual targets via physical modeling, are needed to study the predictions of the VDD model under specific circumstances. In the case of nearby, bright Be star $β$ CMi,… ▽ More The viscous decretion disk (VDD) model is able to explain most of the currently observable properties of the circumstellar disks of Be stars. However, more stringent tests, focusing on reproducing multitechnique observations of individual targets via physical modeling, are needed to study the predictions of the VDD model under specific circumstances. In the case of nearby, bright Be star $β$ CMi, these circumstances are a very stable low-density disk and a late-type (B8Ve) central star. The aim is to test the VDD model thoroughly, exploiting the full diagnostic potential of individual types of observations, in particular, to constrain the poorly known structure of the outer disk if possible, and to test truncation effects caused by a possible binary companion using radio observations. We use the Monte Carlo radiative transfer code HDUST to produce model observables, which we compare with a very large set of multitechnique and multiwavelength observations that include ultraviolet and optical spectra, photometry covering the interval between optical and radio wavelengths, optical polarimetry, and optical and near-IR (spectro)interferometry. Due to the absence of large scale variability, data from different epochs can be combined into a single dataset. A parametric VDD model with radial density exponent of $n$ = 3.5, which is the canonical value for isothermal flaring disks, is found to explain observables typically formed in the inner disk, while observables originating in the more extended parts favor a shallower, $n$ = 3.0, density falloff. Modeling of radio observations allowed for the first determination of the physical extent of a Be disk (35$^{+10}_{-5}$ stellar radii), which might be caused by a binary companion. Finally, polarization data allowed for an indirect measurement of the rotation rate of the star, which was found to be $W \gtrsim 0.98$, i.e., very close to critical. △ Less

Submitted 5 October, 2015; originally announced October 2015.

Comments: 19 pages (35 including online material), 17 figures, 2 online figures, 2 online tables with data

Journal ref: A&A 584, A85 (2015)

Interplay between pulsations and mass loss in the blue supergiant 55 Cygnus = HD 198478

Authors: M. Kraus, M. Haucke, L. S. Cidale, R. O. J. Venero, D. H. Nickeler, P. Nemeth, E. Niemczura, S. Tomic, A. Aret, J. Kubat, B. Kubatova, M. E. Oksala, M. Cure, K. Kaminski, W. Dimitrov, M. Fagas, M. Polinska

Abstract: Blue supergiant stars are known to display photometric and spectroscopic variability that is suggested to be linked to stellar pulsations. Pulsational activity in massive stars strongly depends on the star's evolutionary stage and is assumed to be connected with mass-loss episodes, the appearance of macroturbulent line broadening, and the formation of clumps in the wind. To investigate a possible… ▽ More Blue supergiant stars are known to display photometric and spectroscopic variability that is suggested to be linked to stellar pulsations. Pulsational activity in massive stars strongly depends on the star's evolutionary stage and is assumed to be connected with mass-loss episodes, the appearance of macroturbulent line broadening, and the formation of clumps in the wind. To investigate a possible interplay between pulsations and mass-loss, we carried out an observational campaign of the supergiant 55 Cyg over a period of five years to search for photospheric activity and cyclic mass-loss variability in the stellar wind. We modeled the H, He I, Si II and Si III lines using the nonlocal thermal equilibrium atmosphere code FASTWIND and derived the photospheric and wind parameters. In addition, we searched for variability in the intensity and radial velocity of photospheric lines and performed a moment analysis of the line profiles to derive frequencies and amplitudes of the variations. The Halpha line varies with time in both intensity and shape, displaying various types of profiles: P Cygni, pure emission, almost complete absence, and double or multiple peaked. The star undergoes episodes of variable mass-loss rates that change by a factor of 1.7-2 on different timescales. We also observe changes in the ionization rate of Si II and determine a multiperiodic oscillation in the He I absorption lines, with periods ranging from a few hours to 22.5 days. We interpret the photospheric line variations in terms of oscillations in p-, g-, and strange modes. We suggest that these pulsations can lead to phases of enhanced mass loss. Furthermore, they can mislead the determination of the stellar rotation. We classify the star as a post-red supergiant, belonging to the group of alpha Cyg variables. △ Less

Submitted 7 July, 2015; originally announced July 2015.

Comments: 20 pages, 18 figures, 3 tables, accepted to Astronomy & Astrophysics

Journal ref: A&A 581, A75 (2015)

Fundamental parameters of the close interacting binary HD170582 and its luminous accretion disc

Authors: R. E. Mennickent, G. Djurašević, M. Cabezas, A. Cséki, J. Rosales G., E. Niemczura, I. Araya, M. Curé

Abstract: We present a spectroscopic and photometric study of the Double Period Variable HD170582. Based on the study of the ASAS V-band light curve we determine an improved orbital period of 16.87177 $\pm$ 0.02084 days and a long period of 587 days. We disentangled the light curve into an orbital part, determining ephemerides and revealing orbital ellipsoidal variability with unequal maxima, and a long cyc… ▽ More We present a spectroscopic and photometric study of the Double Period Variable HD170582. Based on the study of the ASAS V-band light curve we determine an improved orbital period of 16.87177 $\pm$ 0.02084 days and a long period of 587 days. We disentangled the light curve into an orbital part, determining ephemerides and revealing orbital ellipsoidal variability with unequal maxima, and a long cycle, showing quasi-sinusoidal changes with amplitude $ΔV$= 0.1 mag. Assuming synchronous rotation for the cool stellar component and semi-detached configuration we find a cool evolved star of $M_{2}$ = 1.9 $\pm$ 0.1 $M_{\odot}$, $T_{2}$ = 8000 $\pm$ 100 $K$ and $R_{2}$ = 15.6 $\pm$ 0.2 $R_{\odot}$, and an early B-type dwarf of $M_{1}$ = 9.0 $\pm$ 0.2 $M_{\odot}$. The B-type star is surrounded by a geometrically and optically thick accretion disc of radial extension 20.8 $\pm$ 0.3 $R_{\odot}$ contributing about 35% to the system luminosity at the $V$ band. Two extended regions located at opposite sides of the disc rim, and hotter than the disc by 67% and 46%, fit the light curve asymmetries. The system is seen under inclination 67.4 $\pm$ 0.4 degree and it is found at a distance of 238 $\pm$ 10 pc. Specially interesting is the double line nature of HeI 5875; two absorption components move in anti-phase during the orbital cycle; they can be associated with the shock regions revealed by the photometry. The radial velocity of one of the HeI 5875 components closely follows the donor radial velocity, suggesting that the line is formed in a wind emerging near the stream-disc interacting region. △ Less

Submitted 19 January, 2015; originally announced January 2015.

Comments: 16 pages, 11 tables and 11 figures, accepted for publication in MNRAS

A method to deconvolve mass ratio distribution from binary stars

Authors: Michel Cure, Diego F. Rial, Alejandra Christen, Julia Cassetti, Henri M. J. Boffin

Abstract: To better understand the evolution of stars in binary systems as well as to constrain the formation of binary stars, it is important to know the binary mass-ratio distribution. However, in most cases, i.e. for single-lined spectroscopic binaries, the mass ratio cannot be measured directly but only derived as the convolution of a function that depends on the mass ratio and the unknown inclination a… ▽ More To better understand the evolution of stars in binary systems as well as to constrain the formation of binary stars, it is important to know the binary mass-ratio distribution. However, in most cases, i.e. for single-lined spectroscopic binaries, the mass ratio cannot be measured directly but only derived as the convolution of a function that depends on the mass ratio and the unknown inclination angle of the orbit on the plane of the sky. We extend our previous method to deconvolve this inverse problem (Cure et al. 2014), i.e., we obtain as an integral the cumulative distribution function (CDF) for the mass ratio distribution. After a suitable transformation of variables it turns out that this problem is the same as the one for rotational velocities $v \sin i$, allowing a close analytic formulation for the CDF. We then apply our method to two real datasets: a sample of Am stars binary systems, and a sample of massive spectroscopic binaries in the Cyg OB2 Association.} {We are able to reproduce the previous results of Boffin (2010) for the sample of Am stars, while we show that the mass ratio distribution of massive stars shows an excess of small mass ratio systems, contrarily to what was claimed by Kobulnicky et al. (2014). Our method proves very robust and deconvolves the distribution from a sample in just a single step. △ Less

Submitted 1 December, 2014; originally announced December 2014.

Comments: Accepted for publication in A&A. 9 pages, 3 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

Analytical Solutions for Radiation-Driven Winds in Massive Stars. I: The Fast Regime

Authors: Ignacio Araya, Michel Cure, Lydia S. Cidale

Abstract: Accurate massloss rate estimates are crucial keys to study wind properties of massive stars and test different evolutionary scenarios. From a theoretical point of view, this implies to solve a complex set of differential equations in which the radiation field and the hydrodynamics are strong coupled. The use of analytical expression to represent the radiation force and the solution of the equation… ▽ More Accurate massloss rate estimates are crucial keys to study wind properties of massive stars and test different evolutionary scenarios. From a theoretical point of view, this implies to solve a complex set of differential equations in which the radiation field and the hydrodynamics are strong coupled. The use of analytical expression to represent the radiation force and the solution of the equation of motion have many advantages over numerical integrations. Therefore, in this work, we present an analytical expression as solution of the equation of motion for radiation driven winds, in terms of the force multipliers parameters. This analytical expression is obtained by employing the line acceleration expression given by Villata (1992) and the methodology proposed by Mueller and Vink (2008). On the other hand, we find useful relationships to determine the parameters for the line acceleration given by Mueller and Vink (2008) in terms of the force multiplier parameters. △ Less

Submitted 3 November, 2014; originally announced November 2014.

Comments: Accepted for publication in ApJ (2014 November 1) 27 pages 9 figs

Journal ref: ApJ, 795 81 (2014)

The yellow hypergiant HR 5171 A: Resolving a massive interacting binary in the common envelope phase

Authors: O. Chesneau, A. Meilland, E. Chapellier, F. Millour, A. M. Van Genderen, Y. Nazé, N. Smith, A. Spang, J. V. Smoker, L. Dessart, S. Kanaan, Ph. Bendjoya, M. W. Feast, J. H. Groh, A. Lobel, N. Nardetto, S. Otero, R. D. Oudmaijer, A. G. Tekola, P. A. Whitelock, C. Arcos, M. Curé, L. Vanzi

Abstract: We initiated long-term optical interferometry monitoring of the diameters of unstable yellow hypergiants (YHG) with the goal of detecting both the long-term evolution of their radius and shorter term formation related to large mass-loss events. We observed HR5171 A with AMBER/VLTI. We also examined archival photometric data in the visual and near-IR spanning more than 60 years, as well as sparse s… ▽ More We initiated long-term optical interferometry monitoring of the diameters of unstable yellow hypergiants (YHG) with the goal of detecting both the long-term evolution of their radius and shorter term formation related to large mass-loss events. We observed HR5171 A with AMBER/VLTI. We also examined archival photometric data in the visual and near-IR spanning more than 60 years, as well as sparse spectroscopic data. HR5171A exhibits a complex appearance. Our AMBER data reveal a surprisingly large star for a YHG R*=1315+/-260Rsun\ (~6.1AU) at the distance of 3.6+/-0.5kpc. The source is surrounded by an extended nebulosity, and these data also show a large level of asymmetry in the brightness distribution of the system, which we attribute to a newly discovered companion star located in front of the primary star. The companion's signature is also detected in the visual photometry, which indicates an orbital period of Porb=1304+/-6d. Modeling the light curve with the NIGHTFALL program provides clear evidence that the system is a contact or possibly over-contact eclipsing binary. A total current system mass of 39^{+40}_{-22} solar mass and a high mass ratio q>10 is inferred for the system. The low-mass companion of HR5171 A is very close to the primary star that is embedded within its dense wind. Tight constraints on the inclination and vsini of the primary are lacking, which prevents us from determining its influence precisely on the mass-loss phenomenon, but the system is probably experiencing a wind Roche-Lobe overflow. Depending on the amount of angular momentum that can be transferred to the stellar envelope, HR5171 A may become a fast-rotating B[e]/Luminous Blue Variable (LBV)/Wolf-Rayet star. In any case, HR5171 A highlights the possible importance of binaries for interpreting the unstable YHGs and for massive star evolution in general. △ Less

Submitted 14 January, 2014; v1 submitted 12 January, 2014; originally announced January 2014.

Comments: Accepted in Astronomy & Astrophysics (Jan 2014)

Journal ref: A&A, 563, A71 (23p, 2014)

A method to deconvolve stellar rotational velocities

Authors: Michel Cure, Diego F. Rial, Alejandra Christen, Julia Cassetti

Abstract: Rotational speed is an important physical parameter of stars and knowing the distribution of stellar rotational velocities is essential for the understanding stellar evolution. However, it cannot be measured directly but the convolution of the rotational speed and the sine of the inclination angle, $v \sin i$. We developed a method to deconvolve this inverse problem and obtain the cumulative distr… ▽ More Rotational speed is an important physical parameter of stars and knowing the distribution of stellar rotational velocities is essential for the understanding stellar evolution. However, it cannot be measured directly but the convolution of the rotational speed and the sine of the inclination angle, $v \sin i$. We developed a method to deconvolve this inverse problem and obtain the cumulative distribution function (CDF) for stellar rotational velocities extending the work of Chandrasekhar & Münch (1950). This method is applied a) to theoretical synthetic data recovering the original velocity distribution with very small error; b) to a sample of about 12.000 field main--sequence stars, corroborating that the velocity distribution function is non--Maxwellian, but is better described by distributions based on the concept of maximum entropy, such as Tsallis or Kaniadakis distribution functions. This is a very robust and novel method that deconvolve the rotational velocity cumulative distribution function from a sample of $v \sin i$ data in just one single step without needing any convergence criteria. △ Less

Submitted 19 March, 2014; v1 submitted 6 January, 2014; originally announced January 2014.

Comments: Accepted in A&A

Radial velocity variations in the young eruptive star EX Lup

Authors: Á. Kóspál, M. Mohler-Fischer, A. Sicilia-Aguilar, P. Ábrahám, M. Curé, Th. Henning, Cs. Kiss, R. Launhardt, A. Moór, A. Müller

Abstract: EX Lup-type objects (EXors) are low-mass pre-main sequence objects characterized by outbursts attributed to highly enhanced disk accretion. The trigger mechanism of EXor outbursts is still debated. One theory requires a close (sub)stellar companion that perturbs the inner disk and triggers the onset of the outburst. Here, we study the radial velocity (RV) variations of EX Lup, the prototype of EXo… ▽ More EX Lup-type objects (EXors) are low-mass pre-main sequence objects characterized by outbursts attributed to highly enhanced disk accretion. The trigger mechanism of EXor outbursts is still debated. One theory requires a close (sub)stellar companion that perturbs the inner disk and triggers the onset of the outburst. Here, we study the radial velocity (RV) variations of EX Lup, the prototype of EXors. We conducted a 5-year RV survey with HARPS and FEROS. We analyzed the activity of EX Lup by checking the bisector and several other activity indicators. We complemented the RV data with photometric monitoring to look for signatures of activity or varying accretion. We found that the RV of EX Lup is periodic, with stable period, semi-amplitude, and phase over at least four years of observations. This period is not present in any of the activity indicators. However, the RV of narrow metallic emission lines suggest the same period, but with an anti-correlating phase. The observed absorption line RVs can be fitted with a Keplerian solution around a 0.6 M_Sun central star with M2 sin i = (14.7 +/- 0.7) M_Jup and eccentricity of e = 0.24. Alternatively, we attempted to model the observations with a cold/hot stellar spot as well, but the spot parameters needed to reproduce the RV semi-amplitude are in contradiction with the photometric variability, making the spot scenario unlikely. We discuss two possibilities to explain the RV data: a geometry with two accretion columns rotating with the star, and a single accretion flow synchronized with the orbital motion of the hypothetical companion. In the companion scenario, the companion's mass would fall into the brown dwarf desert, which, together with the unusually small separation would make EX Lup a unique binary system, with interesting implications on the physical mechanisms responsible for triggering the outburst. △ Less

Submitted 17 November, 2013; originally announced November 2013.

Comments: 10 pages, 9 figures, 3 tables, accepted for publication in A&A