Search

Inner regions of protoplanetary disks host many complex physical processes such as star-disk interactions, magnetic fields, planet formation, and the migration of new planets. To directly study this region requires milli-arcsecond angular resolution, beyond the diffraction limit of the world's largest optical telescopes and even too small for the mm-wave interferometer ALMA. However, we can use infrared interferometers to image the inner astronomical unit. Here, we present new results from the CHARA and VLTI arrays for the young and luminous Herbig Be star HD 190073. We detect a sub-AU cavity surrounded by a ring-like structure that we interpret as the dust destruction front. We model the shape with 6 radial profiles, 3 symmetric and 3 asymmetric, and present a model-free image reconstruction. All the models are consistent with a near face-on disk with inclination $\lesssim 20^\circ$, and we measure an average ring radius of 1.4 $\pm 0.2$ mas (1.14 AU). Around $48\%$ of the total flux comes from the disk with ~$15\%$ of that emission appearing to emerge from inside the inner rim. The cause of emission is still unclear, perhaps due to different dust grain compositions or gas emission. The skewed models and the imaging point to an off-center star, possibly due to binarity. Our image shows a sub-AU structure, which seems to move between the two epochs inconsistently with Keplerian motion and we discuss possible explanations for this apparent change., Comment: 20 pages, 9 figures

We present the visual orbits of four spectroscopic binary stars, HD 61859, HD 89822, HD 109510, and HD 191692, using long baseline interferometry with the CHARA Array. We also obtained new radial velocities from echelle spectra using the APO 3.5 m, CTIO 1.5 m, and Fairborn Observatory 2.0 m telescopes. By combining the astrometric and spectroscopic observations, we solve for the full, three-dimensional orbits and determine the stellar masses to 1-12% uncertainty and distances to 0.4-6% uncertainty. We then estimate the effective temperature and radius of each component star through Doppler tomography and spectral energy distribution analyses. We found masses of 1.4-3.5 Msun, radii of 1.5-4.7 Rsun, and temperatures of 6400-10300K. We then compare the observed stellar parameters to the predictions of the stellar evolution models, but found that only one of our systems fits well with the evolutionary models., Accepted for publication in AJ. arXiv admin note: text overlap with arXiv:2005.00546

Because many classical Be stars may owe their nature to mass and angular-momentum transfer in a close binary, the present masses, temperatures, and radii of their components are of high interest for comparison to stellar evolution models. Kappa Dra is a 61.5-day single-lined binary with a B6 IIIe primary. With the CHARA Array instruments MIRC/MIRC-X and MYSTIC, we detected the secondary at (approximately photospheric) flux ratios of 1.49 +- 0.10% and 1.63 +- 0.09% in the H and K band, respectively. From a large and diverse optical spectroscopic database only the radial velocity curve of the Be star could be extracted. However, employing the parallaxes from Hipparcos and Gaia, which agree within their nominal 1-sigma errors, we could derive the total mass and found component masses of 3.65 +- 0.48 Msun and 0.426 +- 0.043 Msun for the Be star and the companion, respectively. Previous cross-correlation of the observed far-UV spectrum with sdO spectral model templates had not detected a companion belonging to the hot O-type subdwarf (sdO) population known from ~20 earlier-type Be stars. Guided by our full 3D orbital solution, we found a strong cross-correlation signal for a stripped subdwarf B-type companion (far-UV flux ratio of 2.3 +- 0.5%), enabling the first firm characterization of such a star, and making kappa Dra the first mid- to late-type Be star with a directly-observed subdwarf companion., Comment: Accepted to ApJ. Figure sets available on request

