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 present paper is the culminating one of a series aimed to contribute to the understanding of the kinematic structures of the solar neighbourhood (SN), explaining the origin of the Local Arm and relating the moving groups with the spiral-arms resonances in the disk. With a model for the Galactic potential, with the Sun inside the spiral corotation resonance (CR), we integrate the 2D orbits of test particles distributed in birthplaces along the main spiral arms, the Local Arm, and in the axisymmetric disk. A comparison of the resulting U-V plane of the SN with that provided by Gaia DR2 confirms our previous conclusion that the moving groups of Coma Berenices, Pleiades, and Hyades are associated with the CR, and that the Hercules stream is formed by the bulk of high-order inner Lindblad resonances. The kinematic structures result from stellar orbits trapped by the spiral resonances in a timespan of ~ 1 Gyr, indicating the long-living nature of the spiral structure and challenging recent arguments in favor of short-lived structures originated from incomplete phase mixing in the Galactic disk. As a byproduct, our simulations give some insight into the birthplaces of the stars presently located in the SN; the majority of stars of the main moving groups and the Hercules stream were likely born in the Local Arm, while stars of the Sirius group possibly originated from the outer segment of the Sagittarius-Carina arm. We also propose the spiral resonances as the dynamical origin for the diagonal ridges in the Galactic distribution of rotation velocities.

The VISCACHA (VIsible Soar photometry of star Clusters in tApii and Coxi HuguA\footnote{LMC and SMC names in the Tupi-Guarani language spoken by native people in Brazil}) Survey is an ongoing project based on deep and spatiallyresolved photometric observations of Magellanic Cloud star clusters, collected using the SOuthern Astrophysical Research (SOAR) telescope together with the SOAR Adaptive Module Imager. So far we have used $>$300h of telescope time to observe $\sim$150 star clusters, mostly with low mass ($M < 10^4 M_{\odot}$) on the outskirts of the LMC and SMC. With this high-quality data set, we homogeneously determine physical properties using deep colour-magnitude diagrams (ages, metallicities, reddening, distances, mass, luminosity and mass functions) and structural parameters (radial density profiles, sizes) for these clusters which are used as a proxy to investigate the interplay between the Magellanic Clouds and their evolution. We present the VISCACHA survey and its initial results, based on our first two papers. The project's long term goals and expected legacy to the community are also addressed.

We construct self-gravitating razor-thin disks of counterrotating dust around Schwarzschild black holes (BHs) by applying the ``displace, cut, and reflect'' method to known seed solutions representing multiholes. All but one of the sources of the seed solution generate the surrounding annular disk, whereas the remaining BH is kept unaltered and lies at the disk center after the transformation. The disks are infinite in extent, have annular character, and are linearly stable up to the innermost stable circular orbit of the system. Moreover, all energy conditions are satisfied along the disk and the spacetime is asymptotically flat, having finite ADM mass. We also comment on charged disks around extremal Reissner-Nordstr\"om BHs constructed from a Majumdar-Papapetrou $N$-BH seed solution. These simple examples can be extended to more general ``BH + disk'' solutions, obtained by the same method from seeds of the type ``BH + arbitrary axisymmetric source''. A natural follow-up of this work would be to construct disks around Reissner-Nordstr\"om BHs with arbitrary charge-to-mass ratio and around Kerr BHs.

The VISCACHA (VIsible Soar photometry of star Clusters in tApii and Coxi HuguA) Survey is an ongoing project based on deep photometric observations of Magellanic Cloud star clusters, collected using the SOuthern Astrophysical Research (SOAR) telescope together with the SOAR Adaptive Module Imager. Since 2015 more than 200 hours of telescope time were used to observe about 130 stellar clusters, most of them with low mass (M < 10$^4$ M$_\odot$) and/or located in the outermost regions of the Large Magellanic Cloud and the Small Magellanic Cloud. With this high quality data set, we homogeneously determine physical properties from statistical analysis of colour-magnitude diagrams, radial density profiles, luminosity functions and mass functions. Ages, metallicities, reddening, distances, present-day masses, mass function slopes and structural parameters for these clusters are derived and used as a proxy to investigate the interplay between the environment in the Magellanic Clouds and the evolution of such systems. In this first paper we present the VISCACHA Survey and its initial results, concerning the SMC clusters AM3, K37, HW20 and NGC796 and the LMC ones KMHK228, OHSC3, SL576, SL61 and SL897, chosen to compose a representative subset of our cluster sample. The project's long term goals and legacy to the community are also addressed.

