Search

We present H-band interferometric observations of the red supergiant (RSG) AZ Cyg made with the Michigan Infra-Red Combiner (MIRC) at the six-telescope Center for High Angular Resolution Astronomy (CHARA) Array. The observations span 5 years (2011-2016), offering insight into the short and long-term evolution of surface features on RSGs. Using a spectrum of AZ Cyg obtained with SpeX on the NASA InfraRed Telescope Facility (IRTF) and synthetic spectra calculated from spherical MARCS, spherical PHOENIX, and SAtlas model atmospheres, we derive $T_{\text{eff}}$ is between $3972 K$ and $4000 K$ and $\log~g$ between $-0.50$ and $0.00$, depending on the stellar model used. Using fits to the squared visibility and Gaia parallaxes we measure its average radius $R=911^{+57}_{-50}~R_{\odot}$. Reconstructions of the stellar surface using our model-independent imaging codes SQUEEZE and OITOOLS.jl show a complex surface with small bright features that appear to vary on a timescale of less than one year and larger features that persist for more than one year. 1D power spectra of these images suggest a characteristic size of $0.52-0.69~R_{\star}$ for the larger, long lived features. This is close to the values of $0.51-0.53~R_{\star}$ derived from 3D RHD models of stellar surfaces. We conclude that interferometric imaging of this star is in line with predictions of 3D RHD models but that short-term imaging is needed to more stringently test predictions of convection in RSGs., 31 pages, 8 tables, 14 figures, accepted to ApJ

Context. Red supergiant stars possess surface features and extended molecular atmospheres. Photospheric convection may be a crucial factor of the levitation of the outer atmospheric layers. However, the mechanism responsible is still poorly understood. Aims. We image the stellar surface of V602 Carinae (V602 Car) to constrain the morphology and contrast of the surface features and of the extended atmospheric layers. Methods. We observed V602 Car with the Very Large Telescope Interferometer (VLTI) PIONIER instrument (1.53-1.78 $\mathrm{\mu}$m) between May and July 2016, and April and July 2019 with different telescope configurations. We compared the image reconstructions with 81 temporal snapshots of 3D radiative-hydrodynamics (RHD) CO$^5$BOLD simulations in terms of contrast and morphology, using the Structural Similarity Index. Results. The interferometric data are compatible with an overall spherical disk of angular diameter 4.4$\pm$0.2 mas, and an extended molecular layer. In 2016, the reconstructed image reveals a bright arc-like feature toward the northern rim of the photospheric surface. In 2019, an arc-like feature is seen at a different orientation and a new peak of emission is detected on the opposite side. The contrasts of the reconstructed surface images are 11%$\pm$2% and 9%$\pm$2% for 2016 and 2019, respectively. The morphology and contrast of the two images are consistent with 3D RHD simulations, within our achieved spatial resolution and dynamic range. The extended molecular layer contributes 10--13% of the total flux with an angular diameter of 6--8 mas. It is present but not clearly visible in the reconstructed images because it is close to the limits of the achieved dynamic range. The presence of the molecular layer is not reproduced by the 3D RHD simulations. Conclusions..., Comment: 16 pages, 13 figures, accepted for publication in Astronomy and Astrophysics

We present models of the inner region of the circumstellar disk of RY Tau which aim to explain our near-infrared ($K$-band: $2.1\,\mu$m) interferometric observations while remaining consistent with the optical to near-infrared portions of the spectral energy distribution. Our sub-milliarcsecond resolution CHARA Array observations are supplemented with shorter baseline, archival data from PTI, KI and VLTI/GRAVITY and modeled using an axisymmetric Monte Carlo radiative transfer code. The $K$-band visibilities are well-fit by models incorporating a central star illuminating a disk with an inner edge shaped by dust sublimation at $0.210\pm0.005\,$au, assuming a viewing geometry adopted from millimeter interferometry ($65^{\circ}$ inclined with a disk major axis position angle of $23^{\circ}$). This sublimation radius is consistent with that expected of Silicate grains with a maximum size of $0.36-0.40\,\mu$m contributing to the opacity and is an order of magnitude further from the star than the theoretical magnetospheric truncation radius. The visibilities on the longest baselines probed by CHARA indicate that we lack a clear line-of-sight to the stellar photosphere. Instead, our analysis shows that the central star is occulted by the disk surface layers close to the sublimation rim. While we do not see direct evidence of temporal variability in our multi-epoch CHARA observations, we suggest the aperiodic photometric variability of RY~Tau is likely related temporal and/or azimuthal variations in the structure of the disk surface layers., Comment: Accepted for publication in The Astrophysical Journal

We present very detailed images of the photosphere of an AGB star obtained with the PIONIER instrument, installed at the Very Large Telescope Interferometer (VLTI). The images show a well defined stellar disc populated by a few convective patterns. Thanks to the high precision of the observations we are able to derive the contrast and granulation horizontal scale of the convective pattern for the first time in a direct way. Such quantities are then compared with scaling relations between granule size, effective temperature, and surface gravity that are predicted by simulations of stellar surface convection., SCOPUS: ar.k, info:eu-repo/semantics/published

