Search

To accurately characterize the planets a star may be hosting, stellar parameters must first be well-determined. $\tau$ Ceti is a nearby solar analog and often a target for exoplanet searches. Uncertainties in the observed rotational velocities have made constraining $\tau$ Ceti's inclination difficult. For planet candidates from radial velocity (RV) observations, this leads to substantial uncertainties in the planetary masses, as only the minimum mass ($m \sin i$) can be constrained with RV. In this paper, we used new long-baseline optical interferometric data from the CHARA Array with the MIRC-X beam combiner and extreme precision spectroscopic data from the Lowell Discovery Telescope with EXPRES to improve constraints on the stellar parameters of $\tau$ Ceti. Additional archival data were obtained from a Tennessee State University Automatic Photometric Telescope and the Mount Wilson Observatory HK project. These new and archival data sets led to improved stellar parameter determinations, including a limb-darkened angular diameter of $2.019 \pm 0.012$ mas and rotation period of $46 \pm 4$ days. By combining parameters from our data sets, we obtained an estimate for the stellar inclination of $7\pm7^\circ$. This nearly-pole-on orientation has implications for the previously-reported exoplanets. An analysis of the system dynamics suggests that the planetary architecture described by Feng et al. (2017) may not retain long-term stability for low orbital inclinations. Additionally, the inclination of $\tau$ Ceti reveals a misalignment between the inclinations of the stellar rotation axis and the previously-measured debris disk rotation axis ($i_\mathrm{disk} = 35 \pm 10^\circ$)., Comment: 14 pages, 3 figures, 4 tables, 1 appendix, accepted for publication to AJ

V838 Mon is a stellar merger remnant that erupted in 2002 in a luminous red novae event. Although it is well studied in the optical, near infrared and submillimeter regimes, its structure in the mid-infrared wavelengths remains elusive. We observed V838 Mon with the MATISSE (LMN bands) and GRAVITY (K band) instruments at the VLTI and also the MIRCX/MYSTIC (HK bands) instruments at the CHARA array. We geometrically modelled the squared visibilities and the closure phases in each of the bands to obtain constraints on physical parameters. Furthermore, we constructed high resolution images of V838 Mon in the HK bands, using the MIRA and SQUEEZE algorithms to study the immediate surroundings of the star. Lastly, we also modelled the spectral features seen in the K and M bands at various temperatures. The image reconstructions show a bipolar structure that surrounds the central star in the post merger remnant. In the K band, the super resolved images show an extended structure (uniform disk diameter $\sim 1.94$ mas) with a clumpy morphology that is aligned along a north-west position angle (PA) of $-40^\circ$. Whereas in the H band, the extended structure (uniform disk diameter $\sim 1.18$ mas) lies roughly along the same PA. However, the northern lobe is slightly misaligned with respect to the southern lobe, which results in the closure phase deviations. The VLTI and CHARA imaging results show that V838 Mon is surrounded by features that resemble jets that are intrinsically asymmetric. This is also confirmed by the closure phase modelling. Further observations with VLTI can help to determine whether this structure shows any variation over time, and also if such bipolar structures are commonly formed in other stellar merger remnants., Comment: 18 pages, 34 figures

SU Aurigae is a widely studied T Tauri star and here we present original state-of-the-art interferometric observations with better uv and baseline coverage than previous studies. We aim to investigate the characteristics of the circumstellar material around SU Aur, constrain the disk geometry, composition and inner dust rim structure. The MIRC-X instrument at CHARA is a 6 telescope optical beam combiner offering baselines up to 331 m. 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 constrains the physical parameters of the disk. Image reconstruction reveals a highly inclined disk with a slight asymmetry consistent with inclination effects obscuring the inner disk rim through absorption of incident star light on the near-side and thermal re-emission/scattering of the far-side. Geometric models find that the underlying brightness distribution is best modelled as a Gaussian with a FWHM of $1.53\pm0.01 \mathrm{mas}$ at an inclination of $56.9\pm0.4^\circ$ and minor axis position angle of $55.9\pm0.5^\circ$. Radiative transfer modelling shows a flared disk with an inner radius at 0.16 au which implies a grain size of $0.14 \mathrm{\mu m}$ assuming astronomical silicates and a scale height of 9.0 au at 100 au. In agreement with literature, only the dusty disk wind successfully accounts for the NIR excess by introducing dust above the mid-plane. Our results confirm and provide better constraints than previous inner disk studies of SU Aurigae. We confirm the presence of a dusty disk wind in the cicumstellar environment, the strength of which is enhanced by a late infall event which also causes very strong misalignments between the inner and outer disks., Comment: arXiv admin note: substantial text overlap with arXiv:2111.06205, arXiv:1905.11907

