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

We present a detailed visible and near-IR spectro-interferometric analysis of the Be-shell star $\omicron$ Aquarii from quasi-contemporaneous CHARA/VEGA and VLTI/AMBER observations. We measured the stellar radius of $\omicron$ Aquarii as 4.0 $\pm$ 0.3 $\mathrm{R_{\odot}}$. We constrained the disk geometry and kinematics using a kinematic model and a MCMC fitting procedure. The disk sizes in H$\alpha$ and Br$\gamma$ were found to be similar, at $\sim$10-12 $\mathrm{D_{\star}}$, which is uncommon since most results for Be stars show a larger extension in H$\alpha$ than in Br$\gamma$. We found that the inclination angle $i$ derived from H$\alpha$ is significantly lower ($\sim$15 deg) than the one derived from Br$\gamma$. The disk kinematics were found to be near to the Keplerian rotation in Br$\gamma$, but not in H$\alpha$. After analyzing all our data using a grid of HDUST models (BeAtlas), we found a common physical description for the disk in both lines: $\Sigma_{0}$ = 0.12 g cm\textsuperscript{-2} and $m$ = 3.0. The stellar rotational rate was found to be very close ($\sim$96\%) to the critical value. Our analysis of multi-epoch H$\alpha$ profiles and imaging polarimetry indicates that the disk has been stable for at least 20 years. Compared to Br$\gamma$, the data in H$\alpha$ shows a substantially different picture that cannot fully be understood using the current physical models of Be star disks. $\omicron$ Aquarii presents a stable disk, but the measured $m$ is lower than the standard value in the VDD model for steady-state. Such long-term stability can be understood in terms of the high rotational rate for this star, the rate being a main source for the mass injection in the disk. Our results on the stellar rotation and disk stability are consistent with results in the literature showing that late-type Be stars are more likely to be fast rotators and have stable disks.

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

An important tool for the development of the next generation of extremely large telescopes (ELTs) is the Systems Engineering (SE). GMACS is the first-generation multi-object spectrograph working at visible wavelengths for the Giant Magellan Telescope (GMT). The aim is to discuss the application of SE in ground-based astronomy for multi-object spectrographs. For this, it is presented the SE of the GMACS spectrograph, currently on its Conceptual Design phase. SE provide means to assist the management of complex projects, and in the case of GMACS, to ensure its success when in operation, maximizing the scientific potential of GMT.

We present the concept of a new Fabry-Perot instrument called BTFI-2, which is based on the design of another Brazilian instrument for the SOAR Telescope, the Brazilian Tunable Filter Imager (BTFI). BTFI-2 is designed to be mounted on the visitor port of the SOAR Adaptive Module (SAM) facility, on the SOAR telescope, at Cerro Pach\'on, Chile. This optical Fabry-Perot instrument will have a field of view of 3 x 3 arcmin, with 0.12 arcsec per pixel and spectral resolutions of 4500 and 12000, at H-alpha, dictated by the two ICOS Fabry-Perot devices available. The instrument will be unique for the study of centers of normal, interacting and active galaxies and the intergalactic medium, whenever spatial resolution over a large area is required. BTFI-2 will combine the best features of two previous instruments, SAM-FP and BTFI: it will use an Electron Multiplication detector for low and fast scanning, it will be built with the possibility of using a new Fabry-Perot etalon which provides a range of resolutions and it will be light enough to work attached to SAM, and hence the output data cubes will be GLAO-corrected.

The use of Python is noticeably growing among the scientific community, and Astronomy is not an exception. The power of Python consists of being an extremely versatile high-level language, easy to program that combines both traditional programming and data reduction and analysis tools. Here I make a brief introduction to Python, mentioning a few programming practices implemented in the language and some of its useful features on the process of data manipulation. I cover in a little more detail the standard scientific libraries (NumPy and SciPy) for data handling, the graphical library (Matplotlib), and tools for specific use in astronomy (PyFITS and PyRAF). Good programming practices and how they are implemented at the language are also viewed. Python resources and references are mentioned through- out the text for those who wish to go deeper and make use of the power of the language.

Adaptive Optics (AO) is an innovative technique that substantially improves the optical performance of ground-based telescopes. The SOAR Adaptive Module (SAM) is a laser-assisted AO instrument, designed to compensate ground-layer atmospheric turbulence in near-IR and visible wavelengths over a large Field of View. Here we detail our proposal to upgrade SAM, dubbed SAMplus, that is focused on enhancing its performance in visible wavelengths and increasing the instrument reliability. As an illustration, for a seeing of 0.62 arcsec at 500 nm and a typical turbulence profile, current SAM improves the PSF FWHM to 0.40 arcsec, and with the upgrade we expect to deliver images with a FWHM of $\approx0.34$ arcsec -- up to 0.23 arcsec FWHM PSF under good seeing conditions. Such capabilities will be fully integrated with the latest SAM instruments, putting SOAR in an unique position as observatory facility.

An important tool for the development of the next generation of extremely large telescopes (ELTs) is a robust Systems Engineering (SE) methodology. GMACS is a first-generation multi-object spectrograph that will work at visible wavelengths on the Giant Magellan Telescope (GMT). In this paper, we discuss the application of SE to the design of next-generation instruments for ground-based astronomy and present the ongoing development of SE products for the GMACS spectrograph, currently in its Conceptual Design phase. SE provides the means to assist in the management of complex projects, and in the case of GMACS, to ensure its operational success, maximizing the scientific potential of GMT.

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

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

