Context. Be stars are physically complex systems that continue to challenge theory to understand their rapid rotation, complex variability and decretion disks. $\gamma$ Cassiopeiae ($\gamma$ Cas) is one such star but is even more curious because of its unexplained hard thermal X-ray emission. Aims. We aim to examine the optical variability of $\gamma$ Cas and thereby to shed more light on its puzzling behaviour. Methods. Three hundred twenty-one archival H$\alpha$ spectra from 2006 to 2017 are analysed to search for frequencies corresponding to the 203.5 day orbit of the companion. Space photometry from the SMEI satellite from 2003 to 2011 and the BRITE-Constellation of nano-satellites between 2015 and 2019 is investigated in the period range from a couple of hours to a few days. Results. The orbital period of the companion of 203.5 days is confirmed with independent measurements from the structure of the H$\alpha$ line emission. A strong blue/red asymmetry in the amplitude distribution across the H$\alpha$ emission line could hint at a spiral structure in the decretion disk. With the space photometry, the known frequency of 0.82 d$^{-1}$ is confirmed in data from the early 2000s. A higher frequency of 2.48 d$^{-1}$ is present in the data from 2015 to 2019 and possibly also in the early 2000s. A third frequency at 1.25 d$^{-1}$ is proposed to exist in both SMEI and BRITE data. The only explanation covering all three rapid variations seems to be nonradial pulsation. The two higher frequencies are incompatible with rotation.

We investigated 12 unclassified B[e] stars or candidates, 8 from the Galaxy, 2 from the Large Magellanic Cloud (LMC) and 2 from the Small Magellanic Cloud (SMC). Based on the analysis of high-resolution spectroscopic (FEROS) and photometric data, we confirmed the presence of the B[e] phenomenon for all objects of our sample, except for one (IRAS 07455-3143). We derived their effective temperature, spectral type, luminosity class, interstellar extinction and, using the distances from Gaia DR2, we obtained their bolometric magnitude, luminosity and radius. Modeling of the forbidden lines present in the FEROS spectra revealed information about the kinematics and geometry of the circumstellar medium of these objects. In addition, we analyzed the light curves of four stars, finding their most probable periods. The evolutionary stage of 11 stars of our sample is suggested from their position on the HR diagram, taking into account evolutionary tracks of stars with solar, LMC and SMC metallicities. As results, we identified B and B[e] supergiants, B[e] stars probably at the main sequence or close to its end, post-AGB and HAeB[e] candidates, and A[e] stars in the main sequence or in the pre-main sequence. However, our most remarkable results are the identification of the third A[e] supergiant (ARDB\,54, the first one in the LMC), and of an "LBV impostor" in the SMC (LHA 115-N82).

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.

In early-type Be stars, groups of nonradial pulsation (NRP) modes with numerically related frequencies may be instrumental for the release of excess angular momentum through mass-ejection events. Difference and sum/harmonic frequencies often form additional groups. The goal of this study is to find out whether a similar frequency pattern occurs in the cooler third-magnitude B7-8\,IIIe shell star $\nu$ Pup. Time-series analyses are performed of space photometry with BRITE-Constellation (2015, 2016/17, and 2017/18), SMEI (2003--011), and Hipparcos (1989-1993). Two IUE SWP and 27 optical echelle spectra spanning 20 years were retrieved from various archives. The optical spectra exhibit no anomalies or well-defined variabilities. A magnetic field was not detected. All three photometry satellites recorded variability near 0.656 c/d which is resolved into three features separated by ~0.0021 c/d. First harmonics form a second frequency group, also spaced by ~0.0021 c/d. The frequency spacing is very nearly but not exactly equidistant. Variability near 0.0021 c/d was not detected. The long-term frequency stability could be used to derive meaningful constraints on the properties of a putative companion star. The IUE spectra do not reveal the presence of a hot subluminous secondary. $\nu$\,Pup is another Be star exhibiting an NRP variability pattern with long-term constancy and underlining the importance of combination frequencies and frequency groups. The star is a good target for efforts to identify an effectively single Be star.