We aim at detecting and characterizing the main-sequence companions of a sample of known and suspected Galactic binary Cepheids. We used the multi-telescope interferometric combiners MIRC and PIONIER to detect and measure the astrometric positions of the high-contrast companions orbiting 16 bright Galactic Cepheids. We made use of the CANDID algorithm to search for the companions and set detection limits. We also present new high-precision RVs which were used to fit radial pulsation and orbital velocities. We report the detection of the companions orbiting the Cepheids U Aql, BP Cir, and S Mus, possible detections for FF Aql, Y Car, BG Cru, X Sgr, V350 Sgr, and V636 Sco, while no component is detected around U Car, YZ Car, T Mon, R Mus, S Nor, W Sgr and AH Vel. For U Aql and S Mus, we performed a preliminary orbital fit combining astrometry with new high-precision RVs, providing the full set of orbital and pulsation parameters. Assuming the distance, we estimated preliminary masses of M(U Aql) = 4.97+/-0.62Msun and M(S Mus) = 4.63+/-0.99Msun. For YZ Car, W Sgr, V350 Sgr, and V636 Sco, we revised the spectroscopic orbits using new high-precision RVs, while we updated the pulsation parameters for BP Cir, BG Cru, S Nor and AH Vel. Our interferometric observations also provide measurements of the angular diameters, that can be used in a Baade-Wesselink type analysis. We have now several astrometric detections of Cepheid companions. When RVs of the companions will be available, such systems will provide accurate and independent masses and distances. Orbital parallaxes with an accuracy, 29 pages, 11 figures, 23 tables. Accepted for publication in A&A

We recently used archival and newly obtained data from the Navy Precision Optical Interferometer to measure the fundamental properties of 87 stars. The sample consisted of 5 dwarfs, 3 subgiants, 69 giants, 3 bright giants, and 7 supergiants, and spanned a wide range of spectral classes from B to M. We combined our angular diameters with photometric and distance information from the literature to determine each star’s physical radius, effective temperature, bolometric flux, luminosity, mass, and age. Several dozen of the stars have visibility curves sampled down to the first null, where the visibilities drop to zero. Here we present preliminary results showing limb-darkening fits for the five zero crossing stars that have the best coverage of the second lobe.

We present comprehensive models of Herbig Ae star, HD 142666, which aim to simultaneously explain its spectral energy distribution (SED) and near-infrared (NIR) interferometry. Our new sub-milliarcsecond resolution CHARA (CLASSIC and CLIMB) interferometric observations, supplemented with archival shorter baseline data from VLTI/PIONIER and the Keck Interferometer, are modeled using centro-symmetric geometric models and an axisymmetric radiative transfer code. CHARA's 330 m baselines enable us to place strong constraints on the viewing geometry, revealing a disk inclined at 58 degrees from face-on with a 160 degree major axis position angle. Disk models imposing vertical hydrostatic equilibrium provide poor fits to the SED. Models accounting for disk scale height inflation, possibly induced by turbulence associated with magneto-rotational instabilities, and invoking grain growth to >1 micron size in the disk rim are required to simultaneously reproduce the SED and measured visibility profile. However, visibility residuals for our best model fits to the SED indicate the presence of unexplained NIR emission, particularly along the apparent disk minor axis, while closure phase residuals indicate a more centro-symmetric emitting region. In addition, our inferred 58 degree disk inclination is inconsistent with a disk-based origin for the UX Ori-type variability exhibited by HD 142666. Additional complexity, unaccounted for in our models, is clearly present in the NIR-emitting region. We propose the disk is likely inclined toward a more edge-on orientation and/or an optically thick outflow component also contributes to the NIR circumstellar flux., Accepted for publication in The Astrophysical Journal

Nearby, active stars with relatively rapid rotation and large starspot structures offer the opportunity to compare interferometric, spectroscopic, and photometric imaging techniques. In this paper, we image a spotted star with three different methods for the first time. The giant primary star of the RS Canum Venaticorum binary $\sigma$ Geminorum ($\sigma$ Gem) was imaged for two epochs of interferometric, high-resolution spectroscopic, and photometric observations. The light curves from the reconstructions show good agreement with the observed light curves, supported by the longitudinally-consistent spot features on the different maps. However, there is strong disagreement in the spot latitudes across the methods., Comment: Accepted to ApJ, 19 pages, 16 figures, 5 tables

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 one au, and the properties of the dust grains. Methods. The program objects (51 in total) were observed with the H-band (1.6micron) 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 Tsub~1800K. A ring morphology is confirmed for approximately half the resolved objects; these rings are wide delta_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., 41 pages. 23 figures. Four appendices

Context. Red giant branch (RGB) stars are very bright objects in galaxies and are often used as standard candles. Interferometry is the ideal tool to characterize the dynamics and morphology of their atmospheres. Aims. We aim at precisely characterising the surface dynamics of a sample of RGB stars. Methods. We obtained interferometric observations for three RGB stars with the MIRC instrument mounted at the CHARA interfer- ometer. We looked for asymmetries on the stellar surfaces using limb-darkening models. Results. We measured the apparent diameters of HD197989 (Epsilon Cyg) = 4.61+-0.02 mas, HD189276 (HR7633) = 2.95+-0.01 mas, and HD161096 (Beta Oph) = 4.43+-0.01 mas. We detected departures from the centrosymmetric case for all three stars with the tendency of a greater effect for lower logg of the sample. We explored the causes of this signal and conclude that a possible explanation to the interferometric signal is the convection-related and/or the magnetic-related surface activity. However, it is necessary to monitor these stars with new observations, possibly coupled with spectroscopy, in order to firmly establish the cause., Accepted for publication as a Letter in Astronomy and Astrophysics, section 1. Letters to the Editor. The official date of acceptance is 06/03/2017