Context: Stellar evolution models are highly dependent on accurate mass estimates, especially for high-mass stars in the early stages of evolution. The most direct method for obtaining model-independent masses is derivation from the orbit of close binaries. Aims: To derive the first astrometric+RV orbit solution for the single-lined spectroscopic binary MWC 166 A, based on CHARA and VLTI near-infrared interferometry over multiple epochs and ~100 archival radial velocity measurements, and to derive fundamental stellar parameters from this orbit. We also sought to model circumstellar activity in the system from K-band spectral lines. Methods: We geometrically modelled the dust continuum to derive astrometry at 13 epochs and constrain individual stellar parameters. We used the continuum models as a base to examine differential phases, visibilities and closure phases over the Br-$\gamma$ and He-I emission lines. Results: Our orbit solution suggests a period of $367.7\pm0.1$ d, twice as long as found with previous RV orbit fits, subsequently constraining the component masses to $M_1=12.2\pm2.2 M_\odot$ and $M_2=4.9\pm0.5 M_\odot$. The line-emitting gas was found to be localised around the primary and is spatially resolved on scales of ~11 stellar radii, with the spatial displacement between the line wings consistent with a rotating disc. Conclusions: The large radius and stable orientation of the line emission are inconsistent with magnetospheric or boundary-layer accretion, but indicate an ionised inner gas disk around MWC 166 Aa. We observe line variability that could be explained either with generic line variability in a Herbig star disc or V/R variations in a decretion disc. We also constrained the age of the system to ~$(7\pm2)\times10^5$ yr, consistent with the system being comprised of a main-sequence primary and a secondary still contracting towards the main sequence., Comment: 24 pages, 19 figures, 7 tables, 1 appendix. Accepted in A&A

Castor is a system of six stars in which the two brighter objects, Castor A and B, revolve around each other every $\sim$450 yr and are both short-period spectroscopic binaries. They are attended by the more distant Castor C, which is also a binary. Here we report interferometric observations with the CHARA array that spatially resolve the companions in Castor A and B for the first time. We complement these observations with new radial velocity measurements of A and B spanning 30 yr, with the Hipparcos intermediate data, and with existing astrometric observations of the visual AB pair obtained over the past three centuries. We perform a joint orbital solution to solve simultaneously for the three-dimensional orbits of Castor A and B as well as the AB orbit. We find that they are far from being coplanar: the orbit of A is nearly at right angles (92 degrees) relative to the wide orbit, and that of B is inclined about 59 degrees compared to AB. We determine the dynamical masses of the four stars in Castor A and B to a precision better than 1%. We also determine the radii of the primary stars of both subsystems from their angular diameters measured with CHARA, and use them together with stellar evolution models to infer an age for the system of 290 Myr. The new knowledge of the orbits enables us to measure the slow motion of Castor C as well, which may assist future studies of the dynamical evolution of this remarkable sextuple system., Comment: 17 pages in emulateapj format, including figures and tables. Accepted for publication in The Astrophysical Journal

Time series of spectroscopic, speckle-interferometric, and optical long-baseline-interferometric observations confirm that $\nu$ Gem is a hierarchical triple system. It consists of an inner binary composed of two B-type stars and an outer classical Be star. Several photospheric spectral lines of the inner components were disentangled, revealing two stars with very different rotational broadening ($\sim$260 and $\sim$140 kms$^{-1}$, respectively), while the photospheric lines of the Be star remain undetected. From the combined spectroscopic and astrometric orbital solution it is not possible to unambiguously cross-identify the inner astrometric components with the spectroscopic components. In the preferred solution based on modeling of the disentangled line profiles, the inner binary is composed of two stars with nearly identical masses of 3.3 M$_\odot$ and the more rapidly rotating star is the fainter one. These two stars are in a marginally elliptical orbit ($e$ = 0.06) about each other with a period of 53.8 d. The third star also has a mass of 3.3 M$_\odot$ and follows a more eccentric ($e$ = 0.24) orbit with a period of 19.1 yr. The two orbits are co-directional and, at inclinations of 79$^{\circ}$ and 76$^{\circ}$ of the inner and the outer orbit, respectively, about coplanar. No astrometric or spectroscopic evidence could be found that the Be star itself is double. The system appears dynamically stable and not subject to eccentric Lidov-Kozai oscillations. After disentangling, the spectra of the components of the inner binary do not exhibit peculiarities that would be indicative of past interactions. Motivations for a wide range of follow-up studies are suggested., Comment: accepted to ApJ

