Interferometric Observations Research Articles

Optimised use of interferometry, spectroscopy, and stellar atmosphere models for determining the fundamental parameters of stars

Thanks to recent progress in the field of optical interferometry, instrument sensitivities have now reached the level achieved in the domain of new space missions dedicated to exoplanet and stellar studies. Combining interferometry with other observational approaches enables the determination of stellar parameters and helps improve our understanding of stellar physics. In this paper, we aim to demonstrate a new way of using stellar atmosphere models for a joint interpretation of spectroscopic and interferometric observations. Starting from a discrete grid of one-dimensional (1D) stellar atmosphere models, we developed a training algorithm, based on an artificial neural network, capable of estimating the spectrum and intensity profile of a star over a range of wavelengths and viewing angles. A minimisation algorithm based on the trained function allowed for the simultaneous fitting of the observational spectrum and interferometric complex visibilities. As a result, coherent and precise stellar parameters can be extracted. We show the ability of the trained function to match the modelled intensity profiles of stars in the effective temperature range of $4500$ K to $7000$ K and surface gravity range of $3$ to $5$ dex, with a relative precision to the model that is better than $0.05<!PCT!>$. Using simulated interferometric data and actual spectroscopic measurements, we demonstrated the performance of our algorithm on a sample of five benchmark stars. Using this method, we achieved an accuracy within $0.5<!PCT!>$ for the angular diameter, radius, and surface gravity, and within $20$ K for the effective temperature. This paper demonstrates a new method of using interferometric data combined with spectroscopic observations. This approach offers an improved determination of the radius, effective temperature, and surface gravity of stars.

Hidden Companions to Intermediate-mass Stars. Upgraded Multiplicity 2 → 3. XXIII. Discovery of a 0.94M ⊙, 9.2 au Companion to the Ap Star Gamma Piscis Austrini A* *Based on observations collected at the European Southern Observatory, Chile, Program ID 113.26GM.001.

Gamma Piscis Austrini is a nearby ρ = 4.″1 ↔ 284 au visual binary consisting of an A0Vp primary (A) and an F5V secondary (B). Here we report that the primary is actually a closer binary based on a VLTI/GRAVITY interferometric observation. The newly discovered companion Ab has a K band flux ratio 4.8% ± 0.2% and a projected separation ρ = 143.50 mas ↔ 9.2 au relative to Aa. We estimate isochrone masses M Aa = 2.65M ⊙, M Ab = 0.94M ⊙ and M B = 1.20M ⊙. γ PsA will become a triple white dwarf system within a Hubble time, with a low enough hierarchy for possibly interesting dynamical effects with a timescale P KL ∼ 200 kyr.

Probing the physics of star formation (ProPStar)

Context. Turbulence is a key component of molecular cloud structure. It is usually described by a cascade of energy down to the dissipation scale. The power spectrum for subsonic incompressible turbulence is ∝k−5/3, while for supersonic turbulence it is ∝k−2. Aims. We determine the power spectrum in an actively star-forming molecular cloud, from parsec scales down to the expected magnetohydrodynamic (MHD) wave cutoff (dissipation scale). Methods. We analyzed observations of the nearby NGC 1333 star-forming region in three different tracers to cover the different scales from ∼10 pc down to 20 mpc. The largest scales are covered with the low-density gas tracer 13CO (1–0) obtained with a single dish, the intermediate scales are covered with single-dish observations of the C18O (3–2) line, while the smallest scales are covered in H13CO+ (1–0) and HNC (1–0) with a combination of NOEMA interferometer and IRAM 30m single-dish observations. The complementarity of these observations enables us to generate a combined power spectrum covering more than two orders of magnitude in spatial scale. Results. We derive the power spectrum in an active star-forming region spanning more than 2 decades of spatial scales. The power spectrum of the intensity maps shows a single power-law behavior, with an exponent of 2.9 ± 0.1 and no evidence of dissipation. Moreover, there is evidence that the power spectrum of the ions to have more power at smaller scales than the neutrals, which is opposite to the theoretical expectations. Conclusions. We show new possibilities for studying the dissipation of energy at small scales in star-forming regions provided by interferometric observations.

Impact of the Epoch of Reionization sources on the 21-cm bispectrum