Context. T Tauri stars are low-mass young stars whose disks provide the setting for planet formation. Despite this, their structure is poorly understood. We present new infrared interferometric observations of the SU Aurigae circumstellar environment that offer resolution that is three times higher and a better baseline position angle coverage than previous observations. Aims. We aim to investigate the characteristics of the circumstellar material around SU Aur, constrain the disk geometry, composition and inner dust rim structure. Methods. The CHARA array offers unique opportunities for long baseline observations, with baselines up to 331 m. Using the CLIMB three-telescope combiner in the K-band allows us to measure visibilities as well as closure phase. We undertook image reconstruction for model-independent analysis, and fitted geometric models such as Gaussian and ring distributions. Additionally, the fitting of radiative transfer models constrain the physical parameters of the disk. For the first time, a dusty disk wind is introduced to the radiative transfer code TORUS to model protoplanetary disks. Our implementation is motivated by theoretical models of dusty disk winds, where magnetic field lines drive dust above the disk plane close to the sublimation zone. Results. Image reconstruction reveals an inclined disk with slight asymmetry along its minor-axis, likely due to inclination effects obscuring the inner disk rim through absorption of incident star light on the near-side and thermal re-emission and scattering of the far-side. Geometric modelling of a skewed ring finds the inner rim at 0.17 ± 0.02 au with an inclination of 50.9 ± 1.0° and minor axis position angle 60.8 ± 1.2°. Radiative transfer modelling shows a flared disk with an inner radius at 0.18 au which implies a grain size of 0.4 μm assuming astronomical silicates and a scale height of 15.0 at 100 au. Among the tested radiative transfer models, only the dusty disk wind successfully accounts for the K-band excess by introducing dust above the mid-plane.

We present near-infrared and optical emission-line and stellar kinematics of the Seyfert 2 galaxy Mrk 573 using the Near-Infrared Field Spectrograph (NIFS) at Gemini North and Dual Imaging Spectrograph (DIS) at Apache Point Observatory, respectively. By obtaining full kinematic maps of the infrared ionized and molecular gas and stellar kinematics in a 700 x 2100 pc^2 circumnuclear region of Mrk 573, we find that kinematics within the Narrow-Line Region (NLR) are largely due to a combination of both rotation and in situ acceleration of material originating in the host disk. Combining these observations with large-scale, optical long-slit spectroscopy that traces ionized gas emission out to several kpcs, we find that rotation kinematics dominate the majority of the gas. We find that outflowing gas extends to distances less than 1 kpc, suggesting that outflows in Seyfert galaxies may not be powerful enough to evacuate their entire bulges., Comment: 19 pages, 17 figures, submitted for publication in ApJ

We report the discovery of a circumbinary disk around the Herbig Ae/Be system v892 Tau. Our detailed mid-infrared images were made using segment-tilting interferometry on the Keck-1 Telescope and reveal an asymmetric disk inclined at ~60 degs with an inner hole diameter of 250 mas (35 AU), approximately 5X larger than the apparent separation of the binary components. In addition, we report a new measurement along the binary orbit using near-infrared Keck aperture masking, allowing a crude estimate of orbital parameters and the system mass for the first time. The size of the inner hole appears to be consistent with the minimum size prediction from tidal truncation theory, bearing a resemblance to the recently unmasked binary CoKu Tau/4. Our results have motivated a re-analysis of the system spectral energy distribution, concluding the luminosity of this system has been severely underestimated. With further study and monitoring, v892 Tau should prove a powerful testing ground for both predictions of dynamical models for disk-star interactions in young systems with gas-rich disks and for calibrations of pre-main-sequence tracks for intermediate-mass stars., Accepted for publication in Astrophysical Journal Letters

The rapidly rotating Be star phi Persei was spun up by mass and angular momentum transfer from a now stripped-down, hot subdwarf companion. Here we present the first high angular resolution images of phi Persei made possible by new capabilities in longbaseline interferometry at near-IR and visible wavelengths. We observed phi Persei with the MIRC and VEGA instruments of the CHARA Array. Additional MIRC-only observations were performed to track the orbital motion of the companion, and these were fit together with new and existing radial velocity measurements of both stars to derive the complete orbital elements and distance. The hot subdwarf companion is clearly detected in the near-IR data at each epoch of observation with a flux contribution of 1.5% in the H band, and restricted fits indicate that its flux contribution rises to 3.3% in the visible. A new binary orbital solution is determined by combining the astrometric and radial velocity measurements. The derived stellar masses are 9.6+-0.3Msol and 1.2+-0.2Msol for the Be primary and subdwarf secondary, respectively. The inferred distance (186 +- 3 pc), kinematical properties, and evolutionary state are consistent with membership of phi Persei in the alpha Per cluster. From the cluster age we deduce significant constraints on the initial masses and evolutionary mass transfer processes that transformed the phi Persei binary system. The interferometric data place strong constraints on the Be disk elongation, orientation, and kinematics, and the disk angular momentum vector is coaligned with and has the same sense of rotation as the orbital angular momentum vector. The VEGA visible continuum data indicate an elongated shape for the Be star itself, due to the combined effects of rapid rotation, partial obscuration of the photosphere by the circumstellar disk, and flux from the bright inner disk., 16 pages, 6 figures, 1 Annex

To measure the stellar and orbital properties of the metal-poor RS CVn binary o Draconis (o Dra), we directly detect the companion using interferometric observations obtained with the Michigan InfraRed Combiner at Georgia State University's Center for High Angular Resolution Astronomy (CHARA) Array. The H-band flux ratio between the primary and secondary stars is the highest confirmed flux ratio (370 +/- 40) observed with long-baseline optical interferometry. These detections are combined with radial velocity data of both the primary and secondary stars, including new data obtained with the Tillinghast Reflector Echelle Spectrograph on the Tillinghast Reflector at the Fred Lawrence Whipple Observatory and the 2-m Tennessee State University Automated Spectroscopic Telescope at Fairborn Observatory. We determine an orbit from which we find model-independent masses and ages of the components (M_A = 1.35 +\- 0.05 M_Sun, M_B = 0.99 +\- 0.02 M_Sun, system age = 3.0 -\+ 0.5 Gyr). An average of a 23-year light curve of o Dra from the Tennessee State University Automated Photometric Telescope folded over the orbital period newly reveals eclipses and the quasi-sinusoidal signature of ellipsoidal variations. The modeled light curve for our system's stellar and orbital parameters confirm these ellipsoidal variations due to the primary star partially filling its Roche lobe potential, suggesting most of the photometric variations are not due to stellar activity (starspots). Measuring gravity darkening from the average light curve gives a best-fit of beta = 0.07 +\- 0.03, a value consistent with conventional theory for convective envelope stars. The primary star also exhibits an anomalously short rotation period, which, when taken with other system parameters, suggests the star likely engulfed a low-mass companion that had recently spun-up the star., Comment: 14 pages, 13 figures, Accepted to ApJ

