High Angular Resolution Astronomy Research Articles

The Orbit and Dynamical Mass of Polaris: Observations with the CHARA Array

The 30 yr orbit of the Cepheid Polaris has been followed with observations by the Center for High Angular Resolution Astronomy (CHARA) Array from 2016 through 2021. An additional measurement has been made with speckle interferometry at the Apache Point Observatory. Detection of the companion is complicated by its comparative faintness—an extreme flux ratio. Angular diameter measurements appear to show some variation with pulsation phase. Astrometric positions of the companion were measured with a custom grid-based model-fitting procedure and confirmed with the CANDID software. These positions were combined with the extensive radial velocities (RVs) discussed by Torres to fit an orbit. Because of the imbalance of the sizes of the astrometry and RV data sets, several methods of weighting are discussed. The resulting mass of the Cepheid is 5.13 ± 0.28 M ⊙. Because of the comparatively large eccentricity of the orbit (0.63), the mass derived is sensitive to the value found for the eccentricity. The mass combined with the distance shows that the Cepheid is more luminous than predicted for this mass from evolutionary tracks. The identification of surface spots is discussed. This would give credence to the identification of a radial velocity variation with a period of approximately 120 days as a rotation period. Polaris has some unusual properties (rapid period change, a phase jump, variable amplitude, and unusual polarization). However, a pulsation scenario involving pulsation mode, orbital periastron passage, and low pulsation amplitude can explain these characteristics within the framework of pulsation seen in Cepheids.

An Angular Diameter Measurement of β UMa via Stellar Intensity Interferometry with the VERITAS Observatory

We use the Very Energetic Radiation Imaging telescope Array System (VERITAS) imaging air Cherenkov telescope array to obtain the first measured angular diameter of β UMa at visual wavelengths using stellar intensity interferometry (SII) and independently constrain the limb-darkened angular diameter. The age of the Ursa Major moving group has been assessed from the ages of its members, including nuclear member Merak (β UMa), an A1-type subgiant, by comparing effective temperature and luminosity constraints to model stellar evolution tracks. Previous interferometric limb-darkened angular-diameter measurements of β UMa in the near-infrared (Center for High Angular Resolution Astronomy (CHARA) Array, 1.149 ± 0.014 mas) and mid-infrared (Keck Nuller, 1.08 ± 0.07 mas), together with the measured parallax and bolometric flux, have constrained the effective temperature. This paper presents current VERITAS-SII observation and analysis procedures to derive squared visibilities from correlation functions. We fit the resulting squared visibilities to find a limb-darkened angular diameter of 1.07 ± 0.04 (stat) ± 0.05 (sys) mas, using synthetic visibilities from a stellar atmosphere model that provides a good match to the spectrum of β UMa in the optical wave band. The VERITAS-SII limb-darkened angular diameter yields an effective temperature of 9700 ± 200 ± 200 K, consistent with ultraviolet spectrophotometry, and an age of 390 ± 29 ± 32 Myr, using MESA Isochrones and Stellar Tracks. This age is consistent with 408 ± 6 Myr from the CHARA Array angular diameter.

The CHARA Array Interferometric Program on the Multiplicity of Classical Be Stars: New Detections and Orbits of Stripped Subdwarf Companions

Rapid rotation and nonradial pulsations enable Be stars to build decretion disks, where the characteristic line emission forms. A major but unconstrained fraction of Be stars owe their rapid rotation to mass and angular momentum transfer in a binary. The faint, stripped companions can be helium-burning subdwarf OB-type stars (sdOBs), white dwarfs (WDs), or neutron stars. We present optical/near-infrared Center for High Angular Resolution Astronomy (CHARA) interferometry of 37 Be stars selected for spectroscopic indications of low-mass companions. From multiepoch H- and/or K-band interferometry plus radial velocities and parallaxes collected elsewhere, we constructed 3D orbits and derived flux ratios and absolute dynamical masses of both components for six objects, quadrupling the number of anchor points for evolutionary models. In addition, a new wider companion was identified for the known Be + sdO binary 59 Cyg, while auxiliary Very Large Telescope Interferometer/GRAVITY spectrointerferometry confirmed circumstellar matter around the sdO companion to HR 2142. On the other hand, we failed to detect any companion to the six Be stars with γ Cas–like X-ray emission, with sdOB and main-sequence companions of the expected spectroscopic mass being ruled out for the X-ray-prototypical stars γ Cas and π Aqr, leaving elusive WDs as the most likely companions, as well as a likely explanation of the X-rays. No low-mass main-sequence close companions were identified for the other stars.

Orbits and dynamical masses for the active Hyades multiple system HD 284163