The morphology of the 21-cm signal emitted by the neutral hydrogen present in the intergalactic medium (IGM) during the Epoch of Reionization (EoR) depends both on the properties of the sources of ionizing radiation and on the underlying physical processes within the IGM. Variation in the morphology of the IGM 21-cm signal due to the different sources of the EoR is expected to have a significant impact on the 21-cm bispectrum, which is one of the crucial observable statistics that can evaluate the non-Gaussianity present in the signal and which can be estimated from radio interferometric observations of the EoR. Here we present the 21-cm bispectrum for different reionization scenarios assuming different simulated models for the sources of reionization. We also demonstrate how well the 21-cm bispectrum can distinguish between different IGM 21-cm signal morphologies, arising due to the differences in the reionization scenarios, which will help us shed light on the nature of the sources of ionizing photons. Our estimated large-scale bispectrum for all unique k-triangle shapes shows a significant difference in the magnitude and sign across different reionization scenarios. Additionally, our focused analysis of bispectrum for a few specific k-triangle shapes (e.g. squeezed-limit, linear, and shapes in the vicinity of the squeezed-limit) shows that the large scale 21-cm bispectrum can distinguish between reionization scenarios that show inside-out, outside-in and a combination of inside-out and outside-in morphologies. These results highlight the potential of using the 21-cm bispectrum for constraining different reionization scenarios.

CHARA Near-infrared Imaging of the Yellow Hypergiant Star ρ Cassiopeiae: Convection Cells and Circumstellar Envelope

Massive evolved stars such as red supergiants and hypergiants are potential progenitors of Type II supernovae, and they are known for ejecting substantial amounts of matter, up to half their initial mass, during their final evolutionary phases. The rate and mechanism of this mass loss play a crucial role in determining their ultimate fate and the likelihood of their progression to supernovae. However, the exact mechanisms driving this mass ejection have long been a subject of research. Recent observations, such as the Great Dimming of Betelgeuse, have suggested that the activity of large convective cells, combined with pulsation, could be a plausible explanation for such mass-loss events. In this context, we conducted interferometric observations of the famous yellow hypergiant, ρ Cassiopeiae using the CHARA Array in H- and K-band wavelengths. ρ Cas is well known for its recurrent eruptions, characterized by periods of visual dimming (∼1.5–2 mag) followed by recovery. From our observations, we derived the diameter of the limb-darkened disk and found that this star has a radius of 1.04 ± 0.01 mas, or 564–700 R ⊙. We performed image reconstructions with three different image reconstruction software packages, and they unveiled the presence of giant hot and cold spots on the stellar surface. We interpret these prominent hot spots as giant convection cells, suggesting a possible connection to mass ejections from the star’s envelope. Furthermore, we detected spectral CO emission lines in the K band (λ = 2.31–2.38 μm), and the image reconstructions in these spectral lines revealed an extended circumstellar envelope with a radius of 1.45 ± 0.10 mas.

LP349-25AB: Application of Magneto-convection Model to an M Dwarf Binary Where Radio and Astrometric Data have Yielded Masses with 0.7% Precision

Reliable (±(1–2)%) estimates of stellar masses have typically been obtained from eclipsing binaries. In such cases, derived values of radii and T eff are in many cases discrepant relative to standard models, especially in active stars. These discrepancies led to the development of magneto-convection theory. Recently, an independent technique has emerged for determining precise masses of radio-emitting binary components: interferometric observations at multiple points around the orbit. S. Curiel et al. have used this technique to determine the masses of both components of an M8+M9 binary (LP349-25AB) with remarkably high precision: ±0.7%. However, applying standard models suggests that the components are not co-eval, differing in age by ∼30%. But LP349-25 is such a close binary (with a separation of 1.2 au) that non-coevality seems unrealistic. Here, we report on magneto-convective modeling which resolves the apparent non-coevality while simultaneously explaining why the primary component is variable, but the secondary component is not.

Kinematics of the 2023 Kahramanmaraş Earthquake Doublet: Biased Near-Fault Data and Shallow Slip Deficit