Long-baseline interferometry is an important technique to spatially resolve binary or multiple systems in close orbits. By combining several telescopes together and spectrally dispersing the light, it is possible to detect faint components around bright stars. Aims. We provide a rigorous and detailed method to search for high-contrast companions around stars, determine the detection level, and estimate the dynamic range from interferometric observations. We developed the code CANDID (Companion Analysis and Non-Detection in Interferometric Data), a set of Python tools that allows us to search systematically for point-source, high-contrast companions and estimate the detection limit. The search pro- cedure is made on a N x N grid of fit, whose minimum needed resolution is estimated a posteriori. It includes a tool to estimate the detection level of the companion in the number of sigmas. The code CANDID also incorporates a robust method to set a 3{\sigma} detection limit on the flux ratio, which is based on an analytical injection of a fake companion at each point in the grid. We used CANDID to search for the companions around the binary Cepheids V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul. First, we showed that our previous discoveries of the components orbiting V1334 Cyg and AX Cir were detected at > 13 sigmas. The companion around AW Per is detected at more than 15 sigmas with a flux ratio of f = 1.22 +/- 0.30 %. We made a possible detection of the companion orbiting RT Aur with f = 0.22 +/- 0.11 %. It was detected at 3.8{\sigma} using the closure phases only, and so more observations are needed to confirm the detection. We also set the detection limit for possible undetected companions. We found that there is no companion with a spectral type earlier than B7V, A5V, F0V, B9V, A0V, and B9V orbiting V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul, respectively., Comment: Accepted for publication in Astronomy and Astrophysics

The authors acknowledge funding from the NSF through award AST-0807577 to the University of Michigan, as well as funding by the Australian Research Council and the "Lendület Fiatal Kutatói Program" of the Hungarian Academy of Sciences for the initial work on these objects. Two red supergiants (RSGs) of the Per OB1 association, RS Per and T Per, have been observed in the H band using the Michigan Infra-Red Combiner (MIRC) instrument at the CHARA array. The data show clear evidence of a departure from circular symmetry. We present here new techniques specially developed to analyze such cases, based on state-of-the-art statistical frameworks. The stellar surfaces are first modeled as limb-darkened disks based on SATLAS models that fit both MIRC interferometric data and publicly available spectrophotometric data. Bayesian model selection is then used to determine the most probable number of spots. The effective surface temperatures are also determined and give further support to the recently derived hotter temperature scales of RSGs. The stellar surfaces are reconstructed by our model-independent imaging code SQUEEZE, making use of its novel regularizer based on Compressed Sensing theory. We find excellent agreement between the model-selection results and the reconstructions. Our results provide evidence for the presence of near-infrared spots representing about 3%-5% of the stellar flux. Publisher PDF

The Event Horizon Telescope (EHT) is a project to assemble a Very Long Baseline Interferometry (VLBI) network of mm wavelength dishes that can resolve strong field General Relativistic signatures near a supermassive black hole. As planned, the EHT will include enough dishes to enable imaging of the predicted black hole "shadow", a feature caused by severe light bending at the black hole boundary. The center of M87, a giant elliptical galaxy, presents one of the most interesting EHT targets as it exhibits a relativistic jet, offering the additional possibility of studying jet genesis on Schwarzschild radius scales. Fully relativistic models of the M87 jet that fit all existing observational constraints now allow horizon-scale images to be generated. We perform realistic VLBI simulations of M87 model images to examine detectability of the black shadow with the EHT, focusing on a sequence of model images with a changing jet mass load radius. When the jet is launched close to the black hole, the shadow is clearly visible both at 230 and 345 GHz. The EHT array with a resolution of 20-30$\mu$as resolution ($\sim$2-4 Schwarzschild radii) is able to image this feature independent of any theoretical models and we show that imaging methods used to process data from optical interferometers are applicable and effective for EHT data sets. We demonstrate that the EHT is also capable of tracing real-time structural changes on a few Schwarzschild radii scales, such as those implicated by VHE flaring activity of M87. While inclusion of ALMA in the EHT is critical for shadow imaging, generally the array is robust against loss of a station., Comment: 10 pages, 7 figures, accepted for publication in ApJ

Epsilon Aurigae (epsilon Aur) is a visually bright, eclipsing binary star system with a period of 27.1 years. The cause of each 18-month-long eclipse has been a subject of controversy for nearly 190 years because the companion has hitherto been undetectable. The orbital elements imply that the opaque object has roughly the same mass as the visible component, which for much of the last century was thought to be an F-type supergiant star with a mass of approximately 15M[symbol:see text] (M[symbol:see text], mass of the Sun). The high mass-to-luminosity ratio of the hidden object was originally explained by supposing it to be a hyperextended infrared star or, later, a black hole with an accretion disk, although the preferred interpretation was as a disk of opaque material at a temperature of approximately 500 K, tilted to the line of sight and with a central opening. Recent work implies that the system consists of a low-mass (2.2M[symbol:see text]-3.3M[symbol:see text]) visible F-type star, with a disk at 550 K that enshrouds a single B5V-type star. Here we report interferometric images that show the eclipsing body moving in front of the F star. The body is an opaque disk and appears tilted as predicted. Adopting a mass of 5.9M[symbol:see text] for the B star, we derive a mass of approximately (3.6 +/- 0.7)M[symbol:see text] for the F star. The disk mass is dynamically negligible; we estimate it to contain approximately 0.07M[symbol:see text] (M[symbol:see text], mass of the Earth) if it consists purely of dust.