We recently found that regular orbits in axially symmetric galactic disks have their envelopes $Z(R)$ accurately described by the relation $Z(R)\propto[\Sigma_I(R)]^{-1/3}$, if their amplitudes are comparable to the disk thickness, where $\Sigma_I$ is the surface density of the disk (integrated over its whole vertical range). Moreover, the usual adiabatic approximation gives a good description of the orbits' envelopes for low vertical amplitudes. However, these two approaches are apparently disconnected, since their expressions differ qualitatively. Our purpose in this paper is to fill this gap by extending these previous formulae to regular orbits with arbitrary vertical amplitudes inside the disk. We compare existing $Z(R)$ estimates: the razor-thin disk case, the adiabatic approximation (low-amplitude orbits in three-dimensional disks), and the integrated surface-density estimate (high-amplitude orbits in three-dimensional disks) in order to establish a connection between them. The formula presented here links the aforementioned results in an elegant and continuous way, being valid for vertical amplitudes throughout the whole vertical extension of the disk and with an expression which has the same form for all regimes. The advantage of the present formalism is the dependence of $Z(R)$ only on observable quantities, namely the disk's vertically integrated surface density, without the need to obtain the gravitational potential for the system.

The aim of this work is to contribute to the understanding of the stellar velocity distribution in the solar neighborhood (SN). We propose that the structures on the $U-V$ planes, known as the moving groups, can be mainly explained by the spiral arms perturbations. The applied model of the Galactic disk and spiral arms, with the parameters defined by observational data and with pattern speed $\Omega_p=$28.0 km s$^{-1}$ kpc$^{-1}$, is the same that allowed us to explain the origin of the Local Arm and the Sun's orbit trapped inside the corotation resonance (CR). We show that the $U-V$ picture of the SN consists of the main component, associated with the CR, and the inner and outer structures, which we could associate with the Hercules and Sirius streams, respectively. The Coma-Berenices and Hyades-Pleiades groups and the Sun itself belong to the main part. The substructures of Hercules are formed mainly by the nearby 8/1, 12/1, and even 6/1 inner Lindblad resonances, while Sirius is shaped by the bulk of overlapping outer Lindblad resonances, -8/1, -12/1, -16/1, which are stuck to the CR. This richness in resonances only exists near corotation, which should be of the spiral arms, not of the Galactic bar, whose stable corotation zone is far away from the Sun. The model's predictions of the velocity distribution match qualitatively and quantitatively the distribution provided by Gaia DR2.

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.

Observational data indicate that the Milky Way is a barred spiral galaxy. Computation facilities and availability of data from Galactic surveys stimulate the appearance of models of the Galactic structure. More efforts to build dynamical models containing both spiral arms and the central bar/bulge are needed. We expand the study of the stellar dynamics in the Galaxy by adding the bar/bulge component to a model with spiral arms introduced in our previous paper. The model is tested by applying it to the solar neighborhood, where observational data are more precise. We model analytically the potential of the Galaxy to derive the force field in its equatorial plane. The model comprises an axisymmetric disc derived from the observable rotation curve, four arms with Gaussian-shaped groove profiles, and a classical elongated/oblate ellipsoidal bar/bulge structure. The parameters describing the bar/bulge are constrained by observations and the stellar dynamics, and their possible limits are determined. A basic model results in a bar of 2.9 kpc in length, with a mass of the order of a few 10$^9M_\odot$. The size and orientation of the bar are also restricted by the position of masers with VLBI distances. The bar's rotation speed is constrained to $\Omega_{\rm bar}<50$ km s$^{-1}$ kpc$^{-1}$ taking into account the allowed mass range. We conclude that our basic model is compatible with observations and with the dynamical constraints. The model explains simultaneously the bulk of the main moving groups, associated here with the spiral corotation resonance, and the Hercules stream, associated with several inner high-order spiral resonances; in particular, with the 8/1 resonance. From the dynamical constraints on the bar's angular speed, it is unlikely that the bar's OLR lies near the solar circle; moreover, its proximity would compromise the stability of the Local Arm structure.

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.

