Search

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., Comment: Accepted for publication in Astronomy & Astrophysics

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., Comment: Published in Galaxies Journal (33 pages, 18 figures)

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., Comment: Accepted for publication in Astronomy & Astrophysics journal

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., Comment: Accepted in A&A

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., Comment: Accepted for publication in Astronomy & Astrophysics, in 17 February 2022 In this updated version, we corrected information about funding data in Acknowledgements section

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 $\lambda$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., Comment: 14 pages, 10 figures, MNRAS, accepted

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 $\beta$-law. We apply the Lambert-procedure to three O supergiant stars ($\zeta$-Puppis, HD~165763 and $\alpha$-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 $\beta$-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., Comment: Published in the Astrophysical Journal (ApJ), in 13 October 2021

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., Comment: 63 pages, 4+13 Figures, 6 Tables, published at Science

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., Comment: 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

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)

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.

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., Comment: 40 pages, 12 figures, to appear in the Astronomical Journal

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$\alpha$ and R-band continuum. Our SMA sub-millimetre interferometry in 880 $\mu$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$\alpha$ 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., Comment: 17 pages, 14 figures

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)., Comment: 19 pages, 6 figures

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] $\lambda$5577, [OI] $\lambda\lambda$6300,6363, and [CaII] $\lambda\lambda$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., Comment: 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

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.

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 $\Omega$-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 $\Omega$-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 $\Omega$-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 $\Omega$-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., Comment: 12 pages, 14 figures, accepted for publication in MNRAS

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., Comment: 19 pages, 14 figures, Table 1 will be available in CDS only. Accepted for publication in AJ

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., Comment: accepted by A&A

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.

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

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., Comment: Accepted for publication in ApJ (2014 November 1) 27 pages 9 figs

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\"unch (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., Comment: Accepted in A&A

During the post-main sequence evolution massive stars pass through several short-lived phases, in which they experience enhanced mass loss in the form of clumped winds and mass ejection events of unclear origin. The discovery that stars populating the blue, luminous part of the Hertzsprung-Russell diagram can pulsate hence suggests that stellar pulsations might influence or trigger enhanced mass loss and eruptions. We present recent results for two objects in different phases: a B[e] star at the end of the main sequence and a B-type supergiant., Comment: 4 pages, 3 figures, to appear in IAU Symp. 301: Precision Asteroseismology, eds. W. Chaplin, J. Guzik, G. Handler, A. Pigulski

Thanks to the high spatial resolution provided by long baseline interferometry, it is possible to understand the complex circumstellar geometry around stars with the B[e] phenomenon. These stars are composed by objects in different evolutionary stages, like high- and low-mass evolved stars, intermediate-mass pre-main sequence stars and symbiotic objects. However, up to now more than 50% of the confirmed B[e] stars are not well classified, being called unclassified B[e] stars. From instruments like VLTI/AMBER and VLTI/MIDI, we have identified the presence of gaseous and dusty circumstellar disks, which have provided us with some hints related to the nature of these objects. Here, we show our results for two galactic stars with the B[e] phenomenon, HD 50138 and CPD-529243, based on interferometric measurements., Comment: Proceedings of SF2A conference in Nice 2012 - complete list of authors

We calculated the influence of the limb-darkened finite disk correction factor in the theory of radiation-driven winds from massive stars. We solved the 1-D m-CAK hydrodynamical equation of rotating radiation-driven winds for all three known solutions, i.e., fast, \Omega-slow and \delta-slow. We found that for the fast solution, the mass loss rate is increased by a factor \sim 10%, while the terminal velocity is reduced about 10%, when compared with the solution using a finite disk correction factor from a uniformly bright star. For the other two slow solutions the changes are almost negligible. Although, we found that the limb darkening has no effects on the wind momentum luminosity relationship, it would affect the calculation of synthetic line profiles and the derivation of accurate wind parameters., Comment: Accepted for publication in ApJ. 19 pages, 6 figures

V393 Scorpii is a Double Periodic Variable characterized by a relatively stable non-orbital photometric cycle of 253 days. Mennickent et al. argue for the presence of a massive optically thick disc around the more massive B-type component and describe the evolutionary stage of the system. In this paper we analyze the behavior of the main spectroscopic optical lines during the long non-orbital photometric cycle. We study the radial velocity of the donor determining their orbital elements and find a small but significant orbital eccentricity (e = 0.04). The donor spectral features are modeled and removed from the spectrum at every observing epoch using the light-curve model given by Mennickent et al. We find that the line emission is larger during eclipses and mostly comes from a bipolar wind. We find that the long cycle is explained in terms of a modulation of the wind strength; the wind has a larger line and continuum emissivity on the high state. We report the discovery of highly variable chromospheric emission in the donor, as revealed by Doppler maps of the emission lines MgII 4481 and CI 6588. We discuss notable and some novel spectroscopic features like discrete absorption components, especially visible at blue-depressed OI 7773 absorption wings during the second half-cycle, Balmer double emission with V/R-curves showing "Z-type" and "S-type" excursions around secondary and main eclipse, respectively, and H_beta emission wings extending up to +- 2000 km/s. We discuss possible causes for these phenomena and for their modulations with the long cycle., Comment: 19 pages, 22 figures, accepted for publication in MNRAS