We present the results of long-baseline optical interferometry observations using the Precision Astronomical Visual Observations (PAVO) beam combiner at the Center for High Angular Resolution Astronomy (CHARA) Array to measure the angular sizes of three bright Kepler stars: {\theta} Cygni, and both components of the binary system 16 Cygni. Supporting infrared observations were made with the Michigan Infrared Combiner (MIRC) and Classic beam combiner, also at the CHARA Array. We find limb-darkened angular diameters of 0.753+/-0.009 mas for {\theta} Cyg, 0.539+/-0.007 mas for 16 Cyg A and 0.490+/-0.006 mas for 16 Cyg B. The Kepler Mission has observed these stars with outstanding photometric precision, revealing the presence of solar-like oscillations. Due to the brightness of these stars the oscillations have exceptional signal-to-noise, allowing for detailed study through asteroseismology, and are well constrained by other observations. We have combined our interferometric diameters with Hipparcos parallaxes, spectrophotometric bolometric fluxes and the asteroseismic large frequency separation to measure linear radii ({\theta} Cyg: 1.48+/-0.02 Rsun, 16 Cyg A: 1.22+/-0.02 Rsun, 16 Cyg B: 1.12+/-0.02 Rsun), effective temperatures ({\theta} Cyg: 6749+/-44 K, 16 Cyg A: 5839+/-42 K, 16 Cyg B: 5809+/-39 K), and masses ({\theta} Cyg: 1.37+/-0.04 Msun, 16 Cyg A: 1.07+/-0.05 Msun, 16 Cyg B: 1.05+/-0.04 Msun) for each star with very little model dependence. The measurements presented here will provide strong constraints for future stellar modelling efforts., Comment: 9 pages, 4 figures, and 5 tables, accepted for publication in Monthly Notices of the Royal Astronomical Society

We present the first high angular resolution observations in the nearinfrared H-band (1.6 microns) of the Luminous Blue Variable star P Cygni. We obtained six-telescope interferometric observations with the CHARA Array and the MIRC beam combiner. These show that the spatial flux distribution is larger than expected for the stellar photosphere. A two component model for the star (uniform disk) plus a halo (two-dimensional Gaussian) yields an excellent fit of the observations, and we suggest that the halo corresponds to flux emitted from the base of the stellar wind. This wind component contributes about 45% of the H-band flux and has an angular FWHM = 0.96 mas, compared to the predicted stellar diameter of 0.41 mas. We show several images reconstructed from the interferometric visibilities and closure phases, and they indicate a generally spherical geometry for the wind. We also obtained near-infrared spectrophotometry of P Cygni from which we derive the flux excess compared to a purely photospheric spectral energy distribution. The H-band flux excess matches that from the wind flux fraction derived from the two component fits to the interferometry. We find evidence of significant near-infrared flux variability over the period from 2006 to 2010 that appears similar to the variations in the H-alpha emission flux from the wind. Future interferometric observations may be capable of recording the spatial variations associated with temporal changes in the wind structure., 34 pages, 5 tables, 8 figures, accepted to ApJ

We aim at determining the masses of Cepheids in binary systems, as well as their geometric distances and the flux contribution of the companions. The combination of interferometry with spectroscopy will offer a unique and independent estimate of the Cepheid masses. Using long-baseline interferometry at visible and infrared wavelengths, it is possible to spatially resolve binary systems containing a Cepheid down to milliarcsecond separations. Based on the resulting visual orbit and radial velocities, we can then derive the fundamental parameters of these systems, particularly the masses of the components and the geometric distance. We therefore performed interferometric observations of the first-overtone mode Cepheid V1334 Cyg with the CHARA/MIRC combiner. We report the first detection of a Cepheid companion using long-baseline interferometry. We detect the signature of a companion orbiting V1334 Cyg at two epochs. We measure a flux ratio between the companion and the Cepheid f = 3.10+/-0.08%, giving an apparent magnitude mH = 8.47+/-0.15mag. The combination of interferometric and spectroscopic data have enabled the unique determination of the orbital elements: P = 1938.6+/-1.2 days, Tp = 2 443 616.1+/-7.3, a = 8.54+/-0.51mas, i = 124.7+/-1.8{\deg}, e = 0.190+/-0.013, {\omega} = 228.7+/-1.6{\deg}, and {\Omega} = 206.3+/-9.4{\deg}. We derive a minimal distance d ~ 691 pc, a minimum mass for both stars of 3.6 Msol, with a spectral type earlier than B5.5V for the companion star. Our measured flux ratio suggests that radial velocity detection of the companion using spectroscopy is within reach, and would provide an orbital parallax and model-free masses., Comment: Published in A&A

One of the great challenges in understanding stars is measuring their masses. The best methods for measuring stellar masses include binary interaction, asteroseismology and stellar evolution models, but these methods are not ideal for red giant and supergiant stars. In this work, we propose a novel method for inferring stellar masses of evolved red giant and supergiant stars using interferometric and spectrophotometric observations combined with spherical model stellar atmospheres to measure what we call the stellar mass index, defined as the ratio between the stellar radius and mass. The method is based on the correlation between different measurements of angular diameter, used as a proxy for atmospheric extension, and fundamental stellar parameters. For a given star, spectrophotometry measures the Rosseland angular diameter while interferometric observations generally probe a larger limb-darkened angular diameter. The ratio of these two angular diameters is proportional to the relative extension of the stellar atmosphere, which is strongly correlated to the star's effective temperature, radius and mass. We show that these correlations are strong and can lead to precise measurements of stellar masses., 7 pages, 3 figures, accepted for publication in ApJ