ABSTRACT We report near-infrared long-baseline interferometric observations of the Hyades multiple system HD 284163, made with the Center for High Angular Resolution Astronomy array, as well as almost 43 yr of high-resolution spectroscopic monitoring at the Center for Astrophysics. Both types of observations resolve the 2.39 d inner binary, and also an outer companion in a 43.1 yr orbit. Our observations, combined with others from the literature, allow us to solve for the 3D inner and outer orbits, which are found to be at nearly right angles to each other. We determine the dynamical masses of the three stars (good to better than 1.4 per cent for the inner pair), as well as the orbital parallax. The secondary component (0.5245 ± 0.0047 M⊙) is now the lowest mass star with a dynamical mass measurement in the cluster. A comparison of these measurements with current stellar evolution models for the age and metallicity of the Hyades shows good agreement. All three stars display significant levels of chromospheric activity, consistent with the classification of HD 284163 as an RS CVn object. We present evidence that a more distant fourth star is physically associated, making this a hierarchical quadruple system.

Integrated turbulence parameters' estimation from NAOMI adaptive optics telemetry data

Context.Monitoring turbulence parameters is crucial in high-angular resolution astronomy for various purposes, such as optimising adaptive optics systems or fringe trackers. The former systems are present at most modern observatories and will remain significant in the future. This makes them a valuable complementary tool for the estimation of turbulence parameters.Aims.The feasibility of estimating turbulence parameters from low-resolution sensors remains untested. We performed seeing estimates for both simulated and on-sky telemetry data sourced from the new adaptive optics module installed on the four Auxiliary Telescopes of the Very Large Telescope Interferometer.Methods.The seeing estimates were obtained from a modified and optimised algorithm that employs a chi-squared modal fitting approach to the theoretical von Kármán model variances. The algorithm was built to retrieve turbulence parameters while simultaneously estimating and accounting for the remaining and measurement error. A Monte Carlo method was proposed for the estimation of the statistical uncertainty of the algorithm.Results.The algorithm is shown to be able to achieve per-cent accuracy in the estimation of the seeing with a temporal horizon of 20 s on simulated data. A (0.76″ ± 1.2%|stat± 1.2%|sys) median seeing was estimated from on-sky data collected from 2018 to 2020. The spatial distribution of the Auxiliary Telescopes across the Paranal Observatory was found to not play a role in the value of the seeing.

Precise Age for the Binary Star System 12 Com in the Coma Berenices Cluster

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 data set and interferometric measurements from the Palomar Testbed Interferometer and the Center for High Angular Resolution Astronomy array, we measured masses M 1 = 2.64 ± 0.07M ⊙ and M 2 = 2.10 ± 0.03M ⊙. Interferometry also allows us to resolve the giant and measure its size as R 1 = 9.12 ± 0.12 ± 0.01R ⊙. With the measured masses and radii, we find an age of 533 ± 41 ± 42 Myr. For comparison, we measure the radius of 31 Com to be 8.36 ± 0.15 R ⊙. 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.

Advances in Optical/Infrared Interferometry

After decades of fast-paced technical advances, optical/infrared (O/IR) interferometry has seen a revolution in recent years: ▪ The GRAVITY instrument at the Very Large Telescope Interferometer (VLTI) with four 8-m telescopes reaches thousand-times-fainter objects than possible with earlier interferometers, and the Center for High Angular Resolution Astronomy array (CHARA) routinely offers up to 330-m baselines and aperture synthesis with six 1-m telescopes. ▪ The observed objects are fainter than 19 mag, the images have submilliarcsecond resolution, and the astrometry reaches microarcsecond precision. ▪ This led to breakthrough results on the Galactic Center, exoplanets, active galactic nuclei, young stellar objects, and stellar physics. Following a primer in interferometry, we summarize the advances that led to the performance boost of modern interferometers: ▪ Single-mode beam combiners now combine up to six telescopes, and image reconstruction software has advanced over earlier developments for radio interferometry. ▪ With a combination of large telescopes, adaptive optics (AO), fringe tracking, and especially dual-beam interferometry, GRAVITY has boosted the sensitivity by many orders of magnitude. Another order-of-magnitude improvement will come from laser guide star AO. In combination with large separation fringe tracking, O/IR interferometry will then provide complete sky coverage for observations in the Galactic plane and substantial coverage for extragalactic targets.

The Orbits and Dynamical Masses of the Castor System

Castor is a system of six stars in which the two brighter objects, Castor A and B, revolve around each other every ∼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 Center for High Angular Resolution Astronomy (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°) relative to the wide orbit, and that of B is inclined about 59° 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 the CHARA array, 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.

