The presence of a nearby companion alters the evolution of massive stars in binary systems, leading to phenomena such as stellar mergers, x-ray binaries, and gamma-ray bursts. Unambiguous constraints on the fraction of massive stars affected by binary interaction were lacking. We simultaneously measured all relevant binary characteristics in a sample of Galactic massive O stars and quantified the frequency and nature of binary interactions. More than 70% of all massive stars will exchange mass with a companion, leading to a binary merger in one-third of the cases. These numbers greatly exceed previous estimates and imply that binary interaction dominates the evolution of massive stars, with implications for populations of massive stars and their supernovae.

Multiplicity is one of the most fundamental observable properties of massive O-type stars and offers a promising way to discriminate between massive star formation theories. Nevertheless, companions at separations between 1 and 100 milli-arcsec (mas) remain mostly unknown due to intrinsic observational limitations. At a typical distance of 2 kpc, this corresponds to projected physical separations of 2-200 AU. The Southern MAssive Stars at High angular resolution survey (smash+) was designed to fill this gap by providing the first systematic interferometric survey of Galactic massive stars. We observed 117 O-type stars with VLTI/PIONIER and 162 O-type stars with NACO/SAM, respectively probing the separation ranges 1-45 and 30-250 mas and brightness contrasts of ∆H < 4 and ∆H < 5. Taking advantage of NACO's field-of-view, we further uniformly searched for visual companions in an 8 -radius down to ∆H = 8. This paper describes the observations and data analysis, reports the discovery of almost 200 new companions in the separation range from 1 mas to 8 and presents the catalog of detections, including the first resolved measurements of over a dozen known long-period spectroscopic binaries.Excluding known runaway stars for which no companions are detected, 96 objects in our main sample (δ < 0 • ; H < 7.5) were observed both with PIONIER and NACO/SAM. The fraction of these stars with at least one resolved companion within 200 mas is 0.53. Accounting for known but unresolved spectroscopic or eclipsing companions, the multiplicity fraction at separation ρ < 8 increases to f m = 0.91 ± 0.03. The fraction of luminosity class V stars that have a bound companion reaches 100% at 30 mas while their average number of physically connected companions within 8 is f c = 2.2 ± 0.3. This demonstrates that massive stars form nearly exclusively in multiple systems. The nine non-thermal radio emitters observed by smash+ are all resolved, including the newly discovered pairs HD 168112 and CPD−47 • 2963. This lends strong support to the universality of the wind-wind collision scenario to explain the non-thermal emission from O-type stars.

The star S2 orbiting the compact radio source Sgr A* is a precision probe of the gravitational field around the closest massive black hole (candidate). Over the last 2.7 decades we have monitored the star’s radial velocity and motion on the sky, mainly with the SINFONI and NACO adaptive optics (AO) instruments on the ESO VLT, and since 2017, with the four-telescope interferometric beam combiner instrument GRAVITY. In this Letter we report the first detection of the General Relativity (GR) Schwarzschild Precession (SP) in S2’s orbit. Owing to its highly elliptical orbit (e = 0.88), S2’s SP is mainly a kink between the pre-and post-pericentre directions of motion ≈±1 year around pericentre passage, relative to the corresponding Kepler orbit. The superb 2017−2019 astrometry of GRAVITY defines the pericentre passage and outgoing direction. The incoming direction is anchored by 118 NACO-AO measurements of S2’s position in the infrared reference frame, with an additional 75 direct measurements of the S2-Sgr A* separation during bright states (“flares”) of Sgr A*. Our 14-parameter model fits for the distance, central mass, the position and motion of the reference frame of the AO astrometry relative to the mass, the six parameters of the orbit, as well as a dimensionless parameter fSP for the SP (fSP = 0 for Newton and 1 for GR). From data up to the end of 2019 we robustly detect the SP of S2, δϕ ≈ 12′ per orbital period. From posterior fitting and MCMC Bayesian analysis with different weighting schemes and bootstrapping we find fSP = 1.10 ± 0.19. The S2 data are fully consistent with GR. Any extended mass inside S2’s orbit cannot exceed ≈0.1% of the central mass. Any compact third mass inside the central arcsecond must be less than about 1000 M⊙.

The highly elliptical, 16-year-period orbit of the star S2 around the massive black hole candidate Sgr A✻ is a sensitive probe of the gravitational field in the Galactic centre. Near pericentre at 120 AU ≈ 1400 Schwarzschild radii, the star has an orbital speed of ≈7650 km s−1, such that the first-order effects of Special and General Relativity have now become detectable with current capabilities. Over the past 26 years, we have monitored the radial velocity and motion on the sky of S2, mainly with the SINFONI and NACO adaptive optics instruments on the ESO Very Large Telescope, and since 2016 and leading up to the pericentre approach in May 2018, with the four-telescope interferometric beam-combiner instrument GRAVITY. From data up to and including pericentre, we robustly detect the combined gravitational redshift and relativistic transverse Doppler effect for S2 of z = Δλ / λ ≈ 200 km s−1/c with different statistical analysis methods. When parameterising the post-Newtonian contribution from these effects by a factor f , with f = 0 and f = 1 corresponding to the Newtonian and general relativistic limits, respectively, we find from posterior fitting with different weighting schemes f = 0.90 ± 0.09|stat ± 0.15|sys. The S2 data are inconsistent with pure Newtonian dynamics.