Sunspots are cool areas caused by strong surface magnetic fields inhibiting convection. Moreover, strong magnetic fields can alter the average atmospheric structure, degrading our ability to measure stellar masses and ages. Stars more active than the Sun have more and stronger dark spots than in the solar case, including on the rotational pole itself. Doppler imaging, which has so far produced the most detailed images of surface structures on other stars than the Sun, cannot always distinguish the hemisphere in which the starspots are located, especially in the equatorial region and if the data quality is not optimal. This leads to problems in investigating the north-south distribution of starspot active latitudes (those latitudes with more spot activity), which are crucial constraints of dynamo theory. Polar spots, inferred only from Doppler tomography, could plausibly be observational artifacts, casting some doubt on their very existence. Here we report imaging of the old, magnetically-active star $\zeta$ Andromedae using long-baseline infrared interferometry. In our data, a dark polar spot is seen in each of two epochs, while lower-latitude spot structures in both hemispheres do not persist between observations revealing global starspot asymmetries. The north-south symmetry of active latitudes observed on the Sun is absent on $\zeta$ And, which hosts global spot patterns that cannot be produced by solar-type dynamos., Comment: Published in Nature, 33 pages, 8 figures, 5 tables

Optical and infrared interferometers definitively established that the photometric standard Vega (alpha Lyrae) is a rapidly rotating star viewed nearly pole-on. Recent independent spectroscopic analyses could not reconcile the inferred inclination angle with the observed line profiles, preferring a larger inclination. In order to resolve this controversy, we observed Vega using the six-beam Michigan Infrared Combiner on the Center for High Angular Resolution Astronomy Array. With our greater angular resolution and dense (u,v)-coverage, we find Vega is rotating less rapidly and with a smaller gravity darkening coefficient than previous interferometric results. Our models are compatible with low photospheric macroturbulence and also consistent with the possible rotational period of ~0.71 days recently reported based on magnetic field observations. Our updated evolutionary analysis explicitly incorporates rapid rotation, finding Vega to have a mass of 2.15+0.10_-0.15 Msun and an age 700-75+150 Myrs, substantially older than previous estimates with errors dominated by lingering metallicity uncertainties (Z=0.006+0.003-0.002)., Comment: Accepted for publication in ApJ Letters

Using CHARA and VLTI near-infrared spectro-interferometry with hectometric baseline lengths (up to 330m) and with high spectral resolution (up to 12000), we studied the gas distribution and kinematics around two classical Be stars. The combination of high spatial and spectral resolution achieved allows us to constrain the gas velocity field on scales of a few stellar radii and to obtain, for the first time in optical interferometry, a dynamical mass estimate using the position-velocity analysis technique known from radio astronomy. For our first target star, Beta Canis Minoris, we model the H+K-band continuum and Br-gamma line geometry with a near-critical rotating stellar photosphere and a geometrically thin equatorial disk. Testing different disk rotation laws, we find that the disk is in Keplerian rotation (v(r)=r^(-0.5+/-0.1)) and derive the disk position angle (140+/-1.7 deg), inclination (38.5+/-1 deg), and the mass of the central star (3.5+/-0.2 M_sun). As a second target star, we observed the prototypical Be star Zeta Tauri and spatially resolved the Br�� emission as well as nine transitions from the hydrogen Pfund series (Pf14-Pf22). Comparing the spatial origin of the different line transitions, we find that the Brackett (Br-gamma), Pfund (Pf14-17), and Balmer (H-alpha) lines originate from different stellocentric radii (R_cont < R_Pf < R_Br-gamma = R_H-alpha), which we can reproduce with an LTE line radiative transfer computation. Discussing different disk-formation scenarios, we conclude that our constraints are inconsistent with wind compression models predicting a strong outflowing velocity component, but support viscous decretion disk models, where the Keplerian-rotating disk is replenished with material from the near-critical rotating star., 16 pages, 15 figures, accepted by ApJ

Context. The star HD 87643, exhibiting the "B[e] phenomenon", has one of the most extreme infrared excesses for this object class. It harbours a large amount of both hot and cold dust, and is surrounded by an extended reflection nebula. Aims. One of our major goals was to investigate the presence of a companion in HD87643. In addition, the presence of close dusty material was tested through a combination of multi-wavelength high spatial 5Aresolution observations. Methods. We observed HD 87643 with high spatial resolution techniques, using the near-IR AMBER/VLTI interferometer with baselines ranging from 60 m to 130 m and the mid-IR MIDI/VLTI interferometer with baselines ranging from 25 m to 65 m. These observations are complemented by NACO/VLT adaptive-optics-corrected images in the K and L-bands, ESO-2.2m optical Wide-Field Imager large-scale images in the B, V and R-bands, Results. We report the direct detection of a companion to HD 87643 by means of image synthesis using the AMBER/VLTI instrument. The presence of the companion is confirmed by the MIDI and NACO data, although with a lower confidence. The companion is separated by ~ 34 mas with a roughly north-south orientation. The period must be large (several tens of years) and hence the orbital parameters are not determined yet. Binarity with high eccentricity might be the key to interpreting the extreme characteristics of this system, namely a dusty circumstellar envelope around the primary, a compact dust nebulosity around the binary system and a complex extended nebula witnessing past violent ejections.