Abstract Accurate estimate of the shallow slip deficit (SSD) for large strike-slip events is highly dependent on near-fault data. Previous studies have estimated the SSD of the 2023 Kahramanmaraş earthquake doublet, which may vary from a few percent to about 50%. Whether this reduced shallow slip is real or artificial is crucial for understanding the seismic hazards during and following the earthquake doublet. In this study, we inverted for the kinematic slip of this earthquake doublet with refined near-fault Interferometric Synthetic Aperture Radar observations and compared the results with the source model without further data processing. The model that excludes nonphysical data produced only ∼6% and ∼22% SSD for the M7.8 and 7.6 events, respectively, compared to ∼44% and ∼53% SSD of the model using original data. The increased data gap generally leads to overestimated SSD, but when the data coverage is almost complete, our result shows that the SSD is very sensitive to biased near-fault data which may induce significant artificial SSD. Our results suggest that overestimated SSD may be accompanied by slip migration from shallow toward deeper depths, to meet the total moment or energy release constrained by surface displacement data. The M7.8 event increased static Coulomb stress at where the second M7.6 event nucleated and thus may have triggered the following event.

McFine: python-based Monte-Carlo multi-component hyperfine structure fitting

Abstract Modelling complex line emission in the interstellar medium (ISM) is a degenerate, high-dimensional problem. Here, we present McFine, a tool for automated multi-component fitting of emission lines with complex hyperfine structure, in a fully automated way. We use Markov chain Monte Carlo (MCMC) to efficiently explore the complex parameter space, allowing for characterising model denegeracies. This tool allows for both local thermodynamic equilibrium (LTE) and radiative-transfer (RT) models. McFine can fit individual spectra and data cubes, and for cubes encourage spatial coherence between neighbouring pixels. It is also built to fit the minimum number of distinct components, to avoid overfitting. We have carried out tests on synthetic spectra, where in around 90 per cent of cases it fits the correct number of components, otherwise slightly fewer components. Typically, Tex is overestimated and τ underestimated, but accurate within the estimated uncertainties. The velocity and line widths are recovered with extremely high accuracy, however. We verify McFine by applying to a large Atacama Large Millimeter/submillimeter Array (ALMA) N2H+ mosaic of an high-mass star forming region, G316.75-00.00. We find a similar quality of fit to our synthetic tests, aside from in the active regions forming O-stars, where the assumptions of Gaussian line profiles or LTE may break down. To show the general applicability of this code, we fit CO(J = 2-1) observations of NGC 3627, a nearby star-forming galaxy, again obtaining excellent fit quality. McFine provides a fully automated way to analyse rich datasets from interferometric observations, is open source, and pip-installable.

XLSSC 122 caught in the act of growing up

Context. How protoclusters evolved from sparse galaxy overdensities to mature galaxy clusters is still not well understood. In this context, detecting and characterizing the hot intracluster medium (ICM) at high redshifts (z ∼ 2) is key to understanding how the continuous accretion from the filamentary large-scale structure and the mergers along it impact the first phases of cluster formation. Aims. We study the dynamical state and morphology of the z = 1.98 galaxy cluster XLSSC 122 with high-resolution observations (≈5″) of the ICM through the Sunyaev–Zeldovich (SZ) effect. XLSSC 122 is the highest redshift optically confirmed galaxy cluster found in an unbiased, widefield survey. Methods. Via Bayesian forward modeling, we mapped the ICM on scales from the virial radius down to the core of the cluster. To constrain such a broad range of spatial scales, we employed a new technique that jointly forward-models parametric descriptions of the pressure distribution to interferometric ACA and ALMA observations and multiband imaging data from ACT. Results. We detect the SZ effect with 11σ significance in the ALMA+ACA observations and find a flattened inner pressure profile that is consistent with a noncool core classification with a significance of ≥3σ. In contrast to the previous works, we find better agreement between the SZ effect signal and the X-ray emission as well as the cluster member distribution. Further, XLSSC 122 exhibits an excess of SZ flux in the south of the cluster where no X-ray emission is detected. By reconstructing the interferometric observations and modeling in the uv-plane, we obtain a tentative detection of an infalling group or filamentary-like structure in the southeast that is believed to boost and heat up the ICM while the density of the gas is still low. In addition to characterizing the dynamical state of the cluster, we provide an improved SZ mass estimate M500,c = 1.66−0.20+0.23 × 1014 M⊙. Conclusions. Altogether, the observations indicate that we see XLSSC 122 in a dynamic phase of cluster formation while a large reservoir of gas is already thermalized.

