In young multiple stellar systems the gravitational influence of the stars shapes the circumstellar disk, controlling accretion and the material available for planet formation. Our observations of the triple system GW Orionis provide direct evidence for disk tearing, where the gravitational torques due to misaligned disk/orbit planes warp the disk and break the inner disk into precessing rings. We image an eccentric ring that is misaligned with respect to the orbital planes and outer disk, and casts shadows on the strongly warped intermediate disk. Our ring/warp geometry constraints and the fully characterized perturber orbits make the system a potential Rosetta Stone for studying disk hydrodynamics. The ring might offer suitable conditions for planet formation, providing a mechanism for forming wide-separation planets on highly oblique orbits.

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 ($\lambda \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.

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.

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.

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.

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.

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$\gamma$ 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$\gamma$ 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.

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).

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.

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.

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$\alpha$ 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.

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{\AA}$. 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}.

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.

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 $\Omega$-slow, dense wind in the equatorial regions. We use the equatorial mass densities to explore H$\alpha$ 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$\alpha$ 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$, $ \, \alpha \sim 0.6$ and $\, \delta \gtrsim 0.1$, and c) the contribution of the fast wind component to the H$\alpha$ emission line profile is negligible; therefore, the line profiles arise mainly from the equatorial disks of Be stars.

The circumstellar disk density distributions for a sample of 63 Be southern stars from the BeSOS survey were found by modelling their H$\alpha$ 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.

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)

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.

We investigate the possible role of line-driven winds in the circumstellar envelope in B[e] stars, mainly the role of the $\Omega$-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 $\Omega$-slow wind solution is obtained. All this properties are enhanced when the oblate correction factor is included in our calculations.

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$\mu$m, with a characteristic radius of ~170 au. While in stellar properties HD 131835 resembles $\beta$ 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 $\beta$ 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.

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 $\beta$ 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.

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 $\Delta 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.

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.

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.

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.

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.

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.

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.

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.