Some evolved binaries, namely post-asymptotic giant branch binaries, are surrounded by stable and massive circumbinary disks similar to protoplanetary disks found around young stars. Around 10% of these disks are transition disks: they have a large inner cavity in the dust. Previous interferometric measurements and modeling have ruled out the cavity being formed by dust sublimation and suggested that the cavity is due to a massive circumbinary planet that traps the dust in the disk and produces the observed depletion of refractory elements on the surface of the post-AGB star. In this study, we test alternative scenario in which the large cavity could be due to dynamical truncation from the inner binary. We performed near-infrared interferometric observations with the CHARA Array on the archetype of such a transition disk around a post-AGB binary: AC Her. We detect the companion at ten epochs over 4 years and determine the 3-dimensional orbit using these astrometric measurements in combination with the radial velocity time series. This is the first astrometric orbit constructed for a post-AGB binary system. We derive the best-fit orbit with a semi-major axis $2.01 \pm 0.01$ mas ($2.83\pm0.08$ au), inclination $(142.9 \pm 1.1)^\circ$ and longitude of the ascending node $(155.1 \pm 1.8)^\circ$. We find that the theoretical dynamical truncation and dust sublimation radius are at least $\sim3\times$ smaller than the observed inner disk radius ($\sim21.5$ mas or 30 au). This strengthens the hypothesis that the origin of such a cavity is due to the presence of a circumbinary planet., Comment: Accepted to be published in The Astrophysical Journal

We present measurements of the interferometrically-resolved binary star system 12 Com and the single giant star 31 Com in the cluster Coma Berenices. 12 Com is a double-lined spectroscopic binary system consisting of a G7 giant and an A3 dwarf at the cluster turnoff. Using an extensive radial velocity dataset and interferometric measurements from PTI and the CHARA array, we measured masses $M_1 =2.64 \pm 0.07 M_\odot$ and $M_2 =2.10 \pm 0.03 M_\odot$. Interferometry also allows us to resolve the giant, and measure its size as $R_1 = 9.12 \pm 0.12 \pm 0.01 R_\odot$. With the measured masses and radii, we find an age of $533 \pm 41 \pm 42$ Myr. For comparison, we measure the radius of 31 Com to be $8.36 \pm 0.15 R_\odot$. Based on the photometry and radius measurements, 12 Com A is likely the most evolved bright star in the cluster, large enough to be in the red giant phase, but too small to have core helium burning. Simultaneous knowledge of 12 Com A's mass and photometry puts strong constraints on convective core overshooting during the main sequence phase, which in turn reduces systematic uncertainties in the age. Increased precision in measuring this system also improves our knowledge of the progenitor of the cluster white dwarf WD1216+260.

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