Examining the Nature of the Starless Dark Matter Halo Candidate Cloud-9 with Very Large Array Observations

Observations with the Five-hundred-meter Aperture Spherical Telescope recently detected H i 21 cm emission near M94, revealing an intriguing object, Cloud-9, without an optical counterpart. Subsequent analysis suggests that Cloud-9 is consistent with a gas-rich (M H I ≈ 106 M ⊙), starless, dark matter (DM) halo of mass M 200 ≈ 5 × 109 M ⊙. Using the Karl G. Jansky Very Large Array in D-array configuration, we present interferometric observations of Cloud-9, revealing it as a dynamically cold (W 50 ≈ 12 km s−1), nonrotating, and spatially asymmetric system, exhibiting gas compression on one side and a tail-like structure toward the other—features likely originating from ram pressure. Our observations suggest Cloud-9 is consistent with a starless ΛCDM DM halo if the gas is largely isothermal. If interpreted as a faint dwarf, Cloud-9 is similar to Leo T, a nearby gas-rich galaxy that would fall below current optical detection limits at Cloud-9's distance (d ≈ 5 Mpc). Further observations with the Hubble Space Telescope reaching magnitudes m g ≈ 30 would help identify such a galaxy or dramatically lower the current limits on its stellar mass (M gal ≲ 105 M ⊙). Cloud-9 thus stands as the firmest starless DM halo candidate to date or the faintest galaxy known at its distance.

Deep extragalactic H i survey of the COSMOS field with FAST

Abstract We present a deep H i survey at L-band conducted with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) over the COSMOS field. This survey is strategically designed to overlap with the MIGHTEE COSMOS field, aiming to combine the sensitivity of the FAST and high-resolution of the MeerKAT. We observed the field with FAST for 11 hours covering ∼2 square degrees, and reduced the raw data to H i spectral cubes over the frequency range 1310-1420 MHz. The FAST-H i data reach a median 3σ column density of NHI ∼ 2 × 1017 cm−2 over a ∼5 km s−1 channel width, allowing for studies of the distribution of H i gas in various environments, such as in galaxies, the Circum-Galactic Medium (CGM) and Intergalactic Medium (IGM). We visually searched the spectral cubes for H i sources, and found a total of 80 H i detections, of which 56 have been cross-matched with the MIGHTEE-H i catalogue. With the cross-matched sources, we compare their H i masses and find that the total H i mass fraction in the IGM and CGM surrounding the galaxy pairs is statistically higher than the H i fraction surrounding the isolated galaxies by a difference of ∼13 ± 4%, indicating that the CGM and IGM associated with interacting systems are richer in neutral hydrogen compared to those around isolated galaxies in the local Universe. We also describe several FAST-MeerKAT synergy projects, highlighting the full potential of exploiting both single-dish and interferometric observations to study the distribution and evolution of the diffuse H i gas.

Hidden Companions to Intermediate-mass Stars. XXII. Inferring a 0.43M ⊙, 1.65 au Companion in the Astrometric Binary within the Triple System Lambda Muscae* *Based on observations collected at the European Southern Observatory, Chile, Program ID 111.24UP.001.

Lambda Muscae is a nearby A-type star reported as a P = 474 days astrometric binary in Gaia DR3 (and previously in Hipparcos) and which also has a 0.16M ⊙ M dwarf companion at a projected separation 6.4 au discovered by VLT/SPHERE. Based on a K band flux ratio upper limit of 1% from a VLTI/GRAVITY interferometric observation, we estimate a primary mass M Aa = 2.28M ⊙ and infer that the closer companion is a M Ab = 0.43M ⊙ star (almost certainly an M dwarf) orbiting at a = 1.65 au. The triple system might be dynamically interesting given the low hierarchy and the ratio of angular momenta of the inner and outer orbits of about one.

ALMA follow-up of ∼ 3000 red-Herschel galaxies: The nature of extreme submillimeter galaxies