Active Optics—Advances of Cycloid-like Variable Curvature Mirrors for the VLTI Array

Elasticity theory and active optics led us to the discovery of three geometrical configurations of variable curvature mirrors (VCMs) that are either cycloid-like or tulip-like thickness distributions. Cycloid-like VCMs are generated by a uniform load—air pressure—applied over the mirror rear surface, and reacts without any bending moment along its circular contour. This particular VCM configuration is of practical interest because it smoothly generates accurate optical curvatures, varying from plane at rest to spherical curvatures up to f/2.9 over 16-mm aperture under 6.5-bar air pressure. Starting from the thin plate theory of elasticity and modeling with NASTRAN finite element analysis, one shows that 3-D optimizations—using a non-linear static flexural option—provide an accurate cycloid-like thickness distribution. VCM elasticity modeling in quenched stainless steel–chromium substrates allows the obtaining of diffraction-limited optical surfaces: Rayleigh’s criterion is achieved over a zoom range from flat to f/3.6 over 13-mm clear aperture up to 6-bar loading. These VCMs were originally developed and built at the Marseille Observatory in 1975 and implemented as a cat’s-eye mirror of IR Fourier-transform interferometers for laboratory recording of fast events in gas molecular spectroscopy. Later, for high-angular resolution astronomy with the ESO VLTI array—an interferometer made of 8 m Unit Telescopes (UTs) and 1.8 m Auxiliary Telescopes (ATs)—such VCMs were inevitable components to provide in a 3″ co-phased field-of-view since 1998. They were implemented (1) as cat’s eye mirrors of the height delay-lines beam recombination lab and (2) as ATs mirror-pair for output pupil conjugation of the movable x–y baseline. From the ESO-AMU approved convention of making 10 VCM spares up to 2024, the present modeling should provide a diffraction-limited extended field-of-view. It is pure coincidence that present results from modeling with an outer collarette are identical to results from analytic theory without collarette.

Long Term Evolution of Surface Features on the Red Supergiant AZ Cyg

Abstract We present H-band interferometric observations of the red supergiant (RSG) AZ Cyg that were made with the Michigan Infra-Red Combiner (MIRC) at the six-telescope Center for High Angular Resolution Astronomy (CHARA) Array. The observations span 5 yr (2011–2016), which offers 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 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 − 50 + 57 R ⊙ . 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. The 1D power spectra of these images suggest a characteristic size of 0.52–0.69 R ⋆ for the larger, long lived features. This is close to the values of 0.51–0.53 R ⋆ that are 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.

Ultrafast laser inscription of asymmetric integrated waveguide 3 dB couplers for astronomical K-band interferometry at the CHARA array

We present the fabrication and characterization of 3 dB asymmetric directional couplers for the astronomical K-band at wavelengths between 2.0 and 2.4 µm. The couplers were fabricated in commercial Infrasil silica glass using an ultrafast laser operating at 1030 nm. After optimizing the fabrication parameters, the insertion losses of straight single-mode waveguides were measured to be ∼ 1.2 ± 0.5 d B across the full K-band. We investigate the development of asymmetric 3 dB directional couplers by varying the coupler interaction lengths and by varying the width of one of the waveguide cores to detune the propagation constants of the coupled modes. In this manner, we demonstrate that ultrafast laser inscription is capable of fabricating asymmetric 3 dB directional couplers for future applications in K-band stellar interferometry. Finally, we demonstrate that our couplers exhibit an interferometric fringe contrast of > --> 90 % . This technology paves the path for the development of a two-telescope K-band integrated optic beam combiner for interferometry to replace the existing beam combiner (MONA) in Jouvence of the Fiber Linked Unit for Recombination (JouFLU) at the Center for High Angular Resolution Astronomy (CHARA) telescope array.

Compensation of differential dispersion: application to multiband stellar interferometry

ABSTRACT With the aim of pushing the limiting magnitude of interferometric instruments, the need for wide-band detection channels and for a coordinated operation of different instruments has considerably grown in the field of long-baseline interferometry. For this reason, the Center for High Angular Resolution Astronomy (CHARA), an array of six telescopes, requires a new configuration of longitudinal dispersion compensators to keep the fringe contrast above 95 per cent simultaneously in all spectral bands, while preserving the transmission above 85 per cent. In this paper, we propose a new method for defining the longitudinal dispersion compensators (LDC) suited for multiband observations. A literal approximation of the contrast loss resulting from the dispersion residues enables us to define a general criterion for fringe contrast maximization on several bands simultaneously. The optimization of this criterion leads to a simple solution with only two LDC stages per arm and existing differential delay lines, to the glass choice and a simple linear formula for thickness control of all these media. A refined criterion can also take into account glass transmission. After presenting this criterion, we give the optimal solution (medium, configuration) and its expected performance for the planned observing modes on CHARA.