Aims. The X-shooter archive of several thousand telluric star spectra was skimmed for Be and Be-shell stars to derive the stellar fundamental parameters and statistical properties, in particular for the less investigated late type Be stars, and the extension of the Be phenomenon into early A stars. Methods. An adapted version of the BCD method is used, utilizing the Balmer discontinuity parameters to determine effective temperature and surface gravity. This method is optimally suited for late B stars. The projected rotational velocity was obtained by profile fitting to the Mg ii lines of the targets, and the spectra were inspected visually for the presence of peculiar features such as the infrared Ca ii triplet or the presence of a double Balmer discontinuity. The Balmer line equivalent widths were measured, but due to uncertainties in determining the photospheric contribution are useful only in a subsample of Be stars for determining the pure emission contribution. Results. A total of 78 Be stars, mostly late type ones, were identified in the X-shooter telluric standard star archive, out of which 48 had not been reported before. The general trend of late type Be stars having more tenuous disks and being less variable than early type ones is confirmed. The relatively large number (48) of relatively bright (V > 8.5) additional Be stars casts some doubt on the statistics of late type Be stars; they are more common than currently thought: The Be/B star fraction may not strongly depend on spectral subtype.

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

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

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)

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

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

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

Optical long baseline interferometry was recently established as a technique capable of resolving stars and their circumstellar environments at the milliarcsecond (mas) resolution level. This high-resolution opens an entire new window to the study of astrophysical systems, providing information inaccessible by other techniques. Astrophysical disks are observed in a wide variety of systems, from galaxies up to planetary rings, commonly sharing similar physical processes. Two particular disk like systems are studied in the thesis: (i) B He-rich stars that exhibits magnetic fields in order of kG and that trap their winds in structures called magnetospheres; and (ii) Be stars, fast rotating stars that create circumstellar viscous disks. This study uses the interferometric technique to investigate both the photosphere proper and the circumstellar environment of these stars. The objective is to combine interferometry with other observational techniques (such as spectroscopy and polarimetry) to perform a complete and well-constrained physical description of these systems. This description is accompanied by radiative transfer models performed by the HDUST code. Keywords: stars: individual (Achernar), stars: fundamental parameters, techniques: interferometric, circumstellar matter, stars: emission-line, Be, stars: magnetic field

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.

Context. The circumstellar disks ejected by many rapidly rotating B stars (so-called Be stars) offer the rare opportunity of studying the structure and dynamics of gaseous disks at high spectral as well as angular resolution. Aims. This paper explores a newly identified effect in spectro-interferometric phase that can be used for probing the inner regions of gaseous edge-on disks on a scale of a few stellar radii. Methods. The origin of this effect (dubbed central quasi-emission phase signature, CQE-PS) lies in the velocity-dependent line absorption of photospheric radiation by the circumstellar disk. At high spectral and marginal interferometric resolution, photocenter displacements between star and isovelocity regions in the Keplerian disk reveal themselves through small interferometric phase shifts. To investigate the diagnostic potential of this effect, a series of models are presented, based on detailed radiative transfer calculations in a viscous decretion disk. Results. Amplitude and detailed shape of the CQE-PS depend sensitively on disk density and size and on the radial distribution of the material with characteristic shapes in differential phase diagrams. In addition, useful lower limits to the angular size of the central stars can be derived even when the system is almost unresolved. Conclusions. The full power of this diagnostic tool can be expected if it can be applied to observations over a full life-cycle of a disk from first ejection through final dispersal, over a full cycle of disk oscillations, or over a full orbital period in a binary system.

Some massive stars possess strong magnetic fields that confine plasma in the circumstellar environment. These \textit{magnetospheres} have been studied spectroscopically, photometrically and, more recently, interferometrically. Here we report on the first firm detection of a magnetosphere in continuum linear polarization, as a result of monitoring of $\sigma$\,Ori\,E at the Pico dos Dias Observatory. {The non-zero intrinsic polarization indicates an asymmetric structure, whose minor elongation axis is oriented $150\fdg0$ east of the celestial north.} A modulation of the polarization was observed, with a period of half of the rotation period, which supports the theoretical prediction of the presence of two diametrally opposed, co-rotating blobs of gas. A phase lag of -0.085 was detected between the polarization minimum and the primary minimum of the light curve, suggestive of a complex shape of the plasma clouds. We present a preliminary analysis of the data with the Rigidly Rotating Magnetosphere model, which could not reproduce simultaneously the photometric and polarimetric data. A toy model comprising two spherical co-rotating blobs {joined by a thin disk} proved more successful in reproducing the polarization modulation. {With this model we were able to determine that the total scattering mass of the thin disk is similar to the mass of the blobs ($2M_{\rm b}/M_{\rm d}=1.2$) and that the blobs are rotating counterclockwise on the plane of the sky.} This result shows that polarimetry can provide a diagnostic of the geometry of clouds, which will serve as an important constraint for {improving} the Rigidly Rotating Magnetosphere model.

We analyse a feature that is detectable at the differential phases of high-resolution spectro-interferometry of Be stars. The origin of this feature (dubbed CQE-PS, Central Quasi Emission Phase Signature) lies in the differential absorption of photospheric radiation by the circumstellar disk, which is spectroscopically detected as a CQE line profile in shell stars. This phenomenon has great diagnostic potential for Be star disks, revealing properties of these disks on the scale of a few stellar radii.