Open Access
Issue
A&A
Volume 682, February 2024
Article Number A165
Number of page(s) 21
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/202348059
Published online 23 February 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1 Introduction

The study of protoplanetary disks is of crucial importance for our understanding of the initial phases of the formation and evolution of planetary systems. In the last few years, new instruments with higher spatial resolution and sensitivity have revealed the complex structure of protoplanetary disks, showing the presence of spirals, rings, gaps, and shadows from a few tens of au up to ~200 au from the central source (Beuzit et al. 2008, 2019; ALMA Partnership 2015; Andrews et al. 2018; Long et al. 2018; Benisty et al. 2015, 2017; de Boer et al. 2016; Pohl et al. 2017; Avenhaus et al. 2018). These structures are thought to arise from different phenomena, such as the presence of planets, winds, and/or disk instabilities (see e.g. the recent review by Benisty et al. 2023). Most of them have been detected as close as a few au from the central source. However, it is not yet clear whether they extend down to the innermost disk region, within ~1–2 au of the source. Some indirect evidence of the presence of structures in the innermost disk is, however, given by the so-called shadows. These regions are dark areas observed in spectro-polarimetric images obtained with extreme adaptive optics instruments such as SPHERE. They are interpreted as shadows due to misalignments between the inner and outer disk and/or the presence of scale-height variations within the disk, such as warps (Pinilla et al. 2018; Benisty et al. 2018; Bohn et al. 2022). The presence of these shadows thus allows us to place strong, albeit indirect, constraints on the complex structure of the inner disk.

Obtaining direct evidence of the inner disk structure is, however, very difficult. This is due to the small spatial scales involved, meaning optical/IR interferometry is the only technique able to spatially resolve the innermost disk. Despite the complexity of this technique, the advent of new interferometric instruments such as VLTI-GRAVITY and VLTI-MATISSE, which have higher sensitivity and better coverage than first generation interferometers, is allowing us now to fill the gap. For example, recent spectro-interferometric surveys have revealed the presence of asymmetric H- and K-band continuum emission, which can be interpreted as unevenly illuminated rims (Lazareff et al. 2017; GRAVITY Collaboration 2019). Moreover, recent VLTI-GRAVITY and VLTI-MATISSE observations hint at the presence of rotating structures within a few au of the central source (Varga et al. 2021; GRAVITY Collaboration 2021). Whereas geometrical modelling can be used to constrain the general structure of the disk, it cannot provide details on the nature of such asymmetries because geometrical models are too simple to catch the complexity of such structures. Indeed, recent reconstructed images at near-IR wavelengths have shown different complex morphologies in the innerdisk, demonstrating that disks at sub-au spatial scales are not symmetric at all and that substructures are often present (Labdon et al. 2019; Kluska et al. 2020; GRAVITY Collaboration 2021; Hofmann et al. 2022). In this context, the task becomes even more complex and/or difficult if one considers the contribution from the gas component. Indeed, within 1 au of the source, a dust-free disk extends towards the star, and matter accretes onto the source and/or is blown away in the form of winds. Therefore, a complete picture of both gas and dust emission is only possible through direct imaging as the overall picture is too complicated to be reproduced by simple centro-symmetric geometrical models.

So far, only one image in line emission has been reconstructed using near-IR interferometry for a young stellar object (YSO), namely the HI Brγ line imaged around the Herbig Be star MWC 297 (Hone et al. 2017). In this case, the authors find that the HI Brγ line emission in MWC 297 is consistent with having originated at the base of a disk wind, showing for the first time that kinematic effects in the sub-au inner regions of a protoplanetary disk can be directly imaged. It should also be noted that this is the only young star for which the HI Brγ line emission is more extended than the K-band continuum emission, and therefore it is yet to be determined whether the same results can be applied to other sources with more compact Brγ emission (e.g. Caratti o Garatti et al. 2015; Garcia Lopez et al. 2015; Kurosawa et al. 2016). Sources of HI Brγ line emission other than (or complementary to) disk winds are, for instance, the accreting gas and/or the disk itself.

In this paper, we present the first results of an interfer-ometric campaign carried out to obtain images of a sample of intermediate- to high-mass stars obtained with the VLTI-GRAVITY instrument as part of the YSO guaranteed time observation (GTO) programme. Here, we focus on imaging the K-band continuum and HI Brγ line emission around the Herbig star HD 58647. This object is a bright Herbig B9 IV (Mora et al. 2001) star at a distance of 302 pc (Gaia Collaboration 2021, 2016) with a luminosity of about 275 L (Vioque et al. 2018). Strong hydrogen emission lines have been detected and studied in detail by Manoj et al. (2002), Brittain et al. (2007), and Harrington & Kuhn (2009). One of the distinguishing features of these lines is their double-peaked profile, which hints at gas in Keplerian rotation.

Previous optical interferometric studies of the circumstel-lar disk of HD 58647 in the H and K bands (Lazareff et al. 2017; GRAVITY Collaboration 2019) indicated the presence of an elongated continuum emission with a size of ~3–4 milliarc-seconds (mas), an inclination of ~65°, and a position angle (PA) of ~15°. Interestingly, VLTI-AMBER spectro-interferometric observations successfully reproduced the double-peaked Brγ line profile as well as the interferometric observables via a disk-wind model (Kurosawa et al. 2016). However, the lack of sensitivity and the limited uv coverage of VLTI-AMBER interferometric observations did not allow for an image reconstruction of the emission line.

This paper is a step further in studying the innermost region of HD 58647. We do so by simultaneously reconstructing an image of the K-band dust continuum and of HI Brγ line emission using the European Southern Observatory(ESO)-VLTI beam combiner GRAVITY. In particular, we use our data to spatially and spectrally resolve both the dust and the hot gas component in the inner disk region within 1 au of the young star.

The paper is organised as follows: Sects. 2 and 3 present our GRAVITY GTO observations and data reduction. Sections 4 and 5 show the analysis and results of the continuum and the Brγ line emission. In Sect. 7, we present our discussion and in Sect. 8 a summary of our findings.

2 Observations and data reduction

HD 58647 was observed with the ESO-VLTI instrument GRAVITY (Eisenhauer et al. 2011; GRAVITY Collaboration 2017) as part of the YSO GTO. A full log of the observations can be found in Table A.1. Our observations were performed in six different runs, between January 2020 and February 2021. We used the four 1.8 m auxiliary telescopes, resulting in six baselines for each of the three adopted array configurations for a range of baseline sizes of ~11 to 132 m. The data were recorded using both the fringe tracker (FT) and the science channel (SC) detectors covering the K-band range from 1.9 µm to 2.4 µm. The FT data were recorded at low spectral resolution (R ~ 20) at a typical frame rate of 1 kHz (Lacour et al. 2019), allowing the atmosphere effects to be frozen. The SC data were recorded at high spectral resolution (HR; R ~ 4000) using an average DIT of 30 s for a rough total exposure time of 300 s per frame, and each dataset has from 7 to 19 (N in Table A.1) frames (see Table A.1 for more details). The complete u–υ coverage of our observations is shown in Fig. 1, which displays sufficient sampling to attempt image reconstruction.

The data were reduced using the GRAVITY data reduction software v.1.3.0 (Lapeyrere et al. 2014). The transfer function was estimated using a calibrator star (see Table A.1 for a full list of calibrators). The spectrum of the calibrator was also used to correct for the telluric absorption features present in the spectrum. The wavelength calibration of our data was refined by using the many telluric features present in the source spectrum before applying the telluric correction.