Dynamical Surface Imaging of λ Andromedae

Abstract We present temperature maps of RS CVn star λ Andromedae, reconstructed from interferometric data acquired in 2010 and 2011 by the Michigan Infra-Red Combiner instrument at the Center for High Angular Resolution Astronomy Array. To constrain the stellar parameters required for this imaging task, we first modeled the star using our GPU-accelerated code SImulation and Modeling Tool for Optical Interferometry. The stellar surface was then imaged using our open source interferometric imaging code ROTational Image Reconstruction, in the process further refining the estimation of stellar parameters. We report that the measured angular diameter is 2.742 ± 0.010 mas with a limb-darkening coefficient of 0.231 ± 0.024. While our images are consistent with those of prior works, we provide updated physical parameters for λ Andromedae (R ⋆ = 7.78 ± 0.05 R ⊙, M ⋆ = 1.24 ± 0.72 M ⊙, and log L/L ⊙ = 1.46 ± 0.04).

Interferometric Imaging of λ Andromedae: Evidence of Starspots and Rotation

Abstract Presented are the first interferometric images of cool starspots on the chromospherically active giant λ Andromedae. Using the Michigan Infra-Red Combiner coupled to the Center for High Angular Resolution Astronomy Array, 26 interferometric observations were made between 2008 August 17 and 2011 September 24. The photometric time series acquired at Fairborn Observatory spanning 2008 September 20 to 2011 January 20 is also presented. The angular diameter and power-law limb-darkening coefficient of this star are 2.759 ± 0.050 mas and 0.229 ± 0.111, respectively. Starspot properties are obtained from both modeled and SQUEEZE reconstructed images. The images from 2010 through 2011 show anywhere from one to four starspots. The cadence in the data for the 2010 and 2011 data sets is sufficient to measure a stellar rotation period based on apparent starspot motion. This leads to estimates of the rotation period (P 2010 = 61 ± 4.0 days, P 2011 = 54.0 ± 2.4 days) that are consistent with the photometrically determined period of 54.8 days. In addition, the inclination and position angle of the rotation axis are computed for both the 2010 and 2011 data sets; values ( Ψ ¯ = 21.°5, i ¯ = 78.°0) for each are nearly identical between the two years.

A New Look into K-giants’ Chromospheres

Abstract We report the preliminary modeling of archival Center for High Angular Resolution Astronomy (CHARA)/Visible spEctroGraph and polArimeter interferometric data of a K-giant star using the PHOENIX atmosphere code. We find that our preparatory model that includes only the chromospheric contribution closely reproduces the observed infrared Ca ii triplet line profiles of a test star: the K-giant, β Cet. This preliminary work requires the additional modeling of the wind contribution to improve the agreement with observations. We plan to perform a systematic study of K-giants chromospheric emission with multi-wavelength and multi-technique observations and modeling. Our plans include extending the modeling work to include the underlying wind component for a larger set of stars. Stellar Parameters and Images with a Cophased Array, the second-generation instrument at CHARA, will be the ideal instrument to perform such observations and reveal the chromospheric activity of K-giants.

Observational aspect of black hole dark matter

Advances in high angular resolution astronomy make it conceivable that black hole dark matter could be detected via angular deviation effects. Assuming the dark matter in the galaxy is made of solar mass black holes, there is a nontrivial probability that a line-of-sight through the galaxy, leads to [Formula: see text]’s deviations, a value that has been discussed for various astronomical projects. In cosmology, the effects are magnified by an increased density at early times and an opening of angles due to redshift. We stress the interest of the resolution in observations on the CMB and emphasize the distinction between this and the power spectrum. If the dark matter is made of primordial black holes, present at the CMB, random deflections of the CMB photons lead to a limit on the angular resolution, approximately [Formula: see text] rad, with [Formula: see text] the mass of the black holes. Using the resolutions of [Formula: see text] demonstrated in observations of the “acoustic peaks” then implies the limit [Formula: see text]. While this large value seems uninteresting, improved resolutions would lead to significant exclusion limits or conceivably the discovery of primordial black holes.

Influence of the input-stage architecture on the in-laboratory test of a mid-infrared interferometer: application to the ALOHA up-conversion interferometer in the L band