Prompted by peculiar spectroscopic variability observed in SDSS/APOGEE $H$-band spectra, we monitored the Be star HD 55606 using optical spectroscopy and found that it is an exotic double-lined spectroscopic binary (SB2) consisting of a Be star and a hot, compact companion that is probably an OB subdwarf (sdOB) star. Motion of the sdOB star is traced by its impact on the strong He~I lines, observed as radial velocity ($V_{r}$) variable, double-peaked emission profiles with narrow central absorption cores. Weak He II 4686 {\AA} absorption associated with the companion star is detected in most spectra. Use of the emission peaks of low-ionization emission lines to trace the Be star $V_{r}$ and the He I lines to trace the companion star $V_{r}$ yields a circular orbital solution with a 93.8-day period and masses of $M_{\rm Be}=6.2$ $M_{\rm \odot}$ and $M_{\rm sdOB}=0.9$ $M_{\rm \odot}$ in the case of $i=80^{\circ}$. HD 55606 exhibits a variety of phase-locked variability, including the development of shell lines twice per orbit. The shell phases coincide with variation in the double emission peak separations, and both forms of variability are likely caused by a two-armed spiral density perturbation in the Be disk. The intensity ratios of the double emission peaks are also phase-locked, possibly indicating heating by the sdOB star of the side of the Be disk facing it. HD 55606 is a new member of the growing sample of Be+sdOB binaries, in which the Be star's rapid rotation and ability to form a disk can be attributed to past mass transfer.

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.

Observations of 25 Ori much expand the picture derived of other early-type Be stars with BRITE and SMEI. Two instead of one difference frequencies rule the variability: (a) The lower one, 0.0129 c/d, is the frequency of events with full amplitudes of 100-200 mmag which may signal mass loss possibly driven by the higher one, 0.1777 c/d. (b) Much of the entire power spectrum is a tightly woven network of combination frequencies: (i) Below 0.25\,c/d, numerous frequencies are difference frequencies. (ii) Many frequencies above 2.5 c/d can be represented as sum frequencies and in a few cases as harmonics. (iii) Many frequencies between 1.1 and 1.75\,c/d can be portrayed as parents of combination frequencies. The number and fraction of combination frequencies increases steeply with decreasing amplitude and and accuracy of the frequency matching.

The BRITE Constellation of nanosatellites obtained mmag photometry of 28 Cygni for 11 months in 2014-2016. Observations with the Solar Mass Ejection Imager in 2003-2010 and 118 H$\alpha$ line profiles were added. For decades, 28 Cyg has exhibited four large-amplitude frequencies: two closely spaced frequencies of spectroscopically confirmed $g$ modes near 1.5 c/d, one slightly lower exophotospheric (Stefl) frequency, and at 0.05 c/d the difference frequency between the two g modes. This top-level framework is indistinguishable from eta Cen (Paper I), which is also very similar in spectral type, rotation rate, and viewing angle. The Stefl frequency is the only one that does not seem to be affected by the difference frequency. The amplitude of the latter undergoes large variations; around maximum the amount of near-circumstellar matter is increased, and the amplitude of the Stefl frequency grows by some factor. During such brightenings dozens of transient spikes appear in the frequency spectrum, concentrated in three groups. Only eleven frequencies were common to all years of BRITE observations. Be stars seem to be controlled by several coupled clocks, most of which are not very regular on timescales of weeks to months but function for decades. The combination of g modes to the low difference frequency and/or the atmospheric response to it appears significantly nonlinear. Like in eta Cen, the difference-frequency variability seems the main responsible for the modulation of the star-to-disc mass transfer in 28 Cyg. A hierarchical set of difference frequencies may reach the longest timescales known of the Be phenomenon.

We use a smoothed particle hydrodynamics (SPH) code to examine the effects of misaligned binary companions on Be star discs. We systematically vary the degree of misalignment between the disc and the binary orbit, as well as the disc viscosity and orbital period to study their effects on the density in the inner and outer parts of the disc. We find that varying the degree of misalignment, the viscosity, and the orbital period affects both the truncation radius and the density structure of the outer disc, while the inner disc remains mostly unaffected. We also investigate the tilting of the disc in the innermost part of the disc and find the tilt increases with radius until reaching a maximum around 5 stellar radii. The direction of the line of nodes, with respect to the equator of the central star, is found to be offset compared to the orbital line of nodes, and to vary periodically in time, with a period of half a orbital phase. We also compare the scale height of our discs with the analytical scale height of an isothermal disc, which increases with radius as $r^{1.5}$. We find that this formula reproduces the scale height well for both aligned and misaligned systems but underestimates the scale height in regions of the disc where density enhancements develop.