Despite the time gap between the observations, analysis of the dataset revealed no discernible temporal variability, as depicted in Fig. B.1. This allowed us to combine the datasets, which resulted in a more robust and comprehensive analysis.

Table 1

Stellar properties from (Vioque et al. 2018); the stellar luminosity is re-scaled to the Gaia EDR 3 distance (Gaia Collaboration 2021).

thumbnail Fig. 1

plane coverage of our observations of HD 58647. See Table A.1 for a detailed description of the observations.

3 Interferometric observables

Each of our interferometric observations provides us with the following observables: spectrum, six spectrally dispersed visibilities, four closure phases, and six differential phases. Visibilities (V) measure the size of the emitting region, with V = 1 and V = 0 indicating spatially unresolved, or fully resolved emission, respectively. The differential phase measures the photocentre shift of the line with respect to the continuum, whereas closure phases give a measure of the departure from centro-symmetry of the emission. An example of one GRAVITY dataset around the position of the H I Brγ lines is shown in Fig. 2. The full sample can be found in Fig. B.1.

The spectrum of HD 58647 shows bright Brγ line emission (red spectrum of Top panel of Fig. 2). Interestingly, the Brγ line profile is double peaked with the central absorption roughly centred at zero velocity. The blue- and redshifted maxima peak roughly at ±70 km s−1. In order to account for the full H I Brγ emission, the contribution from the intrinsic Brγ photospheric feature was removed by using a photospheric template spectrum of the same effective temperature, surface gravity and υ sin i as HD 58647. A full description of the procedure can be found in Garcia Lopez et al. (2006). The stellar parameters and υ sin i values were taken from Vioque et al. (2018) and Montesinos et al. (2009), respectively (see Table 1). The corrected Brγ spectrum is shown in the top panel of Fig. 2 (blue spectrum).

In the middle-top panel of Fig. 2, we observe an overall increase in the visibilities across the Brγ line. This suggests that the region emitting the Brγ line is more compact compared to the adjacent K-band continuum emission. Notably, the visibilities across the line are spectrally resolved and they exhibit a double-peaked shape, with a dip at lower velocities (υ ~ 0 km s−1), as marked by the red-dashed line.

The differential phases (Fig. 2, middle-bottom panel) show the characteristic S-shaped profile across the Brγ line, suggesting the presence of rotating gas. The S-shaped profile is observed in all of our baselines down to a projected baseline of ~ 10 m. Finally, a closure-phase signature of up to ~ 10° across the line is observed (Fig. 2, bottom panel). The same S-shaped profile is observed in the closure phase at the triangle configurations with the largest baseline lengths indicating departure from centro-asymmetry.

thumbnail Fig. 2

Example of HR GRAVITY interferometric observations of HD 58647 around the H I Brγ line position taken in January 2020 using the (J3-J2-G2-A0) configuration. From top to bottom: HR spectrum (red) and photospheric-corrected spectrum (blue); spectrally dispersed visibility amplitudes; differential phases; and closure phases. Different colours represent different projected baselines and baseline orientations (PAs; visibilities and differential phase panels) or triplets (closure phase panel), as indicated in the middle-top and bottom panel, respectively.

4 Imaging the K-band continuum emission around HD 58647

Interferometric imaging is a powerful tool for obtaining a model-independent characterisation of the brightness distribution of an object in the sky. However, interferometric imaging is a challenging task as this is an ill-posed problem due to the sparsity of the u–υ coverage and the non-linearity between the measured squared visibilities and closure phases. A full description of optical interferometric imaging can be found in Thiébaut & Young (2017) and Sanchez-Bermudez et al. (2018). In short, image reconstruction can be thought as minimising the following expression: x=argmin{ 12χ2(x)+inμiR(x)i },${\bf{x}} = argmin\left\{ {{1 \over 2}{\chi ^2}(x) + \mathop {{{\mathop \sum \nolimits^ }^}}\limits_i^n {\mu _i}R{{(x)}_i}} \right\},$(1)

where the first term is the distance to the data (also known as the negative log-likelihood), the second term is the prior term, and µ is an hyper-parameter expressing the trade-off between the two previous terms. The solution is the sought image that best reproduces our data.

In this paper, the image reconstruction algorithm SQUEEZE (Baron & Kloppenborg 2010) is used to solve Eq. (1). Using stochastic methods, the sought image is calculated as a combination of a large number of flux elements, ranging from 1000 to 10 000, that move randomly within the grid iteratively via simulated annealing algorithms and changing the value of the criterion function in Eq. (1) during the exploration of the posterior distribution. The evolution of the solution as the elements move forms a Markov chain. After the stabilisation of the chain, the posterior distribution is still being explored around the best solution. This allows a mean image to be calculated from the different chains and a standard deviation map that quantifies the significance of each pixel in our mean image to be produced. Our image was reconstructed from the HR SC data binned from the original 1634 to 7 spectral channels in order to increase the S/N of our data.

The resolution of the images is λ/2B, with B our maximum baseline. To avoid the pixellation of the reconstructed images, we set the pixel scale to 0.15 mas.

To initialise the different chains (initial length of 500 iterations), a synthetic image was created from a geometrical model fitting of all our interferometric observables. The geometrical model was also based on the SC HR data binned to seven spectral channels. A full description of the geometrical model is shown in Appendix D.

The final mean images and error maps were created from the results of all the converged chains. We used a combination of two regularizers, l0 norm and total variation (TV), for all the images in the continuum.

The reconstructed images are presented in the left panel of Fig. 3 with 3, 5, and 10 σ contours. As shown in Fig. 3 (right panel), the presented images reproduce our observables well.

The resulting image shows an elongated emission in the direction perpendicular to the elongated beam and it strongly suggests the presence of a thin disk that is oriented in the NE–SW plane with size at 3σ of R ~ 3.5 mas (i.e. ~1.1 au), inclination of i ~ 65° and PA from north to east ~15°, (Table 2). The size of the extended emission is well beyond our element of resolution (i.e. the clean beam, white ellipse at the left-bottom of Fig. 3). The clean beam is extracted from the dirty beam by Fourier transform of the uυ plane. The clean beam size is 2.6 mas × 1.9 mas with a PA ~ 122°.

The continuum emission is roughly symmetric within 10 σ (~1.5 × 0.9 mas from the source, i.e. 1.1 × 0.5 au at the distance of HD 58647), although a slight increase in the elongation with wavelength is observed in the direction perpendicular to our element of resolution. However, the extended emission is clearly asymmetric when exploring larger spatial scales (e.g. 5σ to 3σ contours). These asymmetries cannot be reproduced by simple geometrical modelling, as the one presented in Appendix D namely the models are not able to perfectly reproduce the interferometric observables, in particular, the closure phase signatures (see Fig. D.1).

Table 2

Results from the K-band continuum image reconstruction.

5 Imaging the HI Brγ line emission in HD 58647

Imaging across IR line tracers is even more challenging than obtaining images of the continuum due to the reduced number of spectral channels available to obtain an image. Therefore, very few images across the Brγ line are available in the literature, and the very few available usually involve very bright objects. Of the few examples, imaging across the H I Brγ and Brα, and He I lines has been possible in Eta Car (GRAVITY Collaboration 2018; Weigelt et al. 2021), and across the H I Brγ line in the Herbig B[e] star MWC 297 (Hone et al. 2017) the only YSO, so far.

Most of these images were reconstructed using the IRBIS method (see Hofmann et al. 2014, 2016), in combination with the differential-phase method detailed in Millour et al. (2011), Weigelt et al. (2016) in order to obtain velocity-dispersed images of the total emission (line emitting gas plus continuum) across the line. In some cases, the continuum contribution is later subtracted from the resulting image to increase the contrast of the line emitting region and to better study the gas morphology.