We present updated orbital elements for the Wolf-Rayet (WR) binary WR\,140 (HD\,193793; WC7pd + O5.5fc). The new orbital elements were derived using previously published measurements along with {\color{black}160} new radial velocity measurements across the 2016 periastron passage of WR 140. Additionally, four new measurements of the orbital astrometry were collected with the CHARA Array. With these measurements, we derive stellar masses of $M_{\rm WR} = 10.31\pm0.45 M_\odot$ and $M_{\rm O} = 29.27\pm1.14 M_{\odot}$. We also include a discussion of the evolutionary history of this system from the Binary Population and Spectral Synthesis (BPASS) model grid to show that this WR star likely formed primarily through mass loss in the stellar winds, with only a moderate amount of mass lost or transferred through binary interactions., 10 pages, 5 figures

MIRC-X is a six telescope beam combiner at the CHARA array that works in J and H wavelength bands and provides an angular resolution equivalent to a $B$=331m diameter telescope. The legacy MIRC combiner has delivered outstanding results in the fields of stellar astrophysics and binaries. However, we required higher sensitivity to make ambitious scientific measurements of faint targets such as young stellar objects, binary systems with exoplanets, and active galactic nuclei. For that purpose, MIRC-X is built and is offered to the community since mid-2017. MIRC-X has demonstrated up to two magnitudes of improved faint magnitude sensitivity with the best-case H

We present a new polarimetric mode for the MIRC-X 6-telescope beam combiner at CHARA. Utilizing the extensive u − v coverage afforded by CHARA this mode will be able to resolve and constrain scattered light in environs at milliarcsecond separations of target stars, a largely unexplored parameter space to-date in astronomy. Notably, this upgrade will allow for the investigation of the scattering properties of the inner dust wall at the sublimation radius of Herbig Ae/Be star disks, dust shells surrounding evolved stars, and gas-rich disks around Be stars. Our design adds a series of rotating half-wave plates, achromatic across J- and H-bands, and a polarizing beamsplitter into the MIRC-X beam path. In this work, we also preview on-sky observations, discussing ongoing work calibrating instrumental polarization effects in the CHARA beam path as well as upgrades to the MIRC-X data reduction pipeline.

In this contribution we report on our work to increase the spectral range of the Michigan Infrared Combiner- eXeter (MIRC-X) instrument at the CHARA array to allow for dual H and J band interferometric observations. We comment on the key science drivers behind this project and the methods of characterisation and correction of instrumental birefringence and dispersion. In addition, we report on the first results from on-sky commissioning in November 2019.

Young stars are surrounded by a circumstellar disk of gas and dust, within which planet formation can occur. Gravitational forces in multiple star systems can disrupt the disk. Theoretical models predict that if the disk is misaligned with the orbital plane of the stars, the disk should warp and break into precessing rings, a phenomenon known as disk tearing. We present observations of the triple star system GWOrionis, finding evidence for disk tearing. Our images show an eccentric ring that is misaligned with the orbital planes and the outer disk. The ring casts shadows on a strongly warped intermediate region of the disk. If planets can form within the warped disk, disk tearing could provide a mechanism for forming wide-separation planets on oblique orbits., 63 pages, 4+13 Figures, 6 Tables, published at Science