Abstract We present the analysis of over 3000 red-Herschel sources ($S_{\mathrm{250\, \mu m}}<S_{\mathrm{350\, \mu m}}<S_{\mathrm{500\, \mu m} }$) using public data from the ALMA archive and the Herschel-ATLAS survey. This represents the largest sample of red-Herschel sources with interferometric follow-up observations to date. The high ALMA angular resolution and sensitivity (θFWHM ∼1 arcsecond; σ1.3 mm ∼ 0.17 mJy beam−1) allow us to classify the sample into individual sources, multiple systems, and potential lenses and/or close mergers. Interestingly, even at this high angular resolution, 73 per cent of our detections are single systems, suggesting that most of these galaxies are isolated and/or post-merger galaxies. For the remaining detections, 20 per cent are classified as multiple systems, 5 per cent as lenses and/or mergers, and 2 per cent as low-z galaxies or Active Galactic Nuclei. Combining the Herschel/SPIRE and ALMA photometry, these galaxies are found to be extreme and massive systems with a median star formation rate of ∼1500 M⊙ yr−1 and molecular gas mass of Mgas ∼ 1011 M⊙. The median redshift of individual sources is z ≈ 2.8, while the likely lensed systems are at z ≈ 3.3, with redshift distributions extending to z ∼ 6. Our results suggest a common star-formation mode for extreme galaxies across cosmic time, likely triggered by close interactions or disk-instabilities, and with short depletion times consistent with the starburst-type population. Moreover, all galaxies with S1.3mm ≥ 13 mJy are gravitationally amplified which, similar to the established $S_{500\mathrm{ \, \mu m}}>100$ mJy threshold, can be used as a simple criterion to identify gravitationally lensed galaxies.

Hidden Companions to Intermediate-mass Stars. XXI. Discovery of a 0.5 M ⊙, 1.5 au Companion to Sadachbia = Gamma Aquarii* *Based on observations collected at the European Southern Observatory, Chile, Program ID 113.26GM.001.

Sadachbia = Gamma Aquarii is a nearby A0V star that was suspected to be a spectroscopic binary. Here we report on the direct detection of a companion with a K band flux ratio of 1.2% at a projected separation of 38.2 mas based on a VLTI/GRAVITY near-infrared interferometric observation. The implied mass and separation of the companion are 0.46 M ⊙ and 1.5 au based on the Gaia DR3 distance (or 0.55 M ⊙ and 1.9 au if the discrepant Hipparcos distance is used instead). The companion explains the X-ray emission of γ Aqr and the orbital motion with a period of about one year likely explains the discrepant Gaia and Hipparcos parallaxes. Further observations are necessary to measure the binary semimajor axis but it is likely that the system will undergo common envelope evolution. The primary is partially resolved with an angular diameter of 0.52 mas.

Orbits and Dynamical Masses for Six Binary Systems in the Hyades Cluster

We report long-baseline interferometric observations with the CHARA Array that resolve six previously known double-lined spectroscopic binary systems in the Hyades cluster, with orbital periods ranging from 3 to 358 days: HD 27483, HD 283882, HD 26874, HD 27149, HD 30676, and HD 28545. We combine those observations with new and existing radial-velocity measurements, to infer the dynamical masses for the components as well as the orbital parallaxes. For most stars, the masses are determined to be better than 1%. Our work significantly increases the number of systems with mass determinations in the cluster. We find that, while current models of stellar evolution for the age and metallicity of the Hyades are able to reproduce the overall shape of the empirical mass–luminosity relation, they overestimate the V-band fluxes by about 0.1 mag between 0.5 and 1.4 M ⊙. The disagreement is smaller in H, and near zero in K, and depends somewhat on the model. We also make use of the TESS light curves to estimate rotation periods for our targets, and detect numerous flares in one of them (HD 283882), estimating an average flaring rate of 0.44 events per day.

Multiline study of the radial extent of SiO, CS, and SiS in asymptotic giant branch envelopes