In the case of HD 58647, we opted to adopt a slightly different approach to recover the brightness distribution across the Brγ line. Instead of removing the continuum contribution a posteriori, after the image is reconstructed, we reconstructed the velocity-dispersed image from continuum-subtracted visibilities and differential phases. It should be noted that, even if continuum-subtracted, the displacements obtained from the observed differential phases are with respect to the continuum, assumed to be centro-symmetric. Therefore, to take into account possible asymmetries in the brightness distribution of the continuum and line emission, the Brγ line velocity-dispersed image has to be reconstructed from absolute displacements, that is, the Fourier phase across the Brγ line has to be retrieved from the data. In principle this cannot be done, as optical interferometry does not allow us to obtain a one-to-one match between visibilities and phases. However, in our case, the Fourier phase can be retrieved making use of the continuum image presented in Fig. 3, and subtracting this contribution from the observed continuum-subtracted differential phases across the line (see Appendix E in Weigelt et al. 2007 for a detailed description of the interferometric observables).

In this way, the linear combination of the Fourier phases from the continuum images and the continuum-corrected differential phases allow us to derive the absolute phases across the Brγ line. This along with the continuum-subtracted visibilities is then used to compute continuum-subtracted Brγ line velocity-dispersed images. A description on how the continuum-subtracted Brγ line visibilities and differential phases have been retrieved can be found in Appendix E.

To obtain the iso-velocity images across the Brγ line the software SQUEEZE was used in the same fashion as for the continuum image reconstruction, including the same grid parameters and regularizers. The resulting images across the Brγ line are shown in Fig. 4.

The images were only reconstructed across channels where the flux of the emission line is 10% larger than the continuum flux. This resulted in images across seven spectral channels with velocities ranging from −109 km s−1 to 110 km s−1. The right panel in Fig. 4 shows a comparison between the observed continuum-subtracted visibilities and absolute phases, and those extracted from the retrieved Brγ line images. As shown in the middle and right panels of Fig. 4, our images recover the observed quantities well within the errors.

Our recovered velocity-dispersed images show marginally resolved Brγ line emission at all spectral channels with an average size of the line emitting region at 5 σ level of ~2.3 mas ± 0.2 mas (i.e. 0.7 au ± 0.06,au or ~31.3 R*±2.7 R*). No significant variations of the size are observed as a function of the spectral channel.

However, each velocity-dispersed image shows photocen-tre shifts from the north-east direction towards the south-west from the blueshifted channels (2.1653 µm ≤ λ ≤2.1658 µm; −109 km s−1υ ≤ −72 km s−1) to the redshifted channels (2.1666 µm ≤ λ ≤ 2.1668 µm; 73 km s−1υ ≤ 110 km s−1), going through the centre of the image at the two central velocity channels (2.1661 µm≤ λ ≤ 2.1663 µm; −36 km s−1υ < 37 km s−1), suggesting gas in rotation.

To further probe the kinematic structure of the system, we computed the first moment map (M1) from the recovered velocity-dispersed Brγ line images. In doing so, only emission detected above 5 σ was considered. In addition, the new intensity maps were scaled to the flux of the observed GRAVITY spectrum. The result is shown in Fig. 5. The resulting map supports our findings from the Brγ line images shown in Fig. 4: the blueshifted emission is shifted towards the north-east, whereas the redshifted emission is shifted towards the south-west. Given the elliptical shape of the emission, and assuming that the Brγ line is emitted in a disk-like structure, inclination and PA of the Brγ line emission can be derived. Assuming a circular structure at zero inclination, an inclination of 55° ± 1° is derived.

This value is very similar to the one obtained from the continuum. In the same fashion, the major axis of the disk can be estimated by measuring the location of the blueshifted and red-shifted peaks and measuring the angle between them. We find a value of ~13°, very similar to the value of ~15° obtained from the K-band continuum image. Furthermore, the PA obtained from the H I Brγ line velocity-dispersed images, is also in agreement with the astrometric displacements computed from the continuum-corrected differential phase signatures (see a full description of this procedure in Appendix E). The Brγ line astrometric displacements (see Fig. F.1) clearly align along a straight line with PA=15° ± 5°, and show again a clear shift of the blue-and redshifted emission towards the north-east and south-west, respectively.

thumbnail Fig. 3

Results of the image reconstruction. Left: K-band continuum reconstructed images of HD 58647 at different wavelengths (reported in the labels). North is up, and east is to the left. The white hollow ellipse at the lower-left corner of the bottom panel represents the size of the clean beam. Contours represent 3, 5, and 10 σ pixel significance levels. Middle and right: squared visibilities (|V|2; middle panel) and closure phases (T3ϕ; right panel) as a function of the spatial frequency. The observed data with their corresponding error bars are represented as black dots, while the synthetic observables extracted from the reconstructed images are shown in colour.

6 H I Brγ line Keplerian disk modelling

The results from our Brγ line reconstructed images, the S-shaped differential phases, and the Brγ line astrometric displacements points towards the presence of rotating gas emitting in H I Brγ line located in the inner gaseous disk. To give more quantitative information about the hot gas kinematics, a simple Brγ Keplerian disk model was computed and compared with our results.

The Keplerian disk velocity field was derived assuming a M* = 3.87 M (Vioque et al. 2018). The model was created to match the angular and spectral resolution of our GRAVITY observations. The inclination, PA, and size of the disk model are set as free parameters.

Figure 6 shows a comparison between the observed continuum-corrected Brγ line differential phases and those obtained from our best fitting model. Our best fitting model is obtained for a Brγ line disk with an inclination and PA of ~52°4+7$\~52^\circ _{ - 4}^{ + 7}$ and ~14° ± 8°, respectively, and a radius of the emitting region of ~1.1 mas ±0.3 mas (i.e. ~0.3 au). The PA is in agreement with that found from the Brγ line image and K-band continuum reconstructed images, indicating no major misalignment between the continuum and line emitting region. As seen in Figs. 6 and F.2, the model is able to reproduce the general trend observed in our data. However, there is a slight mismatch between the model and the observed values (χv22)$\left( {\chi _v^2 \sim 2} \right)$. On the one hand, lower spatial frequencies are sampled better than the highest ones, hinting to a more complex structure at small spatial scales. On the other hand, the best match between observed and synthetic values is obtained for intermediate velocities (e.g. ~−73 km s−1, ~37 km s−1, and ~74 km s−1). Finally, the red-shifted velocities (but for the highest velocity channel) are roughly fitted better than blueshifted channel velocities, possibly hinting to the presence of an additional blueshifted velocity field. This might be also supported by the asymmetric double peaked Brγ line profile showing brighter emission at blueshifted velocities.

7 Discussion

7.1 K-band continuum emission