Context. Asymptotic giant branch stars are cool luminous evolved stars that are well observable across the Galaxy and populating Gaia data. They have complex stellar surface dynamics Aims. On the AGB star CL Lac, it has been shown that the convection-related variability accounts for a substantial part of the Gaia DR2 parallax error. We observed this star with the MIRC-X beam combiner installed at the CHARA interferometer to detect the presence of stellar surface inhomogeneities. Methods. We performed the reconstruction of aperture synthesis images from the interferometric observations at different wavelengths. Then, we used 3D radiative hydrodynamics simulations of stellar convection with CO5BOLD and the post-processing radiative transfer code Optim3D to compute intensity maps in the spectral channels of MIRC-X observations. Then, we determined the stellar radius and compared the 3D synthetic maps to the reconstructed ones focusing on matching the intensity contrast, the morphology of stellar surface structures, and the photocentre position at two different spectral channels, 1.52 and 1.70 micron, simultaneously. Results. We measured the apparent diameter of CL Lac at two wavelengths and recovered the radius using a Gaia parallax. In addition to this, the reconstructed images are characterised by the presence of a brighter area that largely affects the position of the photocentre. The comparison with 3D simulation shows good agreement with the observations both in terms of contrast and surface structure morphology, meaning that our model is adequate for explaining the observed inhomogenities. Conclusions. This work confirms the presence of convection-related surface structures on an AGB star of Gaia DR2. Our result will help us to take a step forward in exploiting Gaia measurement uncertainties to extract the fundamental properties of AGB stars using appropriate RHD simulations., Comment: Accepted for publication on Astronomy and Astrophysics

Ground-based optical long-baseline interferometry has the power to measure the orbits of close binary systems at ~10 micro-arcsecond precision. This precision makes it possible to detect "wobbles" in the binary motion due to the gravitational pull from additional short period companions. We started the ARrangement for Micro-Arcsecond Differential Astrometry (ARMADA) survey with the MIRC-X instrument at the CHARA array for the purpose of detecting giant planets and stellar companions orbiting individual stars in binary systems. We describe our observations for the survey, and introduce the wavelength calibration scheme that delivers precision at the tens of micro-arcseconds level for, 25 pages, 15 figures, accepted for publication in AJ

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

Nanoparticles grown in reactive plasmas and nanodusty plasmas gain high interest from basic science and technology. One of the great challenges of nanodusty plasmas is the in-situ diagnostic of the nanoparticle size and refractive index. The analysis of scattered light by means of the Mie solution of the Maxwell equations was proposed and used as an in-situ size diagnostic during the past two decades. Today, imaging ellipsometry techniques and the investigation of dense, i. e. optically thick nanoparticle clouds demand for analysis methods to take multiple scattering into account. We present the first 3D Monte-Carlo polarized radiative transfer simulations of the scattered light in a dense nanodusty plasma. This technique extends the existing diagnostic methods for the in-situ analysis of the properties of nanoparticles to systems where multiple scattering can not be neglected., 5 pages, 5 figures

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.

Hot exozodiacal dust emission was detected in recent surveys around two dozen main sequence stars at distances of less than $1\,\text{au}$ using H and K band interferometry. Due to the high contrast as well as the small angular distance between the circumstellar dust and the star, direct observation of this dust component is challenging. An alternative way to explore the hot exozodiacal dust is provided by mid-infrared interferometry. We analyze the L, M and N band interferometric signature of this emission in order to find stronger constraints for the properties and the origin of the hot exozodiacal dust. Considering the parameters of nine debris disc systems derived previously, we model the discs in each of these bands. We find that the M band possesses the best conditions to detect hot dust emission, closely followed by L and N bands. The hot dust in three systems - HD 22484 (10 Tau), HD 102647 ($\beta$ Leo) and HD 177724 ($\zeta$ Aql) - shows a strong signal in the visibility functions which may even allow one to constrain the dust location. In particular, observations in the mid-infrared could help to determine whether the dust piles up at the sublimation radius or is located at radii up to $1\,\text{au}$. In addition, we explore observations of the hot exozodiacal dust with the upcoming mid-infrared interferometer MATISSE at the Very Large Telescope Interferometer (VLTI)., Comment: 6 pages, 4 figures, 3 tables, accepted for publication in MNRAS