The study of the multiplicity of massive stars gives hints on their formation processes and their evolutionary paths, which are still not fully understood. Large separation binaries (>50 milliseconds of arc, mas) can be probed by adaptive-optics-assisted direct imaging and sparse aperture masking, while close binaries can be resolved by photometry and spectroscopy. However, optical long baseline interferometry is mandatory to establish the multiplicity of Galactic massive stars at the separation gap between 1 and 50 mas. In this paper, we aim to demonstrate the capability of the new interferometric instrument MIRC-X, located at the CHARA Array, to study the multiplicity of O-type stars and therefore probe the full range of separation for more than 120 massive stars (H<7.5 mag). We initiated a pilot survey of bright O-type stars (H<6.5mag) observable with MIRC-X. We observed 29 O-type stars, including two systems in average atmospheric conditions around a magnitude of H=7.5 mag. We systematically reduced the obtained data with the public reduction pipeline of the instrument. We analyzed the reduced data using the dedicated python software CANDID to detect companions. Out of these 29 systems, we resolved 19 companions in 17 different systems with angular separations between ~0.5 and 50 mas. This results in a multiplicity fraction fm=17/29=0.59+/-0.09, and an average number of companions fc=19/29=0.66+/-0.13. Those results are in agreement with the results of the SMASH+ survey in the Southern Hemisphere. Thirteen of these companions have been resolved for the first time, including the companion responsible for the nonthermal emission in Cyg OB2-5 A and the confirmation of the candidate companion of HD 47129 suggested by SMASH+. A large survey on more than 120 northern O-type stars (H<7.5) is possible with MIRC-X and will be fruitful., Comment: 15 pages, 9 figures, 5 tables, accepted in A&A

We present preliminary results from our long-baseline interferometry (LBI) survey to constrain the multiplicity properties of intermediate-mass A-type stars within 80pc. Previous multiplicity studies of nearby stars exhibit orbital separation distributions well-fitted with a log-normal with peaks > 15au, increasing with primary mass. The A-star multiplicity survey of De Rosa et al. (2014), sensitive beyond 30au but incomplete below 100 au, found a log-normal peak around 390au. Radial velocity surveys of slowly-rotating, chemically peculiar Am stars identified a significant number of very close companions with periods $\leq$ 5 days, ~ 0.1au, a result similar to surveys of O- and B-type primaries. With the improved performance of LBI techniques, we can probe these close separations for normal A-type stars where other surveys are incomplete. Our initial sample consists of 27 A-type primaries with estimated masses between 1.44-2.49M$_{\odot}$ and ages 10-790Myr, which we observed with the MIRC-X instrument at the CHARA Array. We use the open source software CANDID to detect five companions, three of which are new, and derive a companion frequency of 0.19$^{+0.11}_{-0.06}$ over mass ratios 0.25-1.0 and projected separations 0.288-5.481 au. We find a probability of 10$^{-6}$ that our results are consistent with extrapolations based on previous models of the A-star companion population, over mass ratios and separations sampled. Our results show the need to explore these very close separations to inform our understanding of stellar formation and evolution processes., Comment: 14 pages, 3 figures, Accepted to the Astrophysical Journal on Nov. 2, 2022

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

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., Comment: Accepted for publication in AJ. arXiv admin note: text overlap with arXiv:2005.00546

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

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

We present the complete sample of protoplanetary disks from the Gemini- Large Imaging with GPI Herbig/T-tauri Survey (Gemini-LIGHTS) which observed bright Herbig Ae/Be stars and T-Tauri stars in near-infrared polarized light to search for signatures of disk evolution and ongoing planet formation. The 44 targets were chosen based on their near- and mid-infrared colors, with roughly equal numbers of transitional, pre-transitional, and full disks. Our approach explicitly did not favor well-known, "famous" disks or those observed by ALMA, resulting in a less-biased sample suitable to probe the major stages of disk evolution during planet formation. Our optimized data reduction allowed polarized flux as low as 0.002% of the stellar light to be detected, and we report polarized scattered light around 80% of our targets. We detected point-like companions for 47% of the targets, including 3 brown dwarfs (2 confirmed, 1 new), and a new super-Jupiter mass candidate around V1295 Aql. We searched for correlations between the polarized flux and system parameters, finding a few clear trends: presence of a companion drastically reduces the polarized flux levels, far-IR excess correlates with polarized flux for non-binary systems, and systems hosting disks with ring structures have stellar masses $<$ 3 Msun. Our sample also included four hot, dusty "FS CMa" systems and we detected large-scale ($>100$ au) scattered light around each, signs of extreme youth for these enigmatic systems. Science-ready images are publicly available through multiple distribution channels using a new FITS file standard jointly developed with members of the VLT/SPHERE team., Comment: 51 pages, 31 figures, 7 tables, accepted to AJ