Context. About 2/3 of the Be stars present the so-called V/R variations, a phenomenon characterized by the quasi-cyclic variation in the ratio between the violet and red emission peaks of the H i emission lines. These variations are generally explained by global oscillations in the circumstellar disk forming a one-armed spiral density pattern that precesses around the star with a period of a few years. Aims. This paper presents self-consistent models of polarimetric, photometric, spectrophotometric, and interferometric observations of the classical Be star ζ Tauri. The primary goal is to conduct a critical quantitative test of the global oscillation scenario. Methods. Detailed three-dimensional, NLTE radiative transfer calculations were carried out using the radiative transfer code HDUST. The most up-to-date research on Be stars was used as input for the code in order to include a physically realistic description for the central star and the circumstellar disk. The model adopts a rotationally deformed, gravity darkened central star, surrounded by a disk whose unperturbed state is given by a steady-state viscous decretion disk model. It is further assumed that this disk is in vertical hydrostatic equilibrium. Results. By adopting a viscous decretion disk model for ζ Tauri and a rigorous solution of the radiative transfer, a very good fit of the time-average properties of the disk was obtained. This provides strong theoretical evidence that the viscous decretion disk model is the mechanism responsible for disk formation. The global oscillation model successfully fitted spatially resolved VLTI/AMBER observations and the temporal V/R variations in the Hα and Brγ lines. This result convincingly demonstrates that the oscillation pattern in the disk is a one-armed spiral. Possible model shortcomings, as well as suggestions for future improvements, are also discussed.

Context. It is now generally accepted that the near-infrared excess of Herbig AeBe stars originates in the dust of a circumstellar disk. Aims. The aims of this article are to infer the radial and vertical structure of these disks at scales of order 1 au, and the properties of the dust grains. Methods. The program objects (51 in total) were observed with the H-band (1.6µm) PIONIER/VLTI interferometer. The largest baselines allowed us to resolve (at least partially) structures of a few tenths of an au at typical distances of a few hundred parsecs. Dedicated UBVRIJHK photometric measurements were also obtained. Spectral and 2D geometrical parameters are extracted via fits of a few simple models: ellipsoids and broadened rings with azimuthal modulation. Model bias is mitigated by parallel fits of physical disk models. Sample statistics were evaluated against similar statistics for the physical disk models to infer properties of the sample objects as a group. Results. We find that dust at the inner rim of the disk has a sublimation temperature T sub ≈ 1800K. A ring morphology is confirmed for approximately half the resolved objects; these rings are wide δ r/r ≥ 0.5. A wide ring favors a rim that, on the star-facing side, looks more like a knife edge than a doughnut. The data are also compatible with a the combination of a narrow ring and an inner disk of unspecified nature inside the dust sublimation radius. The disk inner part has a thickness z/r ≈ 0.2, flaring to z/r ≈ 0.5 in the outer part. We confirm the known luminosity-radius relation; a simple physical model is consistent with both the mean luminosity-radius relation and the ring relative width; however, a significant spread around the mean relation is present. In some of the objects we find a halo component, fully resolved at the shortest interferometer spacing, that is related to the HAeBe class.

We measured the radii of 7 low and very low-mass stars using long baseline interferometry with the VLTI interferometer and its VINCI and AMBER near-infrared recombiners. We use these new data, together with literature measurements, to examine the luminosityradius and mass-radius relations for K and M dwarfs. The precision of the new interferometric radii now competes with what can be obtained for double-lined eclipsing binaries. Interferometry provides access to much less active stars, as well as to stars with much better measured distances and luminosities, and therefore complements the information obtained from eclipsing systems. The radii of magnetically quiet late-K to M dwarfs match the predictions of stellar evolution models very well, providing direct confirmation that magnetic activity explains the discrepancy that was recently found for magnetically active eclipsing systems. The radii of the early K dwarfs are reproduced well for a mixing length parameter that approaches the solar value, as qualitatively expected.

Context. For over a decade, the structure of the inner cavity in the transition disk of TW Hydrae has been a subject of debate. Modeling the disk with data obtained at different wavelengths has led to a variety of proposed disk structures. Rather than being inconsistent, the individual models might point to the different faces of physical processes going on in disks, such as dust growth and planet formation. Aims. Our aim is to investigate the structure of the transition disk again and to find to what extent we can reconcile apparent model differences. Methods. A large set of high-angular-resolution data was collected from near-infrared to centimeter wavelengths. We investigated the existing disk models and established a new self-consistent radiative-transfer model. A genetic fitting algorithm was used to automatize the parameter fitting, and uncertainties were investigated in a Bayesian framework. Results. Simple disk models with a vertical inner rim and a radially homogeneous dust composition from small to large grains cannot reproduce the combined data set. Two modifications are applied to this simple disk model: (1) the inner rim is smoothed by exponentially decreasing the surface density in the inner ∼3 AU, and (2) the largest grains (>100 µm) are concentrated towards the inner disk region. Both properties can be linked to fundamental processes that determine the evolution of protoplanetary disks: the shaping by a possible companion and the different regimes of dust-grain growth, respectively. Conclusions. The full interferometric data set from near-infrared to centimeter wavelengths requires a revision of existing models for the TW Hya disk. We present a new model that incorporates the characteristic structures of previous models but deviates in two key aspects: it does not have a sharp edge at 4 AU, and the surface density of large grains differs from that of smaller grains. This is the first successful radiative-transfer-based model for a full set of interferometric data.