Bulge globular clusters (GCs) with metallicities [Fe/H] <~ -1.0 and blue horizontal branches are candidates to harbor the oldest populations in the Galaxy. Based on the analysis of HST proper-motion-cleaned color-magnitude diagrams in filters F435W and F625W, we determine physical parameters for the old bulge globular clusters NGC 6522 and NGC 6626 (M28), both with well-defined blue horizontal branches. We compare these results with similar data for the inner halo cluster NGC 6362. These clusters have similar metallicities (-1.3 <= [Fe/H] <= -1.0) obtained from high resolution spectroscopy. We derive ages, distance moduli, and reddening values by means of statistical comparisons between observed and synthetic fiducial lines employing likelihood statistics and the Markov chain Monte Carlo method. The synthetic fiducial lines were generated using alpha-enhanced BaSTI and Dartmouth stellar evolutionary models, adopting both canonical (Y ~ 0.25) and enhanced (Y~ 0.30-0.33) helium abundances. RR Lyrae stars were employed to determine the HB magnitude level, providing an independent indicator to constrain the apparent distance modulus and the helium enhancement. The shape of the observed fiducial line could be compatible with some helium enhancement for NGC 6522 and NGC 6626, but the average magnitudes of RR Lyrae stars tend to rule out this hypothesis. Assuming canonical helium abundances, BaSTI and Dartmouth models indicate that all three clusters are coeval, with ages between ~ 12.5 and 13.0 Gyr. The present study also reveals that NGC 6522 has at least two stellar populations, since its CMD shows a significantly wide subgiant branch compatible with 14% +/- 2% and 86% +/- 5% for first and second generations, respectively.

We apply the viscous decretion disc (VDD) model to interpret the infrared disc continuum emission of 80 Be stars observed in different epochs. In this way, we determined 169 specific disc structures, namely their density scale, $\rho_0$, and exponent, $n$. We found that the $n$ values range mainly between $1.5$ and $3.5$, and $\rho_0$ varies between $10^{-12}$ and $10^{-10}\,\mathrm{g\,cm^{-3}}$, with a peak close to the lower value. Our large sample also allowed us to firmly establish that the discs around early-type stars are denser than in late-type stars. Additionally, we estimated the disc mass decretion rates and found that they range between $10^{-12}$ and $10^{-9}\,\mathrm{M_{\odot}\,yr^{-1}}$. These values are compatible with recent stellar evolution models of fast-rotating stars. One of the main findings of this work is a correlation between the $\rho_0$ and $n$ values. In order to find out whether these relations can be traced back to the evolution of discs or have some other origin, we used the VDD model to calculate temporal sequences under different assumptions for the time profile of the disc mass injection. The results support the hypothesis that the observed distribution of disc properties is due to a common evolutionary path. In particular, our results suggest that the timescale for disc growth, during which the disc is being actively fed by mass injection episodes, is shorter than the timescale for disc dissipation, when the disc is no longer fed by the star and dissipates as a result of the viscous diffusion of the disc material.

We investigated the properties of $80$ Be star disks, using the time dependent viscous decretion disk (VDD) model to interpret their mid-infrared continuum fluxes. Besides estimating the typical parameter ranges for such disks, this approach allowed us to relate the disk structure to its evolutionary status. In this contribution, we present a brief summary of the recently published results, and discuss some properties and time-scales of the VDD model solution.

The Local arm of the Milky Way, a short spiral feature near the Sun whose existence is known for decades, was recently observed in detail with different tracers. Many efforts have been dedicated to elaborate plausible hypotheses concerning the origin of the main spiral arms of the Galaxy; however, up to now, no specific mechanism for the origin of the Local arm was proposed. Here we explain, for the first time, the Local arm as an outcome of the spiral corotation resonance, which traps arm tracers and the Sun inside it. We show that the majority of maser sources belonging to the Local arm, together with the Sun, evolve inside the corotation resonance, never crossing the main spiral arms but instead oscillating in the region between them. This peculiar behavior of the Sun could have numerous consequences to our understanding of the local kinematics of stars, the Galactic Habitable Zone, and the Solar System evolution.

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)