We observe the brightest member of the Praesepe cluster, Epsilon Cancri, to precisely measure the characteristics of the stars in this binary system, en route to a new measurement of the cluster's age. We present spectroscopic radial velocity measurements and interferometric observations of the sky-projected orbit to derive the masses, which we find to be M_1/M_sun = 2.420 +/- 0.008 and M_2/M_sun = 2.226 +/- 0.004. We place limits on the color-magnitude positions of the stars by using spectroscopic and interferometric luminosity ratios while trying to reproduce the spectral energy distribution of Epsilon Cancri. We re-examine the cluster membership of stars at the bright end of the color-magnitude diagram using Gaia data and literature radial velocity information. The binary star data are consistent with an age of 637 +/- 19 Myr, as determined from MIST model isochrones. The masses and luminosities of the stars appear to select models with the most commonly used amount of convective core overshooting., Comment: 46 pages, 15 figures

We present multi-instrument observations of the disc around the Herbig~Ae star, HD~145718, employing geometric and Monte Carlo radiative transfer models to explore the disc orientation, the vertical and radial extent of the near infrared (NIR) scattering surface, and the properties of the dust in the disc surface and sublimation rim. The disc appears inclined at $67-71^{\circ}$, with position angle, PA\,$=-1.0-0.6^{\circ}$, consistent with previous estimates. The NIR scattering surface extends out to $\sim75\,$au and we infer an aspect ratio, $h_{\rm{scat}}(r)/r\sim0.24$ in $J$-band; $\sim0.22$ in $H$-band. Our GPI images and VLTI+CHARA NIR interferometry suggest that the disc surface layers are populated by grains $\gtrsim \lambda/2\pi$ in size, indicating these grains are aerodynamically supported against settling and/or the density of smaller grains is relatively low. We demonstrate that our geometric analysis provides a reasonable assessment of the height of the NIR scattering surface at the outer edge of the disc and, if the inclination can be independently constrained, has the potential to probe the flaring exponent of the scattering surface in similarly inclined ($i\gtrsim70^{\circ}$) discs. In re-evaluating HD~145718's stellar properties, we found that the object's dimming events - previously characterised as UX~Or and dipper variability - are consistent with dust occultation by grains larger, on average, than found in the ISM. This occulting dust likely originates close to the inferred dust sublimation radius at $0.17\,$au., Comment: 18 pages, 16 figures; Accepted for publication in MNRAS

The distortions of absorption line profiles caused by photospheric brightness variations on the surfaces of cool, main-sequence stars can mimic or overwhelm radial velocity (RV) shifts due to the presence of exoplanets. The latest generation of precision RV spectrographs aims to detect velocity amplitudes $\lesssim 10$ cm s$^{-1}$, but requires mitigation of stellar signals. Statistical techniques are being developed to differentiate between Keplerian and activity-related velocity perturbations. Two important challenges, however, are the interpretability of the stellar activity component as RV models become more sophisticated, and ensuring the lowest-amplitude Keplerian signatures are not inadvertently accounted for in flexible models of stellar activity. For the K2V exoplanet host $\epsilon$ Eridani, we separately use ground-based photometry to constrain Gaussian processes for modeling RVs and TESS photometry with a light-curve inversion algorithm to reconstruct the stellar surface. From the reconstructions of TESS photometry, we produce an activity model, which reduces the rms scatter in RVs obtained with EXPRES from 4.72 m s$^{-1}$ to 1.98 m s$^{-1}$. We present a pilot study using the CHARA Array and MIRC-X beam combiner to directly image the starspots seen in the TESS photometry. With the limited phase coverage, our spot detections are marginal with current data but a future dedicated observing campaign should allow for imaging, as well as the stellar inclination and orientation with respect to its debris disk to be definitely determined. This work shows that stellar surface maps obtained with high cadence, time-series photometric and interferometric data can provide the constraints needed to accurately reduce RV scatter., Comment: 24 pages, 12 figures, 5 tables, accepted for publication in AJ

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