Figure 3 shows our continuum K-band reconstructed image (left panel) along with a comparison with our interferometric observations (middle and right panels). From this image we extracted the K-band continuum size, inclination, and PA. The inclination ~64° and PA ~ 14° derived from the image are similar to the ones found from the geometrical model presented in Appendix D and the ones reported in Kurosawa et al. (2016); Lazareff et al. (2017); GRAVITY Collaboration (2019). Despite the similarities between the final image and the simple 2D Gaussian model, the image is able to better recover the interferometric observables than the geometrical model (see Figs. 3 and D.3), especially at the largest spatial scales. This might be due to the presence of extended emission in our image (the 5-σ contour represents ~90% of the flux) that cannot be properly accounted for using simple single-component Gaussian and/or Lorentzian brightness distributions. This, in turn, allows us to provide a better estimate of the size of the K-band continuum emission. Moreover, image reconstruction allows us to reproduce the small asymmetries in the continuum as well as to directly derive the Brγ brightness distribution. Indeed, our recovered image shows that the size of the K-band continuum emission extends as far as R ~ 3.2 mas at ~90% flux level (5σ contour). This is a factor of ~1.5 more extended than the value estimated from the geometrical modelling presented in Appendix D, and derived in previous studies (e.g. Lazareff et al. 2017; GRAVITY Collaboration 2019). At high flux contributions (~70%, ~10σ contour), the size from the geometrical modelling and our image roughly reconcile to each other. Our image shows thus K-band continuum emission from R ~ 3.2 mas (i.e. ~0.97 au) down to our resolution limit of ~1.9 mas (i.e. ≤0.6 au) as estimated from the clean beam (see Sect. 4 and Fig. 3).

This range of K-band continuum emission is in agreement with hydrodynamic (HD) and magneto-hydrodynamic models of the dust silicate inner rim (see e.g. Flock et al. 2016, 2017). These models include a constant gas inflow that replenish the inner disk of small dust particles and gas, and they take into account the absorption of stellar radiation by the inner gaseous disk interior to the dust sublimation front, backwarming by IR radiation, as well as accretion stress and heating. Under these conditions, the rim has a triangular shape that extends over a relative high spatial extent, showing a tip that lays in the disk mid-plane and points towards the star. The extent of the inner rim triangular shape, both in height and radially, mainly depends on the stellar luminosity. This determines the average location of the dust sublimation radius (e.g. Dullemond & Monnier 2010; Isella & Natta 2005) and the mass accretion rate, which affects the location of the innermost tip of the triangular rim (Rrimin$R_{{\rm{rim}}}^{{\rm{in}}}$; Flock et al. 2016, 2017). In general, the higher the stellar luminosity the further the location of the rim with respect to the star is. In the same way, the higher the accretion rate the closer Rrimin$R_{{\rm{rim}}}^{{\rm{in}}}$ is located and the higher the scale-height of the rim is. In our case, the luminosity of HD 58647 is log L* =2.4 and the mass accretion rate ranges between M˙acc~3.5×107Myr1106Myr1${\dot M_{{\rm{acc}}}}\~3.5 \times {10^{ - 7}}{M_ \odot }{\rm{y}}{{\rm{r}}^{ - 1}} - {10^{ - 6}}{M_ \odot }{\rm{y}}{{\rm{r}}^{ - 1}}$ (Grant et al. 2022; Brittain et al. 2007). The combination of these values will on the one hand bring Rrimin$R_{{\rm{rim}}}^{{\rm{in}}}$ closer to the star, and on the other hand, move further away the location of the sublimation front (intended here as the location of the local maximum in aspect ratio of the IR disk photosphere; see Rrimout$R_{{\rm{rim}}}^{{\rm{out}}}$; Flock et al. 2016).

The measured Rrimout~1$R_{{\rm{rim}}}^{{\rm{out}}}\~1$ au is indeed in general agreement with the location of the rim in classical rim models including an optically thin inner disk in between the star and the rim location with dust temperatures between 1000 K and 1500 K and with or without backwarming (see e.g. Fig. 5 in Koumpia et al. 2021). On the other hand, the high accretion rates might explain the presence of K-band continuum emission within Rrimout$R_{{\rm{rim}}}^{{\rm{out}}}$ and down to at least 0.6 au. With these high accretion rates Rrimin$R_{{\rm{rim}}}^{{\rm{in}}}$ will move very close to the stellar surface. Alternatively, if the highest available value of M˙acc =106Myr1${\dot M_{{\rm{acc}}}} = {10^{ - 6}}{M_ \odot }{\rm{y}}{{\rm{r}}^{ - 1}}$ is considered, the inner gaseous disk might be optically thick (Muzerolle et al. 2004). If this is the case, dust particles might be screened from the stellar radiation and survive very close to the stellar surface, or the gas itself might even be responsible for part of the K-band continuum emission. Therefore, our images seem to support the presence of a silicate triangular shape rim as the one described in (Flock et al. 2016) extending from Rrimin0.06$R_{{\rm{rim}}}^{{\rm{in}}} \le 0.06$ au up to Rrimout~1$R_{{\rm{rim}}}^{{\rm{out}}}\~1$ au.

thumbnail Fig. 4

Results of the image reconstruction. Left: H I Brγ line reconstructed images of HD58647. North is up, and east is to the left. Contours represent 3, 5, and 10σ. Middle and right: comparison of the observed continuum-subtracted visibility amplitudes and absolute phases (black circles with error bars) with the synthetic observables extracted from the reconstructed image.

thumbnail Fig. 5

Brγ line first moment map (M1) obtained from the line velocity-dispersed images. For comparison, the 1, 3, and 5σ K-band continuum contours are overplotted in cyan, black, and yellow. The dashed red and blue lines illustrate the orientation of the semi-major axis of the continuum and the Brγline emission, respectively.

thumbnail Fig. 6

Results of the Keplerian disk model fitting. Left: velocity maps from our Brγ line Keplerian disk model. From top to bottom, we show velocity maps at: ~ −110 km s−1, ~ −73 km s−1, ~ −36km s−1, ~0 km s−1, ~37 km s−1, ~74 km s−1, and ~110 km s−1. See the main text for details. Right: comparison between the observed continuum-corrected Brγ line absolute phases (black dots) and those obtained from the Brγ line Keplerian disk model shown in left panel (coloured dots).

7.2 H I Brγ line emission

We have obtained the first H I Brγ line image reconstruction of the hot gas component around the Herbig Ae star HD 58647 and one of the few reconstructed images of the gas component around a YSO (Sect. 5). Our image reconstruction method retrieves the Fourier phases from the continuum image, allowing us to obtain a more accurate estimate of the photocentre shifts of the Brγ line. It should be recalled that the photocentre shifts obtained from the Brγ line reconstructed image are ‘absolute’, that is, they no longer are a function of the photo-centre of the line with respect to the continuum. As the Brγ line is only marginally resolved, we cannot provide information about variation of the size as a function of the spectral channel, and we can only give an average size of the Brγ line emitting region of ~2 mas, that is, ~30 R*. However, our images are sensitive to changes in the photocentre shift of the line with velocity. These changes are supported by the first moment map, as well as the astrometric displacements of the line with respect to the continuum photocentre (Figs. 5 and F.1, respectively). These measurements clearly show a shift between the blue- and red-shifted Brγ line emission pointing out to the presence of rotating gas around HD 58647.

The hot gas component as traced by the Brγ line emission would partially overlap with a fraction of the K-band continuum emission as measured in Sect. 4. This is in contrast with previous results that reported the size of the Brγ line emitting region smaller than that of the continuum (Kurosawa et al. 2016; GRAVITY Collaboration 2023) . It should be noted however, that, as already discussed in the previous section, our image reconstruction allows us to obtain a better estimate of the full extension of the emission than simple geometrical models do. That said, our limited baseline length does not allow us to spatially resolve the innermost disk regions and we cannot infer whether the line emission extends down to the stellar surface (as suggested by a boundary layer scenario) or is truncated by a magnetosphere.