ABSTRACT In the framework of the Astronomical Light Optical Hybrid Analysis (ALOHA) laboratory mid-infrared (MIR) up-conversion fibred interferometer in the L band, we report on the influence of the input-stage architecture. Using an amplitude division set-up in the visible or near-infrared is a straightforward choice in most cases. In the MIR context, the results are slightly different and we show that a wavefront division set-up is needed. These in-laboratory principle experiments allow us to measure a reliable 88 per cent instrumental contrast with high flux and to obtain fringes from faint sources at 3.5 μm with a spectral bandwith of 37 nm converted to 817 nm. An equivalent limiting L-band magnitude around 3.9, equivalent to 3.0 fW nm−1, could be demonstrated on 1 m class telescopes. This opens the possibility of planning future on-sky tests at the Center for High Angular Resolution Astronomy (CHARA) array and of predicting the performance attained.

Benchmarking Substellar Evolutionary Models Using New Age Estimates for HD 4747 B and HD 19467 B

Abstract Constraining substellar evolutionary models (SSEMs) is particularly difficult due to a degeneracy between the mass, age, and luminosity of a brown dwarf. In cases where a brown dwarf is found as a directly imaged companion to a star, as in HD 4747 and HD 19467, the mass, age, and luminosity of the brown dwarf are determined independently, making them ideal objects to use to benchmark SSEMs. Using the Center for High Angular Resolution Astronomy Array, we measured the angular diameters and calculated the radii of the host stars HD 4747 A and HD 19467 A. After fitting their parameters to the Dartmouth Stellar Evolution Database, MESA Isochrones and Stellar Tracks, and Yonsei-Yale isochronal models, we adopt age estimates of Gyr for HD 4747 A and Gyr for HD 19467 A. Assuming the brown dwarf companions HD 4747 B and HD 19467 B have the same ages as their host stars, we show that many of the SSEMs under-predict bolometric luminosities by ∼0.75 dex for HD 4747 B and ∼0.5 dex for HD 19467 B. The discrepancies in luminosity correspond to over-predictions of the masses by ∼12% for HD 4747 B and ∼30% for HD 19467 B. We also show that SSEMs that take into account the effect of clouds reduce the under-prediction of luminosity to ∼0.6 dex and the over-prediction of mass to ∼8% for HD 4747 B, an L/T transition object that is cool enough to begin forming clouds. One possible explanation for the remaining discrepancies is missing physics in the models, such as the inclusion of metallicity effects.

Angular Sizes, Radii, and Effective Temperatures of B-type Stars from Optical Interferometry with the CHARA Array

Abstract We present interferometric observations of 25 spectral type-B stars that were made with the Precision Astronomical Visible Observations and the CLassic Interferometry with Multiple Baselines beam combiners at the Center for High Angular Resolution Astronomy Array (CHARA). The observations provide the angular sizes of these stars with an average error of 6%. The stars range in size from 1.09 mas for β Tau down to 0.20 mas for 32 Ori. We collected ultraviolet to infrared spectrophotometry and derived temperatures, angular diameters, and reddening estimates that best fit the spectra, as well as solutions with the angular size fixed by the interferometric measurements. There is generally good agreement between the observed and spectral fit angular diameters, indicating that the fluxes predicted from model atmospheres are reliable. On the other hand, the temperatures derived from angular diameters and fluxes tend to be larger (by ≈4%) than those from published results based on analysis of the line spectrum. This discrepancy may in part be attributed to unexplored atmospheric parameters or the existence of unknown companions. The physical radii of the stars are calculated from the angular diameters and Gaia DR2 parallaxes, and the target stars are placed in the Hertzsprung–Russell diagram for comparison with evolutionary tracks.

Angular Sizes and Effective Temperatures of O-type Stars from Optical Interferometry with the CHARA Array

Abstract We present interferometric observations of six O-type stars that were made with the Precision Astronomical Visible Observations beam combiner at the Center for High Angular Resolution Astronomy (CHARA) Array. The observations include multiple brackets for three targets, λ Ori A, ζ Oph, and 10 Lac, but there are only preliminary, single observations of the other three stars, ξ Per, α Cam, and ζ Ori A. The stellar angular diameters range from 0.55 mas for ζ Ori A down to 0.11 mas for 10 Lac, the smallest star yet resolved with the CHARA Array. The rotational oblateness of the rapidly rotating star ζ Oph is directly measured for the first time. We assembled ultraviolet to infrared flux measurements for these stars, and then derived angular diameters and reddening estimates using model atmospheres and an effective temperature set by published results from analysis of the line spectrum. The model-based angular diameters are in good agreement with those observed. We also present estimates for the effective temperatures of these stars, derived by setting the interferometric angular size and fitting the spectrophotometry.