Dust evolution in protoplanetary disks from small dust grains to pebbles is key to the planet formation process. The gas in protoplanetary disks should influence the vertical distribution of small dust grains ($\sim$1 $\mu m$) in the disk.Utilizing archival near-infrared polarized light and millimeter observations, we can measure the scale height and the flare parameter $\beta$ of the small dust grain scattering surface and $^{12}$CO gas emission surface for three protoplanetary disks IM Lup, HD 163296, and HD 97048 (CU Cha). For two systems, IM Lup and HD 163296, the $^{12}$CO gas and small dust grains at small radii from the star have similar heights but at larger radii ($>$100 au) the dust grain scattering surface height is lower than the $^{12}$CO gas emission surface height. In the case of HD 97048, the small dust grain scattering surface has similar heights to the $^{12}$CO gas emission surface at all radii. We ran a protoplanetary disk radiative transfer model of a generic protoplanetary disk with TORUS and showed that there is no difference between the observed scattering surface and $^{12}$CO emission surface. We also performed analytical modeling of the system and found that gas-to-dust ratios larger than 100 could explain the observed difference in IM Lup and HD 163296. This is the first direct comparison of observations of gas and small dust grain heights distribution in protoplanetary disks. Future observations of gas emission and near-infrared scattered light instruments are needed to look for similar trends in other protoplanetary disks., Comment: Accepted ApJ, 13 pages, 3 figures, 4 tables

The great conjunction of 21 December 2020 saw Jupiter and Saturn appear together in the sky, separated by just a tenth of a degree (equivalent to a distance five times smaller than the diameter of the full Moon). This provided a potential once-in-a-lifetime opportunity to view the solar system's two biggest planets - and up to five of their moons - through a telescope eyepiece at the same time. Moreover, this was the first such opportunity, ever; previous observable conjunctions at similarly close separations took place before the development of the telescope in the early 1600s. Our team of scientists from the University of Exeter's Astrophysics Group and Exeter Science Centre worked with local social enterprises to develop a series of promotional and concurrent events to tie in with our live telescope broadcast of Jupiter and Saturn, to celebrate this spectacular celestial event. We hoped not only to inform and educate the public about great conjunctions, and the solar system more generally, but also to bring some light relief in what had been a rather difficult year., Comment: Published in Astronomy & Geophysics of the RAS, April 2021 issue. More information about the project (including links to all digital media) is available at http://jupitersaturn2020.org/

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., Comment: 10 pages, 5 figures

We present the first visual orbit for the nitrogen-rich Wolf-Rayet binary, WR 133 (WN5o + O9I) based on observations made with the CHARA Array and the MIRC-X combiner. This orbit represents the first visual orbit for a WN star and only the third Wolf-Rayet star with a visual orbit. The orbit has a period of 112.8 d, a moderate eccentricity of 0.36, and a separation of $a$= 0.79 mas on the sky. We combine the visual orbit with an SB2 orbit and Gaia parallax to find that the derived masses of the component stars are $M_{\rm WR}$ = $9.3\pm1.6 M_\odot$ and $M_{\rm O}$ = $22.6\pm 3.2 M_\odot$, with the large errors owing to the nearly face-on geometry of the system combined with errors in the spectroscopic parameters. We also derive an orbital parallax that is identical to the {\it Gaia}-determined distance. We present a preliminary spectral analysis and atmosphere models of the component stars, and find the mass-loss rate in agreement with polarization variability and our orbit. However, the derived masses are low compared to the spectral types and spectral model. Given the close binary nature, we suspect that WR 133 should have formed through binary interactions, and represents an ideal target for testing evolutionary models given its membership in the cluster NGC 6871., Comment: Accepted to ApJ Letters

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 <= 8. Here we present a review of the instrument and present early science results, and highlight some of our ongoing science programs.

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., Comment: 9 Pages, 6 figures. Proceedings of SPIE conference

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 <0.2 arcsecond binaries. We test our instrument performance on a known triple system kappa Peg, and show that our survey is delivering a factor of 10 better precision than previous similar surveys. We present astrometric detections of tertiary components to two B-type binaries: a 30-day companion to alpha Del, and a 50-day companion to nu Gem. We also collected radial velocity data for alpha Del with the Tennessee State University Automated Spectroscopic Telescope at Fairborn Observatory. We are able to measure the orbits and masses of all three components in these systems. We find that the previously published RV orbit for the inner pair of nu Gem is not consistent with our visual orbit. The precision achieved for these orbits suggests that our ARMADA survey will be successful at discovering new compact triple systems to A/B-type binary systems, leading to better statistics of hierarchical system architectures and formation history., Comment: 25 pages, 15 figures, accepted for publication in AJ