Taking advantage of the photocentre displacements, and the Keplerian modelling shown in Sect. 7, we see indications that the line emission is extending down to at least ~1.1 mas, or ~15 R*. The extension of a possible magnetosphere around HD 58647 is expected to be smaller than the corotation radius in order to allow accretion to proceed. The corotation radius can be estimated as Rcor = (GM*)1/3(P*/2π)2/3 = 2.3 R*, with P* = 2πR*eq= 2.1 day. On the other hand, Järvinen et al. (2019) measured a weak longitudinal magnetic field in HD 58647 of B𝓏 ~209 G. This small magnetic field, along with the fast rotation of the star (υ sini=114 km s−1; Montesinos et al. 2009) would place the location of a possible magnetosphere very close to the stellar surface. Therefore, we could hypothesise that the extension of the accretion funnel could be from as close as 1.3 R*–1.5 R* (Kurosawa et al. 2016; Garcia Lopez et al. 2015; Caratti o Garatti et al. 2015) up to the location of the corotation, that is, ~2.3 R*. Therefore, the presence of the magnetosphere could not explain the full extension of the Brγ line emission.

Kurosawa et al. (2016) were, however, successful in reproducing VLTI-AMBER high resolution (R = 12000) observations of HD 58647 with a model including emission from a compact magnetosphere and disk wind. The compact magnetosphere would produce a small absorption centred at zero velocity, helping to better reproduce the total line flux at zero velocity, whereas the disk wind would contribute for most of the emission and it will extend up to 23.5 R*. This value is close to our measured size of the Brγ line emission of ~30 R*. If this is the case, the fact that the Brγ line profile is double peaked, and that our Keplerian disk model presented in Sect. 7 roughly reproduces the line differential phase signatures, indicates that the wind is emitted very close to the disk surface where the poloidal component is not yet important. Alternatively, we cannot discard that the surface layers of the disk itself could be at the origin of at least part of the Brγ line emission. Variations in scale-height on the origin of the Brγ line emission could also reproduce the observed asymmetric Brγ line profile. At the same time, this could also explain the small disagreement between our Keplerian disk model and the measured line photocentre shifts, as tiny variations on the depth at which the Brγ line is emitted will produce displacements from the Keplerian velocity as measured from the disk mid-plane (see e.g. Backs et al. 2023). Unfortunately, due to the limited spectral resolution, our observations do not allow us to distinguish between bound and unbound gas. On the other hand, our measured size of the Brγ line emission seems to disagree with the presence of a photo-evaporative wind as the origin of most of the emission. According to photo-evaporative models, and taken into account the mass of the star, even for the case of extreme UV photo-evaporative wind the bulk of the wind emission would be located much further out than our measurements (i.e. REUV > 5 au, that is, REUV > 200 R*).

8 Summary

In this paper, we present GRAVITY observations of the dusty and gaseous disk around the Herbig star HD 58647. Our dataset from 2020 and 2021 includes three different telescope configurations (small, medium, and large) that cover various baseline lengths and orientations (~11–132 m). These observations have allowed us to marginally resolve the target in the K band at HR. The acquired data show no variability in terms of the wavelength-dependent observables, enabling us to combine the entire dataset. Image reconstruction has proven useful in studying the morphology of the innermost region of HD 58647, including both the continuum and HI Br-gamma emission. Based on the analysis of the recovered images, we draw the following conclusions:

  • The geometry of the marginally resolved continuum K-band emission at a significance level of 3σ is similar to our result from the geometrical model. The orientation of the major axis of the disk has a PA of ~14° and an inclination of 64°;

  • Although the size of the K-band continuum emission at the 3σ level is larger than that of the geometrical model, the emission enclosed within the contour represents 90% of the total flux, which is in good agreement with the model. This model predicts that the flux fractions contributed by the unresolved star and the circumstellar disk are approximately 91%;

  • The measured extent of the emission from our continuum images aligns with the location predicted by classical rim models. These models assume the presence of a thin inner disk between the star and the rim, with dust temperatures ranging from 1000 K to 1500 K.

We also present a novel technique for recovering the brightness distribution across the Brγline and reconstructing the Brγ-dispersed images. Our findings from the analysis of the velocity-dispersed images across the Brγ line are as follows:

  • There is no significant misalignment observed between the hot gaseous disk, traced by the Brγ line emission, and the continuum disk. The orientation of the Brγ-emitting region is similar to that obtained from the astrometric displacement calculated using the differential phases and the continuum images;

  • The velocity field of HD 58647 has been recovered thanks to the reconstructed images across the Brγ line. The first moment map shows blueshifted emission towards the northeast, while the redshifted emission originates in the southwest. Both emissions are oriented along the major axis of the disk;

  • The origin of the emission is still not entirely understood. However, when compared with a disk in Keplerian rotation, we find residuals; this suggests the presence of an unresolved complex structure in the hot gaseous disk of HD 58647. This is compatible with a disk wind emitted close to the disk surface; in other words, we may be detecting the base of a disk wind.

Acknowledgements