MIRC-X is a new beam combination instrument at the CHARA array that enables 6-telescope interferometric imaging on object classes that until now have been out of reach for milliarcsecond-resolution imaging. As part of an instrumentation effort lead by the University of Exeter and University of Michigan, we equipped the MIRC instrument with an ultra-low read-noise detector system and extended the wavelength range to the J and H-band. The first phase of the MIRC-X commissioning was successfully completed in June 2017. In 2018 we will commission polarisation control to improve the visibility calibration and implement a 'cross-talk resiliant' mode that will minimise visibility cross-talk and enable exoplanet searches using precision closure phases. Here we outline our key science drivers and give an overview about our commissioning timeline. We comment on operational aspects, such as remote observing, and the prospects of co-phased parallel operations with the upcoming MYSTIC combiner., Comment: 13 pages, 10 figures, presented at SPIE Astronomical Telescopes + Instrumentation 2018, Austin, Texas, USA

Context. FU Orionis is the archetypal FUor star, a subclass of young stellar objects (YSOs) that undergo rapid brightening events, often gaining between four and six magnitudes on timescales of days. This brightening is often associated with a massive increase in accretion, which is one of the most ubiquitous processes in astrophysics for bodies ranging from planets and stars to super-massive black holes. We present multi-band interferometric observations of the FU Ori circumstellar environment, including the first J-band interferometric observations of a YSO. Aims. We investigate the morphology and temperature gradient of the innermost regions of the accretion disk around FU Orionis. We aim to characterise the heating mechanisms of the disk and comment on potential outburst-triggering processes. Methods. Recent upgrades to the MIRC-X instrument at the CHARA array have allowed for the first dual-band J and H observations of YSOs. Using baselines up to 331 m, we present high-angular-resolution data of a YSO covering the near-infrared bands J, H, and K. The unprecedented spectral range of the data allowed us to apply temperature gradient models to the innermost regions of FU Ori. Results. We spatially resolved the innermost astronomical unit of the disk and determine the exponent of the temperature gradient of the inner disk to T ∝ r−0.74 ± 0.02. This agrees with theoretical works that predict T ∝ r−0.75 for actively accreting, steady-state disks, which is a value only obtainable through viscous heating within the disk. We found a disk that extends down to the stellar surface at 0.015 ± 0.007 au, where the temperature is found to be 5800 ± 700 K. We found a disk inclined at 32 ± 4° with a minor-axis position angle of 34 ± 11°. Conclusions. We demonstrate that J-band interferometric observations of YSOs are feasible with the MIRC-X instrument at CHARA. The temperature gradient power-law derived for the inner disk is consistent with theoretical predictions for steady-state, optically thick, viciously heated accretion disks.

We started a survey with CHARA/MIRC-X and VLTI/GRAVITY to search for low mass companions orbiting individual components of intermediate mass binary systems. With the incredible precision of these instruments, we can detect astrometric "wobbles" from companions down to a few tens of micro-arcseconds. This allows us to detect any previously unseen triple systems in our list of binaries. We present the orbits of 12 companions around early F to B-type binaries, 9 of which are new detections and 3 of which are first astrometric detections of known RV companions. The masses of these newly detected components range from 0.45-1.3 solar masses. Our orbits constrain these systems to a high astrometric precision, with median residuals to the orbital fit of 20-50 micro-arcseconds in most cases. For 7 of these systems we include newly obtained radial velocity data, which help us to identify the system configuration and to solve for masses of individual components in some cases. Although additional RV measurements are needed to break degeneracy in the mutual inclination, we find that the majority of these inner triples are not well-aligned with the wide binary orbit. This hints that higher mass triples are more misaligned compared to solar and lower mass triples, though a thorough study of survey biases is needed. We show that the ARMADA survey is extremely successful at uncovering previously unseen companions in binaries. This method will be used in upcoming papers to constrain companion demographics in intermediate mass binary systems down to the planetary mass regime., Comment: 25 pages, 14 figures, accepted for publication in AJ