Context. FU Orionis is the archetypal FUor star, a subclass of young stellar object (YSO) that undergo rapid brightening events, often gaining 4-6 magnitudes on timescales of days. This brightening is often associated with a massive increase in accretion; one of the most ubiquitous processes in astrophysics 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 inner-most 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 allowed 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 allows us to apply temperature gradient models to the innermost regions of FU Ori. Results. We spatially resolve the innermost astronomical unit of the disk and determine the exponent of the temperature gradient of the inner disk to $T=r^{-0.74\pm0.02}$. This agrees with theoretical work that predicts $T = r^{-0.75}$ for actively accreting, steady state disks, a value only obtainable through viscous heating within the disk. We find a disk which extends down to the stellar surface at $0.015\pm0.007$ au where the temperature is found to be $5800\pm700$ K indicating boundary layer accretion. We find a disk inclined at $32\pm4^\circ$ with a minor-axis position angle of $34\pm11^\circ$., Comment: 10 pages, 6 figures

Spectral energy distributions (SEDs) are useful primary and complementary tools in the analysis of observations of young stars. However, the process of collating, inspecting, and flux-converting archival photometry and spectroscopy to build spectral energy distributions for young stars can be time-consuming. Here, I present SEDBYS (Spectral Energy Distribution Builder for Young Stars), a python-based repository of command-line tools built to (i) query online photometric and spectroscopic catalogs and a distributed database of archival photometry, (ii) use a look-up table of zero points to flux-convert the acquired data, (iii) enable visual inspection of the SED and (iv) handle book-keeping to collate references in bibTeX format. The code is distributed via git and is equipped with additional tools to enable users to add existing or forthcoming catalogs to the list of sources queried, ensuring the longevity of SEDBYS as a tool for the star formation community., Comment: Accepted for publication in Software X (Received 29 October 2020, Accepted 11 March 2021, Available online 3 April 2021.)

MIRC-X (Michigan InfraRed Combiner-eXeter) is a new highly-sensitive six-telescope interferometric imager installed at the CHARA Array that provides an angular resolution equivalent of up to a 330 m diameter baseline telescope in J and H band wavelengths ($\tfrac{\lambda}{2B}\sim0.6$ milli-arcseconds). We upgraded the original MIRC (Michigan InfraRed Combiner) instrument to improve sensitivity and wavelength coverage in two phases. First, a revolutionary sub-electron noise and fast-frame rate C-RED ONE camera based on a SAPHIRA detector was installed. Second, a new-generation beam combiner was designed and commissioned to (i) maximize sensitivity, (ii) extend the wavelength coverage to J-band, and (iii) enable polarization observations. A low-latency and fast-frame rate control software enables high-efficiency observations and fringe tracking for the forthcoming instruments at CHARA Array. Since mid-2017, MIRC-X has been offered to the community and has demonstrated best-case H-band sensitivity down to 8.2 correlated magnitude. MIRC-X uses single-mode fibers to coherently combine light of six telescopes simultaneously with an image-plane combination scheme and delivers a visibility precision better than 1%, and closure phase precision better than $1^\circ$. MIRC-X aims at (i) imaging protoplanetary disks, (ii) detecting exoplanets with precise astrometry, and (iii) imaging stellar surfaces and star-spots at an unprecedented angular resolution in the near-infrared. In this paper, we present the instrument design, installation, operation, and on-sky results, and demonstrate the imaging and astrometric capability of MIRC-X on the binary system $\iota$ Peg. The purpose of this paper is to provide a solid reference for studies based on MIRC-X data and to inspire future instruments in optical interferometry., Comment: 31 pages, 29 figures, accepted for publication in The Astronomical Journal