The sum of squared epicyclic frequencies of nearly circular motion ($\omega_r^2+\omega_\theta^2$) in axially symmetric configurations of Newtonian gravity is known to depend both on the matter density and on the angular velocity profile of circular orbits. It was recently found that this sum goes to zero at the photon orbits of Schwarzschild and Kerr spacetimes. However, these are the only relativistic configurations for which such result exists in the literature. Here, we extend the above formalism in order to describe the analogous relation for geodesic motion in arbitrary static, axially symmetric, asymptotically flat solutions of general relativity. The sum of squared epicyclic frequencies is found to vanish at photon radii of vacuum solutions. In the presence of matter, we obtain that $\omega_r^2+\omega_\theta^2>0$ for perturbed timelike circular geodesics on the equatorial plane if the strong energy condition holds for the matter-energy fluid of spacetime; in vacuum, the allowed region for timelike circular geodesic motion is characterized by the inequality above. The results presented here may be of use to shed light on general issues concerning the stability of circular orbits once they approach photon radii, mainly the ones corresponding to stable photon motion.

During the last few decades, there has been a growing interest in exact solutions of Einstein equations describing razor-thin disks. Despite the progress in the area, the analytical study of geodesic motion crossing the disk plane in these systems is not yet so developed. In the present work, we propose a definite vertical stability criterion for circular equatorial timelike geodesics in static, axially symmetric thin disks, possibly surrounded by other structures preserving axial symmetry. It turns out that the strong energy condition for the disk stress-energy content is sufficient for vertical stability of these orbits. Moreover, adiabatic invariance of the vertical action variable gives us an approximate third integral of motion for oblique orbits which deviate slightly from the equatorial plane. Such new approximate third integral certainly points to a better understanding of the analytical properties of these orbits. The results presented here, derived for static spacetimes, may be a starting point to study the motion around rotating, stationary razor-thin disks. Our results also allow us to conjecture that the strong energy condition should be sufficient to assure transversal stability of periodic orbits for any singular timelike hypersurface, provided it is invariant under the geodesic flow.

Context: Resonances in the stellar orbital motion under perturbations from spiral arms structure play an important role in the evolution of the disks of spiral galaxies. The epicyclic approximation allows the determination of the corresponding resonant radii on the equatorial plane (for nearly circular orbits), but is not suitable in general. Aims: We expand the study of resonant orbits by analysing stellar motions perturbed by spiral arms with Gaussian-shaped profiles without any restriction on the stellar orbital configurations, and we expand the concept of Lindblad (epicyclic) resonances for orbits with large radial excursions. Methods: We define a representative plane of initial conditions, which covers the whole phase space of the system. Dynamical maps on representative planes are constructed numerically, in order to characterize the phase-space structure and identify the precise location of resonances. The study is complemented by the construction of dynamical power spectra, which provide the identification of fundamental oscillatory patterns in the stellar motion. Results: Our approach allows a precise description of the resonance chains in the whole phase space, giving a broader view of the dynamics of the system when compared to the classical epicyclic approach, even for objects in retrograde motion. The analysis of the solar neighbourhood shows that, depending on the current azimuthal phase of the Sun with respect to the spiral arms, a star with solar kinematic parameters may evolve either inside the stable co-rotation resonance or in a chaotic zone. Conclusions: Our approach contributes to quantifying the domains of resonant orbits and the degree of chaos in the whole Galactic phase-space structure. It may serve as a starting point to apply these techniques to the investigation of clumps in the distribution of stars in the Galaxy, such as kinematic moving groups.

We present an analytical simple formula for an approximated third integral of motion associated with nearly equatorial orbits in the Galaxy: $I_{3}=Z\Sigma_{I}^{1/3}$, where $Z(R)$ is the vertical amplitude of the orbit at galactocentric distance $R$ and $\Sigma_{I}(R)$ is the integrated dynamical surface mass density of the disk, a quantity which has recently become measurable. We also suggest that this relation is valid for disk-crossing orbits in a wide variety of axially symmetric galactic models, which range from razor-thin disks to disks with non-negligible thickness, whether or not the system includes bulges and halos. We apply our formalism to a Miyamoto-Nagai model and to a realistic model for the Milky Way. In both cases, the results provide fits for the shape of nearly equatorial orbits which are better than the corresponding curves obtained by the usual adiabatic approximation when the orbits have vertical amplitudes comparable to the disk's scale height. We also discuss the role of this approximate third integral of motion in modified theories of gravity.

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.