Circumstellar envelopes around asymptotic giant branch (AGB) stars contain a rich diversity of molecules, whose spatial distribution is regulated by different chemical and physical processes. In the outer circumstellar layers, all molecules are efficiently destroyed due to interactions with interstellar ultraviolet photons. Here we aim to carry out a coherent and uniform characterization of the radial extent of three molecules (SiO, CS, and SiS) in envelopes around O- and C-rich AGB stars, and to study their dependence on mass-loss rate. To this end, we observed a reduced sample of seven M-type and seven C-type AGB envelopes in multiple lines of SiO, CS, and SiS with the Yebes 40 m and IRAM 30 m telescopes. The selected sources cover a wide range of mass-loss rates, from ~10−7 M⊙ yr−1 to a few times 10−5 M⊙ yr−1, and the observed lines cover a wide range of upper-level energies, from 2 K to 130 K. We carried out excitation and radiative transfer calculations over a wide parameter space in order to characterize the abundance and radial extent of each molecule. A χ2 analysis indicates that the abundance is usually well constrained while the radial extent is in some cases more difficult to constrain. Our results indicate that the radial extent of the molecules considered here increases with increasing envelope density, in agreement with previous observational findings. At high envelope densities of Ṁ/υ∞ > 10−6 M⊙ yr−1 km−1 s, SiO, CS, and SiS show a similar radial extent, while at low envelope densities of Ṁ/υ∞ < 10−7 M⊙ yr−1 km−1 s, differences in radial extent appear among the three molecules, in agreement with theoretical expectations based on destruction due to photodissociation. At low envelope densities, we find a sequence of increasing radial extent, SiS → CS → SiO. We also find a tentative dependence of the radial extent on the chemical type (O- or C-rich) of the star for SiO and CS. Interferometric observations and further investigation of the photodissociation of SiO, CS, and SiS should provide clarification of the situation in regards to the relative photodissociation radii of SiO, CS, and SiS in AGB envelopes and their dependence on envelope density and C/O ratio.

Deciphering the imprint of active galactic nucleus feedback in Seyfert galaxies. Nuclear-scale molecular gas deficits