A crucial diagnostic that can tell us about processes involved in the formation and dynamical evolutionof planetary systems is the angle between the rotation axis of a star and a planet's orbital angular momentum vector ("spin-orbit" alignment or "obliquity"). Here we present the first spin-orbit alignment measurement for a wide-separation exoplanetary system, namely on the directly-imaged planet $\beta$ Pictoris b. We use VLTI/GRAVITY spectro-interferometry with an astrometric accuracy of 1 $\mu$as (microarcsecond) in the Br$\gamma$ photospheric absorption line to measure the photocenter displacement associated with the stellar rotation. Taking inclination constraints from astroseismology into account, we constrain the 3-dimensional orientation of the stellar spin axis and find that $\beta$ Pic b orbits its host star on a prograde orbit. The angular momentum vectors of the stellar photosphere, the planet, and the outer debris disk are well-aligned with mutual inclinations $<3\pm5^{\circ}$, which indicates that $\beta$ Pic b formed in a system without significant primordial misalignments. Our results demonstrate the potential of infrared interferometry to measure the spin-orbit alignment for wide-separation planetary systems, probing a highly complementary regime to the parameter space accessible with the Rossiter-McLaughlin effect. If the low obliquity is confirmed by measurements on a larger sample of wide-separation planets, it would lend support to theories that explain the obliquity in Hot Jupiter systems with dynamical scattering and the Kozai-Lidov mechanism., Comment: 6 pages, 4 figures, accepted for ApJL

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

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., Comment: 63 pages, 4+13 Figures, 6 Tables, published at Science

Debris disks are extrasolar analogs to our own Kuiper Belt and they are detected around at least 17% of nearby Sun-like stars. The morphology and dynamics of a disk encode information about its history, as well as that of any exoplanets within the system. We used ALMA to obtain 1.3 mm observations of the debris disk around the nearby F5V star HD 170773. We image the face-on ring and determine its fundamental parameters by forward-modeling the interferometric visibilities through a Markov Chain Monte Carlo approach. Using a symmetric Gaussian surface density profile, we find a 71 $\pm$ 4 au wide belt with a radius of 193$^{+2}_{-3}$ au, a relatively large radius compared to most other millimeter-resolved belts around late A / early F type stars. This makes HD 170773 part of a group of four disks around A and F stars with radii larger than expected from the recently reported planetesimal belt radius - stellar luminosity relation. Two of these systems are known to host directly imaged giant planets, which may point to a connection between large belts and the presence of long-period giant planets. We also set upper limits on the presence of CO and CN gas in the system, which imply that the exocomets that constitute this belt have CO and HCN ice mass fractions of <77% and <3%, respectively, consistent with Solar System comets and other exocometary belts., Comment: 14 pages, 6 figures, accepted to ApJ

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 3 x higher resolution and better baseline position angle coverage over previous observations. We investigate the characteristics of circumstellar material around SU Aur, constrain the disk geometry, composition and inner dust rim structure. The CHARA array offers opportunities for long baseline observations, with baselines up to 331 m. Using the CLIMB 3-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 geometric modeling. Additionally, the fitting of radiative transfer models constrains 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 dusty disk winds, where magnetic field lines drive dust above the disk plane close to the sublimation zone. 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/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 degrees and minor axis position angle 60.8+/-1.2 degrees. Radiative transfer modelling shows a flared disk with an inner radius at 0.18 au which implies a grain size of 0.4 um and a scale height of 15.0 au 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., Comment: Accepted for publication in Astronomy \& Astrophysics