We present the first resolved images of the eclipsing binary Beta Lyrae, obtained with the CHARA Array interferometer and the MIRC combiner in the H band. The images clearly show the mass donor and the thick disk surrounding the mass gainer at all six epochs of observation. The donor is brighter and generally appears elongated in the images, the first direct detection of photospheric tidal distortion due to Roche-lobe filling. We also confirm expectations that the disk component is more elongated than the donor and is relatively fainter at this wavelength. Image analysis and model fitting for each epoch were used for calculating the first astrometric orbital solution for Beta Lyrae, yielding precise values for the orbital inclination and position angle. The derived semi-ma jor axis also allows us to estimate the distance of Beta Lyrae; however, systematic differences between the models and the images limit the accuracy of our distance estimate to about 15%. To address these issues, we will need a more physical, self-consistent model to account for all epochs as well as the multi-wavelength information from the eclipsing light curves., 13 Pages, 3 Figures; Accepted by ApJL

We present the work developed within the science team of the Very Large Telescope Interferometer Spectro-Imager (VSI) during the Phase A studies. VSI aims at delivering ~ 1 milliarcsecond resolution data cubes in the near-infrared, with several spectral resolutions up to 12 000, by combining up to 8 VLTI telescopes. In the design of an instrument, the science case plays a central role by supporting the instrument construction decision, defining the top-level requirements and balancing design options. The overall science philosophy of VSI was that of a general user instrument serving a broad community. The science team addressed themes which included several areas of astrophysics and illustrated specific modes of operation of the instrument: a) YSO disks and winds; b) Multiplicity of young stars; c) Exoplanets; d) Debris disks; e) Stellar surface imaging; f) The environments of evolved stars; g) AGN tori; h) AGN's Broad Line Region; i) Supermassive black-holes; and j) Microlensing. The main conclusions can be summarized as follows: a) The accessible targets and related science are extremely sensitive to the instrument limiting magnitude; the instrument should be optimized for sensitivity and have its own fringe tracker. b) Most of the science cases are readily achievable with on-axis fringe tracking, off-axis fringe tracking enabling extra science. c) In most targets (YSOs, evolved stars and AGNs), the interpretation and analysis of circumstellar/nuclear dust morphology requires direct access to the gas via spectral resolved studies of emission lines, requiring at least a spectral resolution of 2 500. d) To routinely deliver images at the required sensitivity, the number of telescopes in determinant, with 6 telescopes being favored. e) The factorial increase in the number of closure phases and visibilities, gained in a single observation, makes massive surveys of parameters and related science for the first time possible. f) High dynamic range imaging and very high dynamic range differential closure phase are possible allowing the study of debris disks and characterization of pegasides. g) Spectro-imaging in the near-infrared is highly complementary to ALMA, adaptive optics and interferometric imaging in the thermal infrared.

Fringe tracking in interferometers is typically analyzed with the implicit assumption that there is a single phase associated with each telescope in the array. If the telescopes have apertures significantly larger than r0 and only partial adaptive optics correction, then the phase measured by a fringe sensor may differ significantly from the "piston" component of the aperture phase. In some cases, speckle noise will cause "branch points" in the measured phase as a function of time, causing large and sudden jumps in the phase. We present simulations showing these effects in order to understand their implications for the design of fringe tracking algorithms., 9 pages, to be published in Proc. SPIE conference 7013 "Optical and Infrared Interferometry", Schoeller, Danchi, and Delplancke (eds.)

We report on a total of 106 nights of optical interferometric observations of the $\epsilon$ Aurigae system taken during the last 14 years by four beam combiners at three different interferometric facilities. This long sequence of data provides an ideal assessment of the system prior to, during, and after the recent 2009-2011 eclipse. We have reconstructed model-independent images from the 10 in-eclipse epochs which show that a disk-like object is indeed responsible for the eclipse. Using new 3D, time-dependent modeling software, we derive the properties of the F-star (diameter, limb darkening), determine previously unknown orbital elements ($\Omega$, $i$), and access the global structures of the optically thick portion of the eclipsing disk using both geometric models and approximations of astrophysically relevant density distributions. These models may be useful in future hydrodynamical modeling of the system. Lastly, we address several outstanding research questions including mid-eclipse brightening, possible shrinking of the F-type primary, and any warps or sub-features within the disk., Comment: 105 pages, 57 figures. This is an author-created, un-copyedited version of an article accepted for publication in Astrophysical Journal Supplement Series. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it

To measure the properties of both components of the RS CVn binary sigma Geminorum (sigma Gem), we directly detect the faint companion, measure the orbit, obtain model-independent masses and evolutionary histories, detect ellipsoidal variations of the primary caused by the gravity of the companion, and measure gravity darkening. We detect the companion with interferometric observations obtained with the Michigan InfraRed Combiner (MIRC) at Georgia State University's Center for High Angular Resolution Astronomy (CHARA) Array with a primary-to-secondary H-band flux ratio of 270+/-70. A radial velocity curve of the companion was obtained with spectra from the Tillinghast Reflector Echelle Spectrograph (TRES) on the 1.5-m Tillinghast Reflector at Fred Lawrence Whipple Observatory (FLWO). We additionally use new observations from the Tennessee State University Automated Spectroscopic and Photometric Telescopes (AST and APT, respectively). From our orbit, we determine model-independent masses of the components (M_1 = 1.28 +/- 0.07 M_Sun, M_2 = 0.73 +/- 0.03 M_Sun), and estimate a system age of 5 -/+ 1 Gyr. An average of the 27-year APT light curve of sigma Gem folded over the orbital period (P = 19.6027 +/- 0.0005 days) reveals a quasi-sinusoidal signature, which has previously been attributed to active longitudes 180 deg apart on the surface of sigma Gem. With the component masses, diameters, and orbit, we find that the predicted light curve for ellipsoidal variations due to the primary star partially filling its Roche lobe potential matches well with the observed average light curve, offering a compelling alternative explanation to the active longitudes hypothesis. Measuring gravity darkening from the light curve gives beta < 0.1, a value slightly lower than that expected from recent theory., Accepted to ApJ, 11 pages, 6 figures, 8 tables