We consider the three-dimensional bounded motion of a test particle around razor-thin disk configurations, by focusing on the adiabatic invariance of the vertical action associated with disk-crossing orbits. We find that it leads to an approximate third integral of motion predicting envelopes of the form $Z(R)\propto[\Sigma(R)]^{-1/3}$, where $R$ is the radial galactocentric coordinate, $Z$ is the z-amplitude (vertical amplitude) of the orbit and $\Sigma$ represents the surface mass density of the thin disk. This third integral, which was previously formulated for the case of flattened 3D configurations, is tested for a variety of trajectories in different thin-disk models.

At super-Eddington rates accretion flows onto black holes have been described as slim (aspect ratio $H/R \lesssim 1$) or thick (H/R >1) discs, also known as tori or (Polish) doughnuts. The relation between the two descriptions has never been established, but it was commonly believed that at sufficiently high accretion rates slim discs inflate, becoming thick. We wish to establish under what conditions slim accretion flows become thick. We use analytical equations, numerical 1+1 schemes, and numerical radiative MHD codes to describe and compare various accretion flow models at very high accretion rates.We find that the dominant effect of advection at high accretion rates precludes slim discs becoming thick. At super-Eddington rates accretion flows around black holes can always be considered slim rather than thick.

In the last decade, the viscous decretion disk model has emerged as the new paradigm for Be star disks. In this contribution, we propose a simple analytical model to estimate the continuum infrared excess arising from these circumstellar disks, in the light of the currently accepted scenario. We demonstrate that the disk can be satisfactorily described by a two component system: an inner optically thick region, which we call the pseudo-photosphere, and a diffuse outer part. In particular, a direct connexion between the disk brightness profile and the thermal structure is derived, and then confronted to realistic numerical simulations. This result quantifies how the non-isothermality of the disk ultimately affects both infrared measured fluxes and visibilities.

Possible astrophysical consequences of the Ho\v{r}ava quantum gravity theory have been recently studied by several authors. They usually employ the Kehagias-Sfetsos (KS) spacetime which is a spherically symmetric vacuum solution of a specific version of Ho\v{r}ava's gravity. The KS metric has several unusual geometrical properties that in the present article we examine by means of the often used technique of embedding diagrams. We pay particular attention to the transition between naked singularity and black-hole states, which is possible along some particular sequences of the KS metrics.

We investigate the continuum emission of viscous decretion discs around Be stars in this paper. The results obtained from non-LTE (local thermodynamic equilibrium) radiative transfer models show two regimes in the disc surface brightness profile: an inner optically thick region, which behaves as a pseudo-photosphere with a wavelength-dependent size, and an optically thin tenuous outer part, which contributes with about a third of the total flux. The isophotal shape of the surface brightness is well described by elliptical contours with an axial ratio $b/a=\cos i$ for inclinations $i<75^{\circ}$. Based on these properties, a semi-analytical model was developed to describe the continuum emission of gaseous discs. It provides fluxes and spectral slopes at the infrared within an accuracy of $10\%$ and $5\%$, respectively, when compared to the numerical results. The model indicates that the infrared spectral slope is mainly determined by both the density radial slope and the disc flaring exponent, being practically independent of disc inclination and base density. As a first application, the density structure of 15 Be stars was investigated, based on the infrared flux excess, and the results compared to previous determinations in the literature. Our results indicate that the decretion rates are in the range of $10^{-12}$ to $10^{-9}\,{\rm M_{\odot}\,yr^{-1}}$, which is at least two orders of magnitude smaller than the previous outflowing disc model predictions.

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 spatially and spectrally resolved Br-gamma emission around the planet-hosting, transitional Herbig Ae/Be star HD 100546. Aiming to gain insight into the physical origin of the line in possible relation to accretion processes, we carried out Br-gamma spectro-interferometry using AMBER/VLTI from three different baselines achieving spatial and spectral resolutions of 2-4 mas and 12000. The Br-gamma visibility is larger than that of the continuum for all baselines. Differential phases reveal a shift between the photocentre of the Br-gamma line -displaced 0.6 mas (0.06 au at 100 pc) NE from the star- and that of the K-band continuum emission -displaced 0.3 mas NE from the star. The photocentres of the redshifted and blueshifted components of the Br-gamma line are located NW and SE from the photocentre of the peak line emission, respectively. Moreover, the photocentre of the fastest velocity bins within the spectral line tends to be closer to that of the peak emission than the photocentre of the slowest velocity bins. Our results are consistent with a Br-gamma emitting region inside the dust inner rim (<0.25 au) and extending very close to the central star, with a Keplerian, disc-like structure rotating counter-clockwise, and most probably flared (25 deg). Even though the main contribution to the Br-gamma line does not come from gas magnetically channelled on to the star, accretion on to HD 100546 could be magnetospheric, implying a mass accretion rate of a few 10^-7 Msun/yr. This value indicates that the observed gas has to be replenished on time-scales of a few months to years, perhaps by planet-induced flows from the outer to the inner disc as has been reported for similar systems.