We study the distribution of cold molecular gas in the circumunuclear disks (CND; $r pc) of a sample of 64 nearby L $=7-45 Mpc) disk galaxies ---including 45 active galactic nuclei (AGN) and 19 nonAGN--- for which high-spatial-resolution (median value $ pc) multiline CO interferometer observations have been obtained at millimeter wavelengths with the Atacama Large Millimeter Array (ALMA) and/or Plateau de Bure Interferometer (PdBI). We decipher whether or not the concentration and normalized radial distribution of cold molecular gas change as a function of X-ray luminosity in the 2--10 keV range X $) in order to analyze the imprint left by AGN feedback. We also look for similar trends in the concentration and normalized radial distribution of the hot molecular gas and in the hot-to-cold-molecular gas mass ratio in a subset of 35 galaxies using near-infrared (NIR) integral field spectroscopy data obtained for the H$_2$ 1-0 S(1) line. We find a significant turnover in the distribution of the cold molecular gas concentration as a function of X-ray luminosity with a breakpoint that divides the sample into two branches: (1) the `AGN build-up branch'' ($L_ X pm0.3 $erg s$^ $) and (2) the ``AGN feedback branch'' ($L_ X pm0.3 $erg s$^ $) . Lower-luminosity AGN and nonAGN of the AGN build-up branch show high cold molecular gas concentrations and centrally peaked radial profiles on nuclear ($r pc) scales. Higher-luminosity AGN of the AGN feedback branch show a sharp decrease in the concentration of molecular gas and flat or inverted radial profiles. The cold molecular gas concentration index ($CCI$) --- defined as the ratio of surface densities at $r pc ($ gas $) and $r pc $( gas $), namely $CCI gas gas $)--- spans a 0.63 dex range, equivalent to a factor simeq 4-5, between the galaxies lying at the high end of the AGN build-up branch and the galaxies showing the most extreme nuclear-scale molecular gas deficits in the AGN feedback branch. The concentration and radial distributions of the hot molecular gas in our sample follow qualitatively similar but less extreme trends as a function of X-ray luminosity. As a result, we find higher values of hot-to-cold molecular gas mass ratios on nuclear scales in the highest luminosity AGN sources of the AGN feedback branch. These observations confirm ---with a three times larger sample--- previous evidence found in the context of the Galaxy Activity Torus and Outflow Survey (GATOS) that the imprint of AGN feedback on the CND-scale distribution of molecular gas is more extreme in higher luminosity Seyfert galaxies of the local Universe.

An impressionist view of V Hydrae

Context. Asymptotic giant branch (AGB) stars enrich the interstellar medium through their mass loss. The mechanism(s) shaping the circumstellar environment of mass-losing stars is not clearly understood so far. Aims. Our purpose is to study the effect of binary companions located within the first 10 stellar radii from the primary AGB star. In this work, we target the mass-losing carbon star V Hydrae (V Hya) and search for signatures of its companion in the dust-forming region of the atmosphere. Methods. The star was observed in the L and N bands with the VLTI/MATISSE instrument at low spectral resolution. We reconstructed images of the photosphere and surroundings of V Hya using the two bands and compared our interferometric observables with VLTI/MIDI and VISIR archival data. To constrain the dust properties, we used the 1D radiative transfer code DUSTY to model the spectral energy distribution. Results. The star is dominated by dust emission in the L- and N-bands. The MATISSE reconstructed images show asymmetric and elongated structures in both infrared bands. In the L band, we detected an elongated shape of approximately 15 mas that likely is of photospheric origin. In the N band, we found a 20 mas extension northeast from the star and perpendicular to the L-band elongated axis. The position angle and the size of the N-band extension match the prediction of the companion position at the MATISSE epoch. By comparing MATISSE N-band with MIDI data, we deduce that the elongation axis in the N-band has rotated since the previous interferometric measurements 13 yr ago, supporting the idea that the particle enhancement is related to the dusty clump moving along with the companion. The VISIR image confirms the presence of a large-scale dusty circumstellar envelope surrounding V Hya. Conclusions. The MATISSE images unveil the presence of a dust enhancement at the position of the companion. This opens new doors for further analyses of the binary interaction with an AGB component.

Trajectory tracking PID passivity-based control of spacecraft formation flying around Sun-Earth L2 point in the port-Hamiltonian framework

Spacecraft formation flying for interferometric observations around the 2nd Lagrange point in the Sun-Earth system (SEL2) is currently focal point in deep space exploration research, which demands high precision in relative position control of spacecraft. This paper proposes the relative motion dynamics of satellite formations around the L2 point in the port-Hamiltonian framework, and establishes a high-precision nonlinear dynamics model. PID passivity-based control (PID-PBC) is widely used in engineering. However, existing methods of PID-PBC cannot address the trajectory tracking issues in port-Hamiltonian systems. Utilizing the contraction properties of the port-Hamiltonian system, this paper proposes the trajectory tracking PID-PBC (tPID-PBC) approach, effectively resolving trajectory tracking issues for formation dynamics around L2 point in port-Hamiltonian framework. The paper details explicit solutions of the Partial Differential Equations (PDE) for the tPID-PBC method and its controller structure, and references the trajectories of interferometric observation formations, to verify the method’s effectiveness through numerical simulation. The presented control approach is applicable across generic port-Hamiltonian systems, offering substantial theoretical value.

Catalogue of dual-field interferometric binary calibrators

Context. Dual-field interferometric observations with VLTI/GRAVITY sometimes require the use of a binary calibrator. This is a binary star whose individual components remain unresolved by the interferometer, with a separation between 400 and 2000 mas for observations with the Unit Telescopes (UTs), or 1200–3000 mas for the Auxiliary Telescopes (ATs). The separation vector also needs to be predictable to within 10 mas for a proper pointing of the instrument. Aims. No list of properly vetted calibrators was available so far for dual-field observations with VLTI/GRAVITY on the UTs. Our objective is to compile such a list and make it available to the community. Methods. We identified a list of candidates from the Washington Double Star (WDS) catalogue, all with appropriate separations and brightness, scattered over the southern sky. We observed them as part of a dedicated calibration programme, determined whether these objects were true binaries (excluding higher multiplicities resolved interferometrically, but unseen by imaging), and extracted measurements of the separation vectors. We combined these new measurements with those available in the WDS to determine updated orbital parameters for all our vetted calibrators. Results. We compiled a list of 13 vetted binary calibrators for observations with VLTI/GRAVITY on the UTs, and we provide orbital estimates and astrometric predictions for each of them. We show that our list guarantees that there are always two binary calibrators at least at an airmass < 2 in the sky over the Paranal observatory at any point in time. Conclusions. Any principal investigator wishing to use the dual-field mode of VLTI/GRAVITY with the UTs can now refer to this list to select an appropriate calibrator. We encourage the use of whereistheplanet to predict the astrometry of these calibrators, which seamlessly integrates with p2Gravity for VLTI/GRAVITY dual-field observing material preparation.