We present thirteen epochs of near-infrared (0.8-5 micron) spectroscopic observations of the pre-transitional, "gapped" disk system in SAO 206462 (=HD 135344B). In all, six gas emission lines (including Br gamma, Pa beta, and the 0.8446 micron line of O I) along with continuum measurements made near the standard J, H, K, and L photometric bands were measured. A mass accretion rate of approximately 2 x 10^-8 solar masses per year was derived from the Br gamma and Pa beta lines. However, the fluxes of these lines varied by a factor of over two during the course of a few months. The continuum also varied, but by only ~30%, and even decreased at a time when the gas emission was increasing. The H I line at 1.083 microns was also found to vary in a manner inconsistent with that of either the hydrogen lines or the dust. Both the gas and dust variabilities indicate significant changes in the region of the inner gas and the inner dust belt that may be common to many young disk systems. If planets are responsible for defining the inner edge of the gap, they could interact with the material on time scales commensurate with what is observed for the variations in the dust, while other disk instabilities (thermal, magnetorotational) would operate there on longer time scales than we observe for the inner dust belt. For SAO 206462, the orbital period would likely be 1-3 years. If the changes are being induced in the disk material closer to the star than the gap, a variety of mechanisms (disk instabilities, interactions via planets) might be responsible for the changes seen. The He I feature is most likely due to a wind whose orientation changes with respect to the observer on time scales of a day or less. To further constrain the origin of the gas and dust emission will require multiple spectroscopic and interferometric observations on both shorter and longer time scales that have been sampled so far., Comment: 42 pages, 10 figures

The line driven- and rotation modulated-wind theory predicts an alternative slow solution, besides from the standard m-CAK solution, when the rotational velocity is close to the critical velocity. We study the behaviour of the winds of A-type supergiants (Asg) and show that under particular conditions, e.g., when the \delta line-force parameter is about 0.25, the slow solution could exist over the whole star, even for the cases when the rotational speed is slow or zero. We discuss density and velocity profiles as well as possible observational counterparts., Comment: 2 pages, to appear in the proceedings of the IAUS 272 on "Active OB stars: structure, evolution, mass loss and critical limits" (Paris, July 19-23, 2010), Cambridge University Press. Editors C. Neiner, G. Wade, G. Meynet and G. Peters

We have incorporated the oblate distortion of the shape of the star due to the stellar rotation, which modifies the finite disk correction factor (f_D) in the m-CAK hydrodynamical model. We implement a simplified version for the f_D allowing us to solve numerically the non-linear m- CAK momentum equation.We solve this model for a classical Be star in the polar and equatorial directions. The star's oblateness modifies the polar wind, which is now much faster than the spherical one, mainly because the wind receives radiation from a larger (than the spherical) stellar surface. In the equatorial direction we obtain slow solutions, which are even slower and denser than the spherical ones. For the case when the stellar rotational velocity is about the critical velocity, the most remarkable result of our calculations is that the density contrast between the equatorial density and the polar one, is about 100. This result could explain a long-standing problem on Be stars., Comment: 2 pages, to appear in the proceedings of the IAUS 272 on "Active OB stars: structure, evolution, mass loss and critical limits" (Paris, July 19-23, 2010), Cambridge University Press. Editors C. Neiner, G. Wade, G. Meynet and G. Peters

We present a general method for solving the non--linear differential equation of monotonically increasing steady--state radiation driven winds. We graphically identify all the singular points before transforming the momentum equation to a system of differential equations with all the gradients explicitly give. This permits a topological classification of all singular points and to calculate the maximum and minimum mass--loss of the wind. We use our method to analyse for the first time the topology of the non--rotating frozen in ionisation m--CAK wind, with the inclusion of the finite disk correction factor and find up to 4 singular points, three of the x--type and one attractor--type. The only singular point (and solution passing through) that satisfies the boundary condition at the stellar surface is the standard m--CAK singular point., Comment: Accepted for publication in Astronomische Nachrichten / Astronomical Notes

The effects of rapid rotation and bi--stability upon the density contrast between the equatorial and polar directions of a B[e] supergiant are investigated. Based on a new slow solution for different high rotational radiation--driven winds and the fact that bi--stability allows a change in the line--force parameters ($\alpha$, $k$, and $\delta$), the equatorial densities are about $10^2$--$10^3$ times higher than the polar ones. These values are in qualitative agreement with the observations. This calculation also permits to obtain the aperture angle of the disk., Comment: 6 pages, 2 figures. To appear in the proceedings "Stars with the B[e] phenomenon"

The topological analysis from Bjorkman (1995) for the standard model that describes the winds from hot stars by Castor, Abbott & Klein (1975) has been extended to include the effect of stellar rotation and changes in the ionization of the wind. The differential equation for the momentum of the wind is non--linear and transcendental for the velocity gradient. Due to this non--linearity the number of solutions that this equation possess is not known. After a change of variables and the introduction of a new physically meaningless independent variable, we manage to replace the non--linear momentum differential equation by a system of differential equations where all the derivatives are {\it{explicitely}} given. We then use this system of equations to study the topology of the rotating--CAK model. For the particular case when the wind is frozen in ionization ($\delta=0$) only one physical solution is found, the standard CAK solution, with a X--type singular point. For the more general case ($\delta \neq 0$), besides the standard CAK singular point, we find a second singular point which is focal--type (or attractor). We find also, that the wind does not adopt the maximal mass--loss rate but almost the minimal., Comment: 11 pages, Astronomy and Astrophysics Accepted

The theory of radiation driven wind including stellar rotation is re-examined. After a suitable change of variables, a new equation for the mass loss rate is derived analytically. The solution of this equation remains within 1% confidence when compared with numerical solutions. Also, a non-linear equation for the position of the critical (singular) point is obtained. This equation shows the existence of an additional critical point, besides the standard m--CAK critical point. For a stellar rotation velocity larger than aprox. 0.7 - 0.8 V_{breakup}, there exists only one critical point, located away from the star's surface. Numerical solutions crossing through this new critical point, are attained. In these cases, the wind has a very low terminal velocity and therefore a higher density wind. Disk formation in Be stars is discussed in the frame of this new line driven stellar wind solution., Comment: 40 pages, 15 figures, Accepted for publication in ApJ