The highly eccentric Be binary system {delta} Sco reached periastron during early 2011 July, when the distance between the primary and secondary was a few times the size of the primary disk in the H band. This opened a window of opportunity to study how the gaseous disks around Be stars respond to gravitational disturbance. We first refine the binary parameters with the best orbital phase coverage data from the Navy Precision Optical Interferometer. Then we present the first imaging results of the disk after the periastron, based on seven nights of five telescope observations with the MIRC combiner at the CHARA array. We found that the disk was inclined 27.{sup 0}6 {+-} 6.{sup 0}0 from the plane of the sky, had a half-light radius of 0.49 mas (2.2 stellar radii), and consistently contributed 71.4% {+-} 2.7% of the total flux in the H band from night to night, suggesting no ongoing transfer of material into the disk during the periastron. The new estimation of the periastron passage is UT 2011 July 3 07:00 {+-} 4:30. Re-analysis of archival VLTI-AMBER interferometry data allowed us to determine the rotation direction of the primary disk, constraining it to be inclined either {approx}119more » Degree-Sign or {approx}171 Degree-Sign relative to the orbital plane of the binary system. We also detect inner disk asymmetries that could be explained by spot-like emission with a few percent of the disk total flux moving in Keplerian orbits, although we lack sufficient angular resolution to be sure of this interpretation and cannot yet rule out spiral density waves or other more complicated geometries.« less

Wolf-Rayet stars represent one of the final stages of massive stellar evolution. Relatively little is known about this short-lived phase and we currently lack reliable mass, distance, and binarity determinations for a representative sample. Here we report the first visual orbit for WR 140(=HD193793), a WC7+O5 binary system known for its periodic dust production episodes triggered by intense colliding winds near periastron passage. The IOTA and CHARA interferometers resolved the pair of stars in each year from 2003--2009, covering most of the highly-eccentric, 7.9 year orbit. Combining our results with the recent improved double-line spectroscopic orbit of Fahed et al. (2011), we find the WR 140 system is located at a distance of 1.67 +/- 0.03 kpc, composed of a WR star with M_WR = 14.9 +/- 0.5 Msun and an O star with M_O = 35.9 +/- 1.3 Msun. Our precision orbit yields key parameters with uncertainties times 6 smaller than previous work and paves the way for detailed modeling of the system. Our newly measured flux ratios at the near-infrared H and Ks bands allow an SED decomposition and analysis of the component evolutionary states., Comment: Complete OIFITS dataset included via Data Conservancy Project

The close environment of Herbig stars starts to be revealed step by step and it appears to be quite complex. Many physical phenomena interplay: the dust sublimation causing a puffed-up inner rim, a dusty halo, a dusty wind or an inner gaseous component. To investigate more deeply these regions, getting images at the first Astronomical Unit scale is necessary. This has become possible with near infrared instruments on the VLTI. We have developed a new imaging method adapted to young stellar objects where we process separately the stellar component from the rest of the image to reveal the environment by using the spectral differences between these two components. We present the result of this method on the first imaging survey of Herbig stars carried out by PIONIER on the VLTI.

Context. The disks around some Herbig Be stars have been observed to be more compact than the expected dust sublimation radius for such objects. Highly refractory dust grains and optically thick gas emission have been proposed as possible explanations for this phenomenon. Aims. Previously, the “undersized” Herbig Be star MWC 147 was observed with interferometry, and the results indicated a presence of a compact gaseous disk based on the measured wavelength-dependence of near-infrared or mid-infrared visibilities. Our aim is to search for direct evidence for the presence of hot gas inside of the expected dust sublimation radius of MWC 147. Methods. By combining VLTI/AMBER spectro-interferometry (R = 12 000) with CRIRES spectroscopy (R = 100 000) we can both spectrally and spatially resolve the Brγ line-emitting gas around MWC 147. Additionally, using CHARA/CLIMB enables us to achieve baseline lengths up to 330 m, offering ~2 times higher angular resolution (and a better position angle coverage) than has previously been achieved with interferometry for MWC 147. To model the continuum we fit our AMBER and CLIMB measurements with a geometric model of an inclined Gaussian distribution as well as a ring model. We fit our high-resolution spectra and spectro-interferometric data with a kinematic model of a disk in Keplerian rotation. Results. Our interferometric visibility modelling of MWC 147 indicates the presence of a compact continuum disk with a close to face-on orientation. We model the continuum with an inclined Gaussian and a ring with a radius of 0.60 mas (0.39 au), which is well within the expected dust sublimation radius of 1.52 au. We detect no significant change in the measured visibilities across the Brγ line, indicating that the line-emitting gas is located in the same region as the continuum-emitting disk. Using our differential phase data we construct photocentre displacement vectors across the Brγ line, revealing a velocity profile consistent with a rotating disk. We fit our AMBER spectro-interferometry data with a kinematic model of a disk in Keplerian rotation, where both the line-emitting and continuum-emitting components of the disk originate from the same compact region close to the central star. The presence of line-emitting gas in the same region as the K-band continuum supports the interpretation that the K-band continuum traces an optically thick gas disk. Conclusions. Our spatially and spectrally resolved observations of MWC 147 reveal that the K-band continuum and Brγ emission both originate from a similar region which is 3.9 times more compact than the expected dust sublimation radius for the star; Brγ is emitted from the accretion disk or disk wind region and exhibits a rotational velocity profile. We conclude that we detect the presence of a compact, gaseous accretion disk in Keplerian rotation around MWC 147.