The disks around some Herbig Be stars have been observed to be more compact than the expected dust sublimation radius for such objects, with highly refractory dust grains and optically thick gas emission having been proposed as possible explanations for this phenomenon. Our aim is to search for direct evidence for the presence of hot gas inside of the expected dust sublimation radius of MWC147. By combining VLTI/AMBER spectro-interferometry (R=12,000) with CRIRES spectroscopy (R=100,000) we can both spectrally and spatially resolve the Br-gamma line-emitting gas around MWC147. Our interferometric visibility modelling of MWC147 indicates the presence of a compact continuum disk with a close to face-on orientation. We model the continuum with an inclined Gaussian, as well as a ring with a radius of 0.60mas (0.39au) which is well within the expected dust sublimation radius of 1.52au. We detect no significant change in the measured visibilities across the Br-gamma line, indicating that the line-emitting gas is located in the same region as the continuum-emitting disk. We fit our AMBER spectro-interferometry data with a kinematic model of a disk in Keplerian rotation, with both the line-emitting and continuum-emitting components of the disk originating 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. Our spatially and spectrally resolved observations of MWC147 reveal that the K-band continuum and Br-gamma emission both originate from a similar compact region, with Br-gamma emitted from the accretion disk or disk wind region and exhibiting a rotational velocity profile. We conclude that we detect the presence of a compact, gaseous accretion disk in Keplerian rotation around MWC147.

The discovery of close in, giant planets (hot Jupiters) with orbital angular momentum vectors misaligned with respect to the rotation axis of their host stars presents problems for planet formation theories in which planets form in discs with angular momentum vectors aligned with that of the star. Violent, high eccentricity migration mechanisms purported to elevate planetary orbits above the natal disc plane predict populations of proto-hot Jupiters which have not been observed with Kepler. Alternative theories invoking primordial star-disc misalignments have recently received more attention. Here, the relative alignment between stars and their protoplanetary discs is assessed for the first time for a sample of 20 pre-main-sequence stars. Recently published rotation periods derived from high quality, long duration, high cadence K2 light curves for members of the $\rho$ Ophiuchus and Upper Scorpius star forming regions are matched with high angular resolution observations of spatially resolved discs and projected rotational velocities to determine stellar rotation axis inclination angles which are then compared to the disc inclinations. Ten of the fifteen systems for which the stellar inclination could be estimated are consistent with star-disc alignment while five systems indicate potential misalignments between the star and its disc. The potential for chance misalignment of aligned systems due to projection effects and characteristic measurement uncertainties is also investigated. While the observed frequency of apparent star-disc misalignments could be reproduced by a simulated test population in which 100% of systems are truly aligned, the distribution of the scale of inferred misalignment angles could not., Comment: Accepted for publication in MNRAS

The physical processes occurring within the inner few astronomical units of proto-planetary disks surrounding Herbig Ae stars are crucial to setting the environment in which the outer planet-forming disk evolves and put critical constraints on the processes of accretion and planet migration. We present the most complete published sample of high angular resolution H- and K-band observations of the stars HD 163296 and HD 190073, including 30 previously unpublished nights of observations of the former and 45 nights of the latter with the CHARA long-baseline interferometer, in addition to archival VLTI data. We confirm previous observations suggesting significant near-infrared emission originates within the putative dust evaporation front of HD 163296 and show this is the case for HD 190073 as well. The H- and K-band sizes are the same within $(3 \pm 3)\%$ for HD 163296 and within $(6 \pm 10)\%$ for HD 190073. The radial surface brightness profiles for both disks are remarkably Gaussian-like with little or no sign of the sharp edge expected for a dust evaporation front. Coupled with spectral energy distribution analysis, our direct measurements of the stellar flux component at H and K bands suggest that HD 190073 is much younger (<400 kyr) and more massive (~5.6 M$_\odot$) than previously thought, mainly as a consequence of the new Gaia distance (891 pc)., Comment: 19 pages, 6 figures