The authors extend their gratitude to the anonymous referee for their insightful and valuable feedback. This material is based upon works supported by Science Foundation Ireland under Grant No. 18/SIRG/5597. A.C.G. has been supported by PRIN-INAF MAIN-STREAM 2017 “Proto-planetary disks seen through the eyes of new generation instruments” and from PRIN-INAF 2019 “Spectroscopically tracing the disk dispersal evolution (STRADE)”. J.S.-B. acknowledges the support received from the UNAM PAPIIT project IA 105023; and from the CONAHCyT “Ciencia de Frontera” project CF-2019/263975. This research has made use of the NASA Astrophysics Data System, CDS Astronomical Data bases SIMBAD and VIZIER (http://cdsweb.u-strasbg.fr/) and of the Jean-Marie Mariotti Center Aspro and SearchCal service (http://www.jmmc.fr/aspro).

Appendix A Observation log

The log of the observations is given in Table A.1.

Table A.1

HD58647 GRAVITY observations. UD diameters were obtained from JMMC SearchCal (Bonneau et al. 2006, 2011)

Appendix B Observational data

Figure B.1 shows the VLTI-GRAVITY observations of HD5867 around the wavelength of the Brγ line. Each figure shows the results of merged data for different dates, as displayed in the caption. There are three panels, showing the visibility amplitude (left), differential phase (middle) and differential closure phase (right). The top panels show the photospheric-corrected spectrum. The length and PA of the projected baselines are given in the left panels.

thumbnail Fig. B.1

Merged spectrum, visibilities, differential phases, and closure phases of the January 27, 2020, GRAVITY interferometric data.

thumbnail Fig. B.2

Same as Fig. B.1, but for January 28, 2020.

thumbnail Fig. B.3

Same as Fig. B.1, but for February 4, 2020.

thumbnail Fig. B.4

Same as Fig. B.1, but for December 23, 2020.

thumbnail Fig. B.5

Same as Fig. B.1, but for February 10, 2021.

thumbnail Fig. B.6

Same as Fig. B.1, but for February 21, 2021.

Appendix C Image reconstruction parameters

Table C.1

HD 58647 image reconstruction parameters.

Table C.1 shows the parameters of the image reconstructions of HD 58547.

Appendix D Geometrical modelling of the K-band continuum

We assumed that the K-band continuum detected by GRAVITY has two main components the stellar and circumstellar contribution. In order to take into account possible over-resolved emission, an over-resolved component (halo) has been also taken into account. The total complex visibility at the spatial frequencies (u, υ) can be written as described in (Lazareff et al. 2017; GRAVITY Collaboration 2019): V(u,υ)=F(λ0/λ)ksV+Fc(λ0/λ)kcVc(u,v)×e2πj(Δxu+Δyυ)(F+Fh)(λ0/λ)ks+Fc(λ0/λ)kc,$V(u,\upsilon ) = {{{F_ \star }{{\left( {{\lambda _0}/\lambda } \right)}^{ks}}{V_ \star } + {F_c}{{\left( {{\lambda _0}/\lambda } \right)}^{kc}}{V_c}\left( {u',v'} \right) \times {{\rm{e}}^{ - 2\pi j(\Delta xu + \Delta y\upsilon )}}} \over {\left( {{F_ \star } + {F_h}} \right){{\left( {{\lambda _0}/\lambda } \right)}^{ks}} + {F_c}{{\left( {{\lambda _0}/\lambda } \right)}^{kc}}}},$(D.1)

where F, Fc, Fh are the different flux contributions to the total flux from the star, the circumstellar environment and the halo, respectively. The total flux is normalised at the reference wavelength λ0. At the distance of HD 58647 the stellar component is unresolved, and therefore V = 1. Finally, ks and kc are the spectral indices of the stellar and the circumstellar components defined as ki=dlogFvdlogv.${k_i} = {{d\log {F_v}} \over {d\log v}}.$(D.2)

As the stellar photosphere can be approximated by a black-body, the spectral index value of the star ks is calculated from the reported surface temperature in Table 1 and considered as a fix parameter in the fitting. In order to derive the size, inclination and PA of the circumstellar emission, a 2D Gaussian brightness distribution was assumed: Vc(u,υ)=exp((πau2+υ2)24ln2),${V_c}\left( {u',\upsilon '} \right) = \exp \left( { - {{{{\left( {\pi a\sqrt {{{u'}^2} + {{\upsilon '}^2}} } \right)}^2}} \over {4\ln 2}}} \right),$(D.3)

where u′, υ′ are the spatial frequencies corrected for the PA of the semi-major axis of the disk and its inclination, i, and u=ucosPA+υsinPAυ=(usinPA+υcosPA)cosi.$\matrix{ {} \hfill & {u' = u\cos PA + \upsilon \sin PA} \hfill \cr {} \hfill & {\upsilon ' = ( - u\sin PA + \upsilon \cos PA)\cos i} \hfill \cr } .$(D.4)

In this model, the point-like source is centred in the phase reference point and the 2D Gaussian is displaced by ∆x and ∆y in right ascension and declination to account for possible displacements from centro-symmetry.

Table D.1

K-band continuum model-fitting results with 1σ error bars.

In order to increase the signal-to-noise ratio (S/N), we re-binned the SC HR observables from 1634 to 7 spectral channels. The differential observable across the Brγ line were bracketed to keep only the continuum contribution. A Levenberg-Marquardt algorithm implemented in lmfit (Newville et al. 2018) was used for the minimisation. The squared visibilities and the closure phases were fitted simultaneously at all spectral channels. The results of our best fitting model are reported in Table D.1. A comparison between the observations and the best-fitting geometrical model is displayed in Fig. D.1.

thumbnail Fig. D.1

Comparison between the observed squared visibilities (|V|2) and closure phase signatures (T3ϕ) of HD 58647 (black dots and error bars) with the results from a 2D Gaussian brightness distribution (coloured dots).

Our geometric model-fitting gives a size (FWHM) of 3.69±0.02 mas (i.e. 1.114±0.006 au at a distance of 302 pc), and an inclination of i=64.71º±0.12º, with a disk major-axis PA (N to E) of 13.67º±0.14º. We note that the error bars on the derived parameters resulting from the model fitting represent only the formal errors arising from the fitting procedures, and do not account for other sources of error.

The derived size is in agreement with GRAVITY Collaboration (2019). As, expected our continuum K-band size (FWHM~ 3.7mas) is larger than the H-band size (FWHM~3 mas), which is consistent with a cooler emission in the K-band than in the H-band. In addition, our modelling suggests the presence of more extended emission as traced by the over-resolved component in our model. The flux ratio associated with this component is of ~9% of the total flux. Again this is in agreement with previous H-band results (Lazareff et al. 2017), although the H-band over-resolved component has a much smaller contribution to the total flux than the one reported here (flux ratio of 2% vs the 9% observed in our data). The origin of this over-resolved component is still under debate and it might be attribute to the presence of scattered light (e.g. Pinte et al. 2008; Benisty et al. 2010). In order to obtain more information about this component a better u-v coverage as well as additional information from short baselines are needed in order to constrain its nature and morphology.

thumbnail Fig. D.2

Brγ line model-fitting results (red dots) as a function of the spatial frequency at three different spectral channels. For comparison, the continuum-corrected Brγ line visibilities are overplotted in black.

Finally, it should be pointed out that although our geometrical model is roughly reproducing the general trend of the visibilities and closure phase signatures with spatial frequency, it is poorly recovering their exact values at some spatial frequencies, especially at the shortest baselines (see Fig. D.1). This is reflected by the observed residuals (χv24.5)$\left( {\chi _v^2 \sim 4.5} \right)$ and can be explained by the fact that our model is too simple to reproduce the complexity of the structure at the largest spatial scales probed by those baselines. To fully investigate the nature of the departures from centro-symmetry, image reconstruction is needed (see Sect. 4).

Appendix E Brγ line continuum-corrected visibilities and differential phases

The HR differential visibilities show an increase across the Brγ line at all baselines indicating that the line emitting region is more compact than the continuum. However, the observed differential observables are contaminated by the continuum emission and thus this contribution needs to be subtracted. In order to estimate the size of the Brγ line emitting region we took the contribution from the continuum emission into account and calculated the continuum-corrected visibilities (see Weigelt et al. (2007) for more details): | FlineVline |2=| FtotVtot |2+| FcVc |22FtotVtotFcVccosϕ,${\left| {{F_{{\rm{line}}}}{V_{{\rm{line}}}}} \right|^2} = {\left| {{F_{{\rm{tot}}}}{V_{{\rm{tot}}}}} \right|^2} + {\left| {{F_c}{V_c}} \right|^2} - 2{F_{{\rm{tot}}}}{V_{{\rm{tot}}}}{F_c}{V_c}\cos \phi ,$(E.1)

where Ftot is the total observed flux, Vtot the observed visibility, Fline the flux line, Vc and Fc the continuum visibility and the continuum flux inside the Brγ line, and ϕ describes the observed differential phases.

We used a polynomial interpolation to calculate the flux of the continuum across the Brγ line assuming that the flux of the continuum outside the line is equal to the level of the continuum inside the line. The visibility amplitude of the continuum across the line was also estimated from a polynomial interpolation. The continuum channels nearby the Brγ line were used to perform the interpolation.

We calculated the line visibility for seven spectral channels where the flux of the line is 10% higher than that of the continuum flux and the errors are estimated by propagating Eq. E.1.

We used the 1σ errors as estimated by the data reduction software for the visibility amplitudes and phases. We assumed that the Brγ line emitting region has the following components: an unresolved component with a flux ratio described by fs, and over-resolved component fh, and finally a 2D Gaussian (similar to Eq. D.3) with visibility V associated with a flux ratio fc such as fs + fh + fc = 1 describing the resolved component. The total visibility is then given by Vamp(u,υ)=fs+(1fsfh)V(u,υ).${V_{amp}}(u,\upsilon ) = {f_s} + \left( {1 - {f_s} - {f_h}} \right)V\left( {u',\upsilon '} \right).$(E.2)

We chose to fit the visibilities across the line at blue (channels at ~ −72 km s−1; 2.1655), red (channels at ~0 km s−1; 2.1661 µm), and zero (channels at ~ 109 km s−1; 2.1668 µm) velocities. The model-fitting parameters are listed in Table E.1. The blue-and redshifted line emitting regions have a similarsizes, inclinations and PA. The results of our fitting show that about 50% of the flux is coming from an unresolved component with negligible contribution from an over-resolved component. However, the results shown in Table E.1 should be taken with caution. Unlike the K-band continuum model-fitting, our simple 2D Gaussian brightness distribution does not reproduce the complex geometry of the line emitting region and only the visibility amplitude is modelled (see Fig. D.2 for a comparison between the model-fitting results and the observations). In order to further probe the structure of the system image reconstruction is needed (see Sect. 5).

Table E.1

Brγ line model-fitting results.

Appendix F Continuum-corrected differential phases

The observed differential phases (see Fig. 2 bottom panel) show deviation from zero and allow us to study the displacement of the photocentre of the emission across the Brγ line with respect to the continuum. In order to analyse any displacement of the photocentre of the emission, the differential phase needs to be corrected for the continuum emission: sinϕline=sinϕ| FtotVtot || FlineVline |,$\sin {\phi _{{\rm{line}}}} = \sin \phi {{\left| {{F_{{\rm{tot}}}}{V_{{\rm{tot}}}}} \right|} \over {\left| {{F_{{\rm{line}}}}{V_{{\rm{line}}}}} \right|}},$(F.1)

where ϕline is the continuum-corrected differential phase and ϕ is the observed differential phase (see Weigelt et al. 2007, for more details).

By combining all the continuum-corrected differential phases, we can derive a 2D displacement of the photocentre (see Lachaume 2003; Le Bouquin et al. 2009, for additional details): P=ϕlinei2πλBi.$P = - {{\phi _{line}^i} \over {2\pi }} \cdot {\lambda \over {{B_i}}}.$(F.2)

The astrometric displacements vectors projected on the sky are shown in Fig. F.1. The figure shows all displacements lying along a straight line, with a clear measurement of the displacement between the red- and blueshifted velocities. By a linear fit of the displacements, a PA of the line emitting region of ~15°±5° is estimated. This value is in agreement with that derived from the geometric model of the continuum disk.

thumbnail Fig. F.1

Astrometric displacements obtained from the Brγ line continuum-corrected differential phases.

thumbnail Fig. F.2

Example of the synthetic differential phases obtained from the Brγ Keplerian disk model (red dots). The observed continuum-corrected Brγ line phases are shown in black.

References

  1. ALMA Partnership (Brogan, C. L., et al.) 2015, ApJ, 808, L3 [Google Scholar]
  2. Andrews, S. M., Huang, J., Pérez, L. M., et al. 2018, ApJ, 869, L41 [NASA ADS] [CrossRef] [Google Scholar]
  3. Avenhaus, H., Quanz, S. P., Garufi, A., et al. 2018, ApJ, 863, 44 [NASA ADS] [CrossRef] [Google Scholar]
  4. Backs, F., Poorta, J., Rab, C., et al. 2023, A&A, 671, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Baron, F., & Kloppenborg, B. 2010, SPIE Conf. Ser., 7734, 77344D [NASA ADS] [Google Scholar]
  6. Benisty, M., Natta, A., Isella, A., et al. 2010, A&A, 511, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Benisty, M., Juhasz, A., Boccaletti, A., et al. 2015, A&A, 578, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Benisty, M., Stolker, T., Pohl, A., et al. 2017, A&A, 597, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Benisty, M., Juhász, A., Facchini, S., et al. 2018, A&A, 619, A171 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Benisty, M., Dominik, C., Follette, K., et al. 2023, Protostars and Planets VII, ASP Conference Series, 534, 605 [NASA ADS] [Google Scholar]
  11. Beuzit, J.-L., Feldt, M., Dohlen, K., et al. 2008, SPIE Conf. Ser., 7014, 701418 [Google Scholar]
  12. Beuzit, J. L., Vigan, A., Mouillet, D., et al. 2019, A&A, 631, A155 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Bohn, A. J., Benisty, M., Perraut, K., et al. 2022, A&A, 658, A183 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Bonneau, D., Clausse, J. M., Delfosse, X., et al. 2006, A&A, 456, 789 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Bonneau, D., Delfosse, X., Mourard, D., et al. 2011, A&A, 535, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Brittain, S. D., Simon, T., Najita, J. R., & Rettig, T. W. 2007, ApJ, 659, 685 [NASA ADS] [CrossRef] [Google Scholar]
  17. Caratti o Garatti, A., Tambovtseva, L. V., Garcia Lopez, R., et al. 2015, A&A, 582, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. de Boer, J., Salter, G., Benisty, M., et al. 2016, A&A, 595, A114 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Dullemond, C. P., & Monnier, J. D. 2010, ARA&A, 48, 205 [NASA ADS] [CrossRef] [Google Scholar]
  20. Eisenhauer, F., Perrin, G., Brandner, W., et al. 2011, The Messenger, 143, 16 [NASA ADS] [Google Scholar]
  21. Flock, M., Fromang, S., Turner, N. J., & Benisty, M. 2016, ApJ, 827, 144 [CrossRef] [Google Scholar]
  22. Flock, M., Fromang, S., Turner, N. J., & Benisty, M. 2017, ApJ, 835, 230 [CrossRef] [Google Scholar]
  23. Gaia Collaboration (Brown, A. G. A., et al.) 2016, A&A, 595, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Gaia Collaboration (Brown, A. G. A., et al.) 2021, A&A, 649, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Garcia Lopez, R., Natta, A., Testi, L., & Habart, E. 2006, A&A, 459, 837 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Garcia Lopez, R., Tambovtseva, L. V., Schertl, D., et al. 2015, A&A, 576, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Grant, S. L., Espaillat, C. C., Brittain, S., Scott-Joseph, C., & Calvet, N. 2022, ApJ, 926, 229 [NASA ADS] [CrossRef] [Google Scholar]
  28. GRAVITY Collaboration (Abuter, R., et al.) 2017, A&A, 602, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. GRAVITY Collaboration (Sanchez-Bermudez, J., et al.) 2018, A&A, 618, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. GRAVITY Collaboration (Perraut, K., et al.) 2019, A&A, 632, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. GRAVITY Collaboration (Sanchez-Bermudez, J., et al.) 2021, A&A, 654, A97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. GRAVITY Collaboration (Garcia-Lopez, R., et al.) 2023, A&A, accepted, https://doi.org/10.1051/0004-6361/202245804 [Google Scholar]
  33. Harrington, D. M., & Kuhn, J. R. 2009, ApJS, 180, 138 [NASA ADS] [CrossRef] [Google Scholar]
  34. Hofmann, K. H., Weigelt, G., & Schertl, D. 2014, A&A, 565, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Hofmann, K.-H., Heininger, M., Schertl, D., et al. 2016, SPIE Conf. Ser., 9907, 99073H [NASA ADS] [Google Scholar]
  36. Hofmann, K. H., Bensberg, A., Schertl, D., et al. 2022, A&A, 658, A81 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Hone, E., Kraus, S., Kreplin, A., et al. 2017, A&A, 607, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Isella, A., & Natta, A. 2005, A&A, 438, 899 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Järvinen, S. P., Carroll, T. A., Hubrig, S., Ilyin, I., & Schöller, M. 2019, MNRAS, 489, 886 [Google Scholar]
  40. Kluska, J., Berger, J. P., Malbet, F., et al. 2020, A&A, 636, A116 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Koumpia, E., de Wit, W. J., Oudmaijer, R. D., et al. 2021, A&A, 654, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Kurosawa, R., Kreplin, A., Weigelt, G., et al. 2016, MNRAS, 457, 2236 [NASA ADS] [CrossRef] [Google Scholar]
  43. Labdon, A., Kraus, S., Davies, C. L., et al. 2019, A&A, 627, A36 [EDP Sciences] [Google Scholar]
  44. Lachaume, R. 2003, A&A, 400, 795 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Lacour, S., Dembet, R., Abuter, R., et al. 2019, A&A, 624, A99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Lapeyrere, V., Kervella, P., Lacour, S., et al. 2014, SPIE Conf. Ser., 9146, 91462D [Google Scholar]
  47. Lazareff, B., Berger, J. P., Kluska, J., et al. 2017, A&A, 599, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Le Bouquin, J. B., Absil, O., Benisty, M., et al. 2009, A&A, 498, L41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Long, F., Pinilla, P., Herczeg, G. J., et al. 2018, ApJ, 869, 17 [Google Scholar]
  50. Manoj, P., Maheswar, G., & Bhatt, H. C. 2002, MNRAS, 334, 419 [NASA ADS] [CrossRef] [Google Scholar]
  51. Millour, F., Meilland, A., Chesneau, O., et al. 2011, A&A, 526, A107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Montesinos, B., Eiroa, C., Mora, A., & Merín, B. 2009, A&A, 495, 901 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Mora, A., Merín, B., Solano, E., et al. 2001, A&A, 378, 116 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Muzerolle, J., D’Alessio, P., Calvet, N., & Hartmann, L. 2004, ApJ, 617, 406 [CrossRef] [Google Scholar]
  55. Newville, M., Otten, R., Nelson, A., et al. 2018, https://doi.org/10.5281/zenodo.1699739 [Google Scholar]
  56. Pinilla, P., Benisty, M., de Boer, J., et al. 2018, ApJ, 868, 85 [NASA ADS] [CrossRef] [Google Scholar]
  57. Pinte, C., Padgett, D. L., Ménard, F., et al. 2008, A&A, 489, 633 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Pohl, A., Benisty, M., Pinilla, P., et al. 2017, ApJ, 850, 52 [Google Scholar]
  59. Sanchez-Bermudez, J., Millour, F., Baron, F., et al. 2018, Exp. Astron., 46, 457 [NASA ADS] [CrossRef] [Google Scholar]
  60. Thiébaut, É., & Young, J. 2017, J. Opt. Soc. Am. A, 34, 904 [Google Scholar]
  61. Varga, J., Hogerheijde, M., van Boekel, R., et al. 2021, A&A, 647, A56 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Vioque, M., Oudmaijer, R. D., Baines, D., Mendigutía, I., & Pérez-Martínez, R. 2018, A&A, 620, A128 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Weigelt, G., Kraus, S., Driebe, T., et al. 2007, A&A, 464, 87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Weigelt, G., Hofmann, K. H., Schertl, D., et al. 2016, A&A, 594, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Weigelt, G., Hofmann, K. H., Schertl, D., et al. 2021, A&A, 652, A140 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

All Tables

Table 1

Stellar properties from (Vioque et al. 2018); the stellar luminosity is re-scaled to the Gaia EDR 3 distance (Gaia Collaboration 2021).

Table 2

Results from the K-band continuum image reconstruction.

Table A.1

HD58647 GRAVITY observations. UD diameters were obtained from JMMC SearchCal (Bonneau et al. 2006, 2011)

Table C.1

HD 58647 image reconstruction parameters.

Table D.1

K-band continuum model-fitting results with 1σ error bars.

Table E.1

Brγ line model-fitting results.

All Figures

thumbnail Fig. 1

plane coverage of our observations of HD 58647. See Table A.1 for a detailed description of the observations.

In the text
thumbnail Fig. 2

Example of HR GRAVITY interferometric observations of HD 58647 around the H I Brγ line position taken in January 2020 using the (J3-J2-G2-A0) configuration. From top to bottom: HR spectrum (red) and photospheric-corrected spectrum (blue); spectrally dispersed visibility amplitudes; differential phases; and closure phases. Different colours represent different projected baselines and baseline orientations (PAs; visibilities and differential phase panels) or triplets (closure phase panel), as indicated in the middle-top and bottom panel, respectively.

In the text
thumbnail Fig. 3

Results of the image reconstruction. Left: K-band continuum reconstructed images of HD 58647 at different wavelengths (reported in the labels). North is up, and east is to the left. The white hollow ellipse at the lower-left corner of the bottom panel represents the size of the clean beam. Contours represent 3, 5, and 10 σ pixel significance levels. Middle and right: squared visibilities (|V|2; middle panel) and closure phases (T3ϕ; right panel) as a function of the spatial frequency. The observed data with their corresponding error bars are represented as black dots, while the synthetic observables extracted from the reconstructed images are shown in colour.

In the text
thumbnail Fig. 4

Results of the image reconstruction. Left: H I Brγ line reconstructed images of HD58647. North is up, and east is to the left. Contours represent 3, 5, and 10σ. Middle and right: comparison of the observed continuum-subtracted visibility amplitudes and absolute phases (black circles with error bars) with the synthetic observables extracted from the reconstructed image.

In the text
thumbnail Fig. 5

Brγ line first moment map (M1) obtained from the line velocity-dispersed images. For comparison, the 1, 3, and 5σ K-band continuum contours are overplotted in cyan, black, and yellow. The dashed red and blue lines illustrate the orientation of the semi-major axis of the continuum and the Brγline emission, respectively.

In the text
thumbnail Fig. 6

Results of the Keplerian disk model fitting. Left: velocity maps from our Brγ line Keplerian disk model. From top to bottom, we show velocity maps at: ~ −110 km s−1, ~ −73 km s−1, ~ −36km s−1, ~0 km s−1, ~37 km s−1, ~74 km s−1, and ~110 km s−1. See the main text for details. Right: comparison between the observed continuum-corrected Brγ line absolute phases (black dots) and those obtained from the Brγ line Keplerian disk model shown in left panel (coloured dots).

In the text
thumbnail Fig. B.1

Merged spectrum, visibilities, differential phases, and closure phases of the January 27, 2020, GRAVITY interferometric data.

In the text
thumbnail Fig. B.2

Same as Fig. B.1, but for January 28, 2020.

In the text
thumbnail Fig. B.3

Same as Fig. B.1, but for February 4, 2020.

In the text
thumbnail Fig. B.4

Same as Fig. B.1, but for December 23, 2020.

In the text
thumbnail Fig. B.5

Same as Fig. B.1, but for February 10, 2021.

In the text
thumbnail Fig. B.6

Same as Fig. B.1, but for February 21, 2021.

In the text
thumbnail Fig. D.1

Comparison between the observed squared visibilities (|V|2) and closure phase signatures (T3ϕ) of HD 58647 (black dots and error bars) with the results from a 2D Gaussian brightness distribution (coloured dots).

In the text
thumbnail Fig. D.2

Brγ line model-fitting results (red dots) as a function of the spatial frequency at three different spectral channels. For comparison, the continuum-corrected Brγ line visibilities are overplotted in black.

In the text
thumbnail Fig. F.1

Astrometric displacements obtained from the Brγ line continuum-corrected differential phases.

In the text
thumbnail Fig. F.2

Example of the synthetic differential phases obtained from the Brγ Keplerian disk model (red dots). The observed continuum-corrected Brγ line phases are shown in black.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.