We present new optical broad-band (UBVRI) aperture polarimetric observations of 53 post-asymptotic giant branch (AGB) stars selected to exhibit a large near-infrared excess. 24 out of the 53 stars (45% of our sample) are presented for the first time. A statistical analysis shows four distinctive groups of polarized post-AGB stars: unpolarized or very lowly polarized (degree of polarization or DoP < 1%), lowly polarized (1% < DoP < 4%), moderately polarized (4% < DoP < 8%) and highly polarized (DoP > 8%). 23 out of the 53 (66%) belong to the first group, 10 (19%) to the second, five (9%) to the third and only three (6%) to the last group. Approximately, 34% of our sample was found to be unpolarized objects, which is close to the percentage of round planetary nebulae. On average, the low and moderate groups show a wavelength-dependent polarization that increases towards shorter wavelength, implying an intrinsic origin of the polarization, which signifies a Rayleigh-like scattering spectrum typical for non-symmetrical envelopes composed principally of small dust grains. The moderately polarized stars exhibit higher K-W3 and W1-W3 colour indices compared with the group of lowly polarized stars suggesting a possible relation between DoP and mass-loss rate. Moreover, they are found to be systematically colder (redder in B-V), which may be associated with the condensation process close to these stars that results in higher degree of polarization. We also provide evidence that multiple scattering in optically thin polar outflows is the mechanism that gives high DoP in post-AGB stars with a bipolar or multi-polar envelopes.

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.

Be stars are surrounded by outflowing circumstellar matter structured in the form of decretion discs. They are often members of binary systems, where it is expected that the decretion disc interacts both radiatively and gravitationally with the companion. In this work we study how various orbital (period, mass ratio and eccentricity) and disc (viscosity) parameters affect the disc structure in coplanar systems. We simulate such binaries with the use of a smoothed particle hydrodynamics code. The main effects of the secondary on the disc are its truncation and the accumulation of material inwards of truncation. We find two cases with respect to the effects of eccentricity: (i) In circular or nearly circular prograde orbits, the disc maintains a rotating, constant in shape, configuration, which is locked to the orbital phase. The disc is smaller in size, more elongated and more massive for low viscosity parameter, small orbital separation and/or high mass ratio. (ii) Highly eccentric orbits are more complex, with the disc structure and total mass strongly dependent on the orbital phase and the distance to the secondary. We also study the effects of binarity in the disc continuum emission. Since the infrared and radio SED are sensitive to the disc size and density slope, the truncation and matter accumulation result in considerable modifications in the emergent spectrum. We conclude that binarity can serve as an explanation for the variability exhibited in observations of Be stars, and that our model can be used to detect invisible companions.

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.

The modeling and analysis generic interface for external numerical codes (MAGIX) is a model optimizer developed under the framework of the coherent set of astrophysical tools for spectroscopy (CATS) project. The MAGIX package provides a framework of an easy interface between existing codes and an iterating engine that attempts to minimize deviations of the model results from available observational data, constraining the values of the model parameters and providing corresponding error estimates. Many models (and, in principle, not only astrophysical models) can be plugged into MAGIX to explore their parameter space and find the set of parameter values that best fits observational/experimental data. MAGIX complies with the data structures and reduction tools of ALMA (Atacama Large Millimeter Array), but can be used with other astronomical and with non-astronomical data.

Molecular counterparts to atomic jets have been detected within 1000 AU of young stars. Reproducing them is a challenge for proposed ejection models. We explore whether molecules may survive in an MHD disk wind invoked to reproduce the kinematics and tentative rotation signatures of atomic jets in T Tauri stars. The coupled ionization, chemical and thermal evolution along dusty flow streamlines is computed for a prescribed MHD disk wind solution, using a method developed for magnetized shocks in the interstellar medium. Irradiation by wind-attenuated coronal X-rays and FUV photons from accretion hot spots is included, with self-shielding of H2 and CO. Disk accretion rates of 5e-6, 1e-6 and 1e-7 solar masses per year are considered, representative of low-mass young protostars (Class 0), evolved protostars (Class I) and very active T Tauri stars (Class II). The disk wind has an onion-like thermo-chemical structure, with streamlines launched from larger radii having lower temperature and ionisation, and higher H2 abundance. The coupling between charged and neutral fluids is sufficient to eject molecules from the disk out to 9 AU. The launch radius beyond which most H2 survives moves outward with evolutionary stage. CO survives in the Class 0 but is significantly photodissociated in the Class I/II. Balance between ambipolar heating and molecular cooling establishes an asymptotic temperature 700-3000 K, with cooler jets at earlier protostellar stages. Endothermic formation of H2O is efficient with abundances up to 1e-4, while CH+ and SH+ can exceed 1e-6 in the Class I/II winds. A centrifugal MHD disk wind launched from beyond 0.2-1 AU can produce molecular jets/winds up to speeds 100 km/s in young low-mass stars. The model predicts a high ratio H2/CO and an increase of molecular launch radius, temperature, and flow width as the source evolves, in agreement with current observed trends.