The study of off-equatorial orbits in razor-thin disks is still in its beginnings. Contrary to what was presented in the literature in recent publications, the vertical stability criterion for equatorial circular orbits cannot be based on the vertical epicyclic frequency, because of the discontinuity in the gravitational field on the equatorial plane. We present a rigorous criterion for the vertical stability of circular orbits in systems composed by a razor-thin disk surrounded by a smooth axially symmetric distribution of matter, the latter representing additional structures such as thick disk, bulge and (dark matter) halo. This criterion is satisfied once the mass surface density of the thin disk is positive. Qualitative and quantitative analyses of nearly equatorial orbits are presented. In particular, the analysis of nearly equatorial orbits allows us to construct an approximate analytical third integral of motion in this region of phase-space, which describes the shape of these orbits in the meridional plane.

We discuss photon and test-particle orbits in the Kehagias-Sfetsos (KS) metric. For any value of the Ho\v{r}ava parameter $\omega$, there are values of the gravitational mass $M$ for which the metric describes a naked singularity, and this is always accompanied by a vacuum "antigravity sphere" on whose surface a test particle can remain at rest (in a zero angular momentum geodesic), and inside which no circular geodesics exist. The observational appearance of an accreting KS naked singularity in a binary system would be that of a quasi-static spherical fluid shell surrounded by an accretion disk, whose properties depend on the value of $M$, but are always very different from accretion disks familiar from the Kerr-metric solutions. The properties of the corresponding circular orbits are qualitatively similar to those of the Reissner-Nordstr\"om naked singularities. When event horizons are present, the orbits outside the Kehagias-Sfetsos black hole are qualitatively similar to those of the Schwarzschild metric.

We make a revision of the stability criteria for equatorial circular orbits, obtained from the epicyclic approximation, which is widely used in Newtonian models for axisymmetric galaxies. We find that, for the case of thin disk models, the criteria of vertical stability must be reformulated, due to the discontinuity in the gravitational field. We show that, for a model characterized by a surface mass density $\Sigma$, the necessary and sufficient condition to have vertically stable orbits is that $\Sigma>0$. On the other hand, the criteria for radial stability is the same as in thick diks, i.e. that the epicyclic frequency is positive.

A sample of 27 sources, catalogued as pre-main sequence stars by the Pico dos Dias Survey (PDS), is analyzed to investigate a possible contamination by post-AGB stars. The far-infrared excess, due to dust present in the circumstellar envelope, is typical for both categories: young stars and objects that have already left the main sequence and are suffering a severe mass-loss. The presence of two known post-AGB stars in our sample inspired us to seek for other very likely or possible post-AGB objects among PDS sources previously suggested to be Herbig Ae/Be stars, by revisiting the observational database of this sample. In a comparative study with well known post-AGBs, several characteristics were evaluated: (i) parameters related to the circumstellar emission; (ii) spatial distribution to verify the background contribution from dark clouds; (iii) spectral features, and (iv) optical and infrared colors. These characteristics suggest that 7 objects of the studied sample are very likely post-AGBs, 5 are possible post-AGBs, 8 are unlikely post-AGBs, and the nature of 